|
|
||||||||
a Dep. of Crop and Soil Sci., Michigan State Univ., East Lansing, MI 48824-1325
b Dep. of Crop Sci., Univ. of Illinois, Urbana, IL 61801
c Agri Business Consultants, 911 Edison Ave., Lansing, MI 48910
* Corresponding author (kravche1{at}msu.edu).
Received for publication December 12, 2002.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: EC, apparent electrical conductivity ECd, apparent electrical conductivity measured at 0- to 90-cm depth ECs, apparent electrical conductivity measured at 0- to 30-cm depth OM, organic matter
| INTRODUCTION |
|---|
|
|
|---|
The other easily obtained intensive field information that has recently gained attention as a potential predictor of crop yield variability is soil EC. Soil EC is of particular interest for site-specific management for three reasons. First, newly developed technology allows one to quickly and inexpensively obtain fast, dense, and accurate soil EC measurements (Kitchen et al., 1999; Sudduth et al., 2001). Second, soil EC is related to several soil properties important to plant growth. These include level of soil compaction; depth to claypan; soil water content; soil sand, silt, and clay contents; soil drainage; total C and N; extractable P; and soil pH (Johnson et al., 2001; Kravchenko et al., 2002; Sheets and Hendrickx, 1995; Sudduth et al., 1998). Third, soil EC provides information about subsoil properties at a range of depths that are important to plant growth. This feature adds to the unique importance of soil EC for site-specific management because neither topographical information nor remote sensing can directly assess subsoil properties.
A number of statistical procedures have been used for studying relationships among crop yield, soil properties, and topography. These include simple linear regression, multiple regression, principal components, and other multivariate techniques (Dieleman et al., 2000; Kravchenko et al., 2002). The main drawback of the traditional regression/multivariate approach is that they ignore data location and the spatial structure of data distributions. From that perspective, methods of characterizing relationships between variables that account for data locations are superior to traditional statistics. In particular, cross-correlograms and standardized cross-variograms allow one to evaluate the strength and the direction of the relationship between crop yield and soil/topography. In addition, they can also be used to examine the spatial aspects of the relationships, such as the range of distances over which the correlation between studied variables exists and the direction in space of the strongest and weakest correlation. By using cross-correlograms, Stein et al. (1997) observed the negative correlation between millet [Pennisetum glaucum (L.) R. Br.] yield and soil cation exchange capacity extending to distances of 30 to 40 m. Cassel et al. (2000) observed significant cross-correlogram values between wheat (Triticum aestivum L.) yield and the depth of the Ap horizon at distances of 25 to 30 m. Kravchenko and Bullock (2002a) observed significant cross-correlogram values between elevation and soybean protein and oil concentrations up to distances of 100 to 150 m. Cross-correlograms were also shown to be informative and useful in delineating the size of potential management zones for site-specific soybean harvesting (Kravchenko and Bullock, 2002b).
The advantage of the cross-correlogram over the correlation coefficient is that the cross-correlogram combines the information about the strength of the correlation with the information about the range of the distances over which the correlation is significant. However, the disadvantage of cross-correlograms is that in contrast to the single number (i.e., the correlation coefficient) used to compare the strength of the relationship between the two variables, comparisons between cross-correlograms are awkward because values at several distances are involved. Hence, when it comes to deciding which soil property is more strongly related to yield, cross-correlograms are not a very convenient comparison tool. Previous studies that used cross-correlograms for yieldtopographysoil relationships were based on descriptive characterizations of the differences in cross-correlogram values observed at certain distances and descriptive comparisons of correlation ranges. More efficient analysis can be achieved using a single parameter that combines both the magnitude of cross-correlogram values, the correlation range, and the cross-correlogram's shape. This parameter can be used to obtain quantitative comparisons of the relationships between crop yield and other variables in their spatial context.
The objectives of this study are (i) to develop a single cross-correlogram parameter that will conveniently combine the information about strength and spatial range of the correlation between variables and (ii) to use the developed parameter to quantify relationships between corn and soybean yield and soil EC and elevation.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Georeferenced EC measurements were taken within each field using a Veris 3100 sensor cart (Veris Technol., Salina, KS) that operates on a principle of electromagnetic induction. The measurements were taken every 3 to 5 m, with a distance between cart passes of about 10 m. Two sets of EC measurements were collected corresponding to depths of approximately 0 to 30 cm (ECs) and 0 to 90 cm (ECd) (Sudduth et al., 1998).
Corn and soybean grain yield were recorded via yield monitors each year during the overall period of 1996 to 2001. The yield measurements were taken on a 1-s interval by grain sensors mounted on a combine, with each site measurement covering an area of about 2 by 5 m. Simultaneously, the site's coordinates were determined by a GPS unit.
Daily precipitation data from the nearest weather-recording station for each location were provided by the Midwestern Climate Center (Illinois State Water Survey, Champaign, IL). The closest weather stations were approximately 10 km away from the WL, WS, and Carr fields; 15 km away from the 210 and F62 fields; and 20 km away from the Cur field.
Comparison of the strength of spatial correlation among yield/ECs, yield/ECd, and yield/elevation for different crops was conducted with PROC MIXED function (SAS Inst., 2000). Statistical model included crop (corn vs. soybean) as a fixed factor and field as a random factor. The importance of the other sources of variability, such as precipitation, soil properties (i.e., OM), and slope, was explored by including them as covariates in the statistical model.
| THEORY |
|---|
|
|
|---|
![]() | [1] |
21-h and
22+h are the head and tail variances of Variable 1 and Variable 2, respectively; and n is the number of data pairs used to calculate the cross-correlogram at each distance h.
The cross-correlogram defines the correlation existing between the values of Variable 1 and values of Variable 2 separated by the distance h. At zero distance, the cross-correlogram is equal to the Pearson correlation coefficient. At any h
0, the cross-correlogram values depend on the direction of h. In this study, cross-correlogram values were averaged over the two opposite directions.
For two positively correlated soil/crop yield variables, a typical cross-correlogram has its maximum value at zero distance (Fig. 1) . As the distance increases, the cross-correlogram value decreases until it becomes statistically insignificant. Further behavior of cross-correlogram depends on the specifics of the spatial patterns of the two variables and the strength of the spatial correlation between them. Cross-correlogram may reach a constant value at a certain distance that can be defined as the range of spatial cross-correlation between the variables, or it may keep decreasing to negative values and then become statistically significant again at larger distances but with opposite sign. For the negatively correlated variables, the cross-correlogram behavior follows a similar pattern but with opposite signs. Occasionally, cross-correlogram values that are insignificant at short distances may become significant at larger distances, but this is not typical. In this study, we will concentrate on cross-correlograms with a typical behavior.
|
The cross-correlogram parameter that conveniently combines all three components is the area under the cross-correlogram curve (for positively correlated variables) or above the cross-correlogram curve (for negatively correlated variables). In this study, we considered not the whole area under the cross-correlogram curve but rather only the area of significant correlation for distances ranging from zero to the distance at which the cross-variogram value becomes statistically insignificant.
The procedure for determining the area under the cross-correlogram curve consisted of three steps. First, the sample cross-correlogram was fitted with a mathematical equation. We used polynomial equations for all cross-correlogram fitting:
![]() | [2] |
(h) is the cross-correlogram at distance h and ß0 through ßk are polynomial coefficients for real-number polynomial orders k. A least-square fitting of the polynomial equations to the data was conducted using TableCurve 2D 5.0 software (SYSTAT Software, 2002). The order of the polynomial used for each specific sample cross-correlogram was determined by examining the polynomials with k ranging from 1 to 5 and selecting the polynomial that produced the highest r2 based on the differences between the equation predictions and the sample cross-correlogram values. In addition to r2, visual assessment of the polynomial equation fitting was performed in each case to assure a correct description of the sample cross-correlogram by the polynomial equation. In the case of comparable performance of several polynomials, the one with fewer parameters was selected. Although in some instances other types of equations were better suited to fit a particular cross-correlogram, for the sake of consistency and easy integration, only polynomial equations were used in this study. An example of an experimental cross-correlogram fitted with a polynomial equation is shown in Fig. 2
.
|
![]() | [3] |
The area under the curve P(h) was calculated for h ranging from 0 distance to a, where a is the distance at which the cross-correlogram value became statistically insignificant (P = 0.05), also called zone of influence (Gajem et al., 1981). In one-dimensional space and with a equal to
, Eq. [3] becomes identical to the integral scale as defined by Bakr et al. (1978) and Warrick et al. (1986).
In the third step, we quantified the significant portion of the area under the cross-correlogram curve. That was achieved by calculating minimum statistically significant cross-correlogram values based on the number of pairs, n, used in cross-correlogram calculation (P = 0.05). The minimum significant cross-correlogram values were fitted with a polynomial equation s(h) similar to Eq. [2]. The area under the s(h) curve of the minimum significant cross-correlogram values was obtained by integrating s(h) as:
![]() | [4] |
Finally, the parameter characterizing statistically significant spatial correlation between the variables, SC, was calculated as:
![]() | [5] |
A schematic representation of the cross-correlogram equation; the level of significance, a; and the areas of P(h) and S(h) are shown in Fig. 2.
The SC parameter is more informative than the regular correlation coefficient because it conveniently summarizes the information about the regular correlation and spatial correlations between the variables in a single characteristic number.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
An observed trend can be partially explained if we recall that in Illinois and Michigan, higher soil EC values are very often associated with sites that have greater clay content, greater water content, and poorer drainage. During wet springs, such sites are characterized by lower soil temperature and poor aeration. Crops in such sites are particularly prone to waterlogging. Hence, in wet springs, these sites provide poor conditions for plant growth compared with sites with lower EC, thus the negative correlations between crop yield and EC. However, in years with less spring precipitation, these sites supply more water for plant growth compared with the areas with coarse soil. In dry years, water availability is often a yield-limiting factor leading to a positive correlation between EC and crop yield data. During years of average spring precipitation, these sites influence yield to a much smaller extent, which results in insignificant correlation coefficients between yield and EC. It is interesting to note that early in the growing season (March), precipitation data can be used to estimate the relationship between EC and yield for the upcoming year and thus provide information as to whether or not soil EC data will be useful in building field management zones for a given year.
Analysis of covariance of the SC with maximum March precipitation as a covariate revealed significant differences between corn and soybean in terms of SC values (Table 3). The least-square estimates of the SC means at the maximum March precipitation of 3 cm are shown in Table 3. The value of 3 cm was selected as an example of a high maximum daily March precipitation value for the studied fields. Our results indicate that in years with high March precipitation, the cumulative spatial correlation between soybean yield and EC was stronger than that between corn yield and EC. That is, when the maximum daily March precipitation was high, the strength and the spatial component of the relationship between yield and EC for soybean were larger than those for corn. An equal slope model adequately fitted the data (Table 3). None of the other weather variables or soil properties were significant covariates for either correlation coefficient or SC.
|
Some of the variation in the range of the significant spatial correlations between crop yield and elevation, a, was related to the precipitation data as well as to the soil characteristics of the fields. The a values for elevation were correlated with the OM content of the studied fields. The regression equation relating a values for elevation (aEl) to maximum March precipitation (MaxMarchP) and soil OM explained approximately 42% of the variability of the elevation values, with levels of significance for OM and maximum daily March precipitation of 0.009 and 0.063, respectively:
![]() |
This observation can be explained if we consider the general patterns of the OM distribution across the studied Michigan and Illinois fields. As a rule, much higher OM contents are observed in depression areas. In this study, fields with relatively large areas occupied by such depressions also had higher average OM contents. As has been discussed, such areas are particularly influential on crop yield in both positive (in years with normal and dry weather) and negative (in wet years) ways. The presence of such areas in a field is a decisive factor in determining whether correlation between yield and elevation will be observed and at what distance such relationship will extend. It might seem that field topographic characteristics should also be able to reflect this relationship with the yield. However, the simple topographical characteristic, used in this study, such as elevation, was not able to quantify this feature. It is likely that an estimate of flow accumulation (Kravchenko and Bullock, 2000) would more accurately quantify this effect. However, flow accumulation is a representative characteristic of the studied field if a whole watershed area draining into the field's depression is included into the topographical analysis. In fields studied in this project, the elevation measurements were limited to the field itself, leaving no opportunity to quantify flow accumulation accurately. Hence, we believe that here, the OM content served as a proxy for those lower-elevation sites with substantial water flow accumulation.
| SUMMARY |
|---|
|
|
|---|
The utility of the proposed parameter lies in its potential application in development of zones for site-specific crop management. For example, if a soil property of management importance is strongly correlated with soil EC, we expect that EC will be useful for delineating the zones for managing that particular soil property. However, the size and location of the zones will depend on the spatial component of the correlation between the soil property and EC as reflected in SC values. If the SC parameter is relatively small, then it would be advisable that the size and shape of the management zones exactly follow the EC map. However, if the SC is large, then the size of effective management zones can be larger than what is suggested by the EC map. Indeed, large negative SC value would indicate that high values of the soil properties are observed not only directly in the areas with low EC but also in vicinity of such areas. Hence, specific zones would include not only the areas of the extremely high/low EC but also the areas of intermediate EC surrounding the extremes. This process should enable efficient optimization of management zones within a field without loss in accuracy.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
A. M. Tarquis, N. R. A. Bird, A. P. Whitmore, M. C. Cartagena, and Y. Pachepsky Multiscale Entropy-based Analysis of Soil Transect Data Vadose Zone J., May 27, 2008; 7(2): 563 - 569. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. S. Grigera, R. A. Drijber, K. M. Eskridge, and B. J. Wienhold Soil Microbial Biomass Relationships with Organic Matter Fractions in a Nebraska Corn Field Mapped using Apparent Electrical Conductivity Soil Sci. Soc. Am. J., August 3, 2006; 70(5): 1480 - 1488. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Crop Science | Vadose Zone Journal | |||
| Journal of Natural Resources and Life Sciences Education |
Soil Science Society of America Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||