SSSAJ Journal of Natural Resources and Life Sciences Education
HOME HELP FEEDBACK SUBSCRIPTIONS ARCHIVE SEARCH TABLE OF CONTENTS
 QUICK SEARCH:   [advanced]


     


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (18)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
Related Collections
Right arrow Water Management
Right arrow Watershed and Landscape Processes
Right arrow Spatial Distribution
Soil Science Society of America Journal 66:235-243 (2002)
© 2002 Soil Science Society of America

DIVISION S-6 - SOIL & WATER MANAGEMENT & CONSERVATION

Quantitative Mapping of Soil Drainage Classes Using Topographical Data and Soil Electrical Conductivity

A. N. Kravchenko, G. A. Bollero, R. A. Omonode and D. G. Bullock*

Dep. of Crop Sciences, 1102 S. Goodwin Ave., Univ. of Illinois, Urbana, IL 61801

* Corresponding author (dbullock{at}uiuc.edu)


    ABSTRACT
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
In this study we applied discriminant analysis and geostatistics to create soil drainage maps using topographical and soil electrical conductivity (EC) data as auxiliary information. Drainage classes were determined on 107 soil cores collected from a 20-ha field in Illinois. Approximately 1500 elevation points and 6500 EC points were collected from the field. Slope, curvature, and flow accumulation were derived. Discriminant analysis and logistic discrimination were applied to study the effect of topography and EC on drainage. Spatial variability of the soil drainage data and its relationship with variability of topography and EC were analyzed using variograms and cross-variograms. Indicator kriging and soft indicator cokriging were used. Soil EC, terrain slope, and distance to a drainageway were selected by a stepwise discriminant procedure as significant predictors of the soil drainage class (P = 0.15). When these variables were used as additional information in predicting soil drainage class using either discriminant analysis or cokriging procedures, they slightly improved overall prediction accuracy, compared with the soil survey map (scale 1:15840) and indicator kriging. Discriminant analysis and cokriging correctly estimated drainage classes for more than 90% of the sites, compared with 85 and 63% correct estimates obtained from indicator kriging and soil survey data, respectively. Indicator kriging and cokriging are exact estimators and thus are more suitable for mapping than is discriminant analysis. We recommend the use of stepwise discriminant procedure to select among the available secondary variables and to create drainage maps using cokriging with the selected variables.

Abbreviations: EC, electrical conductivity • GPS, global positioning system • MWD, moderately well-drained soil • PD, poorly drained soil • SWPD, somewhat poorly drained soil


    INTRODUCTION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
SOIL DRAINAGE is an important soil property affecting plant growth, water flow, and solute transport in soils. Accurate and inexpensive prediction and mapping of soil drainage classes is beneficial for both agricultural and environmental management. Various sources of secondary information, such as topography or easily measured soil properties, can be used to facilitate drainage mapping and prediction, however, efficient quantitative methods are needed to utilize this information. Quantitative approach to predicting soil drainage classes using discriminant analysis has been proposed by Bell et al. (1992) and has been found to be an effective method for soil drainage mapping (Bell et al., 1994). The main idea of the discriminant analysis is that, based on the measured data for a certain categorical variable of interest (i.e., drainage) and for several auxiliary variables, a quantitative decision rule is formulated. Every new observation where auxiliary variable data are available but the variable of interest has not been measured can be placed into the appropriate class of the variable of interest based on the previously determined decision rule. The necessity of finding the general decision rule pertains to the fact that each one of the studied auxiliary variables by itself may not adequately describe the variable of interest. However, the decision rule is calculated based on the measured data from all the auxiliary variables, hence it can have better predictive powers than any single variable.

One of the disadvantages of discriminant analysis approach to mapping soil properties is that the information on the data locations is not used in development of the predictive model. However, this information can be utilized by geostatistical methods of data analysis to determine spatial structure in data distribution across the studied area. Combining the predictive capabilities of discriminant analysis with those of geostatistics may provide more accurate tools for soil drainage mapping and prediction.

Bell et al. (1992)(1994) incorporated a variety of topographical and parent material data to map drainage classes of a relatively large study area (about 14000 ha) in south-central Pennsylvania with diverse lithology and topography. They found that soil parent material was one of the most important factors in predicting soil drainage classes. The factors that were found to be significant in discriminating between the drainage classes within the same parent material included slope gradient, elevation above local stream, and distances to local streams and drainageways. Soil drainage mapping and prediction for relatively small areas (e.g., agricultural fields or small watersheds) are also of great importance, particularly for farming and site-specific field management. However, the factors affecting soil drainage classes within farm fields might be different from the factors acting at regional scales. Parent materials within a single field are not usually as diverse as in a relatively large area, but at the same time they might be highly spatially variable. Unfortunately, very often it is not possible to define parent material of a particular site with a considerable degree of accuracy without taking soil cores—a task that is time-consuming and expensive. Therefore, it would be beneficial to determine other easily measured factors for quantitative prediction of soil drainage on small scales and to estimate potential accuracy of such prediction.

Topography is one of the easily evaluated drainage affecting factors that can greatly facilitate mapping and prediction of soil drainage (Bell et al., 1992; Bell et al., 1994; Troeh, 1964). Another easily measured factor potentially related to drainage is soil EC. It can be measured in the field using recently developed fast and nondestructive methods (Doolittle et al., 1994; Kitchen et al., 1996; Kitchen et al., 1999). Electrical conductivity depends on a number of soil physical properties, including soil salinity, and soil water and clay content (Rhoades et al., 1989). Sheets and Hendrickx (1995) found a linear relationship between the EC measurements and soil profile water content. Williams and Hoey (1987) observed positive correlations between conductivity and soil clay content. Since both soil profile water contents and clay contents are related to the soil drainage properties, it would be reasonable to expect that EC data might be helpful in mapping and predicting soil drainage classes on a field scale basis.

The objectives of this study are (i) to examine the relationships between soil drainage and various topographical factors and to determine the factors that can be most helpful in predicting soil drainage classes on a field scale using the results from a farm field in Central Illinois, (ii) to study possibility of using soil EC data with or without topographical information for predicting soil drainage classes, (iii) to compare accuracy and efficiency of multivariate statistical methods with geostatistical procedures in mapping soil drainage classes on an agricultural field scale.


    MATERIALS AND METHODS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Soil and Topographical Data
The studied area was a square 20-ha central portion of a 259-ha agricultural field located in central Illinois, USA. The field has been in a corn (Zea mays L.)–soybean [Glycine max (L.)] rotation for at least 20 yr. Soil sampling was conducted on a semiregular grid and sample locations are shown in Fig. 1 . A total of 107 soil cores, 100 cm deep, was collected from the studied area.



View larger version (101K):
[in this window]
[in a new window]
 
Fig. 1. Locations of the core samples used to determine soil drainage classes along with the soil survey map and elevation map of the studied area. Each sampling site is labeled according to its drainage class.

 
The drainage class for each sample location was determined based on the visual examination of the cores. Drainage classes were assigned according to the depth to the seasonal high water table, as indicated by the presence of low-chroma (<2) mottles (Fehrenbacher et al., 1984; Soil Survey Staff, 1993). The soil was assigned to the well-drained class (WD) if no mottling was observed in the soil core within a 100-cm depth. Soil with mottling occurring at 50- to 100-cm depth was classified as moderately well-drained (MWD), and soils with mottling occurring within 25 to 50 cm and within top 25 cm were classified as somewhat poorly drained (SWPD), and poorly drained (PD), respectively. Out of 107 soil cores, 48 were classified as WD soils, 38 as MWD soils, 15 as SWPD soils, and 6 as PD soils (Fig. 1). Since the number of PD soil cores was insufficient for meaningful quantitative analysis, they were combined in one class with SWPD soils (SWPD/PD).

Soil survey map (Soil Survey of McLean County, Illinois, 1998) characterizes the soils of the field as the Varna series (fine, illitic, mesic Typic Argiudoll), the Drummer series (fine-silty, mixed, mesic Typic Endoaquoll), and the Chenoa series (fine, illitic, mesic Aquic Argiudoll). According to the soil classification of Illinois, soils of Varna series are described as WD and MWD, Drummer soils are described as PD, and Chenoa soils as SWPD (Fehrenbacher et al., 1984). Since soils of Varna series can belong to both WD and MWD drainage classes, for comparing accuracy in predicting drainage classes of soil survey map with that of quantitative methods we used only two drainage classes: WD/MWD (combined WD and MWD sites) and SWPD/PD. Soil series map unites of the studied area are shown in Fig. 1.

Soil EC data were collected in Fall of 1999 using a Veris 3100 sensor cart (Division of Geoprobe Systems, Salina, KS) which operates on the principle of electromagnetic induction. Georeferenced EC measurements were taken every 3 to 5 m with the distance between cart passes about 10 m, resulting in about 6900 EC measurement points for the studied area. Two sets of EC measurements were collected corresponding to depths of ~0 to 30 cm (shallow EC) and 0 to 90 cm (deep EC) (Sudduth et al., 1998).

Survey grid global positioning system (GPS) (Leica 500 RTK, Leica, Heerbrugg, Switzerland) was used to measure elevation. About 1500 elevation measurements were taken on a semiregular grid with a mean distance between measurements of ~10 m. The map of field elevation ranges is shown in Fig. 1. ArcView Spatial Analyst (Environmental Systems Research Institute, 1996) was used to analyze elevation data and to derive the topographical land features, namely, slope, curvature, and flow accumulation. Detailed description of the derivation procedures is presented by Kravchenko and Bullock (2000). ArcView was also used to calculate the shortest distance from each of the sampling sites to the drainageway that ran in the middle of the studied area (Fig. 1). This distance was used as one of the variables with potential effect on soil drainage class.

Discriminant Analysis
In this study discriminant analysis was used to find a decision rule for separating sites with different drainage classes based on the values of the topographical and EC variables measured at each site. The decision rule was defined based on the measured data, then, the decision rule was applied to predict drainage classes at unsampled sites where measurements of topographical and EC variables were available.

Two combinations of the variables were considered. The first combination included topographical variables, such as elevation, slope, curvature, flow accumulation, and distance to a drainageway. The purpose of this combination was to assist with the first objective of the study, i.e., to determine how accurate drainage class prediction can be if it is based on topographical data only. The second combination of the variables was related to the second objective of the study, namely, analysis of the possibility of using soil EC data with or without topographical information for predicting soil drainage classes. All the topographical and EC variables measured at each site were used in the second combination.

Since meaningful application of discriminant analysis requires data to be normally distributed, we analyzed the normality of all studied variables. Nonnormally distributed variables were transformed using either log-normal or normal score (Goovaerts, 1997) transformations. Stepwise discriminant procedure, STEPDISC (SAS Institute, Cary, NC) was applied to both combinations of variables to select the variables that had significant influence on the soil drainage class. Only variables significant at 0.15 significance level (Bell et al., 1992) were used in further discriminant analysis.

Based on the values of soil and topographical variables combined in a variable measurement vector x, for each data location discriminant analysis calculated the probability that this location belonged to each of the drainage classes, t, called posterior probability, {pi}. Then, the location was classified into the drainage class for which it had the highest posterior probability (Khattree and Naik, 2000):

[1]

The highest posterior probability is equivalent to the smallest squared distance from the studied location to the drainage class t, D2t(x), where k is the total number of drainage classes. Subscript t in Eq. [1] refers to a certain drainage class for which the posterior probability is calculated, subscript j indicates that information on all drainage classes is used to calculate posterior probability for each particular drainage class t. Depending on the variance-covariance matrix of the studied combination of variables, the squared distance is calculated as either a linear discriminant function:

[2]
if all the drainage classes have the same variance-covariance matrixes, or as a quadratic discriminant function if the variance-covariance matrixes are not equal:

[3]
where µt is the mean measurement vector of drainage class t, {sum} is the common variance-covariance matrix in Eq. [2], St is the variance-covariance matrix of drainage class t in Eq. [3], and {pi}t is the prior probability of occurrence of drainage class t in the studied area. In this study, we calculated prior probabilities based on the measured proportions of drainage classes. Based on the number of locations belonging to a certain drainage class, prior probabilities were equal to 0.499, 0.355, and 0.196 for WD, MWD, and SWPD/PD classes, respectively. Discriminant analysis was conducted using DISCRIM procedure (SAS Institute, Cary, NC).

Cross-validation was applied to evaluate accuracy of drainage data prediction and classification by discriminant analysis. For cross-validation each value from the data set was eliminated in turn and, then, estimated using information from the rest of the data (Khattree and Naik, 2000). Posterior probabilities of the three drainage classes obtained for each data point were compared, and the site was assigned a drainage class with the highest posterior probability. Percentage of the correct drainage class estimates was used to compare different combinations of variables with their effectiveness in predicting drainage classes.

Effect of the EC and topographical variables in discrimination between the drainage classes was also studied using logistic discrimination (Khattree and Naik, 2000). Unlike discriminant analysis, logistic discrimination is not sensitive to the departures from normality, hence, no data transformation is necessary. We conducted logistic discrimination for each drainage class separately, using constructed indicator variable i({alpha}) and, then, compared the results from all drainage classes. For each of the three drainage classes, drainage class at each data location {alpha} was coded using the indicator approach as (Goovaerts, 1997; Khattree and Naik, 2000):

[4]
resulting in three new indicator variables. For each drainage class, a new indicator variable was represented by sets of 1's corresponding to the locations belonging to this drainage class and sets of 0's corresponding to all other locations belonging to the other drainage classes.

Value of the indicator variable is equivalent to a probability, {pi}, of finding a certain drainage class at a certain location. For a location with known drainage such probability equals to either 1 (site belongs to a certain drainage class) or to 0 (site belongs to any of the other drainage classes). The response variable in the logistic determination is expressed in terms of probability, {pi}, as a logit function, logit({pi}):

[5]

The response variable is assumed to be a linear function of a number of explanatory variables (x1...xn):

[6]
where the input of each explanatory variable xi to the logit({pi}) is determined by the unknown parameter ßi. In this study, EC and topography variables are used as the explanatory variables and the relationship expressed in Eq. [6] produces a model for discriminating between the drainage classes. Similarly to the discriminant analysis, the posterior probabilities of the three drainage classes obtained using logistic discrimination were compared and the location was assigned a drainage class with the highest posterior probability value. Stepwise variable selection was used to determine the variables with the highest discriminative effect (P = 0.15) on drainage classes using logistic discrimination procedure, LOGISTIC (SAS Institute, Cary, NC).

Geostatistical Analysis
For geostatistical analysis, soil drainage class, t, was treated as a categorical variable. At each data location this variable assumed only one of the three mutually exclusive possible states, corresponding to either WD, MWD, or SWPD/PD drainage classes. Indicator transformation was applied to each drainage class data (Eq. [4]) resulting in three indicator variables.

Spatial variability of each indicator variable was described by a sample variogram:

[7]
where {alpha} and {alpha} + h were sampling locations separated by a distance h, i({alpha}) and i({alpha} + h) were values of the indicator variable i at the corresponding locations. The sample variograms were fitted with variogram models and accuracy of the fitting was examined using cross-validation approach (Goovaerts, 1997). The model that produced the cross-validation results of the highest accuracy was retained for further analysis (Kravchenko and Bullock, 1999). Variable values at unsampled locations were obtained using kriging procedure based on the variogram model parameters. For an unsampled location, the kriging estimate was equivalent to the probability of finding certain drainage class at this location. Such probability assumes a value between 0 and 1 depending on the variable values from the surrounding sampled locations and on the spatial structure of variable distribution reflected in the variogram model parameters. For each drainage class kriging produces a map of probabilities of finding this drainage class at each particular location. The cross-validation results for kriging and cokriging procedures were further used to compare their accuracy with accuracy of discriminant analysis.

Two geostatistical procedures compared in the study included ordinary indicator kriging and soft indicator cokriging (Goovaerts, 1997). Indicator kriging uses the drainage data only, hence, the estimates of the indicator variables depend only on spatial structure of the available drainage data. Cokriging allows us to combine primary drainage data with any available secondary data that are related to drainage, hence, drainage estimates depend not only on distribution of the primary data, but also on spatial variability of the secondary data and their relationship with the primary variable. In case of a strong correlation between primary and secondary variables and when abundantly sampled secondary data complement sparse primary data sampling, cokriging can produce much more accurate and reliable estimation results comparing with ordinary kriging technique. One of the disadvantages of cokriging is that it becomes extremely cumbersome and time-consuming with large number of secondary variables. In this study, we conducted cokriging with one primary variable (indicator transformed drainage class) and two secondary variables. The variables that were found to have a significant effect on drainage class during stepwise discrimination procedure were further used in cokriging estimation. Variograms and crossvariograms for cokriging were modeled using linear model of coregionalization (Deutsch and Journel, 1998). Geostatistical analysis including variogram calculation, cross-validation, kriging, and cokriging was performed using the geostatistical software package GSLIB (Deutsch and Journel, 1998).


    RESULTS AND DISCUSSION
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Spatial aspects of the relationships between drainage classes and topographical and EC variables were examined using variograms and cross-variograms. Sample variograms of the three studied drainage classes along with their respective variogram models are shown in Fig. 2 . Examples of sample cross-variograms between drainage classes and slope are shown in Fig. 3 . Different spatial relationships existed between slope and WD/MWD classes, and slope and SWPD/PD class (Fig. 3). Cross-variogram values for WD and MWD classes increased with distance indicating positive relationships between these drainage classes and slope. For SWPS/PD class, cross-variogram value decreased with distance, indicating negative correlation between this class and the slope. Indeed, the higher the slope, the lower was the probability of occurrence of SWPD/PD soils, while probability of occurrence of WD and MWD soils increased with increasing slope values. Similar results were observed for cross-variograms of drainage classes vs. distance to drainageway. For WD and MWD classes, cross-variograms increased with distance and for SWPD/PD class, they decreased with distance reflecting respective positive and negative relationships between the probabilities of finding WD/MWD and SWPD/PD soils and the distance to the drainageway. Cross-variograms of drainage class vs. deep EC values decreased with distance for WD and MWD drainage classes and increased with distance for SWPD/PD soils, indicating that there were negative relationships between WD and MWD drainage classes and soil EC values and positive relationship between the SWPD/PD class and soil EC.



View larger version (18K):
[in this window]
[in a new window]
 
Fig. 2. Sample variograms and variogram models of indicator transformed data for (a) well-drained sites, (b) moderately well-drained sites, and (c) somewhat poorly and poorly drained sites.

 


View larger version (23K):
[in this window]
[in a new window]
 
Fig. 3. Sample cross-variograms and cross-variogram models for terrain slope and indicator transformed data for (a) well-drained sites, (b) moderately well-drained sites, and (c) somewhat poorly and poorly drained sites.

 
When only topographical variables (i.e., elevation, slope, curvature, flow accumulation, and distance to drainageway) were used in the stepwise discriminant procedure, the variables with significant effect on drainage were slope and distance to drainageway (Table 1). All the other topographical variables were not significant in discriminating between the drainage classes for this field. When all the above listed topographical variables and both shallow and deep EC measurements were used in the stepwise discrimination, only distance to drainageway and deep EC were selected as significant variables. Mean values of the variables that were selected as significant during stepwise discriminant procedure (P = 0.15) are shown in Table 1.


View this table:
[in this window]
[in a new window]
 
Table 1. Mean values of the discriminant variables selected by stepwise discriminant procedure (P = 0.15) for the three soil drainage classes.

 
The SWPD/PD soils were located close to drainageway, at sites with the lowest slope values (Table 1). These sites were also characterized by the highest values of deep EC. One of the factors influencing the relationship between EC and soil drainage classes is soil clay content. Soil EC was found to be positively correlated with clay content (Williams and Hoey, 1987) and, at the same time, higher clay content is one of the characteristic properties of PD and SWPD in this study. Analysis of soil texture revealed that SWPD/PD soils had somewhat higher clay contents than MWD and WD soils. The differences in clay contents were most notable at lower parts of the profile at depths of 60 to 100 cm, where clay contents of WD, MWD, and SWPD/PD soils were equal to 26.3, 28.2, and 31.3%, respectively. Another factor of potential influence on the relationship between EC and soil drainage classes is soil water content. Similar to clay content, soil water content is positively correlated to soil EC and, in the field conditions, PD soils often are characterized by higher water contents comparing with MWD and WD soils. Unfortunately, soil water content was not monitored during the EC measurements, hence, it was not possible to separate influences of soil water content and clay content on EC data in this study. Another factor of potential influence on EC is soil salinity, however, it can be assumed to be negligible for the field used in the study.

Distance to drainageway seemed to be a useful parameter in discriminating between WD and MWD drainage classes in this field. Sites with MWD soils were located generally closer to the drainageway than the sites with WD soils (Table 1). Slopes and deep EC values of WD and MWD soils were similar, hence, these parameters were not effective in separating WD and MWD drainage classes using stepwise discriminant procedure.

Further insight in the factors discriminating between soil drainage classes was achieved using logistic discrimination (Table 2). Since logistic regression procedure was conducted separately for each drainage class, a set of variables was obtained for each drainage class that were significant in discriminating between this class and all the other drainage classes. When only topographical variables were included in the logistic discrimination procedure, the variables significant in discriminating between WD soils and other drainage classes were distance to drainageway and flow accumulation. Well drained soils were characterized by lower flow accumulation values than either MWD or SWPD/PD soils. For SWPD/PD soils significant discriminating variables were distance to drainageway and slope. When shallow, and deep EC measurements were added to the topographical data, then, distance to drainageway, flow accumulation and deep EC were found to be significant for separating WD soils. For SWPD/PD sites, both shallow and deep EC measurements along with the distance to drainageway were significant for separating this drainage class from WD and MWD soils. Slope was a significant variable for discriminating between MWD soils and all the other drainage classes, however, the effectiveness of the logistic discrimination was extremely low for MWD soils (Table 2).


View this table:
[in this window]
[in a new window]
 
Table 2. Results of logistic discrimination between soil drainage classes using topographical and topographical plus soil electrical conductivity (EC) variables. Values in parentheses are percentages of correctly estimated sites for each drainage class.

 
As described previously, to compare accuracy of soil survey map with quantitative methods in predicting soil drainage classes, we combined the drainage classes into two groups: WD/MWD (combined WD and MWD sites) and SWPD/PD. Soil survey map correctly predicted locations of 90% of the SWPD/PD soils and 56% of WD/MWD soils (Table 3). Discriminant analysis with topographical variables correctly classified 93% of WD/MWD soils and 76% of SWPD/PD soils, and discriminant analysis with topography and EC correctly classified 95 and 71% of WD/MWD and SWPD/PD soils, respectively (Table 3). Cokriging methods also correctly classified more than 90% of WD/MWD soils and more than 70% of SWPD/PD soils (Table 3). Soil survey classified part of the higher located areas of the field as SWPD Chenoa series (Fig. 1). However, soil classification conducted based on the soil cores collected in the study placed the soils from these locations in Dana, Catlin, and Graymont series (fine-silty, mixed, superactive mesic Oxyaquic Argiudolls) belonging to WD/MWD drainage class.


View this table:
[in this window]
[in a new window]
 
Table 3. Number of correctly and incorrectly classified well-drained/moderately well-drained (WD/MWD) and somewhat poorly drained/poorly drained (SWPD/PD) drainage classes and the total percentage of correctly identified classes from discriminant analysis, soil survey map, indicator kriging, and indicator cokriging. Values in parentheses are percentages of correctly estimated sites for each drainage class.

 
Comparison of the quantitative method accuracy in predicting soil drainage was conducted using three drainage classes: WD, MWD, and SWPD/PD. Cross-validation results of discriminant analysis with topographical variables alone and topographical plus EC variables are presented in Table 4. Both sets of variables produced similar results in predicting soil drainage classes. The WD and SWPD/PD soils were predicted most accurately, with ~70% of each being correctly estimated based on topography and topography plus EC data. Separation of WD and MWD classes was less successful because of previously reported insignificant differences between slope and deep EC values corresponding to WD and MWD drainage classes. Discriminant analysis with topography and EC slightly improved the estimation accuracy comparing with topography only because of better separation of WD and MWD drainage classes (Table 4).


View this table:
[in this window]
[in a new window]
 
Table 4. Number of correctly and incorrectly classified well-drained (WD), moderately well-drained (MWD), and somewhat poorly drained/poorly drained (SWPD/PD) drainage classes and the total percentage of correctly identified classes form discriminant analysis, indicator kriging, and indicator cokriging. Values in parentheses are percentages of correctly estimated sites for each drainage class.

 
Cross-validation results of indicator kriging were not as accurate as those from discriminant analysis, although indicator kriging was better than soil survey map in estimating WD locations (Table 3). Indicator kriging correctly estimated drainage classes at only 56% of the sites (Table 4) comparing with more than 60% of correct estimations from discriminant analysis. Indicator kriging also was less accurate than discriminant analysis and cokriging for predicting WD/MWD and SWPD/PD classes (Table 3). Better results were obtained when information on spatial variability of drainage class data was supplemented by the topographical and EC data in cokriging procedure. Indicator cokriging with slope and distance to drainageway produced the most accurate drainage class prediction comparing with the other methods (Table 4). It correctly estimated all but one WD site, with less accurate estimations of MWD and SWPD/PD sites (Table 4). Indicator cokriging with deep EC and distance to drainageway correctly estimated 18 out of 21 SWPD/PD sites, and was as accurate as discriminant analysis in classifying WD locations, however, it misclassified most of the MWD sites (Table 4). Higher prediction accuracy provided by cokriging was expected, since the primary variable (soil drainage class) was related to the secondary variables (topographical and EC variables) and the secondary variables were sampled much more intensively comparing with the primary variable. However, cokriging with several secondary variables is a time- and labor-consuming process. Indeed, cokriging with just two secondary variables requires computation and model fitting for three variograms (one for primary and one for each secondary variable) and three cross-variograms (two for primary variable versus each secondary variable and one for the secondary variables). Having more than two secondary variables would increase the time and effort spent in variogram computing and modeling even more. Hence, selecting a limited number of secondary variables with the strongest relationship to the primary variable by means of discriminant analysis would greatly facilitate cokriging.

The most substantial difference between discriminant analysis and kriging and cokriging for creating maps pertains to the fact that kriging and cokriging procedures are exact estimators, while discriminant analysis is not. Hence, each of the sampled locations is estimated correctly in the maps produced by kriging and cokriging, while the values in between the sampled locations are estimated using kriging and cokriging procedure. The maps created based on the discriminant analysis ignore the actual data in sampled locations. Hence, indicator kriging and cokriging are preferable comparing with discriminant analysis for creating maps. Drainage maps obtained from (i) discriminant analysis with deep EC and distance to drainageway, (ii) indicator kriging, and (iii) indicator cokriging with slope and distance to drainageway are shown on Fig. 4a, 4b, and 4c , respectively. The data for the maps were obtained on a 11 by 11 m grid. The grid size was selected in correspondence with sampling density of the secondary data, i.e., elevation and EC.



View larger version (96K):
[in this window]
[in a new window]
 
Fig. 4. Drainage maps obtained from (a) discriminant analysis with deep electrical conductivity and distance to drainageway, (b) indicator kriging based on the measured drainage classes, and (c) indicator cokriging based on the measured drainage classes with slope and distance to drainageway as secondary variable.

 

    CONCLUSIONS
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 
Discriminant analysis and logistic discrimination of the soil drainage classes versus topographical and electrical conductivity data revealed that different factors influenced discrimination of each drainage class. Terrain slope, and shallow and deep EC measurements were among main factors separating SWPD/PD drainage class from WD and MWD drainage classes. Distance to drainageway, flow accumulation, and deep EC measurements were significant for separating WD sites from the MWD and SWPD/PD sites. Distance to drainageway was also found to be useful for separating WD and MWD soils. Discriminant analysis accurately predicted ~70% of sites with WD and SWPD/PD drainage classes. However, much lower accuracy (~45% correct) was achieved in separating sites with MWD soils from WD and SWPD/PD sites.

Accuracy of discriminant analysis (62–64% correct) in predicting soil drainage class was comparable with that of cokriging (58–65% correct), whereas kriging, which does not use secondary information in data estimation, was less accurate (56% correct). However, since kriging and cokriging are exact estimators, they are more effective for mapping purposes than discriminant analysis. Since cokriging procedure becomes time-consuming and cumbersome with a large number of secondary variables, the most advantageous seems to be an approach to drainage class mapping that would combine both discriminant analysis and geostatistics. In this approach the data are, first, analyzed with discriminant analysis and the limited number of variables with significant effect on the drainage classes is selected among all the available topographical and soil variables. Then, the significant variables are used in the cokriging for most accurate mapping of soil drainage classes.

Received for publication March 1, 2001.


    REFERENCES
 TOP
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 CONCLUSIONS
 REFERENCES
 




This article has been cited by other articles:


Home page
Agron. J.Home page
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]


Home page
Vadose Zone JHome page
A. M. Tarquis, N. R. A. Bird, A. P. Whitmore, M. C. Cartagena, and Y. Pachepsky
Multiscale Entropy-based Analysis of Soil Transect Data
Vadose Zone J., May 27, 2008; 7(2): 563 - 569.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
W. K. Jung, N. R. Kitchen, K. A. Sudduth, and S. H. Anderson
Spatial Characteristics of Claypan Soil Properties in an Agricultural Field
Soil Sci. Soc. Am. J., June 21, 2006; 70(4): 1387 - 1397.
[Abstract] [Full Text] [PDF]


Home page
Agron. J.Home page
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]


Home page
Agron. J.Home page
A. N. Kravchenko, K. D. Thelen, D. G. Bullock, and N. R. Miller
Relationship among Crop Grain Yield, Topography, and Soil Electrical Conductivity Studied with Cross-Correlograms
Agron. J., September 1, 2003; 95(5): 1132 - 1139.
[Abstract] [Full Text] [PDF]


This Article
Right arrow Abstract Freely available
Right arrow Figures Only
Right arrow Full Text (PDF) Free
Right arrow Alert me when this article is cited
Right arrow Alert me if a correction is posted
Services
Right arrow Similar articles in this journal
Right arrow Similar articles in ISI Web of Science
Right arrow Alert me to new issues of the journal
Right arrow Download to citation manager
Right arrow reprints & permissions
Citing Articles
Right arrow Citing Articles via HighWire
Right arrow Citing Articles via ISI Web of Science (18)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Kravchenko, A. N.
Right arrow Articles by Bullock, D. G.
Related Collections
Right arrow Water Management
Right arrow Watershed and Landscape Processes
Right arrow Spatial Distribution


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