|
|
||||||||
USDA-ARS, George E. Brown, Jr., Salinity Lab., 450 West Big Springs Rd., Riverside, CA 92507-4617
* Corresponding author (dcorwin{at}ussl.ars.usda.gov)
Received for publication April 23, 2001.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: EC, electrical conductivity ECa, apparent soil electrical conductivity ECe, electrical conductivity of the saturated soil paste extract ECw, electrical conductivity of soil water EM, electromagnetic induction EMavg, the geometric mean of the vertical and horizontal electromagnetic induction readings EMh, electromagnetic induction measurement in the horizontal coil-mode configuration EMv, electromagnetic induction measurement in the vertical coil-mode configuration GIS, geographical information system GPS, global positioning systems NPS, nonpoint source SP, saturation percentage TDR, time domain reflectometry
w, total volumetric water content
| INTRODUCTION |
|---|
|
|
|---|
From a global perspective, irrigated agriculture makes an essential contribution to the food needs of the world. While only 15% of the world's farmland is irrigated, roughly 35 to 40% of the total supply of food and fiber comes from irrigated agriculture (Rhoades and Loveday, 1990). However, vast areas of irrigated land are threatened by salinization. Although accurate worldwide data are not available, it is estimated that roughly half of all existing irrigation systems (totaling about 250 million ha) are affected by salinity and waterlogging (Rhoades and Loveday, 1990).
Salinity within irrigated soils clearly limits productivity in vast areas of the USA and other parts of the world. It is generally accepted that the extent of salt-affected soil is increasing. In spite of the fact that salinity buildup on irrigated lands is responsible for the declining resource base for agriculture, we do not know the exact extent to which soils in our country are salinized, the degree to which productivity is being reduced by salinity, the increasing or decreasing trend in soil salinity development, and the location of contributory sources of salt loading to ground and drainage waters. Suitable soil inventories do not exist and until recently, neither did practical techniques to monitor salinity or assess the impacts of changes in management on soil salinity and salt loading. A means of assessing soil salinity across the landscape is essential to the management of soil salinity. Because of the influence of soil salinity on crop productivity and the dynamic spatio-temporal nature of salinity, real-time measurement and monitoring of the spatial and temporal distribution of soil salinity is a crucial piece of information for precision agriculture applications on irrigated agricultural soils of the arid southwestern USA.
Precision agriculture (or site-specific agriculture) utilizes rapidly evolving electronic information technologies to modify land management in a site-specific manner as conditions change spatially and temporally (van Schilfgaarde, 1999). The intent of precision agriculture is to optimize crop production while minimizing detrimental environmental effects. First conceived in the mid-1980s, the technological pieces needed to bring precision agriculture into its own fell into place in the mid-1990s with the maturation of global positioning systems (GPS) and geographical information systems (GIS). As such, precision agriculture is a technologically driven system (van Schilfgaarde, 1999). The measurement of soil EC is among the technologies that have helped to bring precision agriculture from a concept to a potential tool for addressing the issue of agricultural sustainability.
The determination of total solute concentration (i.e., salinity) through the measurement of electrical conductance has been well established for decades (U.S. Salinity Lab. Staff, 1954). Soil EC measurement is a means of easily quantifying and monitoring soil salinity in irrigated agricultural areas of arid-zone soils. Soil salinity refers to the presence of the major dissolved inorganic solutes in the aqueous phase consisting of soluble and readily dissolvable salts in soil, including charged species (e.g., Na+, K+, Mg2+, Ca2+, Cl-, HCO-3, NO-3, SO2-4 and CO2-3), nonionic solutes, and ions that combine to form ion pairs. Soil salinity is quantified in terms of the total concentration of the soluble salts as measured by the EC of the solution in dS m-1 (U.S. Salinity Lab. Staff, 1954). For pure solutions, it is known that the EC of the pure solution
w is a function of the chemical composition as characterized by Eq. [1]:
![]() | [1] |
is the molar limiting ion conductivity (S m2 mol-1), M is the molar concentration (mol m-3), v is the absolute value of the ion charge, and i denotes the ion species in solution. Marion and Babcock (1976) were among the first to confirm the existence of a relationship between EC and molar concentrations of ions in the soil solution. The relationship of EC of soil water (ECw) to aggregate or individual ions in soil has also been confirmed by recent work of Kachanoski et al. (1992), Vanclooster et al. (1993), Wraith et al. (1993), Mallants et al. (1994), Ward et al. (1994), Heimovaara et al. (1995), Nissen et al. (1998), and Das et al. (1999). Operationally, soil EC is determined for an aqueous extract of a soil sample. Ideally, it would be desirable to determine the concentrations of individual solutes in soil water over the entire range of field water contents with a quick, easy measurement, but practical methods are not currently available to do so. Because of the time, labor, and cost of obtaining soil solution extracts, developments in soil salinity measurement over the past two decades have shifted to EC measurement of the bulk soil, referred to as the apparent soil electrical conductivity (ECa). The ECa measures conductance through not only the soil solution, but also through the solid soil particles and via exchangeable cations that exist at the solidliquid interface of clay minerals.
Apparent soil electrical conductivity has become one of the most reliable and frequently used measurements to characterize field variability for application to precision agriculture due to its ease of measurement and reliability (Rhoades et al., 1999b). The value of spatial measurements of ECa to precision agriculture is widely acknowledged, but ECa is still often misunderstood and misinterpreted. To help clarify misconceptions, a general overview of the application of ECa to precision agriculture is presented. It is the objective of this paper to describe the practical technology and methodology for real-time measurement of ECa, with particular emphasis on spatial ECa measurements applied to precision agriculture. The following areas are discussed: a brief history of the measurement of soil salinity with EC, the basic theories and principles of the soil EC measurement and what it actually measures, an overview of the measurement of soil salinity with various EC measurement techniques and equipment (specifically, electrical resistivity with the Wenner array and EM), examples of spatial EC surveys and their interpretation, applications, and value of spatial measurements of soil EC to precision agriculture, and current and future developments.
| BRIEF HISTORY OF THE MEASUREMENT OF SOIL SALINITY |
|---|
|
|
|---|
The first method, visual crop observation, is quick and economical, but it has the disadvantage that salinity development is detected after crop damage has occurred. For this reason, visual observation is the least desirable method for assessing soil salinity. However, on fields where long-term, field-scale patterns of soil salinity exist, visual observations from one year provide useful information for succeeding years.
The second method, electrical conductance of the soil solution extract, gives a quantitative measure of soil salinity but requires considerable resources for field sampling and laboratory analysis, plus the volume of soil measurement is very small and ill-suited to characterize and map the extreme variability of salinity at field scales and larger (Corwin, 2002a). Customarily, soil salinity has been defined in terms of EC of the saturated soil paste extract (ECe). This is because it is impractical for routine purposes to extract soil water from samples at typical water contents; consequently, soil solution extracts must be made at higher-than-normal water contents. Unfortunately, the partitioning of solutes over the three soil phases (i.e., gas, liquid, and solid) is influenced by the soil/water ratio at which the extract is made so that the ratio needs to be standardized to obtain results that can be applied and interpreted universally. Commonly used extract ratios aside from a saturated soil paste are 1:1 (EC1:1), 1:2 (EC1:2), and 1:5 (EC1:5) soil/water mixtures. However, soil salinity can also be determined from the measurement of ECw. Theoretically, ECw is the best index of soil salinity because this is the soil solution salinity actually experienced by the plant root. Nevertheless, ECw has not been widely used to express soil salinity for various reasons: (i) It varies over the irrigation cycle as the soil water content changes, so it is not single valued (Rhoades, 1978), and (ii) the methods for obtaining soil solution samples for routine ECw analysis are too labor, time, and cost intensive at typical field water contents to be practical for field-scale applications (Rhoades et al., 1999a). For disturbed soil samples, soil solution can be obtained in the laboratory by displacement, compaction, centrifugation, molecular adsorption, and vacuum or pressure extraction methods (Rhoades and Oster, 1986). For undisturbed soil samples, ECw can be determined either in the laboratory on a soil solution sample collected with a soil solution extractor (Fig. 1) or directly in the field using in situ, imbibing-type porous-matrix salinity sensors (Fig. 2) .
|
|
Electrical resistivity (e.g., Wenner array) and EM are both well suited for precision agriculture applications because their volumes of measurement are large, which reduces the influence of local-scale variability. However, electrical resistivity is an invasive technique that requires good contact between the soil and four electrodes inserted into the soil; consequently, it produces less reliable measurements in dry or stony soils than the noninvasive EM measurement. The EM has become the first choice for measuring soil salinity in a geospatial context because (i) measurements can be taken as quickly as one can move from one location to the next, (ii) the large volume of soil measured reduces local-scale variability, and (iii) measurements in relatively dry or stony soils are possible because no contact is necessary between the soil and the EM sensor (Hendrickx et al., 1992). Nevertheless, electrical resistivity has a flexibility that has proven advantageous for field application, i.e., the depth and volume of measurement can be easily changed by changing the spacing between the electrodes. Even though the depth of penetration and volume of measurement by EM can be varied by raising the instrument above the soil surface, it is a complicated matter of determining the subsequent depth and volume of measurement, whereas in the case of the electrical resistivity method, it is a simple calculation to estimate the depth of penetration and the volume of measurement once the interelectrode spacing is adjusted. Both electrical resistivity and EM measure the ECa. Electrical resistivity measures apparent resistivity. In homogeneous mediums, resistivity is the reciprocal of conductivity. Resistivity devices convert measurements of apparent resistivity into measurements of apparent conductivity.
Time domain reflectometry (TDR) was initially adapted for use in measuring water content. Later, Dalton et al. (1984) demonstrated the utility of TDR to also measure ECa. The measurement of ECa with TDR is based on attenuation of applied signal voltage as it traverses the medium of interest with the relative magnitude of energy loss related to ECa (Wraith, 2002). The advantages of TDR for measuring ECa include (i) a relatively noninvasive nature because there is only minor interference with soil processes, (ii) an ability to measure both soil water content and ECa, (iii) an ability to detect small changes in ECa under representative soil conditions, (iv) the capability of obtaining continuous unattended measurements, and (v) no need for a calibration of soil water content measurements in many cases (Wraith, 2002).
| PRINCIPLES OF ECa MEASUREMENT |
|---|
|
|
|---|
Apparent soil electrical conduction in sufficiently moist soils is primarily via salts contained in soil water occupying the larger pores; consequently, measurement of EC of bulk soil is closely related to soil salinity (Rhoades et al., 1999b). However, there is also a contribution by the solid phase to ECa in moist soils primarily via the exchangeable cations associated with clay minerals (Rhoades et al., 1999b). A third pathway exists through soil particles in direct and continuous contact with one another (Rhoades et al., 1999b). These three pathways of current flow contribute to the ECa (Fig. 3) .
|
![]() | [2] |
ws and
wc are the volumetric soil water contents in the soil water pathway (cm3 cm-3) and in the continuous liquid pathway (cm3 cm-3), respectively;
ss and
sc are the volumetric contents of the surface-conductance (cm3 cm-3) and indurated (cm3 cm-3) solid phases of the soil, respectively; ECws and ECwc are the specific ECs of the soil water pathway (dS m-1) and continuous liquid pathway (dS m-1); and ECss and ECsc are the ECs of the surface-conductance (dS m-1) and indurated (dS m-1) solid phases, respectively. Equation [2] was reformulated by Rhoades et al. (1989) into Eq. [3]:
![]() | [3] |
w =
ws +
wc = total volumetric water content (cm3 cm-3) and
sc · ECsc was assumed to be negligible. The following simplifying approximations are also known:
![]() | [4] |
![]() | [5] |
![]() | [6] |
![]() | [7] |
![]() | [8] |
b is the bulk density (Mg m-3), SP is the saturation percentage, ECw is average EC of the soil water assuming equilibrium (i.e., ECw = ECws = ECwc), and ECe is the EC of the saturation extract (dS m-1). By measuring soil ECe, SP, PW, and
b and using Eq. [3] through [8], the ECa can be estimated. Very simply, Eq. [3] through [8] indicate that ECa is a function of the soil physical and chemical properties of (i) soil salinity, (ii) SP, (iii) water content, and (iv) bulk density. The SP and bulk density are both closely associated with clay content.
Another factor influencing ECa is temperature. Electrolytic conductivity increases at a rate of approximately 1.9% per degree centigrade increase in temperature. Customarily, EC is expressed at a reference temperature of 25°C for purposes of comparison. The EC measured at a particular temperature t, ECt, can be adjusted to a reference EC at 25°C, EC25, using the below relation from Handbook 60 (U.S. Salinity Lab. Staff, 1954):
![]() | [9] |
![]() | [10] |
Traditionally, ECe has been the standard measure of salinity used in all salt-tolerance plant studies. As a result, a relation between ECa and ECe was needed to relate ECa back to ECe, which is in turn related to crop yield. Over the past two decades, research has been directed at developing reliable and efficient conversion techniques from ECa back to ECe (Williams and Baker, 1982; McNeill, 1980, 1986; McKenzie et al., 1989; Rhoades, 1992, 1996; Rhoades and Corwin, 1990; Rhoades et al., 1989, 1990, 1991, 1999a; Slavich, 1990; Cook and Walker, 1992; Diaz and Herrero, 1992; Yates et al., 1993; Lesch et al., 1992, 1995a, 1995b, 1998).
| METHODS AND EQUIPMENT FOR GEOREFERENCED MEASUREMENT OF ECa |
|---|
|
|
|---|
Electrical Resistivity
Electrical resistivity methods introduce an electrical current into the soil through current electrodes at the soil surface, and the difference in current flow potential is measured at potential electrodes that are placed in the vicinity of the current flow. These methods were developed in the second decade of the 1900s by Conrad Schlumberger in France and Frank Wenner in the United States for the evaluation of ground electrical resistivity (Burger, 1992; Telford et al., 1990).
The electrode configuration is referred to as a Wenner array when four electrodes are equidistantly spaced in a straight line at the soil surface, with the two outer electrodes serving as the current, or transmission, electrodes and the two inner electrodes serving as the potential, or receiving, electrodes (Fig. 4)
. The depth of penetration of the electrical current and the volume of measurement depend on the interelectrode spacing. The larger the spacing is, the deeper the measurement and larger the volume of measurement. The resistivity,
, measured with the Wenner array is (Burger, 1992):
![]() | [11] |
). Because the ECa is the inverse of R, then Eq. [11] becomes:
![]() | [12] |
|
a3. Other electrode configurations are frequently used, and the equations for these configurations are discussed by Burger (1992), Dobrin (1960), and Telford et al. (1990).
The basic equipment for measuring ECa with the Wenner array technique (Fig. 5)
includes an electrical current source, a resistance meter, four metal electrodes, connecting wire, measuring tape, and a soil thermometer (Rhoades and Halvorson, 1977). The current source can be a hand-cranked generator or battery powered, which should measure from 0.1 to 1000
. The electrodes can be made from any noncorrosive conductive metal (e.g., stainless steel, copper, brass). The connecting wire should be flexible, well-insulated, multistranded, 12- to 18-gauge wire. Detailed instructions concerning electrode construction and information concerning the resistivity meters can be found in Rhoades and Halvorson (1977).
|
|
![]() | [13] |
A mobilized, tractor-mounted version of the fixed-electrode array (Fig. 7) has been developed that georeferences the ECa measurement with a GPS (Carter et al., 1993; Rhoades, 1992, 1993). The fixed-electrode array is a particularly appealing approach for field applications because of the relative ease of measurement and the large volume of soil that can be measured. However, because of the large volume of measurement and the complex localized spatial variability of pore geometry, water content, and salinity, which influence the measurement of ECa, some error is introduced, making it difficult to calibrate ECa, as measured with this equipment, to ECw or ECe as determined for soil samples taken at the sites of the ECa measurements. The error is a consequence of the difference in volumes of measurement between the ECa measurement and the soil core sample used to determine the associated ECw or ECe. The calibration problem is also of concern for EM. The mobile, fixed-electrode array equipment is particularly well suited for collecting detailed maps of the spatial variability of average root-zone soil EC at field scales and larger.
|
Electromagnetic Induction
Apparent soil electrical conductivity can be measured remotely with EM. An EM transmitter coil located at one end of the instrument induces circular eddy-current loops in the soil (Fig. 8)
. The magnitude of these loops is directly proportional to the EC of the soil in the vicinity of that loop. Each current loop generates a secondary electromagnetic field that is proportional to the value of the current flowing within the loop. A fraction of the secondary induced electromagnetic field from each loop is intercepted by the receiver coil of the instrument, and the sum of these signals is amplified and formed into an output voltage, which is related to a depth-weighted bulk soil EC, ECa. The receiver coil measures amplitude and phase of the secondary magnetic field. The amplitude and phase of the secondary field will differ from those of the primary field as a result of soil properties (e.g., clay content, water content, and salinity), spacing of the coils and their orientation, frequency, and distance from the soil surface (Hendrickx and Kachanoski, 2002).
|
|
![]() | [14] |
10-7 H m-1); and s is the intercoil spacing (m).
Mobile EM equipment developed at the Salinity Laboratory (Carter et al., 1993; Rhoades, 1993) is available for the appraisal of soil salinity and other soil properties (e.g., water content and clay content) using an EM-38. The drawback of this early mobile EM equipment was that it produced a point ECa measurement for the vertical (EMv) and horizontal (EMh) coil configurations rather than a continuous-stream reading of ECa as with the mobile fixed-array equipment; consequently, all ECa maps were grid maps. The drawback relates to the inability of the EM-38 meter to measure and record data simultaneously in each dipole orientation. The measurement of soil salinity with EM relies on both horizontal and vertical dipole measurements at each observation point to qualitatively evaluate the salinity distribution profile. However, data can be recorded continuously with the EM-38 meter if surveys are conducted in one dipole orientation. This is commonly done in the Midwest where other towed units that have been developed. In these instances, the concern is not with mapping salinity, but usually with mapping clay, water content, or both. Recently, the mobile EM equipment developed at the Salinity Laboratory has been modified by the addition of a dual-dipole EM-38 unit (Fig. 10)
. The dual-dipole EM-38 meter simultaneously records data in both dipole orientations at time intervals of just a few seconds between readings. A close-up of the sled that houses the dual-dipole EM-38 is provided in Fig. 11
. The mobile EM equipment is suited for the detailed mapping of ECa and correlated soil properties at specified depth intervals through the root zone. The advantages of the mobile dual-dipole EM equipment over the mobile fixed-array resistivity equipment are that (i) the EM technique is noninvasive, so it can be used in dry or stony soils that would not be amenable to the invasive technique of the fixed-array approach due to the need for good electrodesoil contact; (ii) the EM technique provides information for the vertical profiling of ECa rather than just a single integrated ECa although the Veris unit does provide two depths; and (iii) stream tubes can be delineated (see Modeling Soil Salinity in a Geospatial Context section). A disadvantage of the EM approach would be that the ECa is a depth-weighted value that is nonlinear. This depth-weighted nonlinearity is shown in Fig. 12
, which illustrates the cumulative relative contributions of all soil EC, R(z), for a homogeneously conductive material below a normalized depth of z based on Eq. [15] and [16] from McNeill (1980) for vertical and horizontal dipoles, respectively:
![]() | [15] |
![]() | [16] |
|
|
|
| STANDARD OPERATING PROCEDURE FOR A FIELD-SCALE ECa SURVEY FOR PRECISION AGRICULTURE APPLICATIONS |
|---|
|
|
|---|
Mobile ECa Survey and Soil Sample Design
An initial intensive ECa survey with either mobile electrical resistivity or EM equipment is conducted and used to establish soil core sampling locations needed for calibration and understanding the dominant soil properties influencing ECa. Depending on the level of detail desired, from 100 to several thousand spatial measurements of ECa are taken generally in regularly spaced traverses across the field of interest. The use of mobile EM equipment has one slight advantage over the use of mobile resistivity equipment, which is the ability to take measurements on dry and stony soils. Figures 13 and 14
show representative ECa surveys using mobile Wenner array and mobile EM-38 equipment, respectively.
|
|
Stochastic and/or Deterministic Calibration of ECa to ECe
Soil salinity, as conventionally expressed in terms of ECe, can be inferred from ECa by two approaches: a deterministic and a stochastic approach (Rhoades et al., 1999b). The preferred approach will vary with the size of the area to be assessed, availability of equipment, and the specific objectives. In the deterministic approach, either theoretically or empirically determined models are used to convert ECa into salinity (ECe). Deterministic models are static; i.e., all model parameters are considered known, and no ECe data need to be determined. For example, Eq. [3], developed by Rhoades et al. (1989), is a deterministic approach. The deterministic approach is the preferred approach to use when significant, localized variations in soil type exist in the field. This approach typically requires knowledge of additional soil properties (e.g., soil water content, SP, bulk density, temperature, etc.). In the stochastic approach, statistical modeling techniques such as spatial regression or co-kriging are used to directly predict the soil salinity from ECa survey data. In this approach, the models are dynamic; i.e., the model parameters are estimated using soil sample data collected during the survey. This calibration is developed by acquiring soil salinity data (or other soil property data, such SP, texture, bulk density, etc.) from a small percentage of the sensor measurement sites and estimating an appropriate stochastic prediction model for each depth increment using the paired soil sample and ECa measurement data. Then, using the remaining sensor data in conjunction with the established model, the soil salinity levels (or other similarly calibrated properties) are predicted at all of the remaining nonsampled measurement locations. The stochastic and deterministic calibration approaches are described in detail by Lesch et al. (1995a)(1995b, 2000) and incorporated into the ESAP software (Lesch et al., 2000).
Typically, when an ECa survey is performed, the ECa data are used to estimate the values of one or more soil variables influencing ECa. For example, estimates of spatial soil salinity pattern, soil texture, or water-holding capacity may be desired. To facilitate this estimation process, the soil variables are calibrated to the ECa survey information with a stochastic approach. As mentioned in the previous section (Mobile ECa Survey and Soil Sample Design), ESAP optimizes the sample design by selecting sites that optimize estimation of a calibration or prediction model (i.e., an equation that can be used to predict the values of soil variables from the ECa survey data). The prediction model is an ordinary regression model. To use an ordinary regression model to predict the values of a spatial data set, ESAP addresses issues of spatial correlation and collinearity.
Apparent soil electrical conductivity can be converted into estimated soil salinity (i.e., ECe) using a deterministic conversion algorithm. First, raw ECa data are converted into depth-specific ECa. Next, these ECa data are converted into estimated soil salinity data using Eq. [3] through [8], originally developed by Rhoades et al. (1989).
Determination of the Dominant Soil Properties Influencing the ECa Measurement
The fact that ECa is a function of several soil properties (i.e., soil salinity, texture, and water content) is often overlooked in the application of ECa measurements obtained with EM and electrical resistivity to precision agriculture. Precision agriculture studies relating ECa to crop yield have met with inconsistent results due to the fact that a combination of factors influence ECa measurements to varying degrees across units of management, thereby confounding interpretation. These factors include soil salinity, water content, SP, and bulk density. In areas of saline soils, salinity dominates the ECa measurements, and interpretations are often more straightforward. To use spatial measurements of ECa in a precision agriculture context, it is necessary to understand what factors are most significantly influencing the ECa measurements within the field of study. There are two approaches for determining the predominant factors influencing ECa measurement: (i) simple statistical correlation and (ii) wavelet analysis. Even though wavelet analysis is a powerful tool for determining the dominant complex interrelated factors influencing ECa measurement, it requires data on a regular grid or equal-spaced transect. Grid or equal-spaced transect sampling schemes are not as practical for determining spatial distributions of soil salinity from ECa measurements as the stochastic statistical approach developed by Lesch et al. (1995a)(1995b, 2000).
Table 1 is a compilation of correlation data for six field study sites where ECa surveys were performed for the purpose of salinity appraisal. The table shows the variation in the influence of various soil properties on ECa for different field locations. In all cases, the surveys were performed as discussed in the above standard operating procedure. An intensive ECa survey was performed, followed by selection of 6 to 20 sites for soil core sampling. An analysis was performed on the soil cores for various physicochemical properties (e.g., SP, salinity, and water content). A calibration and the correlations in Table 1 were determined using the ECa survey and soil sample data. The predicted correlation column corresponds to the correlation between the predicted ECa using Eq. [3] through [8] and the specified soil property (salinity, texture, and
w) while the observed correlation column corresponds to the correlation between the measured ECa and the specified soil property. The fact that the predicted correlations are in most cases the same as the observed correlations shows that the model by Rhoades et al. (1989) is fairly robust. Only in the case of the Coachella sorghum [Sorghum bicolor (L.) Moench] field does the model appear to show some weakness. Even more convincing evidence of the robustness of the Rhoades et al. (1989) model is the high correlation between the predicted ECa, as reflected by ln(EC'a), and the observed ECa, as reflected by ln(EMavg), in Table 2
. However, the observed correlations in Table 1 are the most significant because this column shows what soil properties are influencing the ECa reading most.
|
|
Coachella Valley Wheat (Triticum aestivum L.) Field
This is an example of a survey where the salinity represented by ln(ECe), the soil texture reflected by the SP, and the
w correlate with the EM data, which is represented as ln(EMavg), where EMavg is the geometric mean of the vertical and horizontal EM readings (i.e.,
). In this field, Eq. [3] through [8] correctly predict the correlations and correctly predict that salinity will correlate best with the averaged EM data. The reason why all three soil properties correlate well with the EM data is because all three soil properties are highly correlated with each other (Table 2).
Coachella Valley Sorghum Field
This field is an example of where only salinity correlates well with the EM data. Equations [3] through [8] correctly predict this condition even though the calculated correlation between the ln(EMavg) and the ln(EC'a), where ECa is calculated from Eq. [3] through [8], is less than ideal (r = 0.85). Note from Table 2 that neither the soil texture nor
w correlate much at all with salinity, with r = -0.10 and r = 0.28, respectively. It is this lack of correlation with salinity and because the texture and water content exhibit minimal sample variation (i.e., sample range for SP is 51.061.1%; sample range for
w is 0.330.41 cm3 cm-3) that they correlate so poorly with the EM data, with r = -0.20 and r = 0.25, respectively (Table 1).
Broadview Water District (Quarter Sections 16-2 and 16-3)
These combined quarter sections display large variability in soil texture (SP ranges from 33.285.0%) and water content (
w ranges from 0.210.39 cm3 cm-3) sample data, with relatively minimal salinity variation (80% of the samples fell below the mean value of 3.65 dS m-1). Equations [3] through [8] correctly predict that both the texture and EM and water content and EM correlation levels will be higher than the salinity and EM correlation level (Table 1). Equations [3] through [8] tend to slightly overestimate the influence of salinity and underestimate the texture and water content influences (Table 1).
Fresno Cotton (Gossypium hirsutum L.) Field
Equations [3] through [8] correctly predict a negative correlation between EM data and water content and also correctly predict the order of the degree of correlation between the EM data and each soil property although the influence of soil water content again appears to be underestimated (Table 1).
Coachella ValleyKohl Ranch Field
This field displays excellent agreement between predicted and observed EMsoil property correlations (Table 1). Salinity correlates very well, water content correlates fairly well, and soil texture exhibits weak negative correlation, indicating that the dominant soil properties influencing the EM reading are salinity and water content. In addition, two secondary properties, sodium adsorption ratio and B, were measured. The fact that these correlated quite well with the EM data suggests the close association of these properties with salinity in this particular field because the EM reading does not directly measure sodium adsorption ratio or B.
Broadview Water District (Quarter Section 10-2)
Equations [3] through [8] do an excellent job of predicting both the sign and magnitude of each EMsoil property correlation (Table 1). The dominant soil property influencing the EM reading is salinity.
From this discussion of six typical ECa surveys for agricultural soils in the arid Southwest, what is known about the interrelationship of soil properties influencing the ECa measurement? First, Eq. [3] through [8] describe how various soil properties (salinity, texture, and water content) interact to influence the ECa signal data and highlight the fact that this interaction can be quite complex. However, there appears to be consistent, high positive correlation between the EC'a (i.e., the calculated ECa using Eq. [3][8]) and EM data (see Table 2), implying that Eq. [3] represents a reasonable model. Therefore, we can get a good idea of how well EM survey data will correlate with a specific soil property by determining how well this same soil property correlates with EC'a data.
Second, it is clear that the innercorrelation structure of the various primary soil properties (ECe, SP, and
w) determines how well each property ultimately correlates with the ECa signal data. However, the variability of each soil property also influences the final correlation estimates because increased variability in any given soil property directly translates into increased variation in the ECa data. Obviously, one may encounter many diverse types of innercorrelation structures and/or different degrees of specific soil property variation, as shown in Tables 1 and 2. Thus, the ultimate correlation between the ECa signal data and any specific soil property may be quite different from field to field. For example, this effect is clearly evident in the ln(EMavg) and SP correlation estimates shown in Table 1 where the observed estimates range from -0.33 to 0.84.
Finally, with respect to ECe data, the best scenario for the prediction of salinity from ECa signal data occurs when the (ECe, SP, and
w) cross-correlation estimates are all positive and high (i.e., near 1) and/or the SP and
w variation is minimal. However, in all of the scenarios shown here, Eq. [3] can be successfully used to explain why certain soil properties correlate well with ECa signal data in some survey scenarios but not in others.
| MODELING SOIL SALINITY IN A GEOSPATIAL CONTEXT |
|---|
|
|
|---|
Soil and ground water quality affected by salinity depend on spatially distributed properties that influence salt transport. The phenomenon of salt transport through the vadose zone (i.e., the zone extending from the soil surface to the ground water table) is affected by the temporal variation in irrigation water quality and the spatial variability of plant water uptake and soil chemical and physical properties. Coupling a GIS to a salt transport model potentially offers a means of dealing with the complex spatial heterogeneity of soils, which influences the intricate biological, chemical, and physical processes of transient-state salt transport in the vadose zone (Corwin, 1996).
Modeling the movement and accumulation of salinity is a spatial problem well suited for the integration of a salt transport model with GIS (Corwin, 1996). A GIS serves as a spatial database to organize, manipulate, and display the complex spatial data used by a deterministic model to describe the regional-scale distribution of soil salinity and salt loading to ground water. Coupling of the spatial datahandling capabilities of a GIS with a one-dimensional salt transport model offers the advantage of utilizing the full information content of spatially distributed data to analyze solute movement on a field scale in three dimensions (Corwin, 1996). As a visualization and analysis tool, GIS is capable of manipulating both spatially referenced input and output parameters of the model.
The basic components of a NPS pollutant model (NPS pollutants, such as salinity and pesticides, are spread over a wide area, as opposed to point source pollutants, which are located at a specific site or point, such as a toxic spill) are model, GIS, and data (Corwin et al., 1997). These components depend on the advanced information technologies of a GPS, GIS, geostatistics, remote sensing, solute transport modeling, neural networks, transfer functions, fuzzy logic, hierarchical theory, and uncertainty analysis (Corwin et al., 1999a).
Published work by Corwin and his colleagues (Corwin and Rhoades, 1988; Corwin et al., 1989, 1995, 1996, 1999b) reflects the current trend in the application of GIS-based models to simulate the transport of NPS pollutants in the vadose zone at field, basin, and regional scales (Corwin, 1996; Corwin et al., 1997, 1999a). These models are specifically designed to account for the spatio-temporal variability of the properties and conditions influencing soil salinity accumulation and transport. Two different GIS-based approaches are described for predicting the areal distribution of salinity. The first approach couples a regression model of salinity development to a GIS of soil salinity development factors (i.e., permeability, leaching fraction, and ground water EC) for the WelltonMohawk Irrigation District near Yuma, AZ, over the study period 19681973 (Corwin and Rhoades, 1988; Corwin et al., 1989). The regression model predicts the composite salinity of the root zone (i.e., top 60 cm.). Areas of low, medium, and high salinization potential are delineated for the entire 44 000-ha irrigation district. Measured salinity data verified that 86% of the predicted salinity categories was accurately predicted. The second approach loosely coupled the one-dimensional, transient-state solute transport model, TETrans, to the GIS ARC/INFO (Corwin et al., 1999b). Slightly less than 2400 ha of the Broadview Water District located on the west side of central California's San Joaquin Valley was used as the test site to evaluate the integrated GIStransport model over the study period 19911996 (Corwin et al., 1999b).
The unique aspect of the Corwin et al. (1999b) approach to landscape-scale modeling of a NPS pollutant (i.e., salinity) in the vadose zone is the delineation of stream tubes from ECa measurements taken on a grid with the mobile EM-38 equipment developed by Carter et al. (1993) and Rhoades (1993). Stream tubes are noninteractive volumes of soil whose physicochemical properties influencing solute transport are relatively homogenous so that solute transport within the column of soil defined by the stream tube can be simulated with a one-dimensional solute transport model. Corwin et al. (1999b) first proposed the use of an intensive EM survey measuring ECa as a means of delineating stream tubes. For field sites where ECa is closely correlated with soil salinity, stream tubes can be delineated based on EMh and EMv measurements of ECa. From the geometric mean of EMh and EMv (i.e.,
), quantiles can be defined. The ratio of EMh to EMv (i.e., EMh/EMv) is determined, and within each quantile, the points are selected where the low and high EMh/EMv ratios exist. These points serve as the centroids of the Thiessen polygons delineating the stream tubes throughout the area of study.
A GIS is used to define the spatial location of the stream tubes with their associated solute transport property attributes. TETrans uses the GIS as a spatial database from which to draw its input data. Simulations are presented over a 5-yr period, 19911996. Display maps show spatial distributions of soil salinity profiles to a depth of 1.2 m, irrigation efficiencies, drainage amounts, and salt loading to ground water over the 2396-ha study area. These maps provide a potential precision agriculture tool for making irrigation management decisions to minimize the environmental impact of salinity on soil and ground water. The first approach is best suited for areas where steady-state conditions are approximated while the second approach can be used under transient-state conditions.
| SUMMARY |
|---|
|
|
|---|
Assessing the impact of salinity at regional and localized scales is a key component to achieving sustainable agriculture. Assessment involves the determination of change in salinity over time, which can be measured in real time or predicted with a model. Real-time measurements reflect activities of the past, whereas model predictions are glimpses into the future based on a simplified set of assumptions. Both means of assessment are valuable; however, the advantage of prediction is that it can be used to alter the occurrence of detrimental conditions before they develop. Forecasting information from model simulations is used in decision-making strategies designed to sustain agriculture. This information permits an alteration in the management strategy before the development of levels of soil salinity that detrimentally impact either agricultural productivity of the soil or quality of the ground water.
The ability to locate the sources of soil salinity within irrigated landscapes and to model the migration of salt through the vadose zone to obtain an estimate of salt loading to drainage and ground water is an essential tool in precision agriculture to combat degradation of our soil and water resources. This information is valuable for selecting alternate crops or alternate irrigation management practices to maintain crop productivity while minimizing the environmental impacts of salinity.
The need for a means of measuring soil salinity in the root zone in a quick, reliable, and cost-effective manner resulted in the development of mobile electrical resistivity and EM techniques to measure ECa. However, the measurement of ECa is complicated by the influence of several soil properties aside from soil salinity, including soil texture, temperature, and water content. For this reason, it is necessary to determine the dominant soil properties influencing the ECa measurement at each field of interest to interpret what information is being conveyed by a map of ECa. From a precision agriculture perspective, when maps of ECa are properly understood, they (i) provide a graphic inventory of the scope of the soil salinity problem, (ii) provide useful spatial information concerning soil texture and water content, (iii) provide information for crop selection, (iv) identify potential areas in need of improved irrigation and drainage management, and (v) provide a means of monitoring spatio-temporal changes in soil properties that potentially influence crop production.
Because precision agriculture is an outgrowth of technological developments such as the ECa measurement, the future of precision agriculture rests on a complete understanding of these technologies. The fruition of the precision agriculture application of ECa maps will likely come from future plant indicator approaches where combinations of data inputs including airborne spectral imagery, ECa measurements with mobile EM, and plant and soil sampling based on response-surface sample design strategies are manipulated, organized, and displayed with GIS, image analysis, and spatial statistical analysis to create maps of soil salinity, soil texture and available water.
| ACKNOWLEDGMENTS |
|---|
| NOTES |
|---|
|
|
|---|
2 Geonics, Limited, Mississauga, ON, Canada. Product identification is provided solely for the benefit of the reader and does not imply the endorsement of the USDA. ![]()
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
D. L. Shaner, R. Khosla, M. K. Brodahl, G. W. Buchleiter, and H. J. Farahani How Well Does Zone Sampling Based on Soil Electrical Conductivity Maps Represent Soil Variability? Agron. J., September 8, 2008; 100(5): 1472 - 1480. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. L. Corwin, S. M. Lesch, J. D. Oster, and S. R. Kaffka Short-Term Sustainability of Drainage Water Reuse: Spatio-Temporal Impacts on Soil Chemical Properties J. Environ. Qual., September 2, 2008; 37(5_Supplement): S-8 - S-24. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Segal, S. A. Bradford, P. Shouse, N. Lazarovitch, and D. Corwin Integration of Hard and Soft Data to Characterize Field-Scale Hydraulic Properties for Flow and Transport Studies Vadose Zone J., August 1, 2008; 7(3): 878 - 889. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Robinson, C. S. Campbell, J. W. Hopmans, B. K. Hornbuckle, S. B. Jones, R. Knight, F. Ogden, J. Selker, and O. Wendroth Soil Moisture Measurement for Ecological and Hydrological Watershed-Scale Observatories: A Review Vadose Zone J., February 25, 2008; 7(1): 358 - 389. [Abstract] [Full Text] [PDF] |
||||
![]() |
S. A. Staggenborg, M. Carignano, and L. Haag Predicting Soil pH and Buffer pH In Situ with a Real-Time Sensor Agron. J., May 11, 2007; 99(3): 854 - 861. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. B. Lobell, J. I. Ortiz-Monasterio, F. C. Gurrola, and L. Valenzuela Identification of Saline Soils with Multiyear Remote Sensing of Crop Yields Soil Sci. Soc. Am. J., April 5, 2007; 71(3): 777 - 783. [Abstract] [Full Text] [PDF] |
||||
![]() |
N. E. Derby, F. X. M. Casey, and D. W. Franzen Comparison of Nitrogen Management Zone Delineation Methods for Corn Grain Yield Agron. J., February 6, 2007; 99(2): 405 - 414. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Abdu, D.A. Robinson, and S.B. Jones Comparing Bulk Soil Electrical Conductivity Determination Using the DUALEM-1S and EM38-DD Electromagnetic Induction Instruments Soil Sci. Soc. Am. J., January 1, 2007; 71(1): 189 - 196. [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] |
||||
![]() |
J. A. Guretzky, K. J. Moore, C. L. Burras, and E. C. Brummer Distribution of Legumes along Gradients of Slope and Soil Electrical Conductivity in Pastures Agron. J., March 1, 2004; 96(2): 547 - 555. [Abstract] [Full Text] [PDF] |
||||
![]() |
D. A. Robinson, I. Lebron, S. M. Lesch, and P. Shouse Minimizing Drift in Electrical Conductivity Measurements in High Temperature Environments using the EM-38 Soil Sci. Soc. Am. J., March 1, 2004; 68(2): 339 - 345. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. R. Schepers, J. F. Shanahan, M. A. Liebig, J. S. Schepers, S. H. Johnson, and A. Luchiari Jr. Appropriateness of Management Zones for Characterizing Spatial Variability of Soil Properties and Irrigated Corn Yields across Years Agron. J., January 1, 2004; 96(1): 195 - 203. [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 | |||