|
|
||||||||
Dep. of Biological and Irrigation Engineering Utah State Univ. Logan, UT 84322
Dep. of Plants, Soils and Biometeorology Utah State Univ. Logan, UT 84322-4820
Currently at: Dep. of Geophysics, Stanford Univ. 397 Panama Mall Stanford, CA 94305
Dep. of Plants, Soils and Biometeorology Utah State Univ. Logan, UT 84322-4820
* Corresponding author (hiruyabdu{at}cc.usu.edu).
| ABSTRACT |
|---|
|
|
|---|
5%) at higher ECa (>100 mS m1). The difference in response of the instruments can be ascribed to the temperature-dependent change in soil ECa due to a 20°C diurnal temperature variation in addition to instrumental drift caused by temperature effects on the processing circuitry. The EM38-DD's real-time display and internal powering are its advantages while the lower priced DUALEM-1S does not need instrument calibration and can store data internally | INTRODUCTION |
|---|
|
|
|---|
The EMI-based ECa measurements with the Geonics EM38 (Geonics Inc., Mississauga, ON) have been used by researchers attempting to infer different properties and characterize a wide range of processes at the field scale for a host of different applications (Hendrickx and Kachanoski, 2002). Doolittle et al. (1994) estimated claypan depth via ECa measurements in central Missouri soils via direct calibration between ECa and topsoil depth above the claypan. Jaynes et al. (1995) correlated EMI-derived ECa measurements to herbicide partition coefficients to predict herbicide application leaching potentials in specific areas of an EMI-surveyed field. Anderson-Cook et al. (2002) exploited the positive correlation between ECa and previous year crop yields to statistically classify four different soils. Sudduth et al. (2001) developed a procedure using ECa measurements to infer topsoil depth in claypan soils. Corwin and Lesch (2003, 2005) outlined standard operating procedures for ECa surveys applied to precision agriculture, specifically surveys that calibrate ECa to the electrical conductivity of saturation extract samples (ECe) for use in salinity studies, and discussed several different applications of field-scale ECa maps. Corwin and Lesch (2003) used ECa measurements to infer ECe to assess salinity effects by comparing ECe measurements within fields to crop yields and chemical analyses. Taylor et al. (2003) used the recently developed DUALEM-2 (DUALEM Inc., Milton, ON) to identify the locations and depths of septic system failure.
The reliability of data collected using EMI instruments depends on the thermal stability of the instrument, while the ECa measurement averaging down the soil profile depends on the configuration of the instrument coils (Wait 1951; 1955). Researchers using the EM38 for field mapping have observed progressive instrument drift affecting the ECa measurements during mapping days (Sudduth et al., 2001; Robinson et al., 2004). Sudduth et al. (2001) investigated accuracy issues in the collection of soil ECa data and recommended a calibration to document and adjust for instrument drift. The EM38 data were collected at four 50-m calibration transects during the day, whereby a maximum of >10% ECa deviation was observed. Robinson et al. (2004), on the other hand, were able to register a drift of 20% in the hot southwest USA, where the EM38 instrument panel temperature reached 60°C.
The depth-weighted response of an EMI instrument depends on coil orientation with respect to the half space and spacing of the coils (Wait, 1955, 1951). Rhoades et al. (1999) conducted several studies to determine the ECadepth distributions for the measurement of salinity profiles and for analyzing saline seeps. The ECadepth relation is calculated by successively raising the EMI instrument and measuring the respective contribution of each soil interval to ECa (Corwin and Rhoades, 1982). Corwin and Rhoades (1990) produced empirical calibration and statistical analysis equations to calculate the ECadepth relationships for different soil types using the EM-38.
The purpose of this study was to characterize the DUALEM-1S instrument and compare it with the Geonics EM38-DD. Characterization of the ECa averaging with depth was conducted on low- and high-ECa soils, and compared with theoretical models. Side-by-side experiments were conducted to compare instrument thermal stability on low- and high-ECa soils.
| THEORETICAL CONSIDERATIONS |
|---|
|
|
|---|
![]() | [1] |
x 107 H m1), s is the interdipole spacing (m), Hs is the secondary magnetic field at the receiver coil (H m1), and Hp is the primary magnetic field at the transmitter coil (H m1). The amplitude and phase of the secondary magnetic field measured by the receiver coil differ from the primary field owing to soil properties and transmitterreceiver spacing and orientation (i.e., horizontally or vertically oriented) with respect to the soil surface (Hendrickx and Kachanoski, 2002).
The transmitters and receivers of EMI devices consist of wound coils that can be treated as magnetic dipoles since the separation between transmitters and receivers is more than several coil diameters. The magnetic field from each dipole penetrates the earth and the vertical dipole has a greater depth of penetration than the horizontal dipole because the vertical field couples more effectively with material down in the earth than does the horizontal field. Figure 1 shows that the magnetic fields of the vertical dipoles entering the soil surface are more dense than the horizontal dipoles.
|
In Fig. 1, the three combinations are illustrated. Horizontalhorizontal (HH): in this combination, both the transmitter and receiver dipoles are oriented parallel to the earth's surface. The bottom unit (horizontal dipole mode) in the EM38-DD uses this combination. Vertical vertical (VV): in this combination the dipoles are oriented perpendicular to the earth's surface. The top unit (vertical dipole mode) in the EM38-DD and the horizontal co-planar mode in the DUALEM-1S use this combination. Verticalhorizontal (VH): the transmitter dipole is vertical, while the receiver dipole is horizontal and its axis intersects the transmitter. The perpendicular mode in the DUALEM-1S uses this combination.
Relative Response
The governing equations for the EMI relative response,
, of the three different orientations are (McNeill, 1980; Wait, 1962)
![]() | [2] |
![]() | [3] |
![]() | [4] |
Figure 2A shows the relative sensitivity of the three coil orientations relative to an increase in depth. Both the HH and VH orientations are sensitive at the surface, and rapidly lose their sensitivity with an increase in depth. The VV orientation is insensitive at the ground surface but the sensitivity rapidly increases with depth, peaking at 0.4 m.
|
![]() | [5] |
Equations for the cumulative response of different orientations from a depth (z) to infinity have been given by McNeill (1980) and Wait (1962). We have adapted the equations so that they give cumulative responses from the surface to a given depth (z) of each orientation (Fig. 2B):
![]() | [6] |
![]() | [7] |
![]() | [8] |
Since the cumulative response is exponential, we need to define a depth beyond which the orientation response is relatively insensitivethe depth of exploration (DOE). For our study, we have defined the DOE to be the depth from which 70% of the cumulative response comes. The DOE for the VV orientation is 1.5 m, while the HH and VH orientations achieve the same response at a DOE of 0.75 and 0.5 m, respectively (Fig. 2B).
In a layered earth, ECa is calculated by summing the conductivity and depth-weighted contribution of each layer. In a system with three distinct layers overlying uniform earth, the bulk electrical conductivity is given by (McNeill and Bosnar, 1999)
![]() | [9] |
1,2,3 are the conductivities of each corresponding layer,
4 is the conductivity of the uniform earth underlying the three layers, and R(Z1,2,3) are the cumulative responses at the bottom of each respective layer. | MATERIALS AND METHODS |
|---|
|
|
|---|
|
Study Sites
Greenville Farm (Millville Series)
The Millville series located at the Utah Agricultural Experiment Station's Greenville Farm receives a mean annual precipitation of 422 mm and the mean annual temperature is 8.6°C. The site has a xeric soil moisture regime and a mesic soil temperature regime. The area surrounding the Millville soil pedon is used for irrigated crops and a distinguishable plow layer can be observed in the A horizon. The pedon contains <1% rock fragments and the texture (silt loam) is uniform with depth. The pH of 8.2 is due to the highly disseminated CaCO3. The parent material of the Millville soil is a fine-textured alluvium (silts and very fine sand) due to lower energy distal fan and overbank flood deposits. The soil is classified as a coarse-silty, carbonatic, mesic Typic Haploxeroll. Soil samples were taken every 0.3 m down to a depth of 1.5 m and were analyzed for water content using oven drying; ECe was measured using the saturation paste extraction method (Soil Survey Staff, 2004; Fig. 3A
).
|
Measurement Response vs. Depth
To study the effect of depth on measurement response, the instruments were lifted from the ground surface to conduct depth sounding. The instruments were raised using a polyvinyl chloride pipe as a guiding support, keeping the instrument parallel to the ground. Five sets of measurements for each instrument were taken at the two study sites by lifting the instruments in increments of 0.15 m up to 1.8 m above the ground. The measurements were obtained directly from the EM38-DD display, while an Allegro CX with HGIS software was used to collect the measurements from the DUALEM-1S.
Measurement Response vs. Temperature
To study the effects of temperature, EM readings from the two instruments were recorded throughout the day at two locations with high and low ECa values. On 1 July 2005, the instruments were placed at the Greenville Farm and on 8 July the instruments were positioned at the Cache Junction Farm. At both locations, the instruments were separated by a distance of 4 m and the measuring instruments were placed midway between them. A CR10 datalogger (Campbell Scientific, Logan, UT) was used to record the readings from the six thermocouples placed on the instruments (a thermocouple each by the receiver dipoles and panel of the EMI instruments) and the surrounding environment (one 0.15 m belowground and another 0.3 m aboveground in the air). The EM38-DD data was acquired using a handheld geographic information system (HGIS, StarPal Inc., Fort Collins, CO) program inside an Allegro CX handheld field computer (Juniper Systems, Logan, UT). The DUALEM-1S data was recorded internally and later downloaded to a computer.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
The results for the Greenville Farm are presented in Fig. 3B and 3C. The broken lines show the fit for the uniform earth and demonstrate a poor fit with the data for both instruments. The application of the three-layer inversion improves the fit for the DUALEM-1S data but is still poor for the EM38 data. The RMSE for the DUALEM-1S over a uniform earth was 0.30 for the VVDLM and 0.27 for the VHDLM orientations, respectively. This improved to a RMSE of 0.09 and 0.05 for VVDLM and VHDLM for the three-layer model. The result of the inversion is presented as the dotted line in Fig. 3A, and, although providing a good fit, it fails to capture the dominant higher conductivity layer over the two lower conductivity layers.
The EM38-DD showed the largest divergence from the models, VVEM38 having a RMSE of 0.79 and HHEM38 having a RMSE of 1.2 for the uniform earth model. Interestingly, the use of the three-layer model didn't lead to a major improvement, with RMSEs of 0.54 and 1.00 for VVEM38 and HHEM38, respectively. The difference between the measured and modeled data for the EM38-DD is due to the difficulty calibrating the EM38-DD at low conductivity values, resulting in poor quality measurements. To achieve a calibration with VVEM38 reading twice HHEM38, the EC of the horizontal had to be raised; this sets a threshold below which EC cannot be measured, resulting in poor data at low conductivities for this soil.
The measurements made in the Cache Junction soil are presented in Fig. 3E and 3F. In this higher conductivity soil, the uniform earth model correlated well with the data from both instruments. The fitting to the DUALEM-1S data gave a RMSE of 1.7 for VVDLM and 2.2 for VHDLM (Fig. 3E). The EM38-DD was better, giving a RMSE of 0.88 for the VVEM38 orientation and 1.4 for the HHEM38 orientation (Fig. 3F). The use of the three-layer model improved the RMSE for the DUALEM-1S to 0.19 for VVDLM and 0.16 for VHDLM and marginally for the EM38-DD to 0.80 for VVEM38 and 1.4 for HHEM38. The results of the inversion at Cache Junction for the DUALEM-1S are presented in Fig. 3D. This time the inversion did a better job of picking out the high and low conductivity layers. The use of these models demonstrates the different abilities of the instruments to measure under different conditions, bringing out the important point that with a small range of ECa response (020 mS m1), the DUALEM-1S does much better than the EM38-DD due to its internal, automatic calibration. The application of these simple inversion models shows that the results leave much to be desired. Inversion of these data could be very useful in vadose zone research for determining depth to conductive layers (salts, water tables, or clay layers), and should form the basis of future research with these instruments.
Measurement Response vs. Temperature
Both instruments were placed on bare soil while ECa and temperature data were collected from 10:30 to 18:30 h. At the low conductivity site (Greenville Farm), the air temperature ranged from 24 to 40°C, while the soil temperature climbed from 18 to 32°C at 0.15 m below the soil surface. The instrumental temperature variations of the DUALEM-1S and the EM38-DD are presented in Fig. 4A
and 4B, respectively. The DUALEM-1S casing reached a maximum temperature of around 40°C while the maximum temperature for VVEM38 panel was >50°C. Figures 4C and 4D show the EMI responses of the coil orientations of the DUALEM-1S and EM38-DD, respectively. The mean response of the VVDLM was 8.05 mS m1 with a SD of 0.87 and the VHDLM had a mean response of 10.4 mS m1 with a SD of 0.35, whereas the mean response of the VVEM38 was 19.9 mS m1 with a SD of 1.7 and the HHEM38 had a mean response of 19.6 mS m1 with a SD of 1.4.
|
It is evident from the statistics that a larger EMI response deviation from the initial reading is observed when the instruments are measuring low ECa values. At the low conductivity site, the EM38-DD had a maximum difference of 23% (6.3 mS m1) for the VVEM38 and 22% (8.3 mS m1) for HHEM38. The maximum difference for the DUALEM was 13% (1.8 mS m1) for VHDLM and a high 42% (4.8 mS m1) for VVDLM. The percentage difference in the EMI readings during the day at the high conductivity site was much smaller than at the low conductivity site. Maximum differences of 5.2% (6.4 mS m1) and 6.6% (8.0 mS m1) were observed for VVEM38 and HHEM38, respectively, while maximum differences of 5.0% (5.7 mS m1) and a low 1.2% (1.1 mS m1) were recorded for VVDLM and VHDLM, respectively.
The difference in response can be ascribed to two things: (i) the change in temperature of the soil during the day, which would tend to increase the ground conductivity, and (ii) instrument drift caused by the inability of the processing circuitry to fully compensate for instrument heating. To differentiate between these two competing factors, we modeled ECa to determine what the response of the instruments should be in the soil. As inputs for the model, we used water content measured from soil samples, porosity, and soil temperature measured at 0.15-m depth. We recognize that by using the temperature measured at 0.15 m, the values would be an upper bound to the actual ECa sensed by the instruments; however, we would expect both instruments to show a trend similar to the modeled ECa.
An extended version of Archie's Law (Friedman, 2005), accounting for unsaturated conditions, was used to model the expected ground conductivity response that the EMI instruments would measure. We adopted Friedman's (2005) simplification, which uses a unit value as the exponent of porosity (n), thus obtaining an equation with only one empirical fitting parameter (d):
![]() | [10] |
v is volumetric water content, and according to Corwin and Lesch (2005):
![]() | [11] |
The modeling was conducted such that each orientation would measure 70% of the cumulative response. The VVEM38 and VVDLM achieved this by measuring down to a depth of 1.5 m while the VHDLM and HHEM38 attained such a response by measuring down to 0.5 and 0.75 m, respectively. The parameter d was chosen to fit the first point of the model to the first data point at the beginning of the experiment. The model, derived experimentally using coarse-textured soils, has limitations in clay soils but demonstrates the expected upward trend in ECa as temperature increases. Table 2 lists the parameters used for modeling the responses of the different coil orientations at the two sites.
|
General Observations
Some of the strong features of the EM38-DD are: a real-time LCD display, a built-in handle, portable internal powering, and years of applied research experience; while a complicated instrumental calibration procedure, exposed control knobs, and a commonly overheating black panel are drawbacks. The DUALEM-1S has avoided some of the problems associated with the EM38-DD by not having any control knobs, having a yellow casing to minimize radiation absorption, and by incorporating an automatic instrument calibration. This automated instrument calibration is a distinct advantage for users working with low conductivity soils, for instance, when used for texture mapping.
Even though the DUALEM-1S can store 50000 records in its internal memory, it lacks a built-in display unit, making an external logging or display unit necessary. The DUALEM-1S also requires an external power source and does not come with a handle for manual measurements. This means the instrument is less suited to single measurements and more suited to applications such as mobile mapping where continuous measurements are made. At the time of testing, the price of the DUALEM-1S was approximately two-thirds the cost of the EM38-DD. Both instruments can be connected to various logging programs running on handheld field computers or laptops for recording georeferenced EMI response while field mapping.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
Received for publication December 6, 2005.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
L. Cockx, M. Van Meirvenne, U. W. A. Vitharana, L. P. C. Verbeke, D. Simpson, T. Saey, and F. M. B. Van Coillie Extracting Topsoil Information from EM38DD Sensor Data using a Neural Network Approach Soil Sci. Soc. Am. J., October 21, 2009; 73(6): 2051 - 2058. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Saey, D. Simpson, H. Vermeersch, L. Cockx, and M. Van Meirvenne Comparing the EM38DD and DUALEM-21S Sensors for Depth-to-Clay Mapping Soil Sci. Soc. Am. J., January 21, 2009; 73(1): 7 - 12. [Abstract] [Full Text] [PDF] |
||||
![]() |
O. R. Harvey and C. L. S. Morgan Predicting Regional-Scale Soil Variability using a Single Calibrated Apparent Soil Electrical Conductivity Model Soil Sci. Soc. Am. J., January 21, 2009; 73(1): 164 - 169. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Robinson, H. Abdu, S. B. Jones, M. Seyfried, I. Lebron, and R. Knight Eco-Geophysical Imaging of Watershed-Scale Soil Patterns Links with Plant Community Spatial Patterns Vadose Zone J., October 29, 2008; 7(4): 1132 - 1138. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Journal of Natural Resources and Life Sciences Education |
Vadose Zone Journal | ||||
| Journal of Plant Registrations | Journal of Environmental Quality |
The Plant Genome | |||