Published online 27 April 2005
Published in Agron J 97:772-782 (2005)
DOI: 10.2134/agronj2004.0287
© 2005 American Society of Agronomy
677 S. Segoe Rd., Madison, WI 53711 USA
Spatial Variability
Spatial and Temporal Variation of Soil Nitrogen Parameters Related to Soil Texture and Corn Yield
H. Shahandeha,*,
A. L. Wrighta,
F. M. Honsa and
R. J. Lascanob
a Dep. of Soil and Crop Sci., Texas A&M Univ., College Station, TX 77843-2474
b Texas A&M Univ. Agric. Res. and Ext. Cent., Rt. 3, Box 219, Lubbock, TX 79403-9803
* Corresponding author (h-shahandeh{at}tamu.edu)
Received for publication November 23, 2004.
 |
ABSTRACT
|
|---|
The spatial variability of soil properties that affect the soil N budget and corn (Zea mays L.) grain yield were studied for 2 yr in south-central Texas to better assess the potential for variable-rate N fertilization. Residual soil NO3N with depth and soil N mineralization (Nmin) potential and their relationships with soil total N, soil organic C, and clay content were characterized. Residual soil NO3N to 60-cm depth was more related to corn yield than NO3N at shallower depths. Residual soil NO3N showed temporal variation with spatial structure existing for NO3N in the first year when NO3N concentrations were high, but this variation was absent in the second year when NO3N concentrations were low. The opposite trend was observed for soil Nmin. Soil total N exhibited temporal persistence. Soil Nmin showed significant correlations with soil clay content in both years. Temporal persistence of soil texture and variation in spatial structure for N parameters would likely result in different strategies for soil N management zones in 2002 vs. 2003. This experiment demonstrated the potential importance of soil texture for modifying fertilizer N recommendations. Texture generally is more easily determined than soil N parameters at the scale and intensity necessary for site-specific N management.
Abbreviations: CV, coefficient of variation Nmin, nitrogen mineralization SOC, soil organic carbon SOM, soil organic matter
 |
INTRODUCTION
|
|---|
NITROGEN FERTILIZATION has historically been an important crop input from production, economic, and environmental standpoints but currently is receiving more scrutiny because of the increased cost of N fertilizer, concerns about environmental degradation, and demands for mandatory nutrient planning (Dinnes et al., 2002; Raun et al., 2002; Schmidt et al., 2002). The current recommended method for determining appropriate amounts of N fertilization for crop production involves soil testing for residual surface soil (0 to 15 cm) nitrate (NO3) concentration and integrating with anticipated crop yield goal (Oberle and Keeney, 1990; Vanotti and Bundy, 1994; Stafford et al., 1996). The prescribed quantity of N is then normally applied uniformly across a field. However, the spatial variation in NO3N accumulated below 15-cm depth and soil N mineralized from soil organic matter (SOM) during the growing season must be assessed to determine their contribution to crop yield variability and to evaluate their necessity for variable-rate N fertilization (Mahmoudjafari et al., 1997; Schmidt et al., 2002; Eghball et al., 2003; Katsvario et al., 2003).
Variable-rate fertilization that considers differences in residual nutrient availability should be more economically and environmentally sustainable than uniform-rate application (Machado et al., 2000; Dinnes et al., 2002; Lopez-Granados et al., 2002; Crowin et al., 2003). To implement variable-rate fertilization or site-specific nutrient management practices, databases that provide spatial distributions of soil properties are needed (Cambardella and Karlen, 1999). However, the spatial variability of residual soil N is highly variable and represents potential problems in estimating the quantity of N available in the field (Parkin et al., 1988; Cabrera and Kissel, 1988).
The value of residual soil NO3 as a guide for N fertilization in precision farming using geostatistical methods has been evaluated in long-term fallow and for corn, winter wheat (Triticum aestivum L.), cornsoybean [Glycine max (L.) Merr.] rotation, and wheatsunflower (Helianthus annus L.) rotation (Goovaerts and Chiang, 1993; Cambardella et al., 1994; Mahmoudjafari et al., 1997; Cambardella and Karlen, 1999; Lopez-Granados et al., 2002). Geostatistical techniques have been used to explore the structure of spatial variation in agricultural soils by mapping areas of environmental concern and developing site-specific fertilizer application maps (Cambardella and Karlen, 1999; Lopez-Granados et al., 2002). Lopez-Granados et al. (2002), for example, found that subsoil NO3N concentrations had a strong spatial dependence. Goovaerts and Chiang (1993) determined that the spatial structure of mineralizable soil N was explained by the amount of soil oxidizable C. Cambardella et al. (1994) found that soil total N was spatially dependent, whereas NO3N and mineralized N were not. However, Bruckler et al. (1997) found that soil NO3N was both spatially and temporally dependent, and Robertson et al. (1988) found that soil Nmin had a high degree of spatial dependence.
Variable-rate fertilizer N recommendations, however, likely should not be based only on a single criterion derived from statistical models. For example, Schmidt et al. (2002) argued that variable-rate N recommendations based only on SOM were too simplistic to reflect variability in N availability within a field. Similarly, Mahmoudjafari et al. (1997) indicated that the mineralization potential of SOM was difficult to estimate because Nmin is a function of both soil temperature and water content, which can spatially vary independently of SOM content. Furthermore, Eghball et al. (2003) showed that spatial variability of corn grain yield was much lower than for soil NO3N, indicating the possible ineffectiveness of using spatial distribution of soil NO3N as the sole criterion for variable-rate N application.
There are other approaches that do not use a single criterion for N recommendations to delineate soil management zones. These approaches are based on the spatial structure of N parameters and related factors, including soil characteristics and cultural practices (Oberle and Keeney, 1990; Sadler et al., 1998; Bakhsh et al., 2000a, 2000b; Frogbrook et al., 2002). For example, soil properties such as clay content with only limited sampling of NO3 might be used to estimate and infer NO3 values. Oberle and Keeney (1990) and Sadler et al. (1998) in their assessments of N fertilizer requirement indicated that soil series and soil mapping information provided a basis for N management evaluation. Tabor et al. (1985) found that leaf petiole NO3N and soil clay content had high correlation and similar spatial structure, indicating that spatial structure and optimum sampling for petiole NO3N could be inferred from the spatial structure of clay for a particular soil mapping unit. Machado et al. (2000) also suggested that variable-rate fertilizer N application should be based on texture and soil NO3. Stafford et al. (1996) found a high correlation between soil series and crop yield, and Cox et al. (2003) found that among the soil factors affecting crop yield variability, clay content could be used as a basis for site-specific management.
The objectives of this study were to determine the temporal and spatial dependence of soil N parameters, such as total N, Nmin potential, and residual soil NO3N with depth, and their relationships with soil texture and corn yield. Strong correlations of soil N parameters with more easily measured soil properties, such as texture, may provide an opportunity to use these parameters to estimate and infer N values for site-specific N management.
 |
MATERIALS AND METHODS
|
|---|
Experimental Site
Research was conducted at the Texas A&M University Agricultural Research Farm in Burleson County near College Station, TX (30°32' N, 96°25' W), in 2002 and 2003. The alluvial soil used for the study is an intergrade of Weswood silt loam (fine-silty, mixed superactive, thermic Udifluventic Haplustepts) and Ships clay (very fine, mixed, active, thermic Chromic Hapluderts) with pH of 7.9 to 8.1. The experimental field was
2.5 ha, previously cropped to cotton (Gossypium hirsutum L.), and was planted with corn in this experiment. The cotton received 125 kg N ha1 as 3200 applied preplant. Sixty-three grids (plots) were superimposed on the field, with each grid being eight 1-m rows wide and 12.2 m long. The experimental layout was nine grids long extending consecutively from west to east and seven grids wide extending from north to south. Eight unsampled rows were situated between sampling grids and ran the entire length of the field.
Soil and Plant Measurements
Four soil cores (5-cm i.d.) were collected near the center of each grid using a tractor-mounted hydraulic soil sampler on 18 and 25 April in 2002 and 2003, respectively. One core was taken directly over each grid center, and three additional cores were taken at 1-m radii from the reference point. Cores for each reference point were sectioned into depths of 0 to 15, 15 to 30, 30 to 60, and 60 to 90 cm and composited with depth. Samples were dried in a forced-draft oven at 50°C and then ground with a flail-type soil grinder (Custom Lab, Orange City, FL). The field used during this study received no N fertilizer so that soil Nmin potential and residual NO3N could be more easily correlated with crop yield.
Corn variety Dekalb 687 (Monsanto, St. Louis, MO) was planted on 21 Feb. 2002 and 28 Feb. 2003 with a Case/IH Early Riser planter (Racine, WI) at a rate of
65000 seed ha1. Bicep herbicide {metolachlor [2-chloro-N-(2-ethyl-6-methylphenyl)-N-(2-methoxy-1-methylethyl)acetamide]/atrazine [6-chloro-N-ethyl-N'-(1-methylethyl)-1,3,5-triazine-2,4-diamine]} was used for weed control, along with in-season cultivation.
Soil Nmin potential, C mineralization, soil organic C (SOC), soil total N, and soil texture were determined on 0- to 15-cm soil samples. Residual NO3N and other plant essential nutrients were determined on soil samples from all depths. Soil C mineralization and Nmin were determined according to Franzluebbers et al. (1994a)( 1994b). Approximately 20 g of oven-dried soil were placed in 50-mL beakers, wetted to 30 J kg1, and incubated at 25°C in airtight containers along with a vial containing 10 mL of 1.0 M KOH and another vial containing water to maintain humidity. Vials of KOH were replaced at 1, 10, and 24 d. Mineralized C as CO2 was determined at each sampling date by titrating KOH with 1.0 M HCl to the phenolphthalein endpoint (Anderson, 1982). Ammonium and NO3N at 0 and 24 d were extracted with 2 M KCl and determined using autoanalyzer techniques (Technicon Industrial Systems, 1977a, 1977b). Initial soil mineral NH4N and NO3N were subtracted from that measured at 24 d to determine net soil Nmin.
Soil organic C was determined using the modified Mebius method (Nelson and Sommers, 1982) while soil total N was determined by autoanalyzer techniques (Technicon Industrial Systems, 1977a) following Kjeldahl digestion (Nelson and Sommers, 1980). Soil particle size distribution was determined on all samples using the procedure of Day (1965), which utilizes hydrometer analysis following dispersion of soil by both chemical and physical means.
Three meters of each of the middle two rows of each grid were hand-harvested for grain yield at maturity and threshed using a stationary plot thresher. Grain moisture was determined by electrical resistance, and yields were calculated at a moisture content of 140 g kg1.
Statistical and Spatial Variability Analyses
Statistical relations between soil Nmin potential and other measured soil properties were calculated by correlation analysis. Correlation and spatial statistics were used to relate soil Nmin potential, surface and profile residual soil NO3N, and other soil characteristics with corn grain yield. Geostatistical methods, variography, and kriging (Isaaks and Srivastava, 1989) were used to map variability of soil and plant parameters. Variograms were constructed to characterize the structure of the spatial variation of measured soil and plant parameters. Spatial variability between samples was measured by calculating the omnidirectional semivariance,
(h):
where N(h) is the number of pairs of observations separated by lag distance h and Z(xi) and Z(xi + h) give the value of the variable Z at two positions separated by h. Geostatistical software (GS+ v5.0, Gamma Design Software, St. Plainwell, MI) was used to analyze the spatial structure of the data and to define the semivariograms. Variograms quantify the spatial variation of a range of properties by measuring the degree of correlation between sampling points at a given distance apart (Webster, 1985).
Spherical, exponential, or linear models were fitted to the semivariograms, and their selections were based on visual best fit and the corresponding coefficient of determination, i.e., r2. The parameters of the modelnugget semivariance, range, and sill, or total semivariancewere calculated. Nugget semivariance is the variance at zero distance and represents field and experimental variability, or random variability, which is undetectable at the sampling scale. Sill is the lag distance between measurements at which one value for a variable does not influence neighboring values. Range is the distance at which values of one variable become spatially independent of another. Furthermore, the ratio of the nugget to sill indicates the degree of randomness in the data's spatial variability. This ratio was used to define three classes of spatial dependence for the measured soil variables (Cambardella et al., 1994), i.e., (i) when the ratio was <0.25, the measured variable was considered strongly spatially dependent; (ii) when the ratio was between 0.25 and 0.75, the soil variable was considered moderately spatially dependent; and (iii) if the ratio was >0.75, or the slope of the semivariogram was
0, the variable was considered random or nonspatially correlated (pure nugget).
Semivariogram models were cross-validated, and the punctual kriging procedure in GS+ v5.0 was used to obtain point estimates of soil parameters at unsampled locations. Kriging was performed every 2.1 m across the entire study area. Counter maps of corn yield distribution and clay content percentage were produced in Surfer (Version 8, Golden Software, Golden, CO) based on grid files created from the kriged values from GS+.
 |
RESULTS AND DISCUSSION
|
|---|
Descriptive Statistics
Univariate statistical analyses of soil and plant parameters for 2002 and 2003 are given in Tables 1 and 2. The parameters measured varied spatially and temporally. The shape of parameter distributions are described by skewness. Skewness coefficients in Tables 1 and 2 showed no shifts in distribution, except for SOC, NO3N concentrations with depth, and initial inorganic N. Grain yield and N measurements were higher, and residual NO3N to 60-cm depth was positively skewed and exhibited nonnormal distributions in 2002. Management and temporal effects seem to be likely reasons for nonnormal distribution. Grain yield and most of the measured parameters, however, had lower values, were less variable, and followed a normal distribution with less skewness in 2003. Soil residual NO3N values exhibiting skewed distributions have been reported to be log normally distributed (Cambardella et al., 1994). However, a logarithmic transformation of residual NO3N values in our study reduced this parameter's skewness, but values were still nonnormally distributed (0.05 significance level). The results in 2003 were similar to those of Parkin et al. (1988) and Cambardella et al. (1994), who showed that soil N measurements tended to become more normally distributed in years following N application.
The coefficient of variation (CV) is a useful statistic for measuring the spatial variability of soil properties. Among the N parameters, Nmin values tended to be the most variable, with CV values of 27 and 46% in 2002 and 2003, respectively (Tables 1 and 2). The range in Nmin for the 63 samples collected across the field in 2002 was from 6 to 35 mg kg1, with a mean value of 20 mg kg1, and from
1 to 33 mg kg1, with a mean value of 16 mg kg1 in 2003. The Nmin variability observed in this experiment was similar to CV values of 24 and 36.4% reported by Goovaerts and Chiang (1993) and by Cambardella et al. (1994), respectively. However, CV values for Nmin in our experiment were higher than those reported by Robertson et al. (1988) and Mahmoudjafari et al. (1997). Differences in variability may result from different incubation methods and whether disturbed or undisturbed soil samples were used. Total N, organic N, and NO3N values with depth were characterized by CVs ranging from 12 to 22% in 2002 and 11 to 16% in 2003. Clay content had a CV of 10%.
Spatial Variability
To characterize spatial variability, spherical, exponential, and linear isotropic models were fitted to calculated semivariograms for soil parameters and yield, and their spatial dependence was categorized into three classes based on the nugget/sill ratio suggested by Cambardella et al. (1994) (Tables 3 and 4). Corn grain yield was moderately spatially distributed both years. The structure of variations for grain yield, total N, and Nmin in 2003 remained the same, with moderate and strong spatial dependence, respectively. However, the spatial dependence for residual NO3N became weak or random in 2003, indicating spatial independence. Others have also shown the spatial dependence of NO3N to be temporally dependent (Bruckler et al., 1997). Other parameters, such as initial soil NH4N and NO3N concentrations and soil C mineralized, were weakly or not spatially correlated in either 2002 or 2003. Soil organic C exhibited moderate spatial dependence in 2002 but was randomly distributed in 2003. In general, less spatial dependence was observed between soil samples in 2003 than 2002.
Semivariograms showed strong spatial dependence for soil NO3N with depth and clay content in 2002, mineralized N in 2003, and total N in both years (Fig. 1). Variograms with spherical and exponential models reached upper bounds, i.e., sills, suggesting that the properties varied in a "patchy" way, resulting in some areas with small values and others with large ones. The range of spatial correlation for each variogram provides an average extent of these patches. The range of spatial dependence for the exponential model is the distance at which the semivariance is 95% of the sill and is estimated as three times the fitted range parameter (Isaaks and Srivastava, 1989) (Tables 3 and 4). The range of spatial dependence for NO3N with depth was 31 to 120 m in 2002. Similarly, Lopez-Granados et al. (2002) found that subsoil NO3N concentrations had a strong spatial dependence. The range of the variogram for NO3N with depth in 2003, however, increased without reaching an upper bound, indicating that the extent of the spatial variation was not realized at the sampling scale employed.

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 1. Variograms for (ad) soil NO3N with depth, (ef) total N, (g) mineralized N, and (h) percentage clay of the surface soil (015 cm) in 2002 and 2003.
|
|
Relationship between Soil Nitrogen Parameters and Texture and Yield
To calculate the spatial variability of crop yield and to help develop variable N rate management plans, researchers have used crop and soil properties like plant height, plant population, and clay content (Cox et al., 2003; Katsvario et al., 2003; Machado et al., 2000; Tabor et al., 1985). Spatial distributions of corn grain yield and surface soil (015 cm) clay content of the research field in 2002 are shown in Fig. 2. In each figure, light shading represents the highest values while darker shading is associated with the lowest values. Grain yields were higher along the southern perimeter of the field, in the transect that extended from southwest to north, and in bands toward the northeast across the field. Lower grain yields tended to be located along a southeast to northwest transect and in pockets along eastwest bands near the central portion of the field. Clay content of the surface 0 to 15 cm followed the opposite trend to yield distribution and tended to be higher in northern and eastern portions and lower in the southwestern quadrant (Fig. 2b).
Corn grain yield and clay contents were correlated against residual soil NO3N with depth and other soil characteristics in 2002 and 2003 (Tables 5 to 7). Total N, SOC, Nmin, soil NO3N values for all depths, and clay content were correlated in 2002 (Table 5). Grain yield was positively correlated with residual soil NO3N. Residual soil NO3N to 60-cm depth had a higher correlation with yield than NO3N to either 15- or 30-cm depths. Total N was significantly related to SOC, Nmin, clay content, and NO3N to 15- cm depth in 2002. Goovaerts and Chiang (1993) reported similar relationships for plots with long-term fallow rotations. Mineralizable N was not related to yield but was highly related to SOC and soil total N.
View this table:
[in this window]
[in a new window]
|
Table 6. Pearson correlation coefficient mean values for corn grain yield and soil N parameters for two transects of high-yielding plots (19) and low-yielding plots (4654) in 2002.
|
|
Corn grain yield was negatively correlated with clay content (Table 5). Soil organic C, total N, and soil N mineralized in 24 d, however, were all positively and significantly related to soil clay content (Table 5). Clay protection of adsorbed organic compounds against microbial oxidation may at least partially explain this result. Soil organic C and total N were significantly correlated with soil N mineralized in 24 d. Grain yield was negatively related to clay content, and since many of the soil characteristics (SOC, total N, etc.) were positively related to clay content, they resulted in negative correlations with grain yield. Johnson et al. (2002) showed a similar result for cotton lint yield and SOC and extractable soil P.
To further illustrate the spatial distribution of corn yield and its relation with soil parameters, two field transects in 2002 with high-yielding plots (1 to 9) and low-yielding plots (46 to 54) (Fig. 2a) were selected. Field transect mean values for yield, NO3N, Nmin, total N, SOC, and clay content with their calculated correlation coefficients are presented in Table 6. Nitrate N and clay content were significantly related to yield and seemed to be responsible for variation in high- and low-yield plots, respectively. Residual NO3N was positively related with yield in both high- and low-yield environments while clay content exhibited a negative relationship with grain yield in low-yielding plots. Total N in low-yielding plots was positively related to NO3N and Nmin while clay content and Nmin were negatively related. Clay protection of adsorbed organic compounds against microbial oxidation may again at least partly explain this result. In high-yielding plots where clay was lower, clay content had a significant positive relationship with Nmin. Mineralized N, however, was not related to corn grain yield in either high- or low-yielding plots (Table 6).
Corn grain yield in 2003 was less related to soil NO3N (Table 7) compared with 2002 (Tables 5 and 6). Part of this difference was because of the wide difference in the quantity of NO3N between years, with N availability becoming more limiting in 2003 (Tables 1 and 2). The rate of N applied before this experiment was 125 kg N ha1 as N solution (3200) for cotton at planting in 2001. Residual NO3N measured in 2002 was sufficient for corn yields of 6271 to 9344 kg ha1 (Table 1) (McFarland et al., 1990). However, residual NO3N was more limited in 2003 when NO3N to a depth of 60 cm decreased by 50% and average yield declined by 2132 kg ha1 (Tables 1 and 2). This temporal dependence of NO3N availability may support the need for consideration of soil N supplied through mineralization (Raun et al., 2002; Mamo et al., 2003). In our study, Nmin was not significantly related to corn grain yield in either year. However, Nmin was significantly related to NO3N to 90-cm depth in 2003 (Table 7), suggesting that a greater proportion of residual NO3N originated from mineralization in 2003. Mineralizable N apparently became more significant to the overall N budget as residual NO3N diminished due to plant uptake, leaching, or denitrification. Clay content also was not related to corn grain yield in 2003, suggesting that a different relationship may occur under limited N supply.
Contour Maps of Spatial Distribution of Variables
After constructing the variograms for soil N parameters (Fig. 1), the corresponding kriged values were plotted to form contour maps of the spatial distribution of residual NO3N, total N, and Nmin and to estimate N values in unsampled locations in 2002 and 2003 (Fig. 3). For example, high soil NO3N concentrations were located close to the western perimeter of the field (Fig. 3a3d) in areas with low clay content. This is particularly evident when Fig. 2b and 3a3d are compared.
A comparison of kriged maps (Fig. 3a3d) of soil NO3N with depth and spatial distribution of corn yield (Fig. 2a) showed a similarity, especially in the western portions of a southnorth band. Tabor et al. (1985) suggested that clay content can be used for estimating NO3 concentrations throughout the growing season. Our results showed that clay content exhibited a similar variation in distribution structure (Fig. 1h) as NO3N in 2002 (Fig. 1a1d), suggesting that spatial structure and an optimum sampling program for soil NO3N can be inferred from the spatial structure of clay contents for particular soil mapping units (Cambardella and Karlen, 1999). However, the large nugget semivariance and the nonspatial dependence of residual soil NO3N in 2003 suggested that a different sampling strategy would have had to have been used to capture soil NO3N variability in 2003 and consequently a different strategy for site-specific N application. For example, Eghball et al. (1997) found that when NO3N was low, variable-rate N application did not significantly reduce residual soil NO3N variability. This result could have been true for our field in 2003 where spatial structure for NO3N was not evident, and its concentration was low, and uniform applications of N across the field would have been more justified. Bruckler et al. (1997) explained that when spatial distribution structure exists, grid sampling is more efficient than random sampling, as might have been the case for applying variable-rate N for corn in 2002.
During our experiment, N fertilizer was not applied, and based on the results presented, the distance and sampling scheme used for soil NO3N was temporally dependent. Bruckler et al. (1997) also reported that spatial dependence of soil NO3N was time dependent. It can be hypothesized that strongly spatially dependent soil variables may be controlled by intrinsic variation in soil characteristics such as was observed for soil texture in 2002 (Cambardella et al., 1994). Machado et al. (2000) indicated that management zones for variable-rate N application should be based on information about soil elevation, texture, and soil nitrate. Cambardella and Karlen (1999) suggested, however, that this may be true only for one or two crops after large fertilizer applications. Extrinsic factors appear to play a larger role in controlling spatial dependence when there are two or more crops after N application and low or no fertilizer input as occurred in 2003.
In contrast to residual NO3N distribution maps, soil total N exhibited temporal persistence, i.e., high and low values generally occurred at the same locations in 2002 and 2003. These distinct zones of high and low total N may provide useful information for variable-rate N management. According to Goovaerts and Chiang (1993), temporal persistence of total N reduces the cost of sampling and analysis for N recommendations. Total N was also correlated with clay content for both 2002 and 2003. These relationships indicated that spatial relatedness may bridge several soil map units, and information derived from one set of measurements for one field may have applicability for fertilizer recommendations at other field sites within the same or similar landscapes. Cox et al. (2003) studied the effect of variability of selected soil chemical and physical properties on soybean yield. Soil fertility varied from year to year and field to field. The authors found that fertility parameters had to be considered along with other soil factors to determine their relationship to yield and reported that areas with higher clay content also had higher yield, suggesting that soil clay content might be used as a basis for site-specific soil management. In our experiment, soil total N was also related to NO3N at 15-cm depth. However, as observed by Schmidt et al. (2002), variable-rate N recommendations based solely on soil total N are too simplistic to reflect variability in N availability within a field.
 |
CONCLUSIONS
|
|---|
Most of the measured variables showed temporal variation except for soil total N. Temporal persistence of soil N parameters may have important implications for N fertilizer recommendations. Spatial structure existed for residual NO3N in the first year when concentrations were higher. Spatial structure was absent, however, the following year when NO3N concentrations were lower. Residual soil NO3N sampled to 60-cm depth was more related to corn yield than shallower sampling. Soil Nmin was not related to residual soil NO3N when NO3 was high in 2002 but was related to NO3 when lower in 2003.
Results showed that clay content exhibited a similar variation in distribution structure as NO3N in 2002, suggesting that spatial structure and an optimum sampling program for soil NO3N can be inferred from the spatial structure of clay contents for particular soil mapping units. However, the nonspatial dependence of residual soil NO3N in 2003 suggested that a different sampling strategy would have had to be used to capture soil N variability in 2003. Uniform application of N fertilizer would have resulted in overfertilization in 2002 in some areas and underapplication in others.
This experiment suggests that variable-rate fertilizer recommendations should not be based on a single criterion. Clay content with limited sampling for NO3N might be used to estimate and infer N requirement. Texture may be more easily determined at the scale and intensity necessary for site-specific prediction. Soil textural influences on spatial variation of soil N parameters influencing crop yield require further research in larger fields.
 |
REFERENCES
|
|---|
- Anderson, J.P.E. 1982. Soil respiration. p. 831866. In A.L. Page et al. (ed.) Methods of soil analysis. Part 2. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Bakhsh, A., D.B. Jaynes, T.S. Colvin, and R.S. Kanvar. 2000a. Spatio-temporal analysis of yield variability for cornsoybean field in Iowa. Trans. ASAE 43:3138.
- Bakhsh, A., D.B. Jaynes, T.S. Colvin, R.S. Kanvar, and U.S. Tim. 2000b. Using soil attributes and GIS for interpretation of spatial variability in yield. Trans. ASAE 43:819828.
- Bruckler, L., A.M. Cockborne, P. Rebault, and C. Claudot. 1997. Spatial and temporal variability of nitrate in irrigated salad crops. Irrig. Sci. 17:5361.[CrossRef]
- Cabrera, M.L., and D.E. Kissel. 1988. Mineralizable nitrogen in disturbed and un-disturbed soil samples. Soil Sci. Soc. Am. J. 52:10101015.[Abstract/Free Full Text]
- Cambardella, C.A., and D.L. Karlen. 1999. Spatial analysis of fertility parameters. Precis. Agric. 1:514.
- Cambardella, C.A., T.B. Moorman, J.M. Novak, T.B. Parkin, D.L. Karlen, R.F. Turco, and A.E. Konopka. 1994. Field scale variability of soil properties in central Iowa soils. Soil Sci. Soc. Am. J. 58:15011511.[Abstract/Free Full Text]
- Cox, M.S., P.D. Gerard, M.C. Wardlaw, and M.J. Abshire. 2003. Variability of selected soil properties and their relationships with soybean yield. Soil Sci. Soc. Am. J. 67:12961302.[Abstract/Free Full Text]
- Crowin, D.L., S.M. Lesch, P.J. Shouse, R. Soppe, and J.E. Ayars. 2003. Identifying soil properties that influence cotton yield using soil sampling directed by apparent soil electrical conductivity. Agron. J. 95:352364.[Abstract/Free Full Text]
- Day, P.R. 1965. Particle fractionation and particle-size analysis. p. 545567. In C.A. Black et al. (ed.) Methods of soil analysis. Part 1. 1st ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Dinnes, D.L., D.L. Karlen, D.B. Jaynes, T.C. Kaspar, J.L. Hatfield, T.S. Colvin, and C.A. Cambardella. 2002. Nitrogen management strategies to reduce nitrate leaching in tile-drained Midwestern soils. Agron. J. 94:153171.[Abstract/Free Full Text]
- Eghball, B., R.B. Ferguson, G.E. Varvel, G.W. Hergert, and C.A. Gotway. 1997. Fractal characterization of spatial and temporal variability in site-specific and long term studies. p. 339348. In M.M. Novak and T.G. Dewey (ed.) Fractal frontiers. World Sci., Singapore.
- Eghball, B., J.S. Schepers, M. Neghaban, and M.R. Schlemmer. 2003. Spatial and temporal variability and corn yield: Multifractal analysis. Agron. J. 95:339346.[Abstract/Free Full Text]
- Franzluebbers, A.J., F.M. Hons, and D.A. Zuberer. 1994a. Long-term changes in soil carbon and nitrogen pools in wheat management systems. Soil Sci. Soc. Am. J. 58:16391645.[Abstract/Free Full Text]
- Franzluebbers, A.J., F.M. Hons, and D.A. Zuberer. 1994b. Seasonal changes in soil microbial biomass and mineralizable C and N in wheat management systems. Soil Biol. Biochem. 26:14691475.[CrossRef]
- Frogbrook, Z.L., M.A. Oliver, M. Salahi, and R.H. Ellis. 2002. Exploring the spatial relations between cereal yield and soil chemical properties and implications for sampling. Soil Use Manage. 18:19.[CrossRef]
- Goovaerts, P., and C.M. Chiang. 1993. Temporal persistence of spatial patterns for mineralizable nitrogen and selected soil properties. Soil Sci. Soc. Am. J. 57:372381.[Abstract/Free Full Text]
- Isaaks, E.H., and R.M. Srivastava. 1989. Applied geostatistics. Oxford Univ. Press, New York.
- Johnson, R.M., R.G. Downer, J.M. Bradow, P.J. Bauer, and J.E. Sadler. 2002. Variability in cotton fiber yield, fiber quality and soil properties in a southeastern Coastal Plain. Agron. J. 94:13051316.[Abstract/Free Full Text]
- Katsvario, T.W., W.J. Cox, and H.M. Van Es. 2003. Spatial growth and nitrogen uptake variability of corn at two nitrogen levels. Agron. J. 95:10001011.[Abstract/Free Full Text]
- Lopez-Granados, F., M. Jurado-Exposito, S. Atenciano, A. Garcia-Ferrer, M.S. de la Orden, and L. Garcia-Torres. 2002. Spatial variability of agricultural soil parameters in southern Spain. Plant Soil 246:97105.[CrossRef]
- Machado, S., E.D. Bynum, Jr., T.L. Archer, R.J. Lascano, L.T. Wilson, J. Bordovsky, E. Segarra, K. Bronson, D.M. Nesmith, and W. Xu. 2000. Spatial and temporal variability of corn grain yield: Site-specific relationships of biotic and abiotic factors. Precis. Agric. 2:359376.
- Mahmoudjafari, M., G.J. Kluitenberg, J.L. Havlin, J.B. Sisson, and A.P. Schwab. 1997. Spatial variability of nitrogen mineralization at the field scale. Soil Sci. Soc. Am. J. 61:12141221.[Abstract/Free Full Text]
- Mamo, M., G.L. Malzer, D.J. Mulla, D.R. Huggins, and J. Strock. 2003. Spatial and temporal variation in economically optimum nitrogen rate for corn. Agron. J. 95:958964.[Abstract/Free Full Text]
- McFarland, M.L., F.M. Hons, and V.A. Saladino. 1990. Effects of furrow diking and tillage on corn grain yield and nitrogen accumulation. Agron. J. 83:382386.
- Nelson, D.W., and L.E. Sommers. 1980. Total nitrogen analysis for soil and plant tissues. J. Assoc. Off. Anal. Chem. 63:770778.
- Nelson, D.W., and L.E. Sommers. 1982. Total carbon, organic carbon, and soil organic matter. p. 539577 In A.L. Page et al. (ed.) Methods of soil analysis. Part 2. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Oberle, S.L., and D.R. Keeney. 1990. Factors influencing corn fertilizer N requirements in the northern U.S. corn belt. J. Prod. Agric. 3:527534.
- Parkin, T.B., J.J. Meisinger, S.T. Chester, J.L. Starr, and J.A. Robinson. 1988. Evaluation of statistical estimation methods for lognormally distributed variables. Soil Sci. Soc. Am. J. 52:323329.[Abstract/Free Full Text]
- Raun, W.R., J.B. Solie, G.V. Johnson, M.L. Stone, R.W. Mullen, K.W. Freeman, W.E. Thomason, and E.V. Lukina. 2002. Improving nitrogen use in cereal grain production with optical sensing and variable rate application. Agron. J. 94:815820.[Abstract/Free Full Text]
- Robertson, G.P., M.A. Hutson, F.C. Evans, and J.M. Tiedje. 1988. Spatial variability in a successional plant community: Patterns of nitrogen availability. Ecology 69:15171524.[CrossRef][ISI]
- Sadler, E.J., W.J. Busscher, P.J. Bauer, and D.L. Karlen. 1998. Spatial scale requirements for precision farming: A case study in the southeastern USA. Agron. J. 90:191197.[Abstract/Free Full Text]
- Schmidt, J.P., A.J. DeJoia, R.B. Ferguson, R.K. Taylor, R.K. Young, and J.L. Havlin. 2002. Corn yield response to nitrogen at multiple in-field locations. Agron. J. 94:798806.[Abstract/Free Full Text]
- Stafford, J.V., B. Ambler, R.M. Lark, and J. Catt. 1996. Mapping and interpreting the yield variation in cereal crops. Comput. Electron. Agric. 14:101119.
- Tabor, J.A., W. Warrick, D.E. Myers, and D.A. Pennington. 1985. Spatial variability of nitrate in irrigated cotton: II. Soil nitrate and correlated variables. Soil Sci. Soc. Am. J. 49:390394.[Abstract/Free Full Text]
- Technicon Industrial Systems. 1977a. Determination of nitrogen in BS digests. Method 33474W/B. Technicon Ind. Syst., Tarrytown, NY.
- Technicon Industrial Systems. 1977b. Nitrate and nitrite in soil extracts. Method 48777A. Technicon Ind. Syst., Tarrytown, NY.
- Vanotti, M.B., and L.G. Bundy. 1994. Corn nitrogen recommendations based on yield response data. J. Prod. Agric. 7:249256.
- Webster, R. 1985. Quantitative spatial analysis of soil in the field. Adv. Soil Sci. 3:170.