|
|
||||||||
a Dep. of Crop Sci., North Carolina State Univ., Raleigh, NC 27695
b Dep. of Soil Sci., North Carolina State Univ., Raleigh, NC 27695
c Plant Sci. Dep., South Dakota State Univ., Brookings, SD 57007
* Corresponding author (Ron_Heiniger{at}ncsu.edu)
Received for publication June 21, 2001.
| ABSTRACT |
|---|
|
|
|---|
Abbreviations: CEC, cation exchange capacity DGPS, differential global positioning system ECa, apparent soil electrical conductivity ECs, the electrical conductivity of the soil particles ECwc, the electrical conductivity of the mobile soil solution associated with large, continuous pores ECws, the electrical conductivity of the soil solution associated with discontinuous pores HM, humic matter PC, principal component
s, the volumetric content of soil particles
w, the total volumetric content of water in the soil
ws, the volumetric soil water content of the small, discontinuous pores
| INTRODUCTION |
|---|
|
|
|---|
Apparent soil electrical conductivity consists of two components: (i) the contribution of the solid soil particles primarily associated with exchangeable cations and (ii) the contribution of the soil solution (Nadler and Frenkel, 1980; Shainberg et al., 1980). Rhoades et al. (1989) described ECa using a two-pathway model, with one path for conductance via the discontinuous soilliquid interface and the other path via continuous, large water-filled pores through the soil solution. They found that apparently, soil structure does not provide enough direct particle-to-particle contact to form a continuous pathway for current flow. Their model identifies five major factors that influence ECa. These include: the electrical conductivity of the soil particles (ECs); the electrical conductivity of the soil solution associated with discontinuous pores (ECws); the electrical conductivity of the mobile soil solution associated with large, continuous pores (ECwc); the volumetric soil water content of the small, discontinuous pores (
ws); the total volumetric content of water in the soil (
w); and the volumetric content of soil particles (
s). Several studies have found that the major factors influencing ECa are
w (Rhoades et al., 1976; Nadler, 1982);
s per unit of soil, influenced primarily by the texture and bulk density of the soil (Rhoades and Corwin, 1990); the CEC of the soil, which influences ECs (Shainberg et al., 1980); and the amount of dissolved salts in the soil solution, which influences ECws and ECwc (Malicki and Walczak, 1999).
The most important factor influencing ECa is the total volumetric water content of the soil (Rhoades et al., 1976). Changes in volumetric soil water content tend to mask differences in the other factors influencing ECa (Nadler, 1982). When the volumetric soil water content is high (near field saturation), the primary path for ECa is through the larger continuous pores, and ECwc is a major influence on electrical conductivity (Rhoades et al., 1989). However, when the volumetric soil water content is low (near wilting point), the primary path is through the soil particlediscontinuous soil pore pathway. In this situation, the ECs and ECws are the major influences on ECa.
Another major influencing factor is soil salinity (Malicki and Walczak, 1999). At high Na concentrations, the conductivity of the soil solution (ECws and ECwc) is greater than that of the bulk soil (Shainberg et al., 1980). However, in nonsaturated, nonsaline soils, the conductivity of the bulk soil (ECs) is greater due to the contribution of adsorbed ions. If we assume that salinity is not a major factor in most productive agricultural fields, it follows that in nonsaturated conditions, changes in soil nutrient levels will most likely influence ECa by changing ECs through differences in the type and number of cations held by the soil particles. Therefore, cations commonly associated with binding sites on the soil particles, such as Ca, Mg, or K, could influence ECa by changing ECs. However, the common assumption is that, in most field situations, the influence that changing levels of soil cations have on ECs is minor compared with the influence associated with changes in soil bulk density and texture (Lund et al., 1999).
Differences in soil nutrient content could also potentially influence ECws and ECwc through changes in the conductivity of the soil solution. Under saturated conditions, changes in nutrient levels could influence ECa by changing ECwc (Rhoades et al., 1989). Therefore, dissolved nutrients such as N and S should have more influence on ECa under saturated conditions than when the soil is close to the wilting point. Unfortunately, large differences in volumetric soil water content across a field will likely mask the small differences that changing nutrient levels will have on ECwc. Based on an analysis of the factors contributing to ECa, it is clear that it will be difficult to directly measure site-specific changes in soil nutrient content using electrical conductivity.
An analysis of the factors that contribute to ECa suggests that ECa might be used along with other measurements of soil properties, such as CEC, texture, pH, and water-holding capacity, to determine nutrient levels in the soil. Studies have suggested that if spatial differences in water-holding capacity and salt concentration were taken into account, there would be a good relationship between soil nutrient levels and ECa (Chang et al., 2001). Others have found that ECa could be used to determine the depth to clay (Doolittle et al., 1994) and that these measurements identified areas where ECa was related to soil nutrient content. Because soil properties such as texture and CEC are relatively static over time, a model that included these measurements along with ECa could be developed to predict soil nutrient content.
Apparent soil electrical conductivity could also be used to indirectly determine field sites where soil nutrient levels differ. Depth of topsoil, soil water content, CEC, and texture are all good indicators of crop productivity (Jaynes et al., 1995; Sudduth et al., 1996). Because ECa is able to directly measure these soil properties, it has the potential to identify management zones with differing productivity and nutrient requirements (Kitchen et al., 1999). Several researchers have reported that ECa was useful in identifying field sites with different soil properties (Sudduth et al., 1999; Franzen and Kitchen, 1999). Sudduth et al. (1996) found that ECa was correlated with depth to clay. Because depth to clay was the primary factor influencing yield levels at different field sites, ECa proved to be a valuable tool in determining site-specific N requirements. These studies indicate that ECa could, indirectly, provide a useful measure of nutrient levels.
The main purpose of this paper is to determine the utility of ECa in mapping soil nutrient levels across a field for the purpose of making variable-rate nutrient applications. The specific objectives of this study are to examine the use of ECa: (i) to directly measure soil pH, texture, CEC, and levels of selected elements (P, K, Ca, Mg, Mn, Cu, and Zn) across a field and (ii) in conjunction with measured soil properties to measure the concentrations of selected elements (P, K, Ca, Mg, Mn, Cu, and Zn) in a field.
| MATERIALS AND METHODS |
|---|
|
|
|---|
|
1999 Data Collection
Nine locations were studied in 1999: four fields in the tidewater region, two in the coastal plain region, and three in the piedmont region (Table 1). At each location, soil samples were collected at two depths (030 cm and 3090 cm) from 30- by 30-m grids. Each composite sample contained six cores collected within a circular area with a diameter of 5 m centered on the sample location recorded using a DGPS receiver. Soil samples were dried, ground to pass through a 2-mm sieve, and analyzed for CEC, HM, pH, active acidity, P, K, Ca, Mg, Mn, Zn, and Cu using the same testing procedures used in 1997. In addition, soil textural analysis was performed on all of the 30- to 90-cm samples and the 0- to 30-cm samples from the coastal plain and piedmont regions using the hydrometer method of Gee and Bauder (1986). The 0- to 30-cm samples from the fields in the tidewater region could not be analyzed because the organic matter content (>15%) of the surface horizon interfered with the testing procedure. On the same day that the soil samples were collected, ECa was measured using an EM38 (Geonics Limited, Mississauga, ON, Canada). When taking readings, the EM38 unit was positioned both horizontally and vertically at ground level and at heights of 30, 60, 90, and 120 cm. At three representative sites in each field, soil samples were collected from depths of 15, 30, 45, 60, 75, 90, 105, and 120 cm and analyzed in the laboratory for electrical conductivity. Using the technique described by Rhoades and Corwin (1981), a calibration analysis was performed relating ECa collected by the EM38 with electrical conductivity by depth as measured in the laboratory. This calibration was then applied to all EM38 measurements to obtain ECa for 0- to 30- and 30- to 90-cm depths.
Statistical Analysis
Correlation analysis (SAS Inst., 1990) was conducted on each field and within fields on each soil series to determine if ECa, CEC, HM, pH, soil type, sand, silt, or clay were related to soil test concentrations of P, K, Ca, Mg, Mn, Zn, or Cu. Significant results with a high Pearson correlation coefficient (>0.70) would indicate situations where the measured soil property could be used to estimate the concentration of that particular element in the soil.
Principal-components analysis was used to examine the relationship among the soil properties measured in this study (ECa, CEC, HM, pH, soil type, percentage sand, percentage silt, and percentage clay) and to determine which soil properties were important influences on soil test concentrations of P, K, Ca, Mg, Mn, Zn, or Cu. To account for soil texture in 1997, a soil type variable was formed by assigning a number to the soil series found in each field. This number was then used in the PC analysis to represent the influence of soil texture on differences in the concentrations of the elements tested.
Due to the colinearity of the independent variables, correlation analysis could not be used to directly relate multiple soil properties to soil test levels of the elements tested. Principal-components analysis puts identified, correlated variables into groups. These groups then become new, independent, random variables that could then be used to identify which soil properties influenced nutrient levels. In this study, the objectives of using the PCstepwise regression analysis were to identify the key soil properties that, in conjunction with ECa, had significant relationships with element concentrations; determine the strength of that relationship; and determine the influence and role of each soil property in the relationship.
The PC groups were identified from the correlation matrix using the FACTOR procedure in SAS (SAS Inst., 1990). Any PC group with an eigenvalue greater than 1 was selected because it explained a significant amount of the variance present in the soil properties at each location. In no case did the PC groups with eigenvalues >1 explain <66% of the cumulative variance. More commonly, these PC groups had a cumulative variance of
75%. The PC groups with eigenvalues >1 were then used in a stepwise regression procedure (SAS Inst., 1990) to determine if there was a significant relationship between the PC group and soil test concentrations of P, K, Ca, Mg, Mn, Zn, or Cu. When PC groups remaining in the regression model accounted for >50% of the variability in the concentration of that particular element, the eigenvectors (loading factors) were examined and the soil properties in the PC groups ranked according to the amount of variability explained by the PC group and the loadings (C|j|). For instance, a soil property that was a component of the PC group that accounted for most of the variability in the regression model and had the highest loading factor in that PC group was ranked first. Soil properties with loading factors <0.5 were not considered key latent variables and were not included in the ranking because they did not substantially influence the relationship between the PC groups and the nutrient concentration being examined. The ranking of the soil properties, strength of the loading factor, and sign (positive or negative) of the loading factor were used to determine the influence and role that each soil property had in explaining the variability in the concentration of each element.
| RESULTS AND DISCUSSION |
|---|
|
|
|---|
|
Because volumetric soil particle density and volumetric soil moisture content have a large influence on ECa and are closely related to soil texture as are CEC, pH, and HM, fields B531, B532, B533, B534, and B535 were divided by soil series, and correlations between ECa and element concentrations were examined on each soil series to determine if there was an improvement in predicting concentrations of particular elements based on ECa. In 40 out of 75 cases, there was an improvement in the strength of the relationship between ECa and the concentration of the element measured (Table 2). Figure 1 shows the relationship between ECa and soil P concentration. While the overall relationship was weak and showed a negative trend, the comparisons on each individual soil series were stronger and had a positive trend, indicating that ECa was responding to increasing ion concentrations in the soil.
|
|
|
|
1997
From the seven fields studied in 1997, there were eight regression models where PC groups accounted for >50% of the variability in the measured nutrient level. Of these eight regression models, ECa was identified as a key latent variable in six of the models (Table 5). In field B533, the PCstepwise regression analysis found that both Mn and Zn were related (R2 = 0.63 and 0.75, respectively) to two PC groups. In PC Group 1 (PC1), pH had the highest loading factor (Table 6), indicating it was negatively related to soil test levels of Mn and Zn. Cation exchange capacity, HM, and ECa had similar loading factors and were positively related to soil test levels of Mn and Zn. Soil type dominated PC Group 2 (PC2), indicating a clear relationship between soil type and soil test levels of Mn and Zn. The change from heavy organic to lighter-textured loam soils across this field and the association of pH, CEC, HM, and ECa with soil series (lower pH and higher CEC, HM, and ECa on the organic Belhaven and Ponzer soil series and higher pH and lower CEC, HM, and ECa on the lighter Wasda and Deloss soil series) clearly influenced the relationship with Mn and Zn.
|
|
In field H08, the PCstepwise regression analysis found that soil test levels of Mg were related to PC1 and PC2 (R2 = 0.68) (Table 5). Principal-Component Group 1 was dominated by ECa (C|j| = 0.912) followed by positive loadings from HM and CEC (Table 6). Principal-Component Group 2 was dominated by pH (C|j| = 0.934), with a much smaller contribution from CEC. The strong influence of ECa in this regression model indicates a situation where ECa along with other measured soil properties could be used successfully to determine element concentrations across a field.
Similarly, the PCstepwise regression analysis found that soil test levels of Mg were related to PC1 in field F1 (Table 5). Principal-Component Group 1 was dominated by the negative loading factors CEC and HM (Table 6), with smaller positive loadings from pH and ECa. Field F1 is located in the rolling hills of the piedmont region of North Carolina and has a clay loam soil texture. Changes in CEC in this field are associated with an increase in clay content and CEC on the sideslopes where runoff has eroded the surface horizon, exposing the B horizon where Mg levels are low. The positive relationship with ECa is primarily the result of higher Mg levels in the toeslope areas of the field where eroded materials accumulate and where soil moisture is generally greater.
19990- to 30-cm Sampling Depth
From the nine fields studied in 1999, there were 20 regression models where PC groups accounted for >50% of the variability in the measured nutrient levels in the 0- to 30-cm depth. Of these 20 regression models, ECa was identified as a key latent variable in 13 of the models (Table 7). The results in 1999 were similar to those found in 1997. In situations where ECa was identified as a key latent variable in describing nutrient concentrations, changes in soil texture (Table 4) and CEC were also key factors in determining nutrient levels in the field.
|
|
|
199931- to 90-cm Sampling Depth
At the 31- to 90-cm sampling depth, there were 25 regression models where PC groups accounted for >50% of the variability in measured nutrient levels. Of these 25 regression models, ECa was identified as a key latent variable in 18 of the models (Table 7). As was the case in the soil samples and ECa measurements taken at the 0- to 30-cm depth, Ca and Mg were the two elements where PC groups comprised of the soil properties measured in this study often accounted for most of the variability in measured levels. In eight of the nine fields, the regression models accounted for >54% of the variability in Ca and Mg levels.
The dominant loading factors in fields B500, MR1, GR1, BOT1, BOT4, Fischer, and Graham were sand, silt, and clay content; CEC; HM; and ECa (Tables 8 and 9). Except for CEC in the Fischer field and HM in the Graham field, CEC, HM, and ECa were positively related to Ca or Mg levels. In fields B500, BOT1, and Graham, sand was negatively related to Ca and Mg levels while silt and clay were positively related. These results indicate that the positive relationship between ECa and percentage clay is behind the positive relationship between ECa and Ca or Mg soil test levels. However, for fields MR1, GR1, BOT4, and Fischer, sand content was a positive loading factor while silt and clay were negatively related to nutrient concentration. The results on these fields indicate that the positive relationship between ECa and either Ca or Mg soil test levels must be influenced by other factors, such as CEC or HM. In the case of GR1, this could be the result of the positive association among salinity from the manure application, ECa, and Mg concentrations.
Other elements in which the PCs accounted for >50% of the variability in the measured samples and ECa was a key loading included Mn on fields B500, BOT1, and Hall and Zn and Cu on MR1 (Table 7). For Mn, CEC, ECa, silt, clay, and pH were positive key loading factors while sand was negatively related to soil test levels (Tables 8 and 9). Although the PCstepwise regression analysis of the data from MR1 indicated that these same key loading factors were related to Zn and Cu levels in this field (accounting for 58 and 51% of the variability, respectively) (Table 7), the relationship between the key loading factors and these elements differed. Apparent soil electrical conductivity, CEC, HM, pH, and sand were positive loading factors while silt and clay were negatively related to Zn and Cu soil test levels (Table 8). Based on the strong influence that ECa and CEC had on these regression models (C|j| = 0.679 and 0.671, respectively) and the negative influence of silt and clay, it is reasonable to assume that ECa was influenced by the relationship between CEC and changes in Zn or Cu concentrations.
| CONCLUSIONS |
|---|
|
|
|---|
Improvements in the correlation between ECa and soil nutrient levels were often found when fields were subdivided by soil series. On field B533, the amount of variability explained by a linear model relating ECa to soil P content increased from 3 to 24% by subdividing the field according to soil series (Table 2 and Fig. 1). The fractions of sand, silt, and clay in the soil sample were often significantly related to ECa (Table 4). Because soil texture influences both volumetric soil moisture content and the volumetric content of soil particles, this would make sense. Therefore, differences in soil texture often influence ECa much more than small differences in element concentrations. By dividing fields into areas of similar soil texture, one source of variation was removed, and this improved the accuracy of using ECa to measure changes in soil concentrations of nutrient elements.
Relationships among Apparent Soil Electrical Conductivity, Other Soil Properties, and Soil Nutrient Content
In the 16 fields studied, the PCstepwise regression analysis was able to identify 37 regression models out of a possible 165 comparisons in which PC groups accounted for
50% of the variability in nutrient levels and where ECa was one of the key loading factors (Tables 5 and 7). Principal-component groups consisting of ECa and other measured soil properties were not able to consistently account for changing levels of most of the elements measured in this study, in particular, P and K. However, as the initial correlation analysis indicated, ECa along with measured soil properties, particularly CEC and soil texture, was able to account for a majority of the variability in soil test levels of Ca and Mg in 13 out of 22 comparisons made at the 0 to 30-cm depth (R2 from 0.520.93) and in 13 out of 18 comparisons made at the 31- to 90-cm depth (R2 from 0.540.88).
In most of these comparisons, ECa, CEC, and clay were all identified as positive loading factors. Under these non-water-saturated conditions,
ws,
s, and ECws would be the most important factors influencing ECa. As clay content increased, there would be an increase in
s and
ws and increased concentration of Ca and Mg cations in the soil solution of the discontinuous pores associated with increases in CEC would have increased ECws. While it is reasonable to assume that the influence of silt and clay content on
s and
ws were important influences on ECa, the relationship between ECa and Ca or Mg soil test levels was not always associated with an increases in silt and clay content. On six fields, clay was negatively related to Ca or Mg soil test levels even though ECa had a positive relationship. This occurred most often in the piedmont where increases in clay content were associated with eroded areas of the field where the concentration of Ca or Mg was low. In these situations, increases in CEC contributed to the increased concentration of Ca and Mg in the soil solution and to increases in ECa.
Unfortunately, even when describing Ca and Mg, the PCstepwise regression analysis was not able to consistently identify models that accounted for most of the variability in soil test levels. Furthermore, the model parameters, PC groups, and key loading factors differed according to the field being studied. This means that a single relationship among ECa, soil properties, and soil element concentrations cannot be developed. This doesn't mean that ECa has no value in determining nutrient levels in the soil. Instead, this study shows that ECa can be valuable tool when used in conjunction with multivariate statistical procedures in identifying soil properties and their relationship to nutrient availability. Because there are close relationships among CEC, soil texture, and ECa, measurements of ECa could be substituted for detailed soil textural sampling to develop a detailed field map showing where element levels are likely to change. This map could then be used to determine nutrient management zones. Franzen and Kitchen (1999) demonstrated the use of ECa to delineate management zones. Apparent soil electrical conductivity could also be used in co-kriging to reduce the density of soil sampling without sacrificing the accuracy of the mapped data.
This study shows that it is unlikely that ECa can be used to directly determine soil nutrient content across a field. However, ECa can be used along with other measured soil properties in a multivariate analysis to describe the key factors influencing changes in nutrient concentrations and to establish nutrient management zones. To use ECa in the process of establishing nutrient management zones, it is important to know and account for site-specific changes in volumetric soil moisture content, soil texture, CEC, and salinity (either from manure applications or on saline soils).
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
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] |
||||
![]() |
M. S. Grigera, R. A. Drijber, K. M. Eskridge, and B. J. Wienhold Soil Microbial Biomass Relationships with Organic Matter Fractions in a Nebraska Corn Field Mapped using Apparent Electrical Conductivity Soil Sci. Soc. Am. J., August 3, 2006; 70(5): 1480 - 1488. [Abstract] [Full Text] [PDF] |
||||
![]() |
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] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| 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 | |||