|
|
||||||||
a USDA-ARS, Cropping Syst. and Water Qual. Res. Unit, 269 Agric. Eng. Bldg., Univ. of Missouri, Columbia, MO 65211
b Dep. of Crop Sci., Univ. of Illinois, 1102 S. Goodwin Ave., Urbana, IL 61801
c Dep. of Agron., Univ. of Missouri, 214 Waters Hall, Columbia, MO 65211
* Corresponding author (SudduthK{at}missouri.edu)
Received for publication July 10, 2001.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: CV, coefficient of variation DGPS, differential global positioning system EC, electrical conductivity ECa, apparent soil electrical conductivity ECa-sh, shallow (030 cm) apparent soil electrical conductivity measured by Veris 3100 ECa-dp, deep (0100 cm) apparent soil electrical conductivity measured by Veris 3100 ECa-em, vertical-mode apparent soil electrical conductivity measured by Geonics EM38 EM, electromagnetic induction GPS, global positioning system TD, topsoil depth
| INTRODUCTION |
|---|
|
|
|---|
Measurements of ECa can be used to provide indirect measures of the soil properties listed above if the contributions of the other soil properties affecting the ECa measurement are known or can be estimated. If the ECa changes due to one soil property are much larger than those attributable to other factors, then ECa can be calibrated as a direct measurement of that dominant factor. Lesch et al. (1995a)(1995b) used this direct-calibration approach to quantify variations in soil salinity within a field where water content, bulk density, and other soil properties were "reasonably homogeneous." Research in Missouri has established direct, within-field calibrations between ECa and the depth of topsoil above a subsoil claypan horizon (Doolittle et al., 1994; Sudduth et al., 1995, 2001; Kitchen et al., 1999).
Mapped ECa measurements have been found to be related to a number of soil properties of interest in precision agriculture, including soil water content (Sheets and Hendrickx, 1995), clay content (Williams and Hoey, 1987), CEC, and exchangeable Ca and Mg (McBride et al., 1990). Because ECa integrates texture and moisture availability, two soil characteristics that affect productivity, it can help to interpret spatial grain yield variations, at least in certain soils (e.g., Sudduth et al., 1995; Jaynes et al., 1993; Kitchen et al., 1999). Other uses of ECa in precision agriculture have included refining the boundaries of soil map units (Fenton and Lauterbach, 1999), interpreting within-field corn rootworm (Diabrotica barberi Smith and Lawrence) distributions (Ellsbury et al., 1999), and creating subfield management zones (Fraisse et al., 2001).
Two types of portable, within-field ECa sensors have been used in agriculturean electrode-based sensor requiring direct contact with the soil and a noncontact electromagnetic induction (EM) sensor. The earliest sensors were of the contact type and included four electrodes inserted into the soil, coupled with an electric current source and resistance meter. Hand-carried four-electrode sensors were initially used in salinity surveys (Rhoades, 1993), and later versions were tractor-mounted for mobile, georeferenced measurement of ECa. The electrode-based sensing concept formed the basis of a commercial product, the Veris 3100 (Veris Technol., Salina, KS). This mobile system (Fig. 1) uses six rolling coulters for electrodes and simultaneously generates shallow (ECa-sh; nominally 030 cm) and deep (ECa-dp; 0100 cm) measurements of ECa (Lund et al., 1999). It includes all necessary components except for the tow vehicle and global positioning system (GPS) receiver and requires no user calibration.
|
|
Using a wheeled cart pulled by an all-terrain vehicle (Fig. 2), an EM38 system is adaptable to a wide variety of data collection conditions. This lightweight system requires little power and makes it possible to collect data under wet or soft soil conditions. Also, it is possible to collect data after a crop has been planted in 76-cm rows, up until the time that the crop is 15 to 20 cm tall. The Veris 3100 is much heavier and requires a tractor or truck to pull it through the field, limiting its use to firmer soil conditions and unplanted fields. The newer Veris 2000XA, which only has four coulters and one measurement depth, can be pulled by a large all-terrain vehicle and can collect data between planted 76-cm crop rows.
Commercial operators are using ECa sensing systems to provide soil variability information to producers. Although many or most of these are coulter-based sensors, the vast majority of research information has been obtained with EM-based sensors. As more use is made of ECa sensing in precision agriculture, it will be important to compare the data obtained with each type of system and to understand how these data are related to soil properties. This study was undertaken to compare ECa data collected on Missouri and Illinois fields with the noncontact Geonics EM38 and the coulter-based Veris 3100 and to relate those data to measured soil properties. Objectives were to (i) interpret differences in ECa sensor data in relation to response curves of the sensors, (ii) document the relationship of ECa data to soil properties, and (iii) investigate the improvement, if any, obtained by combining multiple ECa variables for estimating soil properties.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
Soils of the Illinois fields include the Varna series (fine, illitic, mesic Oxyaquic Argiudolls), Drummer series (fine silty, mixed, superactive, mesic Typic Endoaquolls), and Chenoa series (fine, illitic, mesic Aquic Argiudolls). Surface textures include silt loam and silty clay loam. Drainage classes represented at the Illinois fields range from poorly drained to well drained (Kravchenko et al., 2001).
Apparent Soil Electrical Conductivity Sensors and Response Curves
The EM sensor used in this research (Geonics EM38) has a spacing of 1 m between the transmitting coil located at one end of the instrument and the receiver coil at the other end. Calibration controls and a digital readout of ECa in millisiemens per meter (mS m-1) are included, and an analog data output allows data to be recorded on a data logger or computer. The EM38 was operated in the vertical dipole mode, providing an effective measurement depth of approximately 1.5 m. The vertical-mode ECa measurement from the EM38 by Geonics EM38 (designated by ECa-em in this study) is averaged over a lateral area approximately equal to the measurement depth. The instrument response to soil conductivity varies as a nonlinear function of depth (McNeill, 1992). Sensitivity in the vertical mode is highest at about 0.4 m below the instrument (Fig. 3)
. The ECa measurement is determined by the soil conductivity with depth, as weighted by this instrument response function (McNeill, 1992). The EM38 was combined with a data acquisition computer and differential GPS (DGPS) system for mobile data collection (Fig. 2), as described by Sudduth et al. (2001).
|
If the response curves of Fig. 3 are integrated with respect to depth, differences in the soil volumes measured by the different sensors are readily apparent (Fig. 4) . With ECa-sh, 90% of the response is obtained from the soil above the 30 cm depth. With ECa-dp, 90% of the response is obtained from the soil above the 100 cm depth. With ECa-em, 90% of the response is obtained above 5 m depth while 70% of the response is obtained above about 1.5 m. The curves of Fig. 4 are based on equations that assume a homogeneously conductive soil volume. Actual responses will vary somewhat due to ECa differences between soil layers, with a high-conductivity surface layer reducing the depth of response (Barker, 1989).
|
The Veris 3100 and Geonics EM38 were operated in tandem, taking measurements on transects spaced approximately 10 m apart. Data were recorded on a 1-s interval, corresponding to a 4- to 6-m data spacing. Between 4400 and 11 000 individual ECa measurements were obtained for each field. Data obtained by DGPS were associated with each sensor reading to provide positional information with an accuracy of 1.5 m or better.
Using our previously reported approach (Sudduth et al., 2001), a calibration transect was established in each field to monitor instrument drift during the survey. Data were collected on this transect at least every hour, and raw ECa readings were adjusted based on any change in calibration transect data. As expected, the direct ECasensing approach of the Veris system was much less (<50%) prone to instrument drift than the EM38. We believe that drift compensation would not be a necessary component of Veris ECa surveys although it should be done for EM38 surveys.
Within each field, between 12 and 21 sampling sites were selected to cover the range of ECa values present. These sites were chosen by a soil scientist familiar with the soils in the particular field with the additional goal of including samples from all of the landscape positions and soil map units present. One 4.0-cm-diam. core that was 120 cm long was obtained at each site using a hydraulic soil-coring machine. Cores were examined within the field by a skilled soil scientist and pedogenic horizons identified. Cores were segmented by horizon for laboratory analysis. Soil moisture was determined gravimetrically.
Additionally, samples for each horizon were analyzed at the University of Missouri Soil Characterization Laboratory using methods described by the National Soil Survey Center Staff (1996). Data were obtained for the following properties: sand, silt, and clay fractions (pipette method); CEC (base + Al method); organic C; and saturated paste EC.
Data Analysis
To allow comparison between ECa sensors, a combined data set was created for each field. Each Veris data point was combined with the nearest EM38 data point based on GPS coordinates. If a match was not found within a 2-m radius, that point was removed from the data set. Additional data sets were created to compare across sampling dates on the Missouri fields. Because measurement transect locations were not identical between 1997 and 1999, it was necessary to increase the search radius for these data sets to 3 m. Pearson correlation coefficients (r) were calculated between the various ECa sensors and measurement dates.
In this study, soil property data were obtained by horizon, rather than on an even depth increment. To facilitate comparison across calibration points, a depth-weighted mean was calculated for each soil property at each calibration point. To provide a measure of the variability in each soil property with depth, a depth-weighted coefficient of variation (CV) was also calculated. To account for the fact that the response of each ECa sensor is not constant with depth, three additional sets of data were created by weighting each soil property profile by the sensor response curve (Fig. 3).
Analysis of the relationship between ECa and soil properties was performed for each data source (ECa-em, ECa-sh, and ECa-dp) and profile-weighted soil property, using the 1999 calibration point data. These data were examined for spatial autocorrelation by calculating the Moran coefficient as suggested by Long (1996). No significant autocorrelation was detected in any ECa data. Only 15% of the soil property data sets showed significant (P
0.05) spatial autocorrelation. With this general lack of significant spatial autocorrelation, likely caused by the small number (1219) and spatial dispersion of the calibration points in each field, we conducted a nonspatial analysis between ECa and soil properties. Pearson correlation coefficients were calculated between ECa and soil properties (moisture, clay, silt, sand, organic C, CEC, and saturated paste EC). Regressions were performed to estimate soil properties from (i) each individual ECa measurement, (ii) both Veris 3100 ECa measurements, and (iii) all three ECa measurements. Only parameters statistically significant (P
0.05) were retained in the final regression equations.
Our previous work (Doolittle et al., 1994; Kitchen et al., 1999; Sudduth et al., 2001) established the utility of ECa-em data for estimating TD on claypan soils. In this study, we compared the accuracy of TD estimation by ECa-em and ECa-dp for the Missouri claypan soil fields. Estimations based on ECa-sh were not included because 90% of the theoretical ECa-sh response is within 30 cm of the surface. Therefore, ECa-sh data are unable to estimate TDs greater than approximately 30 cm while the TD on these fields exceeded 100 cm in places. Topsoil depth (depth to the first B horizon) data obtained at calibration points in fields F1 and GV were used to develop linear regression equations for estimating TD as a function of the inverse of ECa (ECa-1). Only those calibration points where TD was <100 cm were used because 90% of the theoretical ECa-dp response is within 100 cm of the surface.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
|
|
When data from both fields of a state were combined, correlations were similar, and better in some cases, than correlations calculated within individual fields. When data were combined for all fields, correlations between the two deeper ECa readings did not decrease, but correlations of ECa-sh to the other ECa readings were much lower (Table 3). These results indicate that although the shallow (ECa-sh) and deep (ECa-dp or ECa-em) ECa data were strongly related within a field or soil association (the fields from each state had similar soils and were located near each other), their relationship was not consistent across different soil associations. Thus, the shallow (ECa-sh) and deeper (ECa-dp and ECa-em) sensors provide unique information, and data from one cannot be inferred from data obtained with the other. However, because the two deeper ECa readings were highly correlated both within and across fields, it appears that little additional information would be gained on these soils by collecting both EM38 and Veris data.
Relationship of Apparent Soil Electrical Conductivity to Measured Soil Properties
A statistical summary of profile-average soil property data measured for the calibration points in each field is shown in Table 4. Analysis of variance indicated that profile-average clay and paste EC were significantly higher for the Missouri fields while sand, organic C, and CEC were significantly higher (P
0.05) for the Illinois fields. Profile CVs of clay, silt, CEC, and paste EC were significantly higher for the Missouri fields while profile CVs of organic C were significantly higher (P
0.05) for the Illinois fields. These higher CVs showed that the claypan soils of the Missouri fields were more layered in terms of the soil properties affecting ECa. To further investigate this layering, mean A-horizon and first B-horizon clay and CEC were calculated for the calibration points on each field. For the Illinois fields, mean clay and CEC for the first B-horizon were within 2% of the means for the A horizons. For the Missouri fields, mean clay was 215% greater and mean CEC 185% greater for the first B horizon compared with the A horizons. This significant layering, combined with differences in response functions (Fig. 3) for the different sensors, explains the nonlinear relationship between data from the different sensors seen on the Missouri fields (Fig. 6)
. The similarity of clay and CEC levels between the A horizon and B horizon for the Illinois fields helps to explain the linear relationship between ECa data obtained from the different sensors on those fields (Fig. 6).
|
|
0.05) correlation coefficients between ECa and profile-weighted soil properties for each field are shown in Table 5. Correlations of ECa with sensor-weighted clay content and sensor-weighted CEC were generally highest and most persistent across all fields and ECa data types. This higher correlation with sensor-weighted data supports our hypothesis that transformation of soil property data by weighting with the sensor response function is an appropriate way to help account for curvilinearity in the functional relationship. Other soil properties that exhibited a significant correlation in most cases were clay, silt, and CEC of the upper soil horizon. Some properties, such as profile-average organic C and CEC were significant on the Missouri fields but not on the Illinois fields. Significant correlations with soil moisture, sand content, and paste EC were observed less frequently.
|
0.05), so regressions were performed for three data sets: (i) Missouri data, (ii) Illinois data, and (iii) all data. Table 6 shows the regression statistics for each analysis. Regressions for some soil properties were more predictive for Missouri fields while others were more predictive for Illinois fields. The most accurate estimates were obtained for clay, silt, and CEC. Estimates of soil moisture, organic C, and paste EC obtained by regression on a single ECa variable were of relatively low accuracy.
|
A second series of regression analyses included multiple ECa data sources for estimating the same soil properties listed above. Stepwise quadratic (plus interaction) analyses included (i) both Veris data setsECa-sh and ECa-dpand (ii) all three ECa data sets (Table 6). In general, this approach provided little, if any, improvement over single-factor ECa regressions for top-layer clay, silt, and CEC, reinforcing our observation that ECa-sh data were a reasonable estimator of these properties. Estimates of single-state, profile-average clay and silt were generally improved by including both Veris ECa data sets and were further improved somewhat by including the ECa-em data. For multistate analyses, estimates were improved when all three ECa variables were allowed to enter the regression but were not improved by including just Veris data.
Estimates of paste EC and profile-average soil moisture were of low accuracy and were not improved by including additional ECa variables. Estimates of organic C for Illinois fields were improved by including additional ECa terms while estimates for Missouri fields were not. Estimates of top-layer soil moisture, available only for Missouri fields, also improved when additional ECa terms were included. For both single ECa and multiple ECa regressions, better estimates of soil properties were obtained within a single state than across both states. For best results, site-specific (or soil-specific) equations relating soil properties to ECa should be used.
Regression equations for estimating claypan TD as a function of ECa yielded standard errors from 6 to 16 cm (Table 7). Comparisons between ECa-em and ECa-dp data were variable between fields. For field F1, ECa-dp data were more predictive of TD while ECa-em data were more predictive of TD on field GV. Variations in the accuracy of TD estimations between years could be explained at least partially by the fact that different calibration points were used between 1997 and 1999. For F1, the 1997 calibration points exhibited a reasonably uniform distribution in TD across the range from 0 to 100 cm. However, in 1999, the calibration-point TDs were clustered between 20 and 50 cm. For GV, a more uniform distribution of calibration points was obtained in 1999. These results point out the importance of properly selecting calibration points for relating ECa data to soil physical properties. One way to remove the subjectivity from this process was proposed by Lesch et al. (1995b), who described an algorithmic approach to the selection of optimized locations for calibrating ECa measurements.
|
| SUMMARY AND CONCLUSIONS |
|---|
|
|
|---|
Within a single field and measurement date, ECa-em and ECa-dp were most highly correlated (r = 0.740.88). Lowest correlations were between ECa-em and ECa-sh. Differences were attributed to differences between the depth-weighted response functions for the three data types, coupled with differences in the degree of soil profile layering between sites. The claypan soils of the Missouri fields exhibited more variation with depth in clay content and CEC, two primary drivers of ECa. Because of this, differences between ECa data types were more pronounced on the Missouri fields.
Correlations of ECa with clay content and CEC were generally highest and most persistent across fields and ECa data types. Significant correlations with silt content and, on the Missouri fields, with organic C were also common. Significant correlations with sand content, paste EC, or soil moisture content were observed only about 10% of the time.
Top-layer clay, silt, and CEC could be estimated reasonably well as a function of a single ECa variable, usually ECa-sh. Profile-average clay and silt were estimated less accurately and often required combination of multiple ECa variables for best results. Organic C and top-layer soil moisture estimates were of variable accuracy while paste EC and profile-average soil moisture and CEC were not accurately estimated with regressions including single or multiple ECa variables.
This study showed that, while qualitatively similar, ECa data obtained with different commercial sensors were quantitatively different. Differences between the two data sources designed to give a deep-profile measurement (ECa-dp and ECa-em) were more pronounced in more layered soils. Both of these measurements were related to profile CEC and clay content while the ECa-sh measurement was strongly related to CEC and clay content in the upper soil horizon(s). With these differences, the selection of an ECa sensing system for a particular application should be based on both practical implementation issues and the intended use of the data.
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
R.N. Lerch, E.J. Sadler, N.R. Kitchen, K.A. Sudduth, R.J. Kremer, D.B. Myers, C. Baffaut, S.H. Anderson, and C.-H. Lin Overview of the Mark Twain Lake/Salt River Basin Conservation Effects Assessment Project Journal of Soil and Water Conservation, November 1, 2008; 63(6): 345 - 359. [Abstract] [PDF] |
||||
![]() |
J. Sawchik and A. P. Mallarino Evaluation of Zone Soil Sampling Approaches for Phosphorus and Potassium Based on Corn and Soybean Response to Fertilization Agron. J., November 6, 2007; 99(6): 1564 - 1578. [Abstract] [Full Text] [PDF] |
||||
![]() |
P. Jiang, S. H. Anderson, N. R. Kitchen, K. A. Sudduth, and E. J. Sadler Estimating Plant-Available Water Capacity for Claypan Landscapes Using Apparent Electrical Conductivity Soil Sci. Soc. Am. J., October 29, 2007; 71(6): 1902 - 1908. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Myers, N. R. Kitchen, K. A. Sudduth, R. E. Sharp, and R. J. Miles Soybean Root Distribution Related to Claypan Soil Properties and Apparent Soil Electrical Conductivity Crop Sci., July 30, 2007; 47(4): 1498 - 1509. [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] |
||||
![]() |
W. K. Jung, N. R. Kitchen, K. A. Sudduth, R. J. Kremer, and P. P. Motavalli Relationship of Apparent Soil Electrical Conductivity to Claypan Soil Properties Soil Sci. Soc. Am. J., May 6, 2005; 69(3): 883 - 892. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. R. Kitchen, S. T. Drummond, E. D. Lund, K. A. Sudduth, and G. W. Buchleiter Soil Electrical Conductivity and Topography Related to Yield for Three Contrasting Soil-Crop Systems Agron. J., May 1, 2003; 95(3): 483 - 495. [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 | |||