Published in J. Environ. Qual. 32:1965-1977 (2003).
© ASA, CSSA, SSSA
677 S. Segoe Rd., Madison, WI 53711 USA
TECHNICAL REPORTS
Atmospheric Pollutants and Trace Gases
Spatial Variability of Nitrous Oxide Emissions and Their Soil-Related Determining Factors in an Agricultural Field
Junta Yanai*,a,
Takuji Sawamotob,
Taku Oea,
Kanako Kusab,
Keisuke Yamakawac,
Kazunori Sakamotoc,
Takahiko Naganawad,
Kazuyuki Inubushic,
Ryusuke Hatanob and
Takashi Kosakia
a Graduate School of Agriculture, Kyoto University, Kitashirakawa Oiwake, Sakyo, Kyoto 606-8502, Japan
b Graduate School of Agriculture, Hokkaido University, Sapporo, Hokkaido 060-8589, Japan
c Faculty of Horticulture, Chiba University, Matsudo, Chiba 271-8510, Japan
d Faculty of Life and Environmental Science, Shimane University, Matsue Shimane 690-8504, Japan
* Corresponding author (yanai{at}kais.kyoto-u.ac.jp).
Received for publication July 14, 2002.
 |
ABSTRACT
|
|---|
To evaluate spatial variability of nitrous oxide (N2O) emissions and to elucidate their determining factors on a field-scale basis, N2O fluxes and various soil properties were evaluated in a 100- x 100-m onion (Allium cepa L.) field. Nitrous oxide fluxes were determined by a closed chamber method from the one-hundred 10- x 10-m plots. Physical (e.g., bulk density and water content), chemical (e.g., total N and pH), and biological (e.g., microbial biomass C and N) properties were determined from surface soil samples (00.1 m) of each plot. Geostatistical analysis was performed to examine spatial variability of both N2O fluxes and soil properties. Multivariate analysis was also conducted to elucidate relationships between soil properties and observed fluxes. Nitrous oxide fluxes were highly variable (average 331 µg N m-2 h-1, CV 217%) and were lognormally distributed. Log-transformed N2O fluxes had moderate spatial dependence with a range of >75 m. High N2O fluxes were observed at sites with relatively low elevation. Multivariate analysis indicated that an organic matter factor and a pH factor of the principal component analysis were the main soil-related determining factors of log-transformed N2O fluxes. By combining multivariate analysis with geostatistics, a map of predicted N2O fluxes closely matched the spatial pattern of measured fluxes. The regression equation based on the soil properties explained 56% of the spatially structured variation of the log-transformed N2O fluxes. Site-specific management to regulate organic matter content and water status of a soil could be a promising means of reducing N2O emissions from agricultural fields.
Abbreviations: MBF, microbial biomass factor OMF, organic matter factor PHF, pH factor SOF, soluble salts factor SSF, soil structure factor WFPS, water-filled pore space
 |
INTRODUCTION
|
|---|
NITROUS OXIDE is a trace gas that is involved in both the greenhouse effect (Yung et al., 1976) and destruction of stratospheric ozone (Johnson, 1971). Total annual emissions are estimated to be 17.7 Tg N (6.736.6 Tg N), about 60% of which is derived from the soil (Ehhalt et al., 2001). Estimates of N2O emissions from agricultural soils vary widely, for example, 2.29 to 3.65 Tg N (Bouwman, 1990), 0.2 to 2.1 Tg N (Eichner, 1990), 0.03 to 3 Tg N (Watson et al., 1992), or 2.4 Tg N (Bouwman, 1996). Accurate evaluation of the N2O emissions from agricultural land on a field scale is urgently required to quantify the amount of total emissions more precisely and to establish practical methods of reducing N2O emissions.
Evaluation of N2O emissions in agricultural fields must take account of spatial variability because it is generally accepted that emissions from within even a small area may vary by orders of magnitude (Bouwman, 1996). This heterogeneity makes predictions highly uncertain because in most cases only weak or no spatial dependencies have been demonstrated with considerable variation in the field (Folorunso and Rolston, 1984, 1985; Ball et al., 1997; Clemens et al., 1999; Roever et al., 1999). It is not certain, however, whether this tendency for weak spatial dependency is universal. Elucidation of the determining factors of N2O emissions is also important because N2O emissions have been reported to be highly variable depending on soil conditions such as organic matter, water and oxygen status, pH, substrate concentration, carbon supply, and temperature as well as fertilization and land management (Sahrawat and Keeney, 1986). In particular, overall effects of soil properties on N2O emission rates in the field should be fully evaluated because in many studies only individual factors have been investigated. Furthermore, understanding of spatial relationships is also required to mitigate N2O emissions in the field through soil management.
The objectives of the present study were to (i) evaluate the spatial variability of N2O fluxes and soil properties in an agricultural field, (ii) determine the possible causeeffect relationship between soil properties and N2O fluxes by multivariate analysis, and (iii) interpret the relationships spatially by combining multivariate analysis with geostatistics.
 |
MATERIALS AND METHODS
|
|---|
Experimental Field
The experiment was performed in an agricultural field located in Mikasa City, Hokkaido prefecture, Japan (43°14' N, 141°50' E) with a mean annual temperature and precipitation of 7.2°C and 1204 mm, respectively (Fig. 1)
. The field examined had an area of 2.1 ha (147 x 146 m) with a slight gradient (0.2%) across the northwest (relatively low) and southeast (relatively high) corners. It has been used to cultivate onion for more than 20 yr after long-term paddy rice (Oryza sativa L.) cultivation. The soil is of alluvial origin (clayey aquic) and classified as a fine, mesic, Mollic Fluvaquent (Soil Survey Staff, 1998). The soil in the plow layer has a texture class of silty clay (12% sand, 51% silt, and 37% clay) and a cation exchange capacity (CEC) of 25.5 cmolc kg-1. Further information on this site can be obtained from Hayashi and Hatano (1999). Fertilization and crop management were conducted homogeneously across the field, and followed standard procedures for this area. On 27 Apr. 1999 chemical fertilizer was applied as ammonium sulfate at rates equivalent to 0.032 kg m-2 of N; onion seedlings were transplanted on May 1; onion tops were cut off by agricultural machinery on 23 August to dry out the aboveground material, and onions were harvested on 6 September.
Gas and Soil Sampling and Their Analyses
Gas and soil samples were collected on 17 Aug. 1999. The weather at the time of sampling was partly cloudy, and no rain had been observed since 6 August. The sampling was performed in August because the highest values of the N2O fluxes of the field were generally observed at this time of the year with relatively high precipitation and temperature according to previous studies (Hatano and Sawamoto, 1997; Sawamoto and Hatano, 2000; Kusa et al., 2002). The central part of the experimental field (100 x 100 m) was divided into one-hundred 10- x 10-m plots. In total, 100 gas samples were collected at the center of the plots, using a closed chamber technique, that is, by covering the soil surface for 15 min with chambers and collecting gas samples at time zero from the open air and after 15 min from the headspace. For the gas sampling, two types of plastic chambers were used with surface areas of 0.0707 and 0.0875 m2. Preliminary experiments suggested that the type of chamber had a negligible effect on the measurement of N2O fluxes (data not shown). Soil temperature at a depth of 0.10 m was measured at the time of gas sampling as a physical property of the soil. After collecting gas samples, soil cores 0.14 m in diameter and 0.10 m in depth were taken from the surface layer, again at the center of each plot, using stainless steel core samplers. Sampling and associated measurements, which were performed from north to south, took a total of 5 h to perform (11001600 h). Air temperature and relative humidity of the area were relatively constant during this period: 30.2, 30.2, 29.7, 30.5, 28.7, and 27.1°C and 71, 68, 71, 64, 76, and 81%, from 1100 till 1600 h at 1-h intervals, respectively.
The N2O concentration of each gas sample was measured using gas chromatography (GC-14B; Shimadzu, Kyoto, Japan) with an electron capture detector. Nitrous oxide flux rate from the soil surface was then calculated based on the differences in N2O concentration between ambient air and the headspace gas sample taken after 15 min.
Each of the soil core samples was subdivided into three subsamples to evaluate physical, chemical, and biological properties of the soil. The first set of subsamples (undisturbed) was used to measure physical properties: bulk density; solid-, liquid-, and gas-phase volumetric percentage; the percentage of water-filled pore space (WFPS) in total pore space; and water content. Three variables, liquid phase percentage, WFPS, and water content, were measured in relation to water status of the soil because water was expected to influence the level of N2O fluxes to a considerable degree. The second set of subsamples was air-dried, sieved to 2 mm, and analyzed for the following chemical properties: total N content; total C content; C to N ratio; inorganic C content; water-soluble NO-3; NH+4 and inorganic N (as NO-3 + NH+4); available P (Bray 2 method); exchangeable Ca, Mg, K, and Na; pH; and electrical conductivity (EC). Total N content and total C content were determined by the dry combustion method (NC-800 analyzer; Sumika Chemical Analysis Service, Osaka, Japan) with C to N ratio calculated. Inorganic C was determined after extracting the soil with 0.5 M K2SO4. The contents of water-soluble NO-3 and NH+4 were determined by the GriessIlosvay method (after reduction of NO-3 to NO-2) and indophenol method, respectively, after shaking 1 h at the soil to water ratio of 1:5 (Mulvaney, 1996). All other analytical procedures were as described by Yanai et al. (2000). Total and readily available fractions of N and C were measured because soil N was assumed to supply substrates, and soil C the energy, for the N2O fluxes. The third set of subsamples (only sieved, not dried) was used to measure biological properties: microbial biomass N and C with the fumigation extraction method (Brookes et al., 1985; Vance et al., 1987). In total, 23 properties (7 physical, 14 chemical, and 2 biological properties) were measured as the soil properties of the field.
Statistical Analysis
Normality and lognormality tests were conducted on the data distributions (chi-square test). In addition to descriptive statistics, to evaluate the determining factors of N2O fluxes, multivariate analysis was performed. Principal component analysis (PCA) of the soil properties was conducted to summarize data and investigate the relationships among the properties (Kosaki and Juo, 1989). Stepwise multiple regression analysis was subsequently performed, using the scores of the extracted principal components as independent variables and N2O flux as a dependent variable. The statistical software SYSTAT 8.0 (SPSS, 1998) was used in the analysis.
Geostatistical Analysis
Geostatistical analysis was performed not only on the original data (i.e., measured N2O fluxes and soil properties), but also on calculated data (i.e., extracted PCA factors and predicted N2O fluxes). In this analysis, a semivariogram was first used to evaluate the spatial variability of the properties (i.e., to describe the average variances of pairs of points at a given distance apart) (Oliver, 1987; Webster and Oliver, 2001). This mirrors the similarity of pairs separated by equal distance. Often, it is found that the semivariance increases with increasing distance between sampling points to a maximum (the sill) at a moderate distance (the range). Points closer together than the range are autocorrelated whereas points further apart are not related to one another. The variation below the scale of investigation and/or due to experimental errors, the nugget variance, is determined as the ordinate intercept. In the analysis, two indices of spatial dependence were employed. One is the Q value [calculated as (sill - nugget)/sill], which indicates the spatial structure at the sampling scale (Goerres et al., 1997), and the other is the range, which indicates the limit of spatial dependence. In the analysis, the semivariogram model with the smallest residual sum of squares was used for the estimation of the semivariogram parameters. Details of the fitting models are available from Webster and Oliver (2001). Maps were computed subsequently using block kriging to evaluate regional patterns of variation rather than local details. The geostatistical software, GS+ Version 5.3 for Windows (Gamma Design Software), was used in the analysis (Robertson, 1998).
 |
RESULTS AND DISCUSSION
|
|---|
Descriptive Statistics of Nitrous Oxide Fluxes and Soil Properties
The average value of the N2O fluxes was 331 µg N m-2 h-1 (Table 1), which was comparable with previously reported values at this site (i.e., 871 µg N m-2 h-1) obtained in the middle of August as the maximum value of the growing period (Hatano and Sawamoto, 1997; Sawamoto and Hatano, 2000; Kusa et al., 2002). This value was relatively high compared with reported values for arable soils, for example, 5.8 to 139 (Goodroad et al., 1984), 41.2 (Ambus and Christensen, 1995), 8.6 (Ball et al., 1997), 0.8 to 11.2 (Roever et al., 1999), and 53.6 µg N m-2 h-1 (Clemens et al., 1999), as reviewed by Mosier (1994) and Bouwman (1996). It was comparable with data after N fertilization (Bremner et al., 1981; Thornton et al., 1996) or after plowing (van der Weerden et al., 2000). This could be due to the high amounts of N fertilizer application (Sahrawat and Keeney, 1986) and relatively moist soil conditions reflecting clayey texture and topography of lowland areas (Sawamoto and Hatano, 2000). Nitrous oxide fluxes displayed high variability with the minimum, maximum, and coefficient of variation (CV) of 6.83 and 5980 µg N m-2 h-1, and 217%, respectively (Table 1). Nitrous oxide data were not normally distributed (i.e., fluxes from several plots appeared to be much higher than the rest). They were, therefore, log-transformed and the descriptive statistics calculated. The log-transformed fluxes had a normal distribution (P < 0.01), as has been reported in previous studies (Ambus and Christensen, 1994; Ball et al., 1997; Roever et al., 1999). This result may be because N2O emissions are related to microbial processes, which are often lognormally distributed (Parkin, 1993).
View this table:
[in this window]
[in a new window]
|
Table 1. Descriptive statistics of N2O flux, log-transformed N2O flux, and physical, chemical, and biological properties of the soil.
|
|
In contrast to N2O fluxes, most of the soil properties displayed a normal distribution. The average values represented a typical surface soil of arable land in this area (Table 1). Physical properties indicated that the soil had a bulk density of 1.04 Mg m-3 with about 60% porosity and WFPS of 47.8%. Chemical properties indicated that the N status of the soil was representative for an alluvial arable soil with a total N content of 2.04 g kg-1 and inorganic N (readily available; NO-3 + NH+4) content of 60.3 mg kg-1. The concentrations of NO-3 and NH+4 were sufficient to provide opportunities for both denitrification and nitrification to occur in this soil. The C status (an indicator of energy sources available to microorganisms) and contents of exchangeable cations and available P were also appreciable. Biological properties further showed that the level of microbial biomass N and C and their ratio to total N and C were relatively high, possibly because the biomass community still reflected the original paddy-field condition, for which such ratios are generally high (Shibahara and Inubushi, 1995). The CV values suggested that the soil properties measured had less variability than the N2O fluxes. In addition, there was a general tendency for the CV values of physical properties (about 10%) to be lower than those of the chemical and biological properties, and CV values of chemical properties representing total and available fractions (1020%) were lower than those of biological properties (3050%) and chemical properties representing readily available fractions (>60%). These tendencies would be generally observed in most arable lands (Roever and Kaiser, 1999) and reflect the overall intrinsic characteristics of each soil property, as well as soil type and land management factors.
Spatial Variability of Nitrous Oxide Fluxes and Soil Properties
The variability of N2O fluxes described in the previous section can be interpreted as the outcome of both spatial and temporal effects, as it took 5 h to perform gas sampling in the 100- x 100-m area. However, measurement of diurnal variation in N2O emission, which was performed at the northwest corner of the same field on another date (one month earlier), suggested that the temporal fluctuation of N2O fluxes during daytime was relatively small (less than two times) and even the diurnal difference was within several times (data not shown) as reported by other researchers (e.g., Blackmer et al., 1982; Thomson et al., 1997). Judging from the fact that the variability of N2O fluxes of this experiment was as much as two orders of magnitude, it was concluded that most of the variation, if not all, was derived from a spatial effect.
Based on this assumption, spatial analysis of N2O fluxes was performed. Figure 2
shows the semivariogram and isarithmic map of the log-transformed N2O fluxes. Table 2 shows the geostatistical parameters, which were estimated from the semivariogram. The semivariogam shows that the semivariance increases with increasing distance between sampling points. The degree of spatial dependence was moderate with the Q value of 0.50. The range was estimated to be more than 75 m, suggesting a considerably long range of spatial dependence. This result was in contrast to almost all previously published studies, which indicated no spatial dependence (Folorunso and Rolston, 1984; Ball et al., 1997; Roever et al., 1999) or weak spatial dependence with ranges of several meters only (Folorunso and Rolston, 1985; Clemens et al., 1999), with one exception reporting strong spatial dependence with a range of up to 48 m in a flat grassland field (Ambus and Christensen, 1994). This discrepancy may be due to the fact that the landscape of this field was relatively flat and thus N2O fluxes were more likely to reflect the gradual and continuous changes of other field conditions, which was rarely the case for ordinary undulating upland fields involved in previous studies. The existence of moderate spatial dependency suggests the possibility for control or regulation of N2O fluxes with site-specific soil and land management. The map also shows that relatively high N2O fluxes were observed at sites with a relatively lower elevation (i.e., in the northwest corner of the field).
View this table:
[in this window]
[in a new window]
|
Table 2. Geostatistical parameters of log-transformed N2O flux and physical, chemical, and biological properties of the soil.
|
|
Figures 3 through 6
show the semivariograms and isarithmic maps of the 23 soil physical, chemical, and biological properties. Their geostatistical parameters are shown in Table 2. The degree of spatial dependence of soil properties varied widely. The Q values were between 0.80 and 1.00 for bulk density, solid phase percentage, water content, soil temperature, total N, total C, and available P, suggesting a highly developed spatial structure. The Q values ranged from 0.50 to 0.80 for liquid phase percentage, gas phase percentage, WFPS, C to N ratio, water-soluble NH+4, exchangeable Ca, and pH, indicating moderate spatial structure. However, water-soluble NO-3; inorganic N; inorganic C; exchangeable Mg, K, and Na; electrical conductivity; and microbial biomass N and C displayed low spatial dependencies with low Q values. It may be concluded that physical properties and chemical properties representing the total or available fractions had strong or moderate spatial dependency, whereas chemical properties representing the water-soluble or readily available fractions and biological properties had weak or no spatial dependency. The ranges or the distance of spatial dependency also varied considerably. The ranges were more than 75 m for bulk density; solid-, liquid-, and gas-phase percentages; WFPS; water content; soil temperature; total C; total N; C to N ratio; and water-soluble NH+4, whereas no spatial range was observed for water-soluble NO-3; inorganic N; inorganic C; exchangeable Mg, K, and Na; electrical conductivity; and microbial biomass N and C. There was a general tendency for the properties with strong spatial dependency to have longer ranges, and vice versa. These results indicate the possibility of regulating N2O fluxes through site-specific management of those soil properties with strong or moderate spatial dependency. Isarithmic maps of soil properties further showed that spatial variability was variable among soil properties, as shown in Fig. 3 through 6, probably reflecting their intrinsic characteristics such as mobility and reactivity in soil.
Multivariate Analysis of the Relationship between Nitrous Oxide Fluxes and Soil Properties
Principal component analysis identified five principal components (PC15) with eigenvalues, or the variances of the components, greater than 1.0. Together these components accounted for 73.3% of the total variance (Table 3). The remaining less-significant components were considered as errors, which included the random components of soil variation and various types of errors in soil sampling and analysis. Based on the component loadings after varimax rotation, which was accomplished to increase the interpretability of the components, the first component showed high loadings for total N, total C, water content, and soil temperature. Since these properties were related to the organic matter status of the soil, this component was referred to as the organic matter factor (OMF). The second component showed high loadings for bulk density; solid-, liquid-, and gas-phase percentages; and WFPS. Since these properties determined soil physical properties and soil structure, this component was referred to as the soil structure factor (SSF). Similarly, the third component showed high loadings for water-soluble NO-3, inorganic N, electrical conductivity, and exchangeable K, and was referred to as the soluble salts factor (SOF). The fourth component, the pH factor (PHF), showed high loadings for the content of exchangeable Ca and pH; and the fifth component, the microbial biomass factor (MBF) showed high loadings for microbial biomass N and C. These factors were considered to correspond to the chemical properties representing total fractions, physical properties, chemical properties representing soluble fractions, chemical properties representing available fractions, and microbiological properties, respectively. Variation of the soil properties was thus summarized into a smaller number of factors, which were independent of each other.
View this table:
[in this window]
[in a new window]
|
Table 3. Component loadings, eigenvalues, and percentage of total variance explained for the first five principal components.
|
|
Figure 7
shows the semivariograms and isarithmic maps of the five PCs and Table 4 indicates their geostatistical parameters. The five PCs showed contrasting spatial dependence; strong or moderate spatial dependency was observed for OMF, SSF, and PHF whereas weak dependency was observed for SOF and MBF. These results reflected the spatial dependency of the original soil properties underlying each PC. Isarithmic maps of the five PCs also showed the similarity of their spatial variability with the original soil properties. This relationship is readily apparent if, for example, total N content, exchangeable Ca, water-soluble NO-3, liquid phase percentage, and microbial biomass C in Fig. 3 through 6 are compared with OMF, PHF, SOF, SSF, and MBF in Fig. 7, respectively.
View this table:
[in this window]
[in a new window]
|
Table 4. Geostatistical parameters of the first five principal components and the predicted and measured values of the log-transformed N2O fluxes.
|
|
Stepwise multiple regression analysis was subsequently performed to obtain a model for predicting N2O fluxes. In this analysis, log-transformed N2O flux was used as a dependent variable, and standardized scores of the five principal components described above were used as independent variables. Since no a priori information was available about the regression model, we assumed the presence of a linear combination of first-degree variables. The most appropriate model obtained (R2 = 0.28) with a significance level of
= 0.15 was:
 | [1] |
Spatial analysis of the residual errors showed that the errors were spatially independent, suggesting that the application of this regression model was statistically appropriate.
The fact that R2 equals 0.28 indicates that the model explains 28% of the total variance of the log-transformed measured N2O fluxes. It should be noted that N2O fluxes showed a slightly weak but significant relationship with soil properties on a field basis, in spite of the fact that N2O fluxes were highly variable (as much as two orders of magnitude), as is often the case for those variables derived from biological activities (Parkin, 1993). As the regression coefficients in the equation generally indicate the magnitude of each factor's contribution to the N2O fluxes, OMF can be identified as the most important factor in this relationship. High N2O fluxes were thus related to high amounts of organic matter, high water content, and low temperatures. The contribution of PHF was of secondary importance; a low pH encouraged higher N2O fluxes in this pH range. The components related to soil structure, soluble salts, and microbial biomass (SSF, SOF, and MBF) were not included in the model, although these factors were generally considered to affect the N2O fluxes. Judging from the fact that N2O fluxes were positively related to organic matter content and wetness, it was assumed that denitrification was the main mechanism of N2O emissions. This was consistent with the results of Christensen et al. (1990) where the spatial distribution of denitrification activity in the soil depended on the soil water content and on the amount and distribution of available organic matter. The fact that the ratio of N2O flux to NO flux was more than 100 also suggested active denitrification (Kusa et al., 2002) although nitrification may have occurred simultaneously based on the data of WFPS (Linn and Doran, 1984). From a spatial perspective, it was reasonable that PCs with relatively strong spatial dependency were extracted to explain the measured N2O fluxes with moderate spatial dependency. As a practical interpretation, it is concluded that OMF variables can be good indicators of the spatial variation of N2O fluxes, and site-specific management to regulate organic matter content as well as water status of soil may become an effective means of reducing N2O fluxes in the field.
Application of Geostatistics to the Interpretation of Determining Factors of Nitrous Oxide Fluxes with Multivariate Analysis
Equation [1] was used to predict N2O fluxes and to prepare a semivariogram and an isarithmic map (Fig. 8)
. The isarithmic map of predicted N2O fluxes was remarkably similar to that of measured N2O fluxes (Fig. 2), although the variation was somewhat underestimated in the predicted map. Nevertheless, the prediction of N2O fluxes from soil properties might be an effective way of estimating N2O fluxes for comparative purposes (e.g., identifying N2O emission "hot spots" in a field).
From the semivariogram, the Q value was calculated to be 1.00, suggesting that the multiple regression model in fact omitted all random components (i.e., all variation expressed in Eq. [1] was spatially structured). This is supported by the fact that the extracted PCs (i.e., OMF and PHF) had strong or moderate spatial dependency (Fig. 7). On the other hand, the semivariogram of the measured fluxes (Fig. 2) indicates that 50% of the total variation of measured fluxes was spatially structured (Table 2), with the remaining 50% originally random or unable to be controlled at this scale. Since Eq. [1] explains 28% of the total measured flux variance, it can be concluded that this equation is actually accounting for as much as 56% of the spatially structured or nonrandom variation of the N2O fluxes (28 out of 50%). This is very high compared with other published reports on this matter. In this way, geostatistical analysis in combination with multivariate analysis enables the relationship between dependent and independent variables to be interpreted in a spatial context. This may in turn open up the possibility of site-specific management of the dependent variables, as discussed for the relationship between soil properties and the yield of rice in a paddy field (Yanai et al., 2001).
 |
CONCLUSIONS
|
|---|
Spatial variability of the N2O fluxes and soil properties was evaluated in an agricultural field, and the important role of soil properties in determining N2O emissions was demonstrated by combining geostatistics with multivariate analysis. These results suggest that site-specific soil management such as the control of organic matter content or water status could be a promising strategy for reducing N2O emissions in the field. Based on the fact that more than half of all N2O emissions for this field were released in August and September (Hatano and Sawamoto, 1997; Sawamoto and Hatano, 2000; Kusa et al., 2002), site and time-specific management could be further established for the efficient reduction of N2O emissions. In this context, investigations of the spatiotemporal variation of the N2O fluxes would lead to a better understanding of N2O emissions. Further studies are to be conducted to investigate the effectiveness of site-specific management (such as the control of organic matter content, water status, or inorganic fertilizer application) on the mitigation of N2O fluxes from agricultural fields.
 |
ACKNOWLEDGMENTS
|
|---|
The authors are very grateful to Mathew Turner, Kyoto University, and three anonymous reviewers for their critical reading of the manuscript. This study was partly supported by a Grant-in-Aid for Scientific Research from the Ministry of Education, Science, Sports, and Culture of Japan (no. 10460030).
 |
REFERENCES
|
|---|
- Ambus, P., and S. Christensen. 1994. Measurement of N2O emission from a fertilized grassland: An analysis of spatial variability. J. Geophys. Res. [Atmos.] 99:1654916555.
- Ambus, P., and S. Christensen. 1995. Spatial and seasonal nitrous oxide and methane fluxes in Danish forest, grassland, and agroecosystem. J. Environ. Qual. 24:9931001.[Abstract/Free Full Text]
- Ball, B.C., G.W. Horgan, H. Clayton, and J.P. Parker. 1997. Spatial variability of nitrous oxide fluxes and controlling soil and topographic properties. J. Environ. Qual. 26:13991409.[Abstract/Free Full Text]
- Blackmer, A.M., S.G. Robbins, and J.M. Bremner. 1982. Diurnal variability of emission of nitrous oxide from soils. Soil Sci. Soc. Am. J. 46:937942.[Abstract/Free Full Text]
- Bouwman, A.F. 1990. Exchange of greenhouse gases between terrestrial ecosystems and the atmosphere. p. 61127. In A.F. Bouwman (ed.) Soils and the greenhouse effect. John Wiley & Sons, London.
- Bouwman, A.F. 1996. Direct emission of nitrous oxide from agricultural soils. Nutr. Cycling Agroecosyst. 46:5370.
- Bremner, J.M., G.A. Breitenbeck, and A.M. Blackmer. 1981. Effect of anhydrous ammonia fertilization on emission of nitrous oxide from soils. J. Environ. Qual. 10:7780.
- Brookes, P.C., A. Landman, G. Pruden, and D.S. Jenkinson. 1985. Chloroform fumigation and the release of soil nitrogen: A rapid direct extraction method to measure microbial biomass nitrogen in soil. Soil Biol. Biochem. 17:837842.
- Christensen, S., S. Simkins, and J.M. Tiedje. 1990. Spatial variation in denitrification: Dependency of activity centers on the soil environment. Soil Sci. Soc. Am. J. 54:16081613.[Abstract/Free Full Text]
- Clemens, J., M.P. Schillinger, H. Goldbach, and B. Huwe. 1999. Spatial variability of N2O emissions and soil parameters of an arable silt loamA field study. Biol. Fertil. Soils 28:403406.
- Ehhalt, D., M. Prather, F. Dentener, R. Derwent, E. Dlugokencky, E. Holland, I. Iaksen, J. Katima, V. Kirchhoff, P. Matson, P. Midgley, and M. Wang. 2001. Atmospheric chemistry and greenhouse gases. p. 239288. In J.T. Houghton et al. (ed.) Climate change 2001: The scientific basis. Cambridge Univ. Press, Cambridge.
- Eichner, M.J. 1990. Nitrous oxide emission from fertilized soils: Summary of available data. J. Environ. Qual. 19:272280.[Abstract/Free Full Text]
- Folorunso, O.A., and D.E. Rolston. 1984. Spatial variability of field-measured denitrification gas fluxes. Soil Sci. Soc. Am. J. 48:12141219.[Abstract/Free Full Text]
- Folorunso, O.A., and D.E. Rolston. 1985. Spatial and spectral relationships between field-measured denitrification gas fluxes and soil properties. Soil Sci. Soc. Am. J. 49:10871093.[Abstract/Free Full Text]
- Goerres, J.H., M.J. Dichiaro, J.B. Lyons, and J.A. Amador. 1997. Spatial and temporal patterns of soil biological activity in a forest and an old field. Soil Biol. Biochem. 30:219230.
- Goodroad, L.L., D.R. Keeney, and L.A. Peterson. 1984. Nitrous oxide emissions from agricultural soils in Wisconsin. J. Environ. Qual. 13:557561.[Abstract/Free Full Text]
- Hatano, R., and T. Sawamoto. 1997. Emission of N2O from a clayey aquic soil cultivated with onion plants. p. 555556. In T. Ando (ed.) Plant nutrition for sustainable food production and environment. Kluwer Academic Publ., Dordrecht, the Netherlands.
- Hayashi, Y., and R. Hatano. 1999. Annual nitrogen leaching to subsurface drainage water from a clayey aquic soil cultivated with onions in Hokkaido, Japan. Soil Sci. Plant Nutr. (Tokyo) 45:451459.
- Johnson, H.S. 1971. Reduction of stratospheric ozone by nitrogen oxide catalysts from SST exhaust. Science (Washington, DC) 173:517522.[Abstract/Free Full Text]
- Kosaki, T., and A.S.R. Juo. 1989. Multivariate approach to grouping soils in small fields. I. Extraction of factors causing soil variation by principal component analysis. Soil Sci. Plant Nutr. (Tokyo) 35:469477.
- Kusa, K., T. Sawamoto, and R. Hatano. 2002. Nitrous oxide emissions for six years from a gray lowland soil cultivated with onions in Hokkaido, Japan. Nutr. Cycling Agroecosyst. (in press).
- Linn, D.M., and J.W. Doran. 1984. Effect of water-filled pore space on carbon dioxide and nitrous oxide production in tilled and nontilled soils. Soil Sci. Soc. Am. J. 48:12671272.[Abstract/Free Full Text]
- Mosier, A.R. 1994. Nitrous oxide emissions from agricultural soils. Fert. Res. 37:191200.
- Mulvaney, R.L. 1996. NitrogenInorganic forms. p. 11231184. In D.L. Sparks (ed.) Methods of soil analysis. Part 3. SSSA Book Ser. 5. SSSA, Madison, WI.
- Oliver, M.A. 1987. Geostatistics and its application to soil science. Soil Use Manage. 3:820.
- Parkin, T.B. 1993. Spatial variability of microbial process in soilA review. J. Environ. Qual. 22:409417.[Abstract/Free Full Text]
- Robertson, G.P. 1998. GS+: Geostatistics for the environmental sciences. Gamma Design Software, Plainwell, MI.
- Roever, M., O. Heinemeyer, J.C. Munch, and E.A. Kaiser. 1999. Spatial heterogeneity within the plough layer: High variability of N2O emission rates. Soil Biol. Biochem. 31:167173.
- Roever, M., and E.A. Kaiser. 1999. Spatial heterogeneity within the plough layer: Low and moderate variability of soil properties. Soil Biol. Biochem. 31:175187.
- Sahrawat, K.L., and D.R. Keeney. 1986. Nitrous oxide emission from soils. p. 103148. In B.A. Stewart (ed.) Advances in soil science. Vol. 4. Springer-Verlag, New York.
- Sawamoto, T., and R. Hatano. 2000. N2O flux from a well-structured gray lowland soil cultivated for onion in Hokkaido, Japan (No. 1): Amount of N2O flux in three years and a comparison with published studies. (In Japanese with English summary.) Nippon Dojo Hiryogaku Zasshi 71:659665.
- Shibahara, F., and K. Inubushi. 1995. Measurements of microbial biomass C and N in paddy soils by the fumigationextraction method. Soil Sci. Plant Nutr. (Tokyo) 41:681689.
- Soil Survey Staff. 1998. Keys to soil taxonomy. 8th ed. USDA Natural Resour. Conserv. Serv., Washington, DC.
- SPSS. 1998. SYSTAT 8.0. Statistics. SPSS, Chicago, IL.
- Thomson, P.E., J.P. Parker, J.R.M. Arah, H. Clayton, and K.A. Smith. 1997. Automated soil monolith-flux chamber system for the study of trace gas fluxes. Soil Sci. Soc. Am. J. 61:13231330.[Abstract/Free Full Text]
- Thornton, F.C., B.R. Bock, and D.D. Tyler. 1996. Soil emission of nitric oxide and nitrous oxide from injected anhydrous ammonium and urea. J. Environ. Qual. 25:13781384.[Abstract/Free Full Text]
- Vance, E.D., P.C. Brookes, and D.S. Jenkinson. 1987. An extraction method for measuring soil microbial biomass C. Soil Biol. Biochem. 19:703707.
- Van der Weerden, T.J., R.R. Sherlock, P.H. Williams, and K.C. Cameron. 2000. Effect of the three contrasting onion (Allium cepa L.) production systems on nitrous oxide emissions from soil. Biol. Fertil. Soils 31:334342.
- Watson, R.T., L.G. Meira Filho, E. Sanhueza, and A. Janetos. 1992. p. 2546. In J.T. Houghton (ed.) Climate change 1992. Cambridge Univ. Press, Cambridge.
- Webster, R., and M.A. Oliver. 2001. Geostatistics for environmental scientists. John Wiley & Sons, Chichester, UK.
- Yanai, J., C.K. Lee, T. Kaho, M. Iida, T. Matsui, M. Umeda, and T. Kosaki. 2001. Geostatistical analysis of soil chemical properties and rice yield in a paddy field and its application to the analysis of yield determining factors. Soil Sci. Plant Nutr. (Tokyo) 47:291301.
- Yanai, J., C.K. Lee, M. Umeda, and T. Kosaki. 2000. Spatial variability of soil chemical properties in a paddy field. Soil Sci. Plant Nutr. (Tokyo) 46:473482.
- Yung, Y.L., W.C. Wang, and A.A. Lacis. 1976. Greenhouse effects due to nitrous oxide. Geophys. Res. Lett. 3:619621.
Related articles in JEQ:
- This Issue in Journal of Environmental Quality
JEQ 2003 32: 1931-1938.
[Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
P. Lu, Y. Su, Z. Niu, and J. Wu
Geostatistical Analysis and Risk Assessment on Soil Total Nitrogen and Total Soil Phosphorus in the Dongting Lake Plain Area, China
J. Environ. Qual.,
May 25, 2007;
36(4):
935 - 942.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
G. Drewitt and J. S. Warland
Continuous Measurements of Belowground Nitrous Oxide Concentrations
Soil Sci. Soc. Am. J.,
January 1, 2007;
71(1):
1 - 7.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
J. Lee, J. Six, A. P. King, C. van Kessel, and D. E. Rolston
Tillage and Field Scale Controls on Greenhouse Gas Emissions
J. Environ. Qual.,
April 3, 2006;
35(3):
714 - 725.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
D. Ginting and B. Eghball
Nitrous Oxide Emission from No-Till Irrigated Corn: Temporal Fluctuation and Wheel Traffic Effects
Soil Sci. Soc. Am. J.,
May 6, 2005;
69(3):
915 - 924.
[Abstract]
[Full Text]
[PDF]
|
 |
|