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 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 Web of Science (22)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
Related Collections
Right arrow Water Management
Right arrow Soil Models
Soil Science Society of America Journal 66:1596-1602 (2002)
© 2002 Soil Science Society of America

DIVISION S-6—SOIL & WATER MANAGEMENT & CONSERVATION

Saturated Hydraulic Conductivity and Its Impact on Simulated Runoff for Claypan Soils

Humberto Blanco-Canqui*,a, Clark J. Gantzera, Stephen H. Andersona, E. E. Albertsb and F. Ghideyc

a Dep. of Soil and Atmospheric Sciences, Univ. of Missouri-Columbia, 302 Anheuser-Busch Natural Resources Building, Columbia, MO 65211
b USDA-ARS, Univ. of Missouri, Columbia,Columbia, MO 65211
c Dep. of Biological Engineering, Univ. of Missouri, Columbia, 269 Agricultural Engineering Building, Columbia, MO 65211

* Corresponding author (hb91d{at}mizzou.edu)


    ABSTRACT
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Saturated hydraulic conductivity (Ksat) is an essential parameter for understanding soil hydrology. This study evaluated the Ksat of in situ monoliths and intact cores and compared the results with other studies for Missouri claypan soils. These Ksat values were used as runoff-model inputs to assess the impact of Ksat variation on simulated runoff. Lateral in situ Ksat of the topsoil was determined on 250 by 500 by 230 mm deep monoliths. These values were compared with the Ksat of 76 by 76 mm diam. intact cores with and without bentonite to seal macropores. Mean (± SD) lateral in situ Ksat was 72 ± 0.7 mm h-1 and mean intact core Ksat without bentonite was 312 ± 58 mm h-1. The mean intact core Ksat without bentonite was significantly larger than the lateral in situ Ksat (P = 0.03). The lateral in situ Ksat was not different from core Ksat with bentonite (71 ± 1.1 mm h-1). The intact core Ksat with bentonite differed from previous studies by 10 times. This was attributed to the variations in soil depth to claypan, macropore presence, and methodology. The impact of using an effective hydraulic conductivity (Keff) computed from measured Ksat on intact cores without bentonite underestimated the Water Erosion Prediction Project (WEPP) simulated runoff by 28% for a measured runoff event of 40 mm. The core Ksat with bentonite was correlated with measured runoff from long-term erosion-runoff plots. A quadratic regression explained 95% of the variability between measured and simulated runoff.

Abbreviations: EC, electrical conductivity • Keff, effective hydraulic conductivity • Ksat, saturated hydraulic conductivity • REV, representative elementary volume • SAR, Na adsorption ratio • WEPP, Water Erosion Prediction Project


    INTRODUCTION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
SATURATED HYDRAULIC CONDUCTIVITY is an essential parameter for understanding soil water movement. It is a fundamental input for modeling runoff, drainage, and movement of solutes in soils (Mallants et al., 1997). While Ksat is widely studied, questions remain about how sample size and boundary conditions influence its determination.

Reports have found that measurements on small samples (<100-mm diam.) tend to give higher Ksat values than do measurements on larger samples (Bagarello and Provenzano, 1996). The values of small samples are also questioned because samples are too small to embody a representative elementary volume (REV) of soil. The REV is a conceptual unit representing the smallest volume of a soil unit (Mallants et al., 1997). Its actual dimensions are ill defined. Bouma (1980) suggests three REV sizes for Ksat determinations: 100 cm3 for sand, 1000 cm3 for silt, and 10000 cm3 for clay soils. As a sample size increases, variability in Ksat values is expected to decrease.

The use of the REV is thought to reduce the sample-size dependence of Ksat, and thus facilitate better measurements (Mallants et al., 1997). Samples based on the REV often reflect the natural boundary conditions (Gupta et al., 1993), and diminish disturbance and compaction of soil during sampling (Vepraskas and Williams, 1995).

Soil texture is generally known to affect Ksat. Clay soils typically have low Ksat values (Bouma, 1980; Jamison and Peters, 1967). This is of interest in the midwest USA because about 4 million ha of claypan soils exist in this region (Jamison et al., 1968). These soils have an argillic horizon 130 to 460 mm deep, with clay contents >450 g kg-1 and are very slowly permeable although published data are limited (Jamison and Peters, 1967).

Because of the argillic horizon, claypan soils may perch water and create lateral flow. A study of claypan hydrology suggests that runoff rates may be equal to rainfall under saturated conditions (Saxton and Whitaker, 1970). Furthermore, studies of runoff and rainfall data from the McCredie rainfall-erosion plots near Kingdom City, MO, indicate that lateral flow known as interflow may be a significant component of the total runoff during springtime when precipitation is usually the most intense and the erosion rates are the highest (Minshall and Jamison, 1965; Ghidey and Alberts, 1998). To date, detailed in situ lateral Ksat studies have not been conducted for Missouri claypan soils because measurements are costly and time-consuming (Blevins et al., 1996). Lateral Ksat measurements are also limited elsewhere (Ahuja and Ross, 1983; Wallach and Zaslavsky, 1991). The need for in situ lateral Ksat determination for Missouri claypan soils has been recognized because of the probability of interflow (Jamison et al., 1968; Wilkinson and Blevins, 1999). Information on in situ lateral Ksat through the horizons above the claypan is important for determining their ability to conduct water laterally and assessing runoff and erosion.

Many have characterized the vertical Ksat for claypan soils (Doll, 1976; Zeng, 1994). However, most of the measurements were made only for the surface horizons (Jamison and Peters, 1967; McGinty, 1989), therefore, studies of Ksat variations with depth are few. Because of their hydrologic attributes, claypan soils probably have quite different effective Ksat values with depth from other Alfisols. The information on Ksat depth distribution would be valuable in explaining the claypan hydrology and for characterization of variability in horizons of low and high permeability required for accurate flow studies.

Because the Ksat values may vary by measurement method (Bouma, 1980; Bagarello and Provenzano, 1996; Mallants et al., 1997), the available Ksat data on these soils need to be studied to determine their consistency and uniformity by method. Data from such measurements should be statistically the same to be used for hydrologic prediction and modeling.

Since knowledge of Ksat is essential for the use of water flow models, it is useful to evaluate the influence of measured Ksat on modeled runoff. One modeling approach for erosion/runoff prediction is the WEPP. This model has been extensively used for runoff prediction since 1995 when it was publicly released by the USDA-ARS. Although Ksat is not the only factor that affects runoff, the WEPP model incorporates the estimated values of hydraulic conductivity as an important soil attribute to predict runoff (Flanagan and Nearing, 1995). The WEPP uses Keff values for surface layers, and internally computes Ksat values for subsurface soil layers. Studies indicate that runoff predictions are sensitive to the initial Keff values (Ghidey et al., 1999).

The objectives of this study are to: (i) measure lateral in situ Ksat of the 0- to 230-mm depth (above the claypan) using 250 mm wide by 500 mm long soil monoliths, (ii) measure the Ksat with and without bentonite of 76-mm diam. soil cores taken at 100-mm intervals to a depth of 2 m, (iii) compare the core Ksat vs. lateral in situ Ksat, (iv) compare previously measured Ksat data for Missouri claypan soils with results of this study, and (v) compare measured runoff vs. WEPP predicted runoff using measured Keff as a model input to illustrate the benefit of using measured Keff values.


    MATERIALS AND METHODS
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Installation of Soil Monoliths for Lateral in situ Ksat Measurements
This study was conducted at the Midwest Research Claypan Farm (McCredie) near Kingdom City, MO. A 10 by 30 m area under long-term (>28 yr) continuous fescue (Festuca arundinacea Schreb.) and blue grass (Poa pratensis L.) was chosen. The soil is a Mexico silt loam (fine, smectitic, mesic, Aeric Vertic Epiaqualf) formed in loess developed over glacial till on a slope of about 3%. The mineralogy of the argillic horizon consists of 38% montmorillonite, 34% quartz, 21% kaolinite, and 7% illite. Selected soil properties are presented in Table 1.


View this table:
[in this window]
[in a new window]
 
Table 1. Soil characteristics and land use of the hydraulic conductivity studies.

 
Three in situ monoliths 250 mm wide by 500 mm long by 230 mm deep (the depth to claypan) were constructed. Depth to the claypan was determined by obtaining 20-mm diam. soil samples using a hand probe. These samples were studied in the field for changes in texture and color to determine the depth. Mean (± SD) depth was 230 ± 7 mm (n > 100). Each monolith site was wet for 12 h and then allowed to drain 24 h to soften the soil sufficiently while reducing puddling during construction. A trench was dug around each monolith to form a rectangular soil block with intact bottom.

The monolith set-up had three compartments: (i) water supply pit, (ii) soil monolith, and (iii) water collection pit. Two 6.3-mm steel plates 500 by 700 mm were installed vertically along the two sides of the monolith, which allowed a 100-mm of additional length of steel plate at each end. The additional length allowed for construction of a water supply pit and a water collection pit. Two plates 250 by 500 mm were installed vertically at each end of the border to form a rectangular box. The steel plates extended from 150 mm above to 350 mm below the soil surface. Silicone caulking was used to waterproof the steel plate seams. A 20-mm discharge hole was made at the lower end of the collection pit. A divider screen was made from a metal screen with geotextile material, which separated the soil monolith from the collection and supply pits. A bentonite-slurry was used to seal the soil–steel plate interfaces. The excavated trench was backfilled with the original native soil.

Measurement of Lateral in situ Ksat
Monoliths were slowly wet for 48 h. The electrical conductivity (EC) of the in situ water used was 0.71 dS m-1, and the Na adsorption ratio (SAR) of the in situ water was 2.39. Once the monoliths were satiated, water was added to the supply pit using a Mariotte bottle for maintaining a constant head and measuring the inflow rate. When the water level rose to the soil surface in the collection pit, excess water flowed through the 20-mm discharge hole. Plastic tubing routed the outflow for measurement.

The monolith lateral Ksat was measured by applying water from the water supply pit and measuring outflow in the collection pit for 5 h. A difference in hydraulic head of 16 mm was measured along the in situ pedon. A polyethylene tent was used to cover the plot throughout the measurement to minimize water loss from evaporation. Time, inflow and outflow volumes, and hydraulic gradient were recorded to facilitate calculation of the lateral Ksat.

The time to steady flow conditions was 48 h. After 12 h of satiation, 18.9% of applied water was moving downward through the soil. Downward movement decreased to 4.8% of total inflow after 24 h. This continued to decrease to 1.5% when the plots were satiated for 48 h. Downward flow through the claypan was obtained by subtracting outflow from the inflow. Measurements were initiated when downward flow was 1.5%.

Laboratory Ksat
One hundred eighty soil cores were taken within 10 m of the in situ study sites, to determine the Ksat distribution with depth, and to facilitate comparison of lateral in situ Ksat with Ksat determined on small intact cores. Nine intact 76-mm diam. soil cores were collected every 100 mm with depth to 2 m using a core sampler (Blake and Hartge, 1986). A replicate area near each monolith was used to collect cores in a vertical orientation when the soil was slightly below field capacity. Samples were transported to the laboratory, and slowly wet from the bottom with tap water using a Mariotte bottle having a supply rate of about 3 mm h-1. The EC of the tap water used was 0.68 dS m-1, and the SAR of the tap water used was 2.34. Cores collected above the 200-mm depth were wet for 24 h. Cores collected at or below the claypan, were wet for 7 d. Measurements for samples with higher Ksat were determined with a constant head, and those with low Ksat were determined with a falling head (Klute and Dirksen, 1986).

Visible macropores (>=1 mm) and interfacial voids located between the soil and the cylinder wall on a set of the cores were plugged using a bentonite-slurry. The reason for using this slurry was to eliminate the free flow of water through these macropores and voids. Elimination of bypass flow in small cores during Ksat determinations is a recommended methodology (Smith and Browning, 1946; Klute, 1965; Fadl, 1979). Blocking of macropores may seem at odds with the goal of estimating in situ Ksat, which measures flow through the naturally occurring macropores. However, a problem arises when small, 76-mm cores are used for Ksat measurements. Macropore continuity in field conditions is intact while this continuity is broken in small cores. These macropores are commonly finite in small cores and are often rapid pathways for bypass flow because of differences in the boundary conditions between in situ and core measurements. The dominant saturated flow in small cores is mainly via these macropores rather than through the soil matrix.

The Ksat values with bentonite were compared with those measured without bentonite to assess the effectiveness of the bentonite. A t-test was used to examine the hypotheses that the lateral in situ Ksat and laboratory Ksat determinations of the topsoil were not different by assuming anisotropic conditions (SAS Institute, 1985). This assumption is well supported by studies, which indicate that Ksat within the plow layer of silt loam soils is not appreciably influenced by core orientation (Dabney and Selim, 1987).

Comparison of Existing Ksat Data for Missouri Claypan Soils
The consistency of available Ksat data for the claypan soils was evaluated by comparing previously collected Ksat data with the results from this study. Data are based on studies of Jamison and Peters (1967), Doll (1976), McGinty (1989), Zeng (1994), and Baer and Anderson (1995). Soil characteristics and land use of the study sites are in Table 1. For data from the current study, a 95% confidence interval of the mean was calculated using the pooled variance of the Ksat with and without bentonite data of each depth separately.

Runoff Prediction Using Existing Ksat Data as Input for the WEPP Model
The study of Ksat influence on runoff prediction was conducted by using the WEPP Hillslope model (Version 98.4) on a single event basis using the input of Keff. The Keff input values for WEPP runoff prediction were determined using the Ksat measured on 76 by 76 mm diam. soil cores. The Keff was calculated as:

[1]
where LT is the total thickness of the 0- to 300-mm depth; L1, L2, and L3 are layer thickness values, and K1, K2, and K3 are the Ksat values for each of the three depth intervals (0–100, 100–200, and 200–300 mm; Jury et al., 1991).

The Keff was computed for the horizons within the upper 0 to 300 mm because this depth has soil that is much more permeable (Ksat {approx} 71 mm h-1) than the underlying very slowly permeable argillic horizon (Ksat {approx} 1.83 µm h-1). Hence, the topsoil Ksat would largely control water flow in saturated conditions. The 0- to 300-mm depth reflects the inherent soil properties of this permeable soil. The best approach for Keff estimation would likely be to evaluate soil properties with depth on a case-by-case basis and allow the soil profile to direct the depth chosen for Keff estimation. However, this approach may be too costly and time-consuming for routine use.

The procedure used to compute Keff is different from that estimated internally by WEPP which predicts Keff based on approximate relationships with soil properties (Zhang et al., 1995). The predicted Keff determined by WEPP is optimized using measured runoff data from a database derived from multiple plots for various soil types. For instance, the Keff for the surface soil of the Mexico claypan soil is 0.34 mm h-1 (Nearing et al., 1996). The WEPP estimate of Keff is useful when measured Ksat data are not available. Because we had measured Ksat from five studies on Missouri claypan soils, we computed Keff for each study to evaluate the effect of Ksat variability on predicted runoff.

The computed Keff was used as an input parameter while using other WEPP input parameters as reported by Ghidey and Alberts (1996) for Missouri claypan soils. The only input parameter that was changed in this study was the Keff for the surface 300-mm depth. Below this depth, the WEPP internally predicted Keff values were used for runoff prediction (Zhang et al., 1995). The WEPP model requires four input files containing information on climate, slope, soil, and crop management to estimate runoff (Flanagan and Nearing, 1995). Ghidey and Alberts parameterized the required input files of WEPP Hillslope Model (Ver. 95.7) using measured runoff and soil data from long-term runoff-erosion plots at (McCredie) Kingdom City, MO.

The WEPP predicted runoff was compared with measured runoff data collected from the natural rainfall erosion plots located at the Midwest Research Claypan Farm (Ghidey and Alberts, 1996). The runoff-erosion plots were managed in no-till corn (Zea mays L.) for an 11-yr period (1983–1993). Only the 11 largest rainfall events were selected for study when runoff was likely to occur. Data from the no-till corn plots were used because these plots had the most protective crop residue (~95% residue cover), and thus Ksat would not be greatly reduced by surface seal from rainfall. Table 2 indicates that prior to the reported dates of largest rainfall event the soil was practically saturated in 1984, 1985, 1988, and 1993, and near saturation in 1983, 1986, 1987, and 1989 through 1992. Based on these data, the use of Keff using measured Ksat for runoff prediction was considered appropriate.


View this table:
[in this window]
[in a new window]
 
Table 2. Rainfall and soil characteristics for 11 large rainfall events from 1983 to 1993.

 
The rate of change in predicted runoff as influenced by the Keff change was quantified by performing a sensitivity analysis of WEPP predicted runoff. The sensitivity index for the WEPP results was computed as described by Lane and Nearing (1989):

[2]
where R1 is the predicted runoff using measured Keff with bentonite, R2 is the predicted runoff using Keff selected from the other studies, R12 is the average of R2 and R1, Keff1 is the Keff with bentonite, Keff2 is the Keff selected from the other studies, and Keff12 is the average of Keff2 and Keff1.


    RESULTS AND DISCUSSION
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Lateral in situ Ksat
In situ Ksat values were calculated by assuming that the argillic horizons were nearly impermeable when saturated. Mean lateral in situ Ksat was 72 ± 0.7 mm h-1. Differences were not significant among sites (P > 0.5). About 98.5% of applied water in the upper end of the monolith moved laterally through the soil layer above the restrictive argillic horizons after 48 h of wetting. A small amount (1.5%) was unaccounted for and likely was downward flow through the claypan equivalent to a Ksat of ~9 µm h-1. Results suggest that the argillic horizons were a barrier directing the vertical flow horizontally above the claypan as lateral flow. A perched water table is thus likely as these soils are ponded for several hours. These results support earlier findings of Jamison and Peters (1967), and Saxton and Whitaker (1970), who reported the occurrence of lateral flow in these soils.

Comparison of Ksat Determined on in situ Monoliths and Intact Cores
The Ksat values for cores without bentonite were significantly larger than for the monoliths (P = 0.033). The mean value was four times more (312 ± 58 mm h-1) than for the monoliths' Ksat values (72 ± 0.7 mm h-1). The difference between the Ksat of the monoliths and the Ksat of intact cores with bentonite (71 ± 1.1 mm h-1) was not significant (P = 0.50). Inspection of the cores showed numerous vertically oriented macropores produced by flora and faunal biological activity (biochannels). Water flow through cores without bentonite was largely governed by flow through macropores extending throughout cores, and resulted in Ksat values that were unrealistically high. Under field conditions, such macropores would be expected to terminate in the subsoil because of the decrease in porosity with depth, soil swelling, and from clay hydration, and are thus much less conductive when satiated than continuous open-ended macropores in cores. Since the Ksat of cores with bentonite is similar to in situ values, the bentonite technique may be used to approximate laboratory Ksat to the Ksat of in situ soils. Small cores may not reflect the in situ Ksat values if continuous macropores that dominate flow in small cores are not eliminated.

The macropore effect on water flow is also a function of the pore orientation. A macropore extending vertically throughout a 76-mm core causes a higher Ksat value compared with values under field conditions. In contrast, a laterally oriented macropore in a core conducts less or no water because free water will not enter the pore, thus reducing the Ksat value (Hillel, 1998). This study found that cores having macropores visible at the exposed surface that were oriented vertically produced four times larger Ksat values compared with in situ measured Ksat, whereas Ksat values measured with bentonite injected to eliminate this effect were not statistically different from in situ measured Ksat values.

The Ksat Profiles with Depth of Intact Cores
Profile plots of Ksat are shown in Fig. 1 . Data show that Ksat with bentonite was significantly lower (P = 0.007) than Ksat without bentonite throughout the profile. Measurement of Ksat on cores without bentonite had higher conductivities even for samples within the claypan with high montmorillonitic clay content. This high clay content is commonly thought to increase swelling and thus close macropores reducing Ksat values. However, the measured data suggest this notion is not correct. The mean Ksat without bentonite (312 ± 58 mm h-1) is four times greater than Ksat with bentonite (71 ± 1.1 mm h-1) for the surface 100 mm of soil. The comparison of Ksat of cores with and without bentonite indicates that ~90% of water flow through cores from the upper 100 mm of the soil can be conducted by the macropores. McGinty (1989) also measured the Ksat of 76-mm diam. cores without bentonite on claypan soils and found high Ksat values for the surface soil (333 mm h-1). This work was done on soil samples collected from no-till sites where some macropores were present and very likely were not closed, and thus conducted water very rapidly.



View larger version (23K):
[in this window]
[in a new window]
 
Fig. 1. Comparison of saturated hydraulic conductivity (Ksat) data from selected studies for Missouri claypan soils. The error bars represent the 95% confidence interval of the mean Ksat value for each depth (n = 9).

 
Differences in the Ksat profile (with bentonite) among the three sampling sites across a depth of 2 m were not significantly different (P = 0.77). A significant variation in Ksat with depth occurred (P = 0.001). The lowest conductivities (2.2–1.8 µm h-1) were between the 550- and 750-mm depth, correlating to soil with weakly developed, compact, and firm structure. This layer corresponds to the region immediately below the claypan (the zone of maximum clay accumulation). The Ksat measured at the 100-mm depth was about 40000 times greater than that found at the 600-mm depth (1.8 x 10-3 mm h-1). An increase is noted from the 600- to 950-mm depth likely because of the textural change from silty clay to silty clay loam (Bohnert, 1967).

Comparison of Ksat Determinations for Missouri Claypan Soils
Figure 1 shows the Ksat measured on selected Missouri claypan soils. The Ksat decreases with depth because of changes in soil density, texture, and structure. The Ksat values with bentonite were nearly 10 times greater than the Ksat values measured by Doll (1976), and Baer and Anderson (1995). The Ksat values were nearly five times less than those reported by McGinty (1989) for the upper 200-mm depth. The Ksat values reported by Zeng (1994) on 76-mm cores were 1.3 times greater than the Ksat with bentonite.

The variation in Ksat presented in Fig. 1 is mainly attributable to (i) the variations in depth to the claypan among the studies, (ii) the presence of conductive macropores, and (iii) the method of Ksat determination. First, the depth to claypan varies between 130 and 370 mm with an average of 250 mm (Jamison and Peters, 1967). Samples taken by previous investigators from different sites at the same depth may also have differed in clay content and bulk density, altering Ksat values (Table 1). The low Ksat values found by Baer and Anderson (1995) for example, may be explained because their samples were collected from severely eroded soil that had exposed the claypan. Secondly, claypan soils often have abundant macropores (Jordan et al., 1997). Small cores may overestimate Ksat values particularly for surface depths with abundant macropores. This is shown in Fig. 1 where the Ksat without bentonite is about four times higher than that the Ksat with bentonite. The Ksat values by McGinty (1989) are higher than the KSat values with bentonite because his measurements were made without bentonite, and had large macroporosity (about 3–5% porosity in the size range of 1- to 2-mm diam.). Thirdly, the Ksat variation may be due to different methods and different aspects of measurement (with vs. without macropores). For example, Jamison and Peters (1967) determined the Ksat with the double tube method, Doll (1976) used the crust method, and the core Ksat in this study was measured with and without bentonite.

Influence of Ksat on Modeled Runoff Prediction
Process-based hydrologic models require input of Ksat. However, model users often have limited access to measured data and thus use published or estimated values. Studies of claypan soils indicate that Ksat values may vary by 100 times due in part to spatial and temporal variability (Fig. 1). This variability in input Ksat has the undesirable effect of producing variable and inaccurate model predictions.

The impact of Ksat variability on runoff was evaluated by performing the WEPP runoff prediction using measured Ksat from selected studies for the Missouri claypan soils. The Keff values for the surface 300-mm depth are: Blanco (with bentonite) = 1.3 mm h-1, Doll = 2.7 mm h-1, Blanco (without bentonite) = 3.4 mm h-1, Baer and Anderson = 5.4 mm h-1, Zeng = 8.2 mm h-1, and McGinty = 183 mm h-1. Prediction results show runoff to vary greatly in response to changes in Keff input. The Keff of the other studies is significantly higher than the Keff with bentonite and underestimated the observed runoff. As expected, higher Keff values produce lower predicted runoff. Figure 2 compares the measured runoff with the predicted runoff using the Keff with and without bentonite (McCredie, MO), and the highest Keff (Novelty, MO) reported by McGinty. The effect on WEPP predicted runoff of using an Keff without bentonite measured at 40-mm runoff was 29 mm versus about 39 mm when using an Keff with bentonite. This indicates that the use of Keff without bentonite underestimated the runoff by 28% at a measured runoff of 40 mm. Use of Keff value calculated from cores with bentonite most closely correlated with the observed runoff (Fig. 2). This is attributed to the fact that Ksat measured with bentonite excluded macropore flow through continuous macropores in small cores and thus better reflected the in situ conditions where the water flow in macropores is reduced when the soil is saturated.



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 2. Comparison of Water Erosion Prediction Project (WEPP) predicted and measured runoff using effective hydraulic conductivity (Keff) with bentonite = 1.3 mm h-1 (McCredie, MO), Keff without bentonite = 3.4 mm h-1 (McCredie, MO), and Keff without bentonite = 183.6 mm h-1 (Novelty, MO) as input values.

 
Figure 3 shows the relationship between WEPP predicted using Keff with bentonite as input versus measured runoff. It was expected that a linear relationship would be found but there was a significant quadratic relationship (r2 = 0.95). The quadratic behavior is probably due to (i) spatial and temporal variability of Keff, (ii) dependence of Keff on rainfall amount and intensity, (iii) variable satiated initial conditions, and (iv) effect of the underlying argillic horizons on runoff. Three points in the WEPP predicted runoff emerge in Fig. 3. The WEPP model (i) overpredicted the runoff in the low range, (ii) underpredicted runoff in the medium range, and (iii) performed well in the high range of measured runoff. These results highlight that selection of Ksat has a great impact on runoff prediction. Previously collected Ksat data cited above all underpredicted runoff when compared with measured values.



View larger version (22K):
[in this window]
[in a new window]
 
Fig. 3. Relationships between Water Erosion Prediction Project (WEPP) predicted runoff and measured runoff, using effective hydraulic conductivity (Keff) as input computed from saturated hydraulic conductivity (Ksat) with bentonite determined on small 76-mm soil cores.

 
This highlights the need for researchers to use caution when using Ksat data as model input without field validation. Model users need to consider both the variability of Ksat data associated with a specific soil location and understand how the method of determination may influence its value. The Ksat data determined on small soil cores with bentonite predicted runoff satisfactorily, indicating that the use of bentonite to plug macropores is advisable. The core Ksat values with bentonite were not significantly different from in situ Ksat values because the bypass flow through the visible pores (>1 mm diam.) in the small cores was eliminated. The Keff values without bentonite were 160% higher than Keff values with bentonite. The use of bentonite was useful to approximate the core Ksat to in situ Ksat.

The sensitivity index for the WEPP results reflects the change in runoff with respect to change in Keff. It was the greatest (0.25) for the highest Keff reported by McGinty indicating that for a 100% increase in Keff, runoff would be increased by 25%. The sensitivity values of other studies were: 0.10 for Zeng, 0.08 for Baer and Anderson, 0.07 for Blanco (without bentonite), 0.05 for Doll, and 0.04 for Blanco (with bentonite). This last sensitivity value was obtained using in situ Ksat values for comparison. Predicted runoff was sensitive to changes in Keff, indicating that Ksat is a critical parameter for obtaining accurate runoff estimates (Fig. 2). Indeed, Flanagan and Nearing (1991) stated that hydraulic conductivity is one of the most sensitive soil input parameters in predicting runoff. Consequently, model users should be cautious in using estimated Ksat without proper evaluation of its accuracy.

As suggested by Kutilec and Nielsen (1994), use of a pedotop-scale model will likely improve Ksat estimation for use as model input by accounting for some of the spatial variability in Ksat (the pedotop-scale consists of a surface area overlying similar soil REV units typically totaling from 100 to 1000 m2 in size). Laboratory Ksat of small cores should only be regarded as a rapid approximation of field conditions rather than a representative measure of the pedotop Ksat. Careful consideration of measurement method, presence of biochannels, natural variations in soil depth, and use of pedotop scaling should be pursued to improve the Ksat estimation.


    ACKNOWLEDGMENTS
 
This study was supported by the Latin American Scholarship Program of American Universities (LASPAU) and the Fulbright Program. Additional support was provided by the Missouri Agric. Exp. Stn. Proj. MO-NRSL0260 and the USDA-ARS Cropping Systems and Water Quality Research Unit. Special thanks are extended to Drs. A. Hjelmfelt, A. Thompson, and C. Baffaut for their critical reviews.


    NOTES
 TOP
 NOTES
 ABSTRACT
 INTRODUCTION
 MATERIALS AND METHODS
 RESULTS AND DISCUSSION
 REFERENCES
 
Contribution of the Missouri Agric. Exp. Stn. Journal. No. 13159

Received for publication July 25, 2001.


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




This article has been cited by other articles:


Home page
Agron. J.Home page
H. Blanco-Canqui and R. Lal
Axle-Load Impacts on Hydraulic Properties and Corn Yield in No-Till Clay and Silt Loam
Agron. J., November 7, 2008; 100(6): 1673 - 1680.
[Abstract] [Full Text] [PDF]


Home page
Journal of Soil and Water ConservationHome page
R.N. Lerch, E.J. Sadler, N.R. Kitchen, K.A. Sudduth, R.J. Kremer, D.B. Myers, C. Baffaut, S.H. Anderson, and C.-H. Lin
Overview of the Mark Twain Lake/Salt River Basin Conservation Effects Assessment Project
Journal of Soil and Water Conservation, November 1, 2008; 63(6): 345 - 359.
[Abstract] [PDF]


Home page
Journal of Soil and Water ConservationHome page
S. Kumar, S.H. Anderson, L.G. Bricknell, R.P. Udawatta, and C.J. Gantzer
Soil hydraulic properties influenced by agroforestry and grass buffers for grazed pasture systems
Journal of Soil and Water Conservation, July 1, 2008; 63(4): 224 - 232.
[Abstract] [PDF]


Home page
J. Environ. Qual.Home page
H. V. Kazemi, S. H. Anderson, K. W. Goyne, and C. J. Gantzer
Atrazine and Alachlor Transport in Claypan Soils as Influenced by Differential Antecedent Soil Water Content
J. Environ. Qual., June 23, 2008; 37(4): 1599 - 1607.
[Abstract] [Full Text] [PDF]


Home page
Agron. J.Home page
P. Jiang, N. R. Kitchen, S. H. Anderson, E. J. Sadler, and K. A. Sudduth
Estimating Plant-Available Water Using the Simple Inverse Yield Model for Claypan Landscapes
Agron. J., May 7, 2008; 100(3): 830 - 836.
[Abstract] [Full Text] [PDF]


Home page
Crop Sci.Home page
D. B. Myers, N. R. Kitchen, K. A. Sudduth, R. E. Sharp, and R. J. Miles
Soybean Root Distribution Related to Claypan Soil Properties and Apparent Soil Electrical Conductivity
Crop Sci., July 30, 2007; 47(4): 1498 - 1509.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
T. Seobi, S. H. Anderson, R. P. Udawatta, and C. J. Gantzer
Influence of Grass and Agroforestry Buffer Strips on Soil Hydraulic Properties for an Albaqualf
Soil Sci. Soc. Am. J., May 6, 2005; 69(3): 893 - 901.
[Abstract] [Full Text] [PDF]


Home page
Agron. J.Home page
K. A. Nelson, P. P. Motavalli, and M. Nathan
Response of No-Till Soybean [Glycine max (L.) Merr.] to Timing of Preplant and Foliar Potassium Applications in a Claypan Soil
Agron. J., April 27, 2005; 97(3): 832 - 838.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
A. Rachman, S. H. Anderson, C. J. Gantzer, and A. L. Thompson
Influence of Stiff-Stemmed Grass Hedge Systems on Infiltration
Soil Sci. Soc. Am. J., November 1, 2004; 68(6): 2000 - 2006.
[Abstract] [Full Text] [PDF]


Home page
J. Environ. Qual.Home page
R. P. Udawatta, P. P. Motavalli, and H. E. Garrett
Phosphorus Loss and Runoff Characteristics in Three Adjacent Agricultural Watersheds with Claypan Soils
J. Environ. Qual., September 1, 2004; 33(5): 1709 - 1719.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
H. Blanco-Canqui, C. J. Gantzer, S. H. Anderson, E. E. Alberts, and A. L. Thompson
Grass Barrier and Vegetative Filter Strip Effectiveness in Reducing Runoff, Sediment, Nitrogen, and Phosphorus Loss
Soil Sci. Soc. Am. J., September 1, 2004; 68(5): 1670 - 1678.
[Abstract] [Full Text] [PDF]


Home page
Soil Sci.Home page
H. Blanco-Canqui, C. J. Gantzer, S. H. Anderson, and E. E. Alberts
Tillage and Crop Influences on Physical Properties for an Epiaqualf
Soil Sci. Soc. Am. J., March 1, 2004; 68(2): 567 - 576.
[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 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 Web of Science (22)
Right arrow Citing Articles via Google Scholar
Google Scholar
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
Right arrow Search for Related Content
PubMed
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
GeoRef
Right arrow GeoRef Citation
Agricola
Right arrow Articles by Blanco-Canqui, H.
Right arrow Articles by Ghidey, F.
Related Collections
Right arrow Water Management
Right arrow Soil Models


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