Agronomy Journal 95:339-346 (2003)
© 2003 American Society of Agronomy
NITROGEN MANAGEMENT
Spatial and Temporal Variability of Soil Nitrate and Corn Yield
Multifractal Analysis
Bahman Eghball*,a,
James S. Schepersa,
Mehrdad Negahbanb and
Michael R. Schlemmera
a USDA-ARS, 121 Keim Hall, Univ. of Nebraska, Lincoln, NE 68583
b Dep. of Eng. Mechanics, Univ. of Nebraska, Lincoln, NE 68583
* Corresponding author (beghball1{at}unl.edu)
Received for publication February 9, 2002.
 |
ABSTRACT
|
|---|
High levels of residual soil NO3N can contaminate ground water by leaching through the soil. Our objective was to reduce the level and spatial variability of residual soil NO3N while maintaining optimum corn (Zea mays L.) production by variable rate N fertilizer application. The experiment was located on a 60-ha sprinkler-irrigated corn field in central Nebraska and included four N management practices: uniform rate, variable rate (VRAT), variable rate at 75% of recommended amount (VRAT @ 75%), and variable rate plus 10% (VRAT + 10%). VRAT @ 75% decreased the amount of residual NO3N in the soil while maintaining similar grain yield to the other treatments, indicating over-application of N with treatments receiving the recommended rate. Increasing the recommended rate by 10% (VRAT + 10%) did not increase corn yield or residual soil NO3N. Based on multifractal spectrum, no consistent pattern of spatial variability of soil NO3N was observed for each treatment across years. Spatial variability in corn grain yield was much lower than that for soil NO3N, indicating noneffectiveness of using soil NO3N spatial distribution for variable rate N application unless some areas in the field are severely N deficient. Variable rate N application did not reduce variability of residual soil NO3N or corn grain yield as compared with uniform N. Multifractal analysis quantitatively characterized the extent and pattern of spatial and temporal variability in corn grain yield and residual soil nitrate.
Abbreviations: adiff, the distance between minimum and maximum a values of each multifractal spectrum CEC, cation exchange capacity VRAT, variable rate VRAT @ 75%, variable rate at 75% of the recommended amount VRAT + 10%, variable rate of the recommended amount plus 10%
 |
INTRODUCTION
|
|---|
RECENT DEVELOPMENTS in agricultural technology have made site-specific fertilizer application a reality. Variable rate (site-specific) N application should provide the plant with the appropriate amount of N while reducing the quantity and variability of residual soil NO3N after harvest. One may also expect to find a more homogeneous yield response across the field following adoption of variable rate N application. By reducing variability and quantity of residual soil NO3N, its leaching and subsequent ground water contamination potential should be reduced. Eghball et al. (1999) found that the extent of variability in residual soil NO3N was significantly reduced following adoption of variable rate N application in a continuous corn system under gravity irrigation. The residual soil NO3N to a depth of 0.9 m was high (avg. 6.8 mg kg-1, max. 12.0 and min. 2.4) across the field before initiation of variable rate N application. After 1-yr variable rate N application, average residual soil NO3N was 5.0 mg kg-1 with a maximum of 7.9 and a minimum of 3.7. In another study where residual soil NO3N was low (avg. 4.0 mg kg-1, max. 7.8 and min. 1.5), variable rate N application did not significantly reduce residual soil NO3N variability (Eghball et al., 1997). Ferguson et al. (2002) found reduction in soil nitrate concentration due to variable rate fertilizer N application in only 3 out of 12 site-years as compared with uniform N application. Machado et al. (2000) indicated that management zones for variable rate fertilizer and water applications should be based on information about soil elevation, texture, and soil nitrate. Spatial dependence of soil NO3N was found to be time dependent in irrigated salad crops (Bruckler et al., 1997).
Fractal analysis can provide insight into the spatial or temporal variability of crop or soil parameters. Fractal analysis has been shown to be useful in a variety of scientific disciplines. The use of fractals for numerical analysis of soil and plant parameters is still a relatively new technique. It has been used for characterizing soil structure (Eghball et al., 1993b; Perfect and Blevins, 1997), soil chemical and physical parameters (Burrough, 1983), root morphology (Eghball et al., 1993a), temporal yield variations (Eghball and Power, 1995; Eghball and Varvel, 1997), and spatial variability of soil and crop yield (Eghball et al., 1997, 1999). Fractal analysis was found to be useful in characterizing soil and plant parameters that was not possible or very difficult to do before. Fractal dimension (D) of a curve can have a value between 1 and 2, giving a quantitative indication of the function's shape or roughness.
Multifractal analysis has been proposed for determination of spatial variability of soil parameters (Folorunso et al., 1994; Kravchenko et al., 1999, 2000). Multifractal parameters were found to reflect many of the major aspects of variability in soil properties, provided a unique quantitative characterization of the data spatial distribution, and multifractal parameters were useful in choosing an appropriate interpolation procedure for mapping soil properties (Kravchenko et al., 1999). Multifractal analysis was used to characterize particle-size distribution of soils with wide range of particle sizes (Posadas et al., 2001). A single fractal dimension might not be sufficient to characterize soil spatial variability because of the heterogeneous nature of soil parameters. A set of fractal dimensions, called a multifractal spectrum, is referred to as multifractal analysis (Frisch and Parisi, 1985). Multifractal analysis needs to be evaluated to determine its usefulness in comparing spatial variability of soils treated with different treatments. The objective of this study was to characterize and compare spatial and temporal variability of residual soil NO3N and corn grain yield in a variable rate N study using multifractal analysis.
 |
MATERIALS AND METHODS
|
|---|
Field Treatments
An experiment was conducted from 1994 to 1997 on a 60-ha center-pivot irrigated corn field located in central Nebraska. The soil types within the field included 40% Blendon loam, 0 to 1% slope (coarse loamy, mixed, superactive, mesic Pachic Haplustolls), 20% Blendon loam, 1 to 3% slope, and 40% Hord silt loam, 0 to 1% slope (fine-silty, mixed, superactive, mesic Cumulic Haplustolls). Growing season rainfalls (1 May31 October) were 418, 448, 528, and 458 mm while average temperature [(maximum + minimum) ÷ 2] were 19.3, 18.1, 18.3, and 18.6°C for 1994, 1995, 1996, and 1997, respectively. Four N management practices were applied to 32-row wide strips (24.4 m) that ran the entire length of the field (780 m). The management practices used were arranged in a randomized complete block design with five replications. Soil samples were collected in the spring to determine NO3N level in the soil, which was utilized to calculate the N application rate for each treatment. The treatments were (i) fixed uniform N rate based on a strip average of soil NO3N and organic matter obtained from grid sampling, (ii) variable rate N applied at 100% of recommended rate determined at each of the grid sample points based on the soil NO3N and organic matter found at that point, (iii) variable rate N applied at 75% of recommended rate with the remainder being applied through fertigation if needed based on chlorophyll meter readings (N was applied only once in 1996 at a rate of 34 kg ha-1 using a high clearance applicator), and (iv) variable rate N plus an additional 10% of the recommended rate. All practices used the University of Nebraska N recommendation equation for corn (Hergert et al., 1995). The NO3N in irrigation water (26 mg L-1) for an expected application of 23 cm water yr-1 was considered when recommended amount of N was calculated. Nitrogen fertilizer was applied as a sidedress application of anhydrous NH3 at growth stages V6 to V9. Anhydrous NH3 was applied with a toolbar-mounted coulter/knife injection unit placed into the furrow midway between plant rows. Nitrogen application rate was set either manually for the uniform N treatment, or adjusted according to field position by a SoilTeq1 Falcon controller for the VRAT treatments. Nitrogen application maps were developed using SoilTeq SGIS software with grid soil sample data. Treatments were applied to the same strips each year. Nitrogen application rates are given in Table 1. The soil nitrate and grain yield maps were generated by kriging the data using GS+ (Gamma Design Software, Plainwell, MI).
Soil samples were collected on a 12.2 by 24.4 m staggered grid to a depth of 0.9 m. Nitrate content of the soil samples was determined by extracting the soil with 2 M KCl and then using a Lachat system (Zellweger Analytics, Milwaukee, WI). A yield monitor linked to the global positioning system (John Deere Green Star) was used for the 1995, 1996, and 1997 harvest. Yield was obtained at each location that soil was sampled by averaging the yield monitor data from a 12.2-m-side rectangle centered around the soil sampling point.
Statistical analysis of the residual NO3N and grain yield was performed using PROC MIXED (Littell et al., 1996) and adjusting the means for the spatial variability of the data. In this analysis, the means for each treatment is adjusted based on spatial variability of the data in the entire field as determined by geostatistics' parameters (nugget, sill, and range). For details see pages 303 to 330 of Littell et al. (1996). Multifractal analysis was performed using a Fortran computer program written by the authors.
Multifractal Analysis
The field was a square of 780 by 780 m divided into five segments (replications). Each segment contained randomly assigned treatment strips. As a result, the distance between strips of each treatment varied throughout the field. As shown in Fig. 1a
for the uniform N treatment, each sampling point was assigned to a rectangular grid. When analyzing a treatment, it was assumed that the entire field received that treatment. This is similar to grid sampling of certain strips in a field to characterize spatial distribution of certain soil property in the entire field. The sides of these rectangular grids were defined by the midpoints between consecutive sampling points, or were defined by the field boundaries. When a sampling point was missing, depending on whether there were readings on both sides of the point or not, either the reading on one side, the other side, or the average of both sides was assigned to the missing value. When consecutive readings were missing, this process was started from the closest nonzero reading and repeated until every rectangle had a number assigned to it. This defined the data grid.

View larger version (35K):
[in this window]
[in a new window]
|
Fig. 1. The grid arrangements for the multifractal analysis. (a) is the sampling grids of the soil NO3N and grain yield for the uniform N treatment (sampling point is indicated by a dot); (b), (c), and (d) are different multifractal cell sizes.
|
|
For the multifractal analysis, a set of cells was constructed. The field was subdivided into 4, 9, 16, and 25 square (Fig. 1b, 1c, 1d for 4, 9, 16, respectively) cells. As a result, the cell hypotenuse size
was 552, 367, 276, and 221 m, respectively. The hypotenuse size was used instead of the side so that the method would be applicable whether square or rectangular cells were used. For each treatment, actual data from the five strips (replications) were used for the multifractal analysis, with the missing data filled as described above.
The multifractal analysis method was as follows: For each treatment, the mass probability function (ui) for multifractal cell i (Fig. 1b, 1c, 1d) with hypotenuse size
was evaluated as
 | [1] |
where
 | [2] |
and
 | [3] |
where F is the initial yield or NO3N value in each data grid (assuming similar value as the sampling point for the entire area of each grid), the integral is a double integral over the area Ai of multifractal cell i. The integral is a double integral over the area Ai of multifractal cell i (data value was multiplied by its sampling grid area and summed for cell i), and n represents the total number of cells of size
.
The distribution of the probability mass function was then analyzed for multifractality using the method of moments (Evertsz and Mandelbrot, 1992). Briefly, a partition function was determined as follows:
 | [4] |
where q is a real number ranging from -
to
. For multifractally distributed measures, the partition function scales with the grid size as
 | [5] |
where
(q) is the mass exponent of order q. The mass exponent for each q value can be obtained by plotting log
q(
) vs. log
. If the probability function µi(
) in the neighborhood of the grid scales with the grid size as µi(
)

, then as
0, the singularity exponent
is a scaling property specific to the cell size. Parameter
is also called a local fractal dimension or a singularity index. The local fractal dimension can be determined by Legendre transformation of the
(q) curve (Evertsz and Mandelbrot, 1992) as
 | [6] |
and be used to determine the multifractal spectrum (discussed next). The number of cells of size
with the same
, N
(
), is related to the cell size as N
(
)
-f(
), where f(
) is a scaling exponent of the cells with common
. Parameter f(
) can be calculated as
 | [7] |
A plot of f(
) vs. a is called a multifractal spectrum. Multifractal spectrum quantitatively characterizes variability of soil or crop parameters with asymmetry to the right and left indicating domination of small and large values, respectively. The width of the multifractal spectrum (adiff) indicates overall variability similar to the nugget effects in geostatistics. For each treatment we calculated multifractal spectrum with q values ranging from -10 to +10 in increments of 0.2. The f(
) spectrum is related to the commonly used generalized multifractal dimension as
 | [8] |
The fractal dimension at q = 0 is the box-counting dimension of the geometric support of the measure being studied, which in our case was Euclidian dimension of a plane (i.e., 2), information fractal dimension, D1, is obtained at q = 1 using the l'Hôpital's rule, and the correlation fractal dimension, D2, is obtained at q = 2. The lower D1 or D2 values indicate domination of long-range variation while higher values indicate domination of short-range variation.
 |
RESULTS AND DISCUSSION
|
|---|
Corn Grain Yield and Residual Soil Nitrate
The total N application rates were not different between uniform and variable rate N application methods as the ratio of VRAT/uniform N was >0.97 (Table 1). The software used to control fertilizer application when this study was conducted was not capable of recording actual N rate applied; therefore, mean N rates for the VRAT treatments are reported. Variable rate N application did not result in less N application. Corn grain yields were influenced by the N application methods in 1995 and 1997 (Table 2). The coefficients of variation for grain yields were small, indicating low overall variability of grain yield in each year. Variable rate @ 75% resulted in less corn yield than VRAT in 1995, while both VRAT @ 75% and VRAT + 10% resulted in less yield than VRAT in 1997 (Table 2). However, because of low variability, the magnitude of the difference was small (max. of 360 kg ha-1). Adding 10% more N to the recommended rate seems to have negatively influenced corn grain yield in 1997. However, reducing the recommended amount by 25% did not result in less yield than other treatments in 2 out of 3 yr, indicating over-application of N using the recommended rate and/or more water applied than the expected 23 cm (26 mg L-1 NO3N in irrigation water).
The VRAT @ 75% treatment resulted in less residual soil NO3N than other treatments indicated the effects of lowering the N application rate on soil residual NO3N (Table 3). Increasing the application rate by 10% in the VRAT + 10% treatment did not result in greater residual soil NO3N than the uniform or variable rate N treatments (Table 3). The CV values were high for residual soil NO3N, indicating large overall variability. The VRAT and Uniform N application methods resulted in similar residual soil NO3N in all years except in 1996, where VRAT resulted in 0.6 mg kg-1 less soil NO3N than uniform N application (Table 3). Variable rate N application was not effective in reducing residual soil NO3N as compared with uniform application unless the application rate was reduced by 25%. Eghball et al. (1999) also found no difference between variable rate and uniform N applications for spatial variability in residual soil NO3N. The soil-based N recommendation used in this study seems to have underestimated N mineralization from soil. Ferguson et al. (2002) also reported that Nebraska N recommendation equation may not be appropriate for variable rate N application. In irrigated systems, a plant-based N management system utilizing corn canopy remote sensing and fertigation may be a viable alternative to soil-based recommendations because of the uncertainty about N mineralization in the soil. Another alternative might be to apply a reduced rate (perhaps 50%) of N fertilizer based on soil testing results and follow up with variable rate N application as an in-season treatment as needed based on remote sensing or crop stress data (Varvel et al., 1997).
Multifractal Analysis
Coefficient of variation provides an indication of variability in the overall data. Multifractal analysis was performed to provide indication of the pattern and nature of spatial and temporal variability in the soil residual NO3N and corn grain yield data. Advantage of multifractal analysis over traditional statistics and geostatistics were tested by Kravchenko et al. (1999) where they exchanged 40 random subsamples of cation exchange capacity (CEC) values in a field with the highest CEC values. In this case, statistical properties of the data set remained intact and only spatial structure of the data distribution was modified. Relocation did not change the statistical properties or the variogram, but multifractal spectra differentiated between the original and the modified data set.
Monofractal analysis performed on the residual soil NO3N and corn grain yield data indicated significant differences among replications of each treatment for fractal dimensions, pointing out the heterogeneity of variability (data not shown). For a discussion on the use of monofractal analysis for characterizing and comparing spatial and temporal variability see Eghball et al. (1999). Because of heterogeneity of variability, multifractal analysis was performed on the combined data from all five replications of each treatment.
Residual Soil Nitrate
If D(q) values decrease for increasing parameter q
0, then the measure is called multifractal (Peitgen et al., 1992, p. 737). Regression lines of log
q(
) vs. log
at different q values were linear for all treatments in all 4 yr with R2 values > 0.99, indicating excellent fit of the models used. The D(q) values were decreasing for increasing q values, indicating the multifractal nature of spatial variability of soil NO3N (Fig. 2)
. The greatest decrease of D(q) with increasing q was observed for VRAT @ 75% and VRAT + 10% treatments in 1995 (Fig. 2). Multifractal spectrums of residual soil NO3N are presented in Fig. 3
. Asymmetry toward the left from a = 2 indicates domination of large or presence of extremely large values in the spatial variability pattern while asymmetry to the right indicates domination of small or presence of extremely small values. In 1994 when initial soil measurements were made before treatment applications, all four treatments had similar variability patterns with VRAT treatment skewed toward the right of a = 2 (Fig. 3). In 1995, there were significant differences among the treatments for spatial variability patterns with VRAT @ 75% and VRAT + 10% skewed to the left, indicating domination of large values in spatial distribution of soil NO3N for these two treatments (Fig. 3). Soil NO3N maps for the four treatments in 1995 are shown in Fig. 4
. The distribution pattern of soil NO3N for the VAR @ 75% and VAR + 10 treatments indicated large areas of low soil NO3N values (2.54.5 mg kg-1) while VRAT and Uniform treatments had a more uniform distribution of the low and medium values (Fig. 4). The very high soil NO3N values were concentrated in the right corners of the field for VRAT @ 75% and VRAT + 10% treatments. These high values skewed the multifractal spectrum to the left for a < 2. In 1996 and 1997, soil NO3N data also showed differences in spatial variability patterns among treatments, but the patterns were not consistent among years.

View larger version (27K):
[in this window]
[in a new window]
|
Fig. 2. Residual soil nitrate D(q) values as a function of increasing q values for four N treatments in 4 yr.
|
|

View larger version (68K):
[in this window]
[in a new window]
|
Fig. 4. Residual soil NO3N distribution for four treatments in 1995. Each map was generated using the data from the strips of that treatment.
|
|
The D1 and D2 values for residual soil NO3N are given in Table 4. Smaller D1 or D2 values indicate domination of long-range variation while higher values indicate domination of short-range variability in the spatial pattern. The differences among treatments for D1 and D2 of residual soil NO3N were significant for all years except 1994, when the samples were collected before initiation of the treatments. The D1 and D2 values were smaller for the VRAT @ 75% and VRAT + 10% treatments in 1995, indicating domination of long-range variability (Fig. 4). This substantiated the patterns in the multifractal spectrums for these treatments (Fig. 3).
View this table:
[in this window]
[in a new window]
|
Table 4. Information (D1) and correlation (D2) fractal dimensions of residual soil nitrate for four treatments across 4 yr.
|
|
The distance between minimum and maximum a values for each multifractal spectrum (adiff) is an important indicator of variability distribution. Kravchenko et al. (1999) found that adiff was highly correlated with the nugget effects when performing geostatistics on the same data. Analysis of variance of adiff values (using year as replication) indicated no significant effect of year (P = 0.35) or treatment (P = 0.98) main effects across years on adiff. The distribution of variability was not significantly different among treatments across years, pointing out the inconsistency of temporal spatial variability for each treatment. However, the multifractal spectrums were significantly different among treatments in each year, indicating different variability pattern for each (Table 5). The differences among treatments became more pronounced with years following initiation of the treatments.
View this table:
[in this window]
[in a new window]
|
Table 5. The distance between a values of the treatments for the f(a) vs. a lines (adiff) for the residual soil nitrate and grain yields (Fig. 3 and 6, respectively) in various years.
|
|
Corn Grain Yield
Regression lines of log
q(
) vs. log
were linear for all treatments in all 3 yr with R2 values > 0.99. However, the D(q) values remained relatively constant for increasing q values (Fig. 5)
. The multifractal spectrum indicated the least variability in 1995 and slight increase in asymmetry to the right in 1996 and 1997 for the treatments (Fig. 6)
. It seems that the spatial variability was increasingly dominated by the average values of corn grain yield each year. Corn grain yield maps for the four treatments in 1997 are shown in Fig. 7
. As indicated in the maps, the spatial variability was not much different among treatments and no strong domination of small or large values was apparent. The adiff values also were not different among treatments (Table 5). It seems that the strong spatial variability of soil NO3N did not influence spatial variability patterns in corn grain yield, indicating adjustment by corn plants for spatial variability in N availability. Adequate N was available throughout the field, even though soil NO3N distribution showed patterns of high or low soil NO3N areas in the field. Eghball et al. (1997) also observed less spatial variability of corn yield while soil NO3N exhibited much stronger spatial variability in another variable rate N study.

View larger version (24K):
[in this window]
[in a new window]
|
Fig. 5. Corn grain yield D(q) values as a function of increasing q values for four N treatments in 3 yr.
|
|

View larger version (60K):
[in this window]
[in a new window]
|
Fig. 7. Corn grain yield distribution for four treatments in 1997. Each map was generated using the data from the strips of that treatment.
|
|
The D1 and D2 values were all 2 for corn grain yield for all treatments across 3 yr, indicating lack of long-range variation in yield spatial variability. Since long-range variability was unimportant in the corn yield data, small short-range variation became dominant when spatial variability of yield was characterized.
 |
CONCLUSIONS
|
|---|
Variable rate N application did not result in greater corn grain yield or less spatial variability of residual soil NO3N or corn yield than uniform N application. When N application rate was reduced by 25% in the VRAT @ 75% treatment, corn grain yield was basically similar to full rate application but residual soil NO3N was significantly reduced, pointing to underestimation of soil N mineralization in the N recommendation equation. Multifractal analysis indicated significantly greater spatial variability for residual soil NO3N than corn grain yield. However, no consistent patterns of spatial variability of soil residual NO3N and corn grain yield were observed across years. It seems that corn grain yield spatial variability is not significantly influenced by soil NO3N distribution, indicating noneffectiveness of using soil NO3N spatial distribution for managing corn yield variability unless some areas in the field are severely N deficient. Increasing the N application rate by 10% over the recommended rate resulted in significant yield reduction, but no difference in residual soil NO3N as compared with the recommended rate. Variable rate N application based on spatial variability of soil NO3N and organic matter was not more effective in terms of corn grain yield and spatial distribution of residual soil NO3N than uniform N application. Multifractal parameters provided quantitative indications of spatial and temporal variability patterns for soil NO3N and corn grain yield.
 |
ACKNOWLEDGMENTS
|
|---|
This project was funded in part through a grant from the USEPA.
 |
NOTES
|
|---|
Joint contribution of the USDA-ARS and the Univ. of Nebraska Agric. Res. Div., Lincoln, NE, as paper no. 13618.
1 Mention or use of a product does not imply endorsement by the USDA-ARS or the University of Nebraska. 
 |
REFERENCES
|
|---|
- Bruckler, L., A.M. Cockborne, P. Rebault, and B. Claudot. 1997. Spatial and temporal variability of nitrate in irrigated salad crops. Irrig. Sci. 17:5361.
- Burrough, P.A. 1983. Multiscale sources of spatial variation in soil: I. The application of fractal concepts to nested levels of soil variation. J. Soil Sci. 34:577597.
- 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 Scientific, Singapore.
- Eghball, B., G.W. Hergert, G.W. Lesoing, and R.B. Ferguson. 1999. Fractal analysis of spatial and temporal variability. Geoderma 88:349362.
- Eghball, B., L.N. Mielke, G.A. Calvo, and W.W. Wilhelm. 1993b. Fractal description of soil fragmentation for various tillage methods and crop sequences. Soil Sci. Soc. Am. J. 55:13371341.
- Eghball, B., and J.F. Power. 1995. Fractal description of temporal yield variability of ten crops in the United States. Agron. J. 87:152156.[Abstract/Free Full Text]
- Eghball, B., J.A. Settimi, J.W. Maranville, and A.N. Parkhurst. 1993a. Fractal analysis for morphological description of corn roots under nitrogen stress. Agron. J. 85:287289.[Abstract/Free Full Text]
- Eghball, B., and G.E. Varvel. 1997. Fractal analysis of temporal yield variability of crop sequences: Implications for site-specific management. Agron. J. 89:851855.[Abstract/Free Full Text]
- Evertsz, C.J.G., and B.B. Mandelbrot. 1992. Multifractal measures (Appendix B). p. 922953. In H.-O. Peitgen et al. (ed.) Chaos and fractals: New frontiers of science. Springer-Verlag, New York.
- Ferguson, R.B., G.W. Hergert, J.S. Schepers, C.A. Gotway, J.E. Cahoon, and T.A. Peterson. 2002. Site-specific nitrogen management of irrigated maize: Yield and soil residual nitrate effects. Soil Sci. Soc. Am. J. 66:544553.[Abstract/Free Full Text]
- Folorunso, O.A., C.E. Puente, D.E. Rolston, and J.E. Pinzon. 1994. Statistical and fractal evaluation of the spatial characteristics of soil surface strength. Soil Sci. Soc. Am. J. 58:284294.[Abstract/Free Full Text]
- Frisch, U., and G. Parisi. 1985. On the singularity structure of fully developed turbulence. p. 8488. In M. Gil et al. (ed.) Turbulence and predictability in geophysical fluid dynamics. North Holland Publ. Co., Amsterdam.
- Hergert, G.W., R.B. Ferguson, and C.A. Shapiro. 1995. Fertilizer suggestions for corn. NebGuide G74174-A. Coop. Ext., Univ. of Nebraska, Lincoln.
- Kravchenko, A.N., C.W. Boast, and D.G. Bullock. 1999. Multifractal analysis of soil spatial variability. Agron. J. 91:10331041.[Abstract/Free Full Text]
- Kravchenko, A.N., D.G. Bullock, and C.W. Boast. 2000. Joint multifractal analysis of crop yield and terrain slope. Agron. J. 92:12791290.[Abstract/Free Full Text]
- Littell, R.C., G.A. Milliken, W.W. Stroup, and R.D. Wolfinger. 1996. SAS system for mixed models. SAS Inst., Cary, NC.
- 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.
- Peitgen, H.-O., H. Jürgens, and D. Saupe. 1992. Chaos and fractals: New frontier of science. Springer-Verlag, New York.
- Perfect, E., and R.L. Blevins. 1997. Fractal characterization of soil aggregation and fragmentation as influenced by tillage treatment. Soil Sci. Soc. Am. J. 61:896900.[Abstract/Free Full Text]
- Posadas, A.N.D., D. Gimenez, M. Nittelli, C.M.P. Vaz, and M. Flury. 2001. Multifractal characterization of soil particle-size distribution. Soil Sci. Soc. Am. J. 65:13611367.[Abstract/Free Full Text]
- Varvel, G.E., J.S. Schepers, and D.D. Francis. 1997. Ability for in-season correction of nitrogen in corn using chlorophyll meters. Soil Sci. Soc. Am. J. 61:12331239.[Abstract/Free Full Text]
This article has been cited by other articles:

|
 |

|
 |
 
Y. Miao, D. J. Mulla, J. A. Hernandez, M. Wiebers, and P. C. Robert
Potential Impact of Precision Nitrogen Management on Corn Yield, Protein Content, and Test Weight
Soil Sci. Soc. Am. J.,
August 9, 2007;
71(5):
1490 - 1499.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
N. Hong, J. G. White, R. Weisz, M. L. Gumpertz, M. G. Duffera, and D. K. Cassel
Groundwater Nitrate Spatial and Temporal Patterns and Correlations: Influence of Natural Controls and Nitrogen Management
Vadose Zone J.,
January 24, 2007;
6(1):
53 - 66.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
W. K. Jung, N. R. Kitchen, K. A. Sudduth, and S. H. Anderson
Spatial Characteristics of Claypan Soil Properties in an Agricultural Field
Soil Sci. Soc. Am. J.,
June 21, 2006;
70(4):
1387 - 1397.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
B. D. Kay, A. A. Mahboubi, E. G. Beauchamp, and R. S. Dharmakeerthi
Integrating Soil and Weather Data to Describe Variability in Plant Available Nitrogen
Soil Sci. Soc. Am. J.,
May 23, 2006;
70(4):
1210 - 1221.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
N. Hong, J. G. White, M. L. Gumpertz, and R. Weisz
Spatial Analysis of Precision Agriculture Treatments in Randomized Complete Blocks: Guidelines for Covariance Model Selection
Agron. J.,
June 17, 2005;
97(4):
1082 - 1096.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
H. Shahandeh, A. L. Wright, F. M. Hons, and R. J. Lascano
Spatial and Temporal Variation of Soil Nitrogen Parameters Related to Soil Texture and Corn Yield
Agron. J.,
April 27, 2005;
97(3):
772 - 782.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
L. Pozdnyakova, D. Gimenez, and P. V. Oudemans
Spatial Analysis of Cranberry Yield at Three Scales
Agron. J.,
January 1, 2005;
97(1):
49 - 57.
[Abstract]
[Full Text]
[PDF]
|
 |
|