Printer Friendly

Measurements of soil carbon dioxide emissions from two maize agroecosystems at harvest under different tillage conditions.

1. 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 (C[O.sub.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.sup.-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 maybe important sources or sinks of atmospheric carbon, with consequent implications and effects on a global scale [5-7].

Tillage, even if it incorporates crop residues, is considered a practice that favours C[O.sub.2] emissions from arable lands since it improves the ventilation of the top-soils inducing a fast biological oxidation of the organic matter [8-10].

Conservative agriculture practices such as the no-tillage and the minimum tillage techniques, on the contrary, are considered less carbon emissive [11-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 content of the deeper layers may decrease, accelerating the decomposition in the surface layer 16].

The exchange of C[O.sub.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 C[O.sub.2] effluxes are required to assess whether a crop management technique is better than another in reducing the C[O.sub.2] emission from soils.

Several approaches have been used to assess the exchange of C[O.sub.2] between the ecosystem and the atmosphere [17-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 C[O.sub.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 AGRIC[O.sub.2]LTURA Project of the Lombardy Region, which was aimed at investigating C[O.sub.2] fluxes from cultivated soil in an intensive agriculture area of the Po Plain (Italy).

Soil C[O.sub.2] fluxes measurements taken in two maize fields managed with different tillage techniques are presented to provide parameters for the soil-atmosphere C[O.sub.2] exchange modeling.

In particular, measurements of C[O.sub.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.

2. Materials and Methods

2.1. 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[degrees]16'45.53" N, 9[degrees]32'59.62"E, 69 m a.s.l.) and the other in the Secugnago village (conventional tillage agriculture site, CT site from now on, 45[degrees]14'24.68"N, 9[degrees]35'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 soils of the two fields were both Hapli-Cutanic Luvisols (IUSS-WRB, 2007) with a coarse texture, a clay content of 10-14%, and a pH of 5.5-5.9 [31].

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.

2.2. 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 x 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.sup.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.sub.2][O.sub.7.sup.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 (CFE) described by Vance et al. [35] and implemented by Tate et al. [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.

2.3. 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 box-shaped lid of transparent Plexiglas (dimensions 40 x 40 x 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 x 35 [cm.sup.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.sup.-1] membrane pump trough a Teflon tube (4 mm inner diameter) to the IRGA (Carbocap 343, Vaisala, SF) which provided one C[O.sub.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 C[O.sub.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 a one day period of cross-measurements in the same field at the end of the measuring campaign.

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 C[O.sub.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 Briiel-Kjiaer photoacoustical device [37].

2.4. Soil C[O.sub.2] Flux Calculation and Data Analysis. The instantaneous flux density of C[O.sub.2] from the soil in each closed chamber, [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], is the quantity of C[O.sub.2] that, flowing through the soil surface A, is able to increase the C[O.sub.2] concentration of the trapped air volume V of a quantity equal to d[C[O.sub.2]] in the time unit dt [22]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (1)

If, as in our case, the C[O.sub.2] concentration is measured as mixing ratio (ppm), the addition of the factor term P/(RT), where P and T are the atmospheric pressure and temperature (Pa and K), and R is the universal gas constant (J [mol.sup.-1] [K.sup.-1]), can be used to convert the C[O.sub.2] concentrations to mass density ([micro]mol [m.sup.-3]), and the expression of the flux as [micro]mol [m.sup.-2] [s.sup.-1]:

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (2)

The flux of the undisturbed soil can be derived by estimating [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], that is, the C[O.sub.2] flux just after the closure of the cambers.

In this work [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] was calculated by linear interpolation of the C[O.sub.2] concentration values versus time during the first 70 seconds of sampling after the chamber's closure in order to get the C[O.sub.2] accumulation velocity (d[C[O.sub.2]]/dt) from the slope of the interpolation line. The selected time window was chosen to ensure the linearity of the first part of the C[O.sub.2] accumulation curve.

In our conditions this simple linear fit guaranteed excellent results even in the (worst) cases of very intense C[O.sub.2] emission rates, when soil is treated with a significant amount of manure.

Two quality controls were employed to check the ill-conditioned samples. The first one regarded the reaching of the full scale value of the IRGA, set to 1000 ppm. When the last C[O.sub.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 C[O.sub.2] accumulation curve in the chosen time window. When the slope of the interpolation line of the first 70 data was not statistically significant (F test with n - 2 = 68 degrees of freedom resulting in a probability P > 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 [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] 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 [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] obtained by all the chambers of each SASSFLUX system in each site (hereafter called C[O.sub.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 C[O.sub.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 (<r) 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 (=[sigma]/[square root of n], with n 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 C[O.sub.2] fluxes in different periods was tested by means of a one-way ANOVA repeated 24 times, one per each hour, on all the available C[O.sub.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 P greater than 0.05 were considered not significant.

2.5. 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.

2.6. 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 [u.sup.*] (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.37 [m.sup.3.sub.water] [m.sup.-3.sub.soil] in the drier months of July and August as a result of the irrigation practice, while in the following September and October months the soil moisture remained quite stable between 0.27 and 0.34 [m.sup.3.sub.water] [m.sup.-3.sub.soil]. The average soil moisture of the two sites during the common measuring period (September-October) was very similar: 0.294 [m.sup.3.sub.water] [m.sup.-3.sub.soil] at the NT site and 0.271 [m.sup.3.sub.water] [m.sup.-3.sub.soil] at the CT site.

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.

3. Results and Discussion

3.1. Soil C[O.sub.2] Fluxes. 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.

The NT site experienced higher C[O.sub.2] emission fluxes than the CT site. Around 4 [micro]mol C[O.sub.2] [m.sup.-2] [s.sup.-1] with peaks of 7 [micro]mol [m.sup.-2] [s.sup.-1] and minimum of 2-3 [micro]mol [m.sup.-2] [s.sup.-1] were recorded at the NT site, and about 1 [micro]mol C[O.sub.2] [m.sup.-2] [s.sup.-1] with peaks of 2 [micro]mol [m.sup.-2] [s.sup.-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 C[O.sub.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 [micro]mol C[O.sub.2] [m.sup.-2] [s.sup.-1] while at CT site was 0.79 [micro]mol [m.sup.-2] [s.sup.-1], 3.5 times lower. This difference is also evidenced by the average daily courses of C[O.sub.2] fluxes of the two sites (Figure 4(b)), which resulted significantly different (P < 0.01) for every hour of the day.

The C[O.sub.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 [micro]mol [m.sup.-2] [s.sup.-1] in the NT site and 0.50 [micro]mol [m.sup.-2] [s.sup.-1] in the CT site in the common period (Figure 4(b)).

The higher C[O.sub.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 C[O.sub.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.

3.2. Potential Role of the Litter. The increase in C[O.sub.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 [+ or -] 0.4kg/[m.sup.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 C[O.sub.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 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 C[O.sub.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.

3.3. 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 C[O.sub.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 [micro]gC/g versus 55.8 [micro]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 C[O.sub.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 C[O.sub.2] effluxes observed in their experiment.

3.4. Influence of Weather. The main driver of the C[O.sub.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 ([R.sup.2] = 0.846 P < 0.001 and [R.sup.2] = 0.724 P < 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 C[O.sub.2] flux decreased at both sites at increasing soil moisture ([R.sup.2] = 0.39, P < 0.001 at NT site, and [R.sup.2] = 0.17, P < 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 ([R.sup.2] = 0.034, P = 0.169 at NT and [R.sup.2] = 0.0103 P = 0.539 at CT site, Figures 5(b) and 5(d)).

Therefore, the decrease of C[O.sub.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.

3.5. Daily Patterns of Soil C[O.sub.2] Fluxes. In general, the diurnal fluctuations of the C[O.sub.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 post-harvest 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 C[O.sub.2] fluxes can be estimated from the observation of the mean diurnal flux decrease (1 [micro]mol C[O.sub.2] [m.sup.-2] [s.sup.-1] at the NT site and 0.8 [micro]mol C[O.sub.2] [m.sup.-2] [s.sup.-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)).

3.6. Influence of Weedings. Weeds can play a possible role in determining the observed diurnal relative minimum of C[O.sub.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 after the 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 canopy to appreciable values after its removal (Figure 8), following dynamics that also other works reported [49-51].

Our results show that despite the C[O.sub.2] fluxes were very different in their absolute values in the two post-harvest months at NT site (Figure 6), the diurnal flux "depression" of about 1 [micro]mol [m.sup.-2] [s.sup.-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 ([u.sup.*]) 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 C[O.sub.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 C[O.sub.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.sup.-3] at the surface layer, a density which is barely different from the organic litter density observed at our site (0.14 [+ or -] 0.042 Mg [m.sup.-3]), has shown how wind speed can deplete the C[O.sub.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 C[O.sub.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 C[O.sub.2] from the litter layer towards the atmosphere than during the night-time, thus leaving less exchangeable C[O.sub.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.

4. Conclusion

In this paper C[O.sub.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 C[O.sub.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 C[O.sub.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 C[O.sub.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 C[O.sub.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.

http://dx.doi.org/10.1155/2014/141345

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgment

This research and this publication were funded by the Catholic University of the Sacred Heart in the frame of its program of promotion and dissemination of the scientific research.

References

[1] A. Fischlin, G. F. Midgley, J. Price et al., "Ecosystems, their properties, goods, and services," in Climate Change 2007: Impacts, Adaptation and Vulnerability, M. L. Parry, O. F. Canziani, J. P. Palutikof, P. J. van der Linden, and C. E. Hanson, Eds., pp. 211-272, Cambridge University Press, Cambridge, UK, 2007.

[2] W. H. Schlesinger and J. A. Andrews, "Soil respiration and the global carbon cycle," Biogeochemistry, vol. 48, no. 1, pp. 7-20, 2000.

[3] A. M. Silva-Olaya, C. E. P. Cerri, N. La Scala Jr., C. T. S. Dias, and C. C. Cerri, "Carbon dioxide emissions under different soil tillage systems in mechanically harvested sugarcane," Environmental Research Letters, vol. 8, no. 1, Article ID 015014, 2013.

[4] B. Bond-Lamberty and A. Thomson, "Temperature-associated increases in the global soil respiration record," Nature, vol. 464, no. 7288, pp. 579-582, 2010.

[5] R. Lal, "Soil carbon sequestration to mitigate climate change," Geoderma, vol. 123, no. 1-2, pp. 1-22, 2004.

[6] R. Lal, "World cropland soils as a source or sink for atmospheric carbon," Advances in Agronomy, vol. 71, pp. 145-191, 2001.

[7] A. R. Mosier, "Soil processes and global change," Biology and Fertility of Soils, vol. 27, no. 3, pp. 221-229, 1998.

[8] R. W. Gesch, D. C. Reicosky, R. A. Gilbert, and D. R. Morris, "Influence of tillage and plant residue management on respiration of a Florida Everglades Histosol," Soil and Tillage Research, vol. 92, no. 1-2, pp. 156-166, 2007

[9] D. Reicosky, Tillage-Induced C[O.sub.2] Emissions and Carbon Sequestration: Effect of Secondary Tillage and Compaction, Springer, 2003.

[10] D. C. Reicosky, R. W. Gesch, S. W Wagner, R. A. Gilbert, C. D. Wente, and D. R. Morris, "Tillage and wind effects on soil C[O.sub.2] concentrations in muck soils," Soil and Tillage Research, vol. 99, no. 2, pp. 221-231, 2008.

[11] N. La Scala, D. Bolonhezi, and G. Pereira, "Short-term soil C[O.sub.2] emission after conventional and reduced tillage of a no-till sugar cane area in Southern Brazil," Soil and Tillage Research, vol. 91, no. 1-2, pp. 244-248, 2006.

[12] M. L. Sanchez, M. I. Ozores, R. Colle et al., "Soil C[O.sub.2] fluxes in cereal land use of the Spanish plateau: influence of conventional and reduced tillage practices," Chemosphere, vol. 47, no. 8, pp. 837-844, 2002.

[13] C. Bayer, J. Mielniczuk, T. J. C. Amado, L. Martin-Neto, and S. V. Fernandes, "Organic matter storage in a sandy clay loam acrisol affected by tillage and cropping systems in southern Brazil," Soil and Tillage Research, vol. 54, no. 1-2, pp. 101-109, 2000.

[14] R. Lal and T. Logan, "Agricultural activities and greenhouse gas emissions from soils of the tropics," in Soil Management and Greenhouse Effect, pp. 293-307, 1995.

[15] J. Six, S. M. Ogle, F. J. Breidt, R. T. Conant, A. R. Mosiers, and K. Paustian, "The potential to mitigate global warming with no-tillage management is only realized when practised in the long term," Global Change Biology, vol. 10, no. 2, pp. 155-160, 2004.

[16] S. Fontaine, G. Bardoux, L. Abbadie, and A. Mariotti, "Carbon input to soil may decrease soil carbon content," Ecology Letters, vol. 7, no. 4, pp. 314-320, 2004.

[17] P. A. Matson and R. C. Harriss, Biogenic Trace Gases: Measuring Emissions from Soil and Water, Blackwell Scientific Publications, Oxford, UK, 2009.

[18] L. S. Jensen, T. Mueller, K. R. Tate, D. J. Ross, J. Magid, and N. E. Nielsen, "Soil surface C[O.sub.2] flux as an index of soil respiration in situ: a comparison of two chamber methods," Soil Biology and Biochemistry, vol. 28, no. 10-11, pp. 1297-1306, 1996.

[19] T. J. Griffis, X. Lee, J. M. Baker, S. D. Sargent, and J. Y. King, "Feasibility of quantifying ecosystem-atmosphere [C.sup.18][O.sup.16]O exchange using laser spectroscopy and the flux-gradient method," Agricultural and Forest Meteorology, vol. 135, no. 1-4, pp. 44-60, 2005.

[20] J. M. Norman, C. J. Kucharik, S. T. Gower et al., "A comparison of six methods for measuring soil-surface carbon dioxide fluxes," Journal of Geophysical Research D: Atmospheres, vol. 102, no. 24, pp. 28771-28777, 1997

[21] P. Rochette, E. Gregorich, and R. Desjardins, "Comparison of static and dynamic closed chambers for measurement of soil respiration under field conditions," Canadian Journal of Soil Science, vol. 72, pp. 605-609, 1992.

[22] A. R. Pedersen, S. O. Petersen, and K. Schelde, "A comprehensive approach to soil-atmosphere trace-gas flux estimation with static chambers," European Journal of Soil Science, vol. 61, no. 6, pp. 888-902, 2010.

[23] R. S. Jassal, T. A. Black, Z. Nesic, and D. Gaumont-Guay, "Using automated non-steady-state chamber systems for making continuous long-term measurements of soil C[O.sub.2] efflux in forest ecosystems," Agricultural and Forest Meteorology, vol. 161, pp. 57-65, 2012.

[24] G. Livingston and G. Hutchinson, "Enclosure-based measurement of trace gas exchange: applications and sources of error," in Biogenic Trace Gases: Measuring Emissions from Soil and Water, pp. 14-51, 1995.

[25] A. Mosier and A. Bouwman, "Gas flux measurement techniques with special reference to techniques suitable for measurements over large ecologically uniform areas," in Soils and the Greenhouse Effect, pp. 289-301, John Wiley & Sons, New York, NY, USA, 1990.

[26] J. Schneider, L. Kutzbach, S. Schulz, and M. Wilmldng, "Overestimation of C[O.sub.2] respiration fluxes by the closed chamber method in low-turbulence nighttime conditions," Journal of Geophysical Research G: Biogeosciences, vol. 114, no. 3, Article ID G03005, 2009.

[27] L. Kutzbach, J. Schneider, T. Sachs et al., "C[O.sub.2] flux determination by closed-chamber methods can be seriously biased by inappropriate application of linear regression," Biogeosciences, vol. 4, no. 6, pp. 1005-1025, 2007.

[28] W. J. Wang, Y. G. Zu, H. M. Wang et al., "Effect of collar insertion on soil respiration in a larch forest measured with a LI-6400 soil C[O.sub.2] flux system," Journal of Forest Research, vol. 10, no. 1, pp. 57-60, 2005.

[29] M. Rayment and P. Jarvis, "An improved open chamber system for measuring soil C[O.sub.2] effluxes in the field," Journal of Geophysical Research: Atmospheres, vol. 102, no. 24, pp. 28779-28784, 1997.

[30] E. A. Davidson, K. Savage, L. V. Verchot, and R. Navarro, "Minimizing artifacts and biases in chamber-based measurements of soil respiration," Agricultural and Forest Meteorology, vol. 113, no. 1-4, pp. 21-37, 2002.

[31] IUSS Working Group, "World reference base for soil resources," World Soil Resources Report 103, FAO, Rome, Italy, 2006.

[32] S. Brenna, A. Rocca, M. Sciaccaluga, M. Valagussa, and R. Scaccabarozzi, "Stock di carbonio organico e fertilita biologica," in Regione Lombardia. Il Ruolo Dellagricoltura Conservativa nel Bilancio del Carbonio, AgriCO2ltura. Quaderni della Ricerca, pp. 53-73, Giugno, Milano, Italy, 2013.

[33] V. Stolbovoy, L. Montanarella, N. Filippi, S. Selvaradjou, P. Panagos, and J. Gallego, "Soil sampling protocol to certify the changes of organic carbon stock in mineral soils of European Union," EUR 21576 EN, Office for Official Publications of the European Communities, Luxembourg, Germany, 2005.

[34] A. Walkley and I. A. Black, "An examination of the Degtjareff method for determining soil organic matter, and a proposed modification of the chromic acid titration method," Soil Science, vol. 37, pp. 29-38, 1934.

[35] E. D. Vance, P. C. Brookes, and D. S. Jenkinson, "An extraction method for measuring soil microbial biomass C," Soil Biology and Biochemistry, vol. 19, no. 6, pp. 703-707, 1987

[36] K. R. Tate, D. J. Ross, and C. W. Feltham, "A direct extraction metho d to estimate soil microbial c: effects of experimental variables and some different calibration procedures," Soil Biology and Biochemistry, vol. 20, no. 3, pp. 329-335, 1988.

[37] G. Gerosa, A. Finco, E. Zaglio, R. Marzuoli, M. Chiazzese, and S. Amaducci, "A novel static automatic sampler for soil flux measurements of carbon dioxide and nitrogen protoxide," in La Chimica Agraria Tra Energia e Ambiente, Proceedings of the 30th Congress of the Italian Society of Agricultural Chemistry, 2012.

[38] G. Gerosa, S. Cieslik, and A. Ballarin-Denti, "Micrometeorological determination of time-integrated stomatal ozone fluxes over wheat: a case study in Northern Italy," Atmospheric Environment, vol. 37, no. 6, pp. 777-788, 2003.

[39] P. Pengthamkeerati, P. P. Motavalli, R. J. Kremer, and S. H. Anderson, "Soil carbon dioxide efflux from a claypan soil affected by surface compaction and applications of poultry litter," Agriculture, Ecosystems and Environment, vol. 109, no. 1-2, pp. 75-86, 2005.

[40] J. Iqbal, R. Hu, S. Lin et al., "C[O.sub.2] emission in a subtropical red paddy soil (Ultisol) as affected by straw and N-fertilizer applications: a case study in Southern China," Agriculture, Ecosystems and Environment, vol. 131, no. 3-4, pp. 292-302, 2009.

[41] L. J. Li, X. Z. Han, M. Y. You, Y. R. Yuan, X. L. Ding, and Y. F. Qiao, "Carbon and nitrogen mineralization patterns of two contrasting crop residues in a Mollisol: Effects of residue type and placement in soils," European Journal of Soil Biology, vol. 54, pp. 1-6, 2013.

[42] K. Oorts, R. Merckx, E. Grehan, J. Labreuche, and B. Nicolardot, "Determinants of annual fluxes of C[O.sub.2] and [N.sub.2]O in long-term no-tillage and conventional tillage systems in northern France," Soil and Tillage Research, vol. 95, no. 1-2, pp. 133-148, 2007

[43] T. N. Maraseni and G. Cockfield, "Does the adoption of zero tillage reduce greenhouse gas emissions? An assessment for the grains industry in Australia," Agricultural Systems, vol. 104, no. 6, pp. 451-458, 2011.

[44] S. De Neve and G. Hofman, "Influence of soil compaction on carbon and nitrogen mineralization of soil organic matter and crop residues," Biology and Fertility of Soils, vol. 30, no. 5-6, pp. 544-549, 2000.

[45] X.-Z. Han, H.-B. Li, and W. R. Horwath, "Temporal variations in soil C[O.sub.2] efflux under different land use types in the black soil zone of northeast china," Pedosphere, vol. 23, no. 5, pp. 636-650, 2013.

[46] A. Sowerby, H. Blum, T. R. Gray, and A. S. Ball, "The decomposition of Lolium perenne in soils exposed to elevated C[O.sub.2]: comparisons of mass loss of litter with soil respiration and soil microbial biomass," Soil Biology and Biochemistry, vol. 32, no. 10, pp. 1359-1366, 2000.

[47] P.-L. Shi, X.-Z. Zhang, Z.-M. Zhong, and H. Ouyang, "Diurnal and seasonal variability of soil C[O.sub.2] efflux in a cropland ecosystem on the Tibetan Plateau," Agricultural and Forest Meteorology, vol. 137, no. 3-4, pp. 220-233, 2006.

[48] Y. Liu, K.-Y. Wan, Y. Tao et al., "Carbon dioxide flux from rice paddy soils in central China: effects of intermittent flooding and draining cycles," PLoS ONE, vol. 8, no. 2, Article ID e56562, 2013.

[49] M. R. Raupach and A. S. Thom, "Turbulence in and above plant canopies," Annual Review of Fluid Mechanics, vol. 13, pp. 97-129, 1981.

[50] K. T. Paw U, Y. Brunet, S. Collineau et al., "On coherent structures in turbulence above and within agricultural plant canopies," Agricultural and Forest Meteorology, vol. 61, no. 1-2, pp. 55-68, 1992.

[51] J. Finnigan, "Turbulence in plant canopies," Annual Review of Fluid Mechanics, vol. 32, pp. 519-571, 2000.

[52] E. S. Takle, W. J. Massman, J. R. Brandle et al., "Influence of high-frequency ambient pressure pumping on carbon dioxide efflux from soil," Agricultural and Forest Meteorology, vol. 124, no. 3-4, pp. 193-206, 2004.

[53] T. O. West and G. Marland, "A synthesis of carbon sequestration, carbon emissions, and net carbon flux in agriculture: comparing tillage practices in the United States," Agriculture, Ecosystems and Environment, vol. 91, no. 1-3, pp. 217-232, 2002.

Gerosa Giacomo, (1) Finco Angelo, (1) Boschetti Fabio, (1) Brenna Stefano, (2) and Marzuoli Riccardo (1)

(1) Department of Mathematics & Physics, Catholic University, via dei Musei 41, 25121 Brescia, Italy

(2) ERSAF, Regional Agency for Services to the Agriculture and Forests of the Lombardy Region, via Pola 12, 20124 Milan, Italy Correspondence should be addressed to Marzuoli Riccardo; riccardo.marzuoli@unicatt.it

Received 13 May 2014; Revised 29 August 2014; Accepted 1 September 2014; Published 28 October 2014

Academic Editor: Ana Iglesias

TABLE 1: Crop management practices in the two sites (CT =
conventional tillage site, NT = no-tillage site). All
dates are referred to year 2012.

                   CT site                    NT site

Plowing            (i) 21 May                 No-tillage
                   30 cm depth tillage,
                     mouldboard plough
                   (ii) 22 May
                   Harrow clods reduction
Species and        Zea mays cv Pioner P0222   Zea mays cv Decalb 6903
  cultivar
Sowing date        23 May                     27 April
                   Seedbed prepared using     Sod-seed
                     harrows/cultivators
Sowing technique   (1-2 passes)
Seed population    78 seed [m.sup.-2]         8.6 seed [m.sup.-2]
Interrow spacing   70 cm                      70 cm
Fertilizations     (i) 15 May                 (i) 27 April
                   Cattle manure: 100         Ammoniumsulphate N/S
                     [m.sup.3] [ha.sup.-1]      21: 24, 150 kg
                   C:N = 9.7, N: 200 kg         [ha.sup.-1]
                     [ha.sup.-1]              Diammoniumphosphate
                   (ii) 23 May                  N/P 18: 46, 60
                   Perphosphate (46% of P),     kg [ha.sup.-1]
                     30 Kg [ha.sup.-1]        (ii) 19 May 2012
                   (iii) 10 June              Urea (46% of N), 175
                   Urea (46% of N), 150 Kg      kg [ha.sup.-1]
                     [ha.sup.-1]
Irrigations        (i) 9 July 2012            (i) 5 July 2012
                   Flooding (2000 [m.sup.3]   Flooding (2000
                     [ha.sup.-1] per            [m.sup.3] [ha.sup.-1]
                     irrigation)                per irrigation)
                   (ii) 6 August 2012         (ii) 27 July 2012
                   Flooding (2000 [m.sup.3]   Flooding (2000
                     [ha.sup.-1] per            [m.sup.3] [ha.sup.-1]
                     irrigation)                per irrigation)
Open field         (i) 21 June 2012           (i) 27 April 2012
  weeding          Du Pont "Codacide" 1.0     Syngenta "Force" 12 kg
                     L [ha.sup.-1]              [ha.sup.-1]
                   Syngenta "Callisto"        (ii) 28 April 2012
                     0.5 L [ha.sup.-1]        Bayer "Merlin Gold"
                   Du Pont "Titus ultra"        1.2 L[ha.sup.-1]
                     50 g [ha.sup.-1]         Dow "Hopper"
                   (ii) 2 October 2012          (glyphosate) 3
                   Dow "Hopper"                 L[ha.sup.-1]
                     (glyphosate), 3.5        (iii) 26 May 2012
                     L [ha.sup.-1]            Du Pont "Codacide"
                                                1.3 L [ha.sup.-1]
                                              Syngenta "Callisto"
                                                0.5 L [ha.sup.-1]
                                              Du Pont "Titus ultra"
                                                42 g[ha.sup.-1]
                                              (iv) 8 October 2012
                                              Dow "Hopper"
                                                (glyphosate), 3
                                                L [ha.sup.-1]
Weedings within    (i) 16 July                (i) 16 July
  SASSFLUX         Dow "Hopper"               Dow "Hopper"
                     (glyphosate) 3             (glyphosate) 3
                     L [ha.sup.-1]              L [ha.sup.-1]
                   (ii) 28 September          (ii) 28 September
                   Dow "Hopper"               Dow "Hopper"
                     (glyphosate) 3             (glyphosate) 3
                     L [ha.sup.-1]              L [ha.sup.-1]
                   (iii) 24 October           (iii) 24 October 2012
                   Dow "Hopper"               Dow "Hopper"
                     (glyphosate) 3             (glyphosate) 3
                     L [ha.sup.-1]              L [ha.sup.-1]
Harvest date       2 September                22 August
Harvest technique  One combined harvester     One combined harvester
                     (thresher and chopper      (thresher and
                     machine for silage)        chopper machine)
                   Two tractors with wagon    Two tractors with
                     for fine                   wagon for fine

TABLE 2: Soil and litter organic carbon content in the two
experimental sites. In both sites, soil type was Luvisol,
sandy loam, subacid, and well drained. CT = conventional
tillage site; NT = no-tillage site.

                                             CT site

Soil type                             Luvisol, sandy loam,
                                   sub acid, and well drained

Bulk density of the soil        1.57 [+ or -] 0.10 g [cm.sup.-3]
Organic carbon of the            10.3 [+ or -] 1.8 g [kg.sup.-1]
  topsoil (0-30 cm)
Carbon stock in the topsoil      48.7 [+ or -] 6.3 t [ha.sup.-1]
Organic carbon in the                  55.8 [+ or -] 16.2
  microbial biomass                   [micro]g C [g.sup.-1]
Litter mass                                 No litter
Organic carbon in the litter                No litter
Carbon stock in the litter                  No litter

                                             NT site

Soil type                             Luvisol, sandy loam,
                                       sub acid, and well
                                             drained
Bulk density of the soil        1.57 [+ or -] 0.07 g [cm.sup.-3]
Organic carbon of the            9.5 [+ or -] 0.2 gC [kg.sup.-1]
  topsoil (0-30 cm)
Carbon stock in the topsoil      44.8 [+ or -] 5.6 t [ha.sup.-1]
Organic carbon in the                  97.4 [+ or -] 41.5
  microbial biomass                   [micro]gC [g.sup.-1]
Litter mass                      1.70 [+ or -] 0.4 kg [m.sup.-2]
Organic carbon in the litter    237.9 [+ or -] 90 gC [kg.sup.-1]
Carbon stock in the litter       4.2 [+ or -] 2.0 t [ha.sup.-1]

TABLE 3: Weather conditions measured at the NT site during
the data measuring campaign.

Month       Mean air temp.      Average of          Average of
             ([degrees]C)    maximum air temp.   minimum air temp.
                               ([degrees]C)        ([degrees]C)

July            22.73              29.20               16.59
August          24.73              31.51               17.82
September       18.57              24.15               14.30
October         12.45              17.12               9.73

Month       Rain   Average soil    Average soil
            (mm)   water content    heat flux *
                       (v/v)       (W/[m.sup.2])

July        12.6       0.26            -0.69
August      24.6       0.25            0.11
September   74.2       0.28            -6.75
October     87.6       0.30            -8.59

* Negative values indicate losses of thermal energy toward
the atmosphere.
COPYRIGHT 2014 Hindawi Limited
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2014 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Gerosa, Giacomo; Finco, Angelo; Boschetti, Fabio; Brenna, Stefano; Marzuoli, Riccardo
Publication:The Scientific World Journal
Article Type:Report
Date:Jan 1, 2014
Words:8581
Previous Article:Response of nitrogen and potassium fertigation to "Waris" almond (Prunus dulcis) under northwestern Himalayan region of India.
Next Article:Molecular cloning and bioinformatics analysis of a new plasma membrane [NA.sup.+]/[H.sup.+] antiporter gene from the halophyte Kosteletzkya virginica.
Topics:

Terms of use | Privacy policy | Copyright © 2021 Farlex, Inc. | Feedback | For webmasters |