Measurements of Soil Carbon Dioxide Emissions from Two Maize Agroecosystems at Harvest under Different Tillage Conditions

In this study a comparison of the soil CO2 fluxes emitted from two maize (Zea mays L.) fields with the same soil type was performed. Each field was treated with a different tillage technique: conventional tillage (30 cm depth ploughing) and no-tillage. Measurements were performed in the Po Valley (Italy) from September to October 2012, covering both pre- and postharvesting conditions, by means of two identical systems based on automatic static soil chambers. Main results show that no-tillage technique caused higher CO2 emissions than conventional tillage (on average 2.78 and 0.79 μmol CO2 m−2 s−1, resp.). This result is likely due to decomposition of the organic litter left on the ground of the no-tillage site and thus to an increased microbial and invertebrate respiration. On the other hand, fuel consumption of conventional tillage technique is greater than no-tillage consumptions. For these reasons this result cannot be taken as general. More investigations are needed to take into account all the emissions related to the field management cycle.


Introduction
Soils represent the main terrestrial reservoir of organic carbon (OC), containing almost three times more carbon than plants biomass [1].
The release of carbon dioxide (CO 2 ) by soil respiration processes, mainly related to the degradation of organic matter in soils, represents the second largest component of the global carbon cycle [2,3] and can play a relevant role in the climate change.
Bond-Lamberty and Thomson estimated in 98 PgC year-1 the current global emission flux from soil over the Earth's land surface, which is 10 times higher than the total anthropogenic emission fluxes from fuel combustion [4].
Depending on the agricultural practices of soil management, soils may be important sources or sinks of atmospheric carbon, with consequent implications and effects on a global scale [5][6][7].
Tillage, even if it incorporates crop residues, is considered a practice that favours CO 2 emissions from arable lands since it improves the ventilation of the top-soils inducing a fast biological oxidation of the organic matter [8][9][10].
Conservative agriculture practices such as the no-tillage and the minimum tillage techniques, on the contrary, are considered less carbon emissive [11][12][13][14], because the mechanical agitation of the soil is avoided or minimized. The new season crop is sown over the previous crop residues in thin furrows, without ploughing. The decomposition of the crop residues lead to an increase of the organic carbon content on the top layer of the soil, favouring its aggregation and stability. However, as observed by Six et al. [15], in a five-year experiment of no-tillage practice in arid climates, the carbon 2 The Scientific World Journal content of the deeper layers may decrease, accelerating the decomposition in the surface layer [16].
The exchange of CO 2 between arable soils and atmosphere is only one of the aspects of the complexity of the carbon budget of an agroecosystem.
Nevertheless, accurate measurements of CO 2 effluxes are required to assess whether a crop management technique is better than another in reducing the CO 2 emission from soils.
Several approaches have been used to assess the exchange of CO 2 between the ecosystem and the atmosphere [17][18][19], each of them with their own advantages and disadvantages [20,21]. The direct methods include the use of static enclosures ("chambers") on the soil surface [22,23] and the analysis of the rate of the CO 2 accumulation (or removal) in the chamber headspace over a certain period of time after enclosing the surface.
The static chamber technique is relatively cheap and easy, and the automation of the process allows the system to perform continuous measurements of soil respiration. On the other hand both static and dynamic chambers methods were criticized for altering soil environment and air-soil pressure gradient thus affecting the fluxes [24] and for removing the turbulent fluctuations [25,26]. Many different sources of error are possible when chambers are used [27] and underestimations of soil respiration may happen [28] when the rigid collar over which the chambers are mounted is inserted in the soil so deeply that it cuts roots. However different solutions based on proper chamber designs, data analyses, and spatial and temporal sampling regimes were presented for correcting most of these problems [29,30].
This work was developed in the context of the AGRICO 2 LTURA Project of the Lombardy Region, which was aimed at investigating CO 2 fluxes from cultivated soil in an intensive agriculture area of the Po Plain (Italy).
Soil CO 2 fluxes measurements taken in two maize fields managed with different tillage techniques are presented to provide parameters for the soil-atmosphere CO 2 exchange modeling.
In particular, measurements of CO 2 fluxes after harvest and for the following two months were analysed evaluating their mean daily patterns and the influence of different environmental factors such as, litter presence, microbial biomass, and wind turbulence. However, it is important to underline that this paper is not intended to provide the carbon budget of the whole crop cycle in the two agroecosystems, but it is referred to a comparison of measurements made only for a short period of time (two months) around harvest.

Sites Description and Agronomical Management.
This study was conducted in the agricultural district of Lodi, Italy, located at the so called "basal level of the plain", just north of the Po river in the Po Valley (Italy).
Two farms located 5 km apart were selected, one located in the San Martino in Strada village (no-Tillage agriculture site, NT site from now on, 45 ∘ 20.64 E, 66 m a.s.l.). The dominant cropping system for both farms consisted of a cereal forage rotation, which is typical for the region due to its extensive livestock production.
In each farm a maize field was chosen for the measurements.
The two maize fields were managed with a wheat-maize rotation and were subjected to different ploughing treatments.
At the NT site, as was the case for the previous ten years of management, a no-tillage agriculture practice was employed. Maize seeds were directly planted into the soil covered by wheat residues and weeds, and two herbicide treatments were performed after maize sowing. N and P fertilizations were performed with chemical fertilizers at the sowing and at the emergence of the 6th-7th leaf.
The field at the CT site experienced a conventional 30 cm depth ploughing followed by a harrow clods reduction before the sowing. Then, two herbicide treatments were performed. Fertilization was made with cattle manure before ploughing and with chemical P fertilizer at sowing.
Irrigation in both fields was performed with flooding until soils reached their field capacity.
After the maize harvest, which occurred at the end of August at the NT site and at the beginning of September at CT site, the maize stalks were left on the soils in both sites.
Detailed information on crop management, including fertilization treatments, irrigations, herbicide applications, mechanical treatments, harvest and relative dates can be found in Table 1.

Soil Measurements and Characteristics.
Soil characteristics of the selected fields were assessed on February 2012 [32] with the Area Frame Randomised Soil Sampling methodology (AFRSS) described by Stolbovoy et al. [33]. In each site one 20 m × 20 m sampling area was randomly selected, and two additional areas were identified at 80 and 40 m, respectively, from the first one, following an "L" shape disposition. In each of these three areas, soil and litter samples were taken from 9 sampling points dislocated over the whole surface with a cross shape design (for detail refer to [33]). In each sampling point the litter was collected from a surface of 0.63 m 2 and the soil samples were taken at three different depths: 0-10, 10-20, and 20-30 cm.
The 9 samples were then merged and mixed into a single sample per area, and the three mixed samples of each field were then sent to the lab for the analysis.
The organic carbon (OC) in each soil and litter sample was assessed by means of the Walkley-Black method [34] which consists of the oxidation of carbon with acidic dichromate (Cr 2 O 7 2− ) followed by the titration of the excess of dichromate with ferrous sulfate. The OC is calculated from the difference between the total dichromate added and the amount of dichromate left unreacted after OC oxidation.
The organic carbon (OC) in the microbial biomass was assessed by the chloroform fumigation extraction technique The Scientific World Journal 3  [36]. Details can be found in the cited references. The two sites presented very similar OC content of the topsoils (0-30 cm), not statistically different between them: 9.5 gC/kg at NT site and 10.3 gC/kg at CT site, on average. As a consequence, the total carbon stocks (defined as the total amount of the OC stored in the first 30 cm of soil depth) were comparable at two sites: 44.8 t/ha at NT site and 48.7 t/ha at CT site.
The organic carbon of the crop residues in the NT site litter was 4.8 t/ha, an amount corresponding to 10.7% of the total carbon stock in the top-soil (0-30 cm depth). Table 2 summaries the OC content in the soil, the litter (where present), and in the microbial biomass of the two sites.

Flux Measurement Technique.
A system of three automatic static chambers (SASSFLUX, Ecometrics srl, Italy) was deployed at the NT site from the middle of July (with green maize) to the end of October 2012, and a second system with four chambers was deployed at the CT site from the beginning of September (with senescent maize) to the end of the measuring campaign (October 2012).
Each system consisted of a central unit where an infrared gas analyzer (IRGA), a pump, 8 solenoid valves and the controlling system were located, and up to 4 measuring chambers deployed to a distance of 5 m from the central unit.
The measuring chambers ( Figure 1) were made of a boxshaped lid of transparent Plexiglas (dimensions 40 × 40 × 10 cm) which is mounted on a steel frame (collar) that is inserted into the soil at a depth of 8 cm, delimiting a measuring soil surface of 35 × 35 cm 2 . The total air volume trapped by each chamber was 16.5 litres taking into account the thickness of the seals.
The chambers were placed below the maize canopy between the rows at a distance of 5 m from the SASSFLUX central unit heading toward the 4 cardinal directions (N excluded at the NT site) and then left until the harvest and the following postharvest period. Measurements restarted one week after the harvest at the NT site and 3 days after the harvest at the CT site. The chambers sequentially closed over their metal collars and the air within was sampled by a 10 L min −1 membrane pump trough a Teflon tube (4 mm inner diameter) to the IRGA (Carbocap 343, Vaisala, SF) which provided one CO 2 concentration measurement every second until the chamber lid opened and the soil surface was reexposed to the atmosphere.
To avoid pressure alterations inside the chambers, the sampled air was redirected to the chamber through a second teflon tube after the IRGA analysis.
A personal computer controlled the solenoid valves system which allowed the IRGA to cyclically analyse the air from the different chambers.
The tubing flushing time before and after each chamber closure, and the duration of each measurement session, was set in order to get 4 cycles of measurement every hour (i.e., each chamber was closed 4 times per hour and the closure time was about 3.5 minutes per cycle, including flushing, lid opening/closure and settling time).
The very short closure time allowed the soil disturbance to be kept as low as possible and to gather a high number of measurement replicates in order to get the hourly variation of CO 2 fluxes. A very similar system is also described by Jassal et al. [23]. The two SASSFLUX systems were calibrated in laboratory before their deployment in the fields and their instrumental drift was checked at the end of the campaign by performing The proper operation of the IRGA and the pneumatic system during the field campaign was also checked daily by comparing the daily courses of the soil level atmospheric CO 2 concentration measured by the two SASSFLUX devices in the two sites.
All of these checks did not reveal any significant difference between the two systems.
Moreover, the SASSFLUX system proved to provide comparable data with other measuring techniques such as vials sampling from the headspace of fixed collars inserted into the soil and subsequent gas-chromatography analysis or sampling made in the same collars with a Brüel-Kjiaer photoacoustical device [37].

Soil CO 2 Flux Calculation and Data Analysis.
The instantaneous flux density of CO 2 from the soil in each closed chamber, Φ CO 2 ( ), is the quantity of CO 2 that, flowing through the soil surface , is able to increase the CO 2 concentration of the trapped air volume of a quantity equal to [CO 2 ] in the time unit [22]: If, as in our case, the CO 2 concentration is measured as mixing ratio (ppm), the addition of the factor term /( ), where P and T are the atmospheric pressure and temperature (Pa and K), and R is the universal gas constant (J mol −1 K −1 ), can The Scientific World Journal 5 be used to convert the CO 2 concentrations to mass density ( mol m −3 ), and the expression of the flux as mol m −2 s −1 : The flux of the undisturbed soil can be derived by estimating Φ CO 2 (0), that is, the CO 2 flux just after the closure of the cambers.
In this work Φ CO 2 (0) was calculated by linear interpolation of the CO 2 concentration values versus time during the first 70 seconds of sampling after the chamber's closure in order to get the CO 2 accumulation velocity ( [CO 2 ]/ ) from the slope of the interpolation line. The selected time window was chosen to ensure the linearity of the first part of the CO 2 accumulation curve.
In our conditions this simple linear fit guaranteed excellent results even in the (worst) cases of very intense CO 2 emission rates, when soil is treated with a significant amount of manure.
Two quality controls were employed to check the illconditioned samples. The first one regarded the reaching of the full scale value of the IRGA, set to 1000 ppm. When the last CO 2 concentration measured during the 70 seconds time window exceeded 999 ppm, the sample was considered "saturated" and then discarded. The second one regarded the linearity of the CO 2 accumulation curve in the chosen time window. When the slope of the interpolation line of the first 70 data was not statistically significant ( test with − 2 = 68 degrees of freedom resulting in a probability > 0.01), the sample was rejected. Significant deviation from linearity generally occurred when the lid failed to close tightly, typically when some plant material fell over the collar seals during the chamber closure.
The Φ CO 2 (0) of the samples which passed the quality check were stored in a common database for the following analysis. For comparison purposes, the arithmetical average of the Φ CO 2 (0) obtained by all the chambers of each SASS-FLUX system in each site (hereafter called CO 2 fluxes for simplicity) was calculated for selected time periods.
With the aim of comparing the diurnal emission patterns of the two sites, all the CO 2 fluxes detected by all the SASSFLUX chambers (three for NT site and four for CT site) in the same hour of different days were averaged for a given period in each site. Whenever fluxes were averaged, the standard deviation ( ) of the averaged values was calculated and presented. When the averaged fluxes were intended to be meant as an estimation of the real (unknown) fluxes of each hour of the day (and not simply to describe the features of the samples), the associated standard error of the mean was calculated (= /√ , with the number of data used for the calculation of the mean) and presented. The statistical significance of the difference between the daily courses of the CO 2 fluxes in different periods was tested by means of a oneway ANOVA repeated 24 times, one per each hour, on all the available CO 2 flux data (all the cambers, all the available data per each hour) grouped by the different periods (the ANOVA factor). Differences with a probability greater than 0.05 were considered not significant.

Data
Capturing. Several interruptions in data sampling, due to technical problems or field care needs, occurred during the measuring campaign. The total data capturing, that is, the ratio between the number of the collected samples and the maximum obtainable with a continuously operative system, resulted 88% in the NT site and 91% in the CT site. However, the samples, which passed the two adopted quality control criteria, were 74% and 66% at NT site, and CT site, respectively.

Meteorological Measurements and Weather Conditions.
Ancillary measurements of air temperature and humidity (HD9008 DeltaOhm, UK) above and below the crop canopy, barometric pressure (PTB101b Vaisala, Finland), ground heat fluxes (HFP01SC, Hukseflux, Netherlands), and soil moisture (CS616 Campbell Sci., USA; EC5 Decagon, USA) were also monitored at the NT site every 30 s and stored as 15 minutes mean with a CR1000 datalogger (Campbell Sci., USA).
The atmospheric turbulence in the two was measured by means of two colocated eddy covariance measuring stations equipped with a sonic anemometer (USA-1, Metek, Germany) mounted on a mast at 3 m above the ground level.
The turbulence parameter * (friction velocity) was calculated from the covariance of the vertical and the horizontal wind speed, by following the procedure reported in Gerosa et al. 2003 [38].
The weather conditions of the measuring campaign are presented in Figure 2 and summarised in Table 3. July and August were sunny and hot with rare strong precipitation events, while September and October were rainy, with quite variable weather in September and more foggy weather in October. The volumetric soil water content of the first 30 cm of soil ranged between 0.19 and 0.  The Scientific World Journal During July-August the soils were in thermal equilibrium with the below canopy atmosphere (Table 3), while they lost energy in September-October following the removal of the canopy with the harvests and the air cooling. Figure 3 shows an example of one week of measurements at both sites in September. The measurement gaps were due to IRGA saturation conditions, particularly frequent at the NT site during the night.

Soil CO 2 Fluxes.
The NT site experienced higher CO 2 emission fluxes than the CT site. Around 4 mol CO 2 m −2 s −1 with peaks of 7 mol m −2 s −1 and minimum of 2-3 mol m −2 s −1 were recorded at the NT site, and about 1 mol CO 2 m −2 s −1 with peaks of 2 mol m −2 s −1 were recorded at the CT site.
Moreover the NT site showed a higher flux variability from chamber to chamber, reflecting a greater spatial variability of the organic matter at the microscale. Conversely, the CO 2 fluxes at the CT site were quite homogeneous, as an expected result of the top soil homogenisation due to the mechanical treatments (see Table 1).
Taking into account only the concomitant measuring time period of both sites (September-October), the average flux at NT site was 2.78 mol CO 2 m −2 s −1 while at CT site was 0.79 mol m −2 s −1 , 3.5 times lower. This difference is also evidenced by the average daily courses of CO 2 fluxes of the two sites (Figure 4(b)), which resulted significantly different ( < 0.01) for every hour of the day.
The CO 2 fluxes were much more intense in the summer months (June-July) than in the following fall months (Figure 4(a), only comparison on NT site was possible) when they decreased by 42% on average. This was surely caused by the general air cooling, but was also the result of the increased soil exposure to the air after the harvest.
The greater flux variability at the NT site is confirmed by the higher standard deviation of the mean flux estimation for each hour of the day: on average 1.32 mol m −2 s −1 in the NT site and 0.50 mol m −2 s −1 in the CT site in the common period (Figure 4(b)).
The higher CO 2 emission fluxes registered at the NT site compared to the CT site was an unexpected result, since the application of no-tillage agriculture practices should have lowered carbon dioxide emissions at the NT site.
This result does not exclude that high CO 2 losses might have occurred at CT site just after the tillage, as observed by other studies [8]. However, since our measurements were made long after the ploughing treatment (4 months later), they might have missed this initial emission peak.

Potential
Role of the Litter. The increase in CO 2 emissions detected at the NT site may find an explanation if we take into account the presence of a thick organic litter (2.02 ± 0.4 kg/m 2 , Table 2) that could have caused a higher decomposition activity than the one experienced at the CT site, and a different distribution of the organic carbon content in the profile of the first 30 cm of soil between the two sites.
The presence of an organic litter in the NT site, completely absent in the CT site, is the direct consequence of the application of the no-tillage agriculture technique in the past 10 years. Residues of the previous crops were left on the field even during the following crop growing season and never buried by ploughing.
The stimulation of the decomposition activity due to the presence of litter has been well described by Pengthamkeerati et al. [39], Iqbal et al. [40], Li et al. [41], Oorts et al. [42], and Maraseni and Cockfield [43]. For example Pengthamkeerati et al. [39] found that claypan soils with the application of a poultry litter, in an incubation experiment, emitted about the double amount of CO 2 than the same soils without litter. Li et al. [41], for example, found that the cumulative carbon emissions from mollisols amended with maize crop residues The Scientific World Journal applied on the soil surface was up to 4 times greater than the not-amended soils. The field trials of Pengthamkeerati et al. [39] confirmed their lab findings, but a great interannual variability was highlighted. These results are well in agreement with those found by De Neve and Hofman [44] with crop residues used as an organic amendment. Morever the stimulation of CO 2 emission due to crop residues is independent on the nitrogen addition, as Iqbal et al. [40] found in red (Ultisol) soils after the application of straws and different amounts of nitrogen fertilisers.

Soil Organic Carbon Profile and Microbial Biomass.
The OC content of the topsoil in the two sites was comparable ( Table 2), but the OC content in the surface layer (0-10 cm) of the NT site (13.0 g/kg) was slightly higher than CT site (11.4 g/kg). However, this difference did not give a statistically significant result and, as a consequence, cannot explain the higher CO 2 fluxes detected at NT site. Conversely, the total carbon content of the microbial biomass in the top soil mineral fraction (0-30 cm) was almost twice higher at the NT site than at CT site, resulting in 97.4 gC/g versus 55.8 gC/g, but once again this difference was not statistically significant.
Nevertheless, it is important to underline that there is a clear relation between increased microbial biomass and increased soil CO 2 effluxes, as for example, shown in a recent work by Han et al. [45]. Moreover Sowerby et al. [46] clearly demonstrated that increased soil microbial biomass, rather than increased microbial activity, was responsible for the rise of soil CO 2 effluxes observed in their experiment.

Influence of Weather.
The main driver of the CO 2 soil fluxes at both sites was temperature, as expected for the microbial metabolism. The correlation of the daily mean fluxes with the daily mean air temperatures, a good proxy of the soil temperature [47,48], was significant at both sites ( 2 = 0.846 < 0.001 and 2 = 0.724 < 0.001, at NT and CT resp.) even if the flux intensities were very different (Figures 5(a) and 5(c)).
On the contrary the CO 2 flux decreased at both sites at increasing soil moisture ( 2 = 0.39, < 0.001 at NT site, and 2 = 0.17, < 0.05 at CT site), following rain events. However, once the effect of the temperature was removed from the data and the residuals were considered for the regression with the soil moisture, this decrease completely disappeared at both sites ( 2 = 0.034, = 0.169 at NT and 2 = 0.0103 = 0.539 at CT site, Figures 5(b) and 5(d)). Therefore, the decrease of CO 2 observed with increasing soil water content is exclusively due to the reduction of air temperature that is linked to rain events, and which also affects the following days.

Daily
Patterns of Soil CO 2 Fluxes. In general, the diurnal fluctuations of the CO 2 fluxes were weak. Nevertheless little differences of the average daily courses in the two periods were detected. In July-August the maximum flux intensity occurred in the central part of the day, while a slight flux minimum appeared in the same part of the day at September-October (Figure 4). The maximum flux observed at noon in July-August reflects the roots activity of the plants during the daylights. Indeed, this characteristic pattern disappeared in the postharvest period as shown in Figure 6, where the flux average daily courses recorded one week before and one week after the harvest is presented. The relative contribution of the root respiration to the diurnal CO 2 fluxes can be estimated from the observation of the mean diurnal flux decrease (1 mol CO 2 m −2 s −1 at the NT site and 0.8 mol CO 2 m −2 s −1 at the CT site).
A relative minimum flux is clearly visible during the light hours after the harvest ( Figure 5(b)) and lasts for the whole second period of measurements (Figure 4(b)).

Influence of Weedings.
Weeds can play a possible role in determining the observed diurnal relative minimum of CO 2 fluxes. In fact, the great increase of the sunlight at ground level, which followed the removal of the canopy with the harvest, could have stimulated the photosynthetic activity and growth of weeds and microalgae within the transparent Plexiglas chambers.
The role of weeds was confirmed with a comparison of the measurements taken three days before and three days after a chemical weeding treatment made at both sites in October ( Figure 6). The diurnal flux minimum at the CT site ( Figure 6(b)) completely disappeared afterthe herbicide treatment, while the relative flux depression was greatly reduced at the NT site after the same herbicide application (Figure 6(a)). Moreover, a slight reduction of the evening and night fluxes was observed at both sites.

3.7.
Influence of Turbulence. The analysis of turbulence at ground level can help us to understand the residual diurnal flux minimum not explained by the weeds' photosynthesis at the NT site, after the herbicide treatment.
The direct exposure to the atmosphere of the two studied soils after the harvest caused an increase of the air turbulence at ground level, which rose from negligible values under The Scientific World Journal the canopy to appreciable values after its removal (Figure 8), following dynamics that also other works reported [49][50][51].
Our results show that despite the CO 2 fluxes were very different in their absolute values in the two postharvest months at NT site (Figure 6), the diurnal flux "depression" of about 1 mol m −2 s −1 shown by the daily courses of September and October was very similar. At the same time, the mean daily patterns of the air turbulence ( * ) in September and October were comparable during the daylight hours ( Figure 7) and were identical during the night-time, thus suggesting a possible role of the air turbulence.
The effect of turbulence on the extraction of CO 2 from the topsoil was well described by Takle et al. [52] who highlighted that the rapid fluctuations in the static pressure fields introduced by wind interactions with terrain may lead literally to pressure pumping of CO 2 at the surface.  Moreover, in the presence of litter on the surface, like for example, in the NT site, the enhanced surface rugosity introduced by the litter can also contribute in increasing this process. In fact, Reicosky et al. [10] working with histosols with a bulk density between 0.23 and 0.27 Mg m −3 at the surface layer, a density which is barely different from the organic litter density observed at our site (0.14 ± 0.042 Mg m −3 ), has shown how wind speed can deplete the CO 2 concentrations in highly porous organic soils.
It is likely that analogous processes occurred also at the NT site, where the litter layer could have acted as a buffer for the CO 2 emitted from the underlying mineral topsoil. During the daylight hours, when the chambers were open, the increased turbulence could have caused a more intense removal of CO 2 from the litter layer towards the atmosphere than during the night-time, thus leaving less exchangeable CO 2 into the buffer. Hence, after the chamber's closure, the flux was lower during the daylight hours than during the night hours when the atmospheric turbulence was negligible (Figure 7). It is worth noticing that the particular SASSFLUX design (Figure 1) does not interfere with the wind action on the topsoil, because there are no emerging parts from the soil when the lid is open. It is likely therefore, that the relative minimum after the harvest is the result of the combination of two processes, photosynthesis/respiration of the weeds and air turbulence increase caused by litter.

Conclusion
In this paper CO 2 emissions from the soil of two maize fields managed with conventional tillage and no-tillage techniques were compared at the time of maize harvest.
The NT site showed CO 2 emissions up to 3.5 times higher than the CT site during the two months following the maize harvest. This result was related to the presence of a conspicuous organic litter on the NT site field's surface due to the accumulation of crop debris and residues in the last 10 years, and thus to an increased decomposition activity of the microbial community. However, the CO 2 emission differences between the two fields slowly declined approaching the fall season, together with the flux intensity which was clearly driven by the air (and soil) temperature.
The spatial flux variability of the NT site field was higher than the CT site, reflecting a more irregular distribution of the organic matter on the soil surface and in the topsoil, which contrasts with the continuous redistribution of the organic matter along the soil profile in the CT site field.
However, results from this case-study do not allow us to get any conclusions on the environmental sustainability of the two different management practices, regarding their impacts on GHG emission and global warming, since they are focused only on one small aspect of the carbon budget of the studied agroecosystems.
Only a careful analysis of the carbon emissions during the whole life-cycle of the agroecosystem would allow a meaningful conclusion. And the simple inclusion in the budget of the CO 2 emitted by the tractors employed for the mechanical treatments at the CT site (around 40 liters of fuel per hour) [53] would have radically changed any conclusion on the impacts on GHG emissions and global warning of the two different management practices.
The presented data can contribute to the development of high-time resolved modules for the CO 2 exchange processes at the soil surface, to be included in more general models aimed at evaluating the total carbon budget of agroecosystems at regional scale.