| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
|
|
|
|||||||
|
|||||||||
a Dep. of Crop and Soil Environmental Sciences, Clemson Univ., Clemson, SC 29634 USA
b Dep. of Agronomy, University of Kentucky, Lexington, KY 40546 USA
vqsnbrr{at}clemson.edu
 |
ABSTRACT |
|---|
 TOP ABSTRACT INTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
), total porosity, the soil water retention curve, and saturated hydraulic conductivity were determined on the same blocks. Volumetric water content increased
as the slope of the water retention curve decreased. Increasing resulted in decreasing v
and thus, because of the linear relationship between D and v
, decreasing D. Structural controls on solute dispersion in this study were mainly indirect, and related to variations in water content produced by differences in pore-size distribution.
Abbreviations: ANOVA, analysis of variance • CDE, convection dispersion equation • MIM, mobile–immobile model • TDR, time-domain reflectometry
|
INTRODUCTION |
|---|
|
TOPABSTRACTINTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
Quisenberry et al. (1994) found that more than one-half of the total drainage occurred through just 12 to 19% of the soil cross-sectional area. For any given leaching event, surface-applied solutes can elute from constantly changing locations in the soil profile. Flow path stability may be related to soil structural stability. Tillage practices influence soil structure and transport processes. Tillage destroys the natural pore structure of surface soils, disrupting macropore continuity, and reducing the extent of bypass flow. Increased infiltration rates in no-till fields have been attributed to the greater number and continuity of macropores in the surface layer compared with moldboard-plowed soils (Edwards et al., 1988; Dick et al., 1989). Reduced macropore flow may promote increased adsorption of agricultural chemicals in soil. This could delay the onset of groundwater contamination from individual storms.
Numerous studies have summarized the results of fitting transport parameters to tracer experiments (e.g., Parker and van Genuchten, 1984; Jury and Sposito, 1985; Kool et al., 1987, Beven et al., 1993). Flow through macropores can be considered part of a continuum of flow velocities that may be quantified using the classical CDE (Parker and van Genuchten, 1984). However, miscible displacement experiments conducted on undisturbed samples often reveal substantial deviations between observed breakthrough curves and those predicted by the CDE (Kay and Elrick, 1967; Nielsen et al., 1986; Jardine et al., 1988). Two-region mobile–immobile models (MIM) were developed to address these discrepancies (Rao et al., 1980a, b; van Genuchten and Dalton, 1986). Stochastic–convective models have also been used to predict the solute movement in structured soils. The most widely used of the latter type is the convective-lognormal transfer function described by Jury (1982).
Relationships between transport parameters and variables such as flow rate were reported by Biggar and Nielsen, (1976), Bowman and Rice (1986), and Dyson and White (1987). However, there is no currently accepted model for the influence of water content on solute dispersion under steady-state unsaturated flow conditions. Several researchers working with simplified glass bead systems have noted that the slope of the relationship between D and the mean v appears to be steeper under unsaturated conditions than under saturated conditions (e.g., De Smedt and Wierenga, 1984; Matsubayashi et al., 1997; Haga et al., 1999). These studies imply that the dispersivity
depends on the degree of saturation. In contrast, the results of Yule and Gardner (1978) suggest that dispersivity is independent of water content and is determined exclusively by pore space geometry. Yule and Gardner (1978) worked with packed beds of Plainfield sand (mixed, mesic Typic Udipsamment). They found a wide range in dispersivity at any given flow rate. The variability in estimates of D increases significantly with undisturbed samples from structured soils (Beven et al., 1993). Given such variability, it is unlikely that the relationship between D and v can be shown to be statistically different under saturated and unsaturated flow conditions for heterogeneous porous media.
Relatively few experimental studies have been conducted in which the solute breakthrough curve and structural characteristics were determined on the same undisturbed samples, thereby permitting the establishment of an empirical relationship between dispersivity and pore-size distribution. Anderson and Bouma (1977) observed greater Cl- dispersion in undisturbed soil samples with subangular blocky structure than in samples with a prismatic structure. Walker and Trudgill (1983) reported significant correlations between solute transport parameters and several pore geometry variables measured by image analysis of soil thin sections. Gist et al. (1990) showed that tracer dispersion in consolidated rocks was a function of the width of the pore-size distribution determined by Hg porosimetry. Recently, Vervoort et al. (1999) related dispersivities determined under saturated conditions to the width of the pore-size distribution, as inferred from the slope of the soil water retention curve (Vervoort et al., 1999).
Our study combined measurements of soil water content, structural properties, and solute breakthrough in order to investigate methods of predicting transport parameters independently. Solute transport was measured through relatively large, undisturbed soil blocks under steady-state, unsaturated flow conditions. The main objectives of the study were to investigate (i) how soil type and land-use management practices affected the solute movement and (ii) how soil water content and structural properties are related to transport parameters.
|
Materials and methods |
|---|
|
TOPABSTRACTINTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
|
10 cm from the bottom of the casing and transported the blocks to a temperature-controlled room for storage at 4°C. All blocks were covered in plastic to ensure that the soil would not dry, crack or pull away from the foam interface during storage. We trimmed the bottom of each individual soil block flush with the wood casing, placed the block on a collection chamber, and caulked it with silicon to make an airtight and waterproof seal. The top of the collection chamber was a metal grid consisting of 144 cells in a 12 by 12 array that collected drainage from the block. The collection cells were 3.05 by 3.05 cm and tapered to a 3-mm-diameter drain hole at the bottom. Nylon screens (Nitex 53-µm mesh, Sefar America, Briarcliff Manor, NY) were placed in the bottom of each cell and the cells were filled with a saturated, 100-µm-diameter glass bead phase barrier placed between the soil and the collection chamber.
The outermost row of cells collected the outflow from the soil–foam interface. All results in this study are based on effluent collected from only the 100 innermost collection cells. Therefore, any possible effects of edge flow were minimized. Plastic trays held 100 plastic centrifuge tubes (50-mL volume) beneath the drain holes of the collection cells to collect soil block drainage. Drainage from the outermost row of cells was collected in a separate tray. The procedures for installing a soil block on this collection chamber are described in detail by Quisenberry et al. (1994). A -2.0-kPa vacuum was applied to the lower boundary of each soil block via the collection chamber (Phillips et al., 1995). Based on the capillary equation (Danielson and Sutherland, 1986), this vacuum drained all pores greater than 0.15-mm equivalent cylindrical diameter.
Transport Experiments
Simulated rainfall was applied to the top of each block at a target rate of 1 cm h-1 (1056 mL h-1). The simulator was a square reservoir, 32 by 32 by 5 cm, constructed of acrylic plastic 0.32 cm thick. It was positioned 20 cm above the soil block. One hundred hypodermic needles, 0.25 mm in diameter, were used for the rainfall application (Quisenberry et al., 1994). The rainfall simulator actually applied fluxes in the range of 0.64 to 1.00 cm h-1; the average value was 0.85 cm h-1. When steady flow was achieved, the influent solution was changed from 0.003 mol L-1 CaSO4 to 0.03 mol L-1 KCl to produce a step increase in solute concentration.
Each experiment lasted 36 h. Fifteen trays of effluent were collected: 11 trays of 50-mL tubes collecting effluent from individual drain holes alternated with four bulk trays (all leachate collected en mass). The bulk trays collected the drainage water after 12, 20, 24, and 32 h. Two subsamples were taken to determine the Cl- concentration. The tube trays were changed at the end of each hour, except Trays 5, 7, and 15, which were run for a 4-h period. The trays with tubes were checked regularly for cells with high flow, and individual tubes were replaced once full. The 50-mL tubes were weighed to determine leachate volume in each cell, and this value was converted into a flux (cm h-1). The Cl- concentration (mg L-1) was measured using a Bio-Tech EL-311 microplate autoreader (Biotech Instruments, Luton, UK) by the automated ferricyanide method (American Public Health Association, 1989). Relative Cl- concentrations (C/C0) were determined as a function of time for each block by dividing the flux averaged effluent Cl- concentration (C) by the concentration of the influent Cl- (C0).
Time-domain reflectometry was used to measure and BEC in situ during each experiment. The TDR method is rapid and nondestructive (Baker and Allmaras, 1990). By measuring changes in BEC with time in response to a change in influent solution concentration, it is possible to determine a resident concentration breakthrough curve at exactly the same location where is measured (Hart and Lowery, 1998). Each soil block was instrumented with 20-cm-long TDR wave guides installed horizontally at the 5-, 15-, and 25-cm depths. The wave guides were connected to a Tektronix 1502C cable tester (Tektronix, Beaverton, OR) via a multiplexer. Data were collected automatically, and the waveforms were analyzed for and BEC with software developed by Wraith et al. (1993). The s within the soil blocks were extremely stable for the course of the 36-h experiment and were therefore averaged over time for each depth. Resident concentration breakthrough curves were calculated from the BEC data by the following relationship:
 | (1) |
Once the breakthrough data were obtained, they were parameterized using the CDE. For the transport of Cl-, the two principal fitting parameters of the CDE are D and v. These parameters are often combined to give the dispersivity (
), defined as
. The dimensionless column Peclet number (van Genuchten and Wierenga, 1986), defined as
(where L is the depth where the breakthrough curves were measured), was also calculated.
The CXTFIT computer program (Parker and van Genuchten, 1984; Toride et al., 1995) was used to estimate the solute transport parameters of the soil blocks. CXTFIT has an inverse capability that can estimate D and v from the observed results of a solute transport experiment. The program uses a nonlinear least squares fitting procedure to optimize the fit between the observed data and the CDE. Flux averaged relative concentrations were used as inputs in the CXTFIT program for the effluent Cl- data. For the TDR resident concentration data, BEC in Eq. [1] was the input, with BECinitial specified and BECmax treated as a fitting parameter. All data were fitted in terms of the equilibrium CDE, with CXTFIT 2.0 in Mode 1 (Mode code = 1 was deterministic equilibrium CDE).
Soil Structural Characterization
After each experiment, the soil blocks were drained overnight, removed from the collection chamber and triplicate soil cores (5.7 cm in diam. by 6 cm long) were taken from the 5-, 15-, and 25-cm depths. These cores were used to determine the bulk density and the soil water retention (drying) curve (Klute, 1986). Tempe cells were used to collect the retention data between tensions of 0 and 3 kPa, and pressure plate apparatuses were used to collect the corresponding data between 3 and 1500 kPa (Klute and Dirksen, 1986). Total porosity was calculated from the bulk density data, assuming a particle density of 2.65 g cm-3. The soil water retention curves were parameterized using the empirical model developed by Campbell (1974). The estimated parameters in this model are 
a, the air-entry value, and b the pore-size distribution index. Coefficients of determination (R2) for the fits ranged from 0.77 to 0.99.
The cores were also used to determine the saturated hydraulic conductivity (Ksat). These measurements were carried out using the constant head method, as described in Klute and Dirksen (1986). Saturated hydraulic conductivity data were collected at each depth for sod Block 3 and tilled Block 3 of the Maury soil and for all four Cecil soil blocks.
Data Analysis
Geostatistical analyses were performed by computing indicator variograms (1 = flowing, 0 = nonflowing) for the individual fluxes for each block. The results indicated little or no spatial structure, so the data were analyzed using conventional statistical procedures, including analysis of variance (ANOVA) and regression analysis (SAS Institute, 1988). The overall statistical design was: two soils (Maury, Cecil), two land-use treatments (conventional-till, pasture), three depths (5, 15, 25 cm), and three replicates. Data for the third replicate of the Cecil soil were treated as missing values. For the effluent transport parameters depth was not a factor. Three determinations were made per replicate for the soil structural properties; these values were averaged prior to performing the ANOVAs. Statistical significance was assessed at the P < 0.01 and P < 0.05 probability levels.
Correlation analyses were performed for each block between the spatial distribution of fluxes in Tray 1 (Hour 1) and the spatial distribution in Tray 15 (Hour 36) for both Maury and Cecil soils.
Linear regression analyses were performed on all of the data in order to evaluate relationships between the soil water contents, structural characteristics and solute transport parameters. Only those relationships that were significant at P < 0.05 or greater will be reported.
|
Results and discussion |
|---|
|
TOPABSTRACTINTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
|
ranged from 0.2 to 0.35 cm3 cm-3 for the Maury blocks and from 0.25 to 0.36 cm3 cm-3 for the Cecil blocks (Table 2). Analysis of variance indicated significant depth and treatment (land use) effects on the water contents (P < 0.05) (Table 3) . The sod-covered blocks had significantly higher than the tilled blocks, and increased with depth for both soils (Maury and Cecil), and both land-use treatments (sod and tilled).
|
a), and pore-size distribution index (b), are summarized in Table 2. The air-entry value greatly increased with increasing depth for the tilled blocks, but increased only slightly with increasing depth for the sod-covered blocks. This trend may be due to loosening of the surface soil and subsurface compaction (i.e., a plow pan) in the tilled blocks compared with the sod-covered blocks. However, because of variability in the data (the coefficient of variation was 95%), the ANOVA for this parameter was not significant (Table 3).
In contrast to a, the b parameter was much more predictable, with a coefficient of variation of 15.2%. Soil, depth, the soil x depth interaction, and the land-use treatment x depth interaction were all significant factors influencing b at P < 0.05 (Table 3). The Maury soil has a significantly higher mean b value than the Cecil soil. The b value increased with increasing depth for the Cecil soil in both treatments, but showed an increase with depth only for the tilled Maury soil (Table 2). The b parameter is an indicator of pore-size distribution; for large values of b, small pores dominate the total porosity, and vice versa.
Values for the saturated hydraulic conductivities ranged from 8 x 10-5 to 8 x 10-3 m s-1 for the Maury soil blocks, and from 4 x 10-4 to 5x10-3 m s-1 for the Cecil soil blocks. Due to missing values an ANOVA could not be performed on these data. However, the data for the Maury soil indicate that the tilled block had higher surface Ksat values than the sod-covered soil. The land-use effect was much less pronounced for the Cecil soil, which generally appeared to have decreasing Ksat values with increasing depth regardless of management history. This trend can probably be attributed to the marked increase in clay content with depth for this soil.
Individual Cell Fluxes
For the Maury sod blocks, an average of 33% of cells collected leachate with a mean flux of 2.67 cm h-1, while the Maury tilled blocks averaged 38% of cells collecting leachate, with a mean flux of 3.48 cm h-1. These results indicate that the sod blocks had a slightly more concentrated spatial flow pattern than the tilled blocks. The opposite trend was observed for the Cecil soil. For the Cecil tilled blocks, the mean percentage of cells collecting drainage water was 40% and the mean flux was 2.23 cm h-1. For the Cecil sod blocks, the mean percentage of cells was 50% with a slightly lower mean flux of 1.96 cm h-1.
The correlation analyses, performed for each block, between the spatial distribution of fluxes in Tray 1 (Hour 1) and the spatial distribution in Tray 15 (Hour 36) provided a measure of the temporal stability of the flow pattern (from the start of the experiment to completion). Although the mean correlation coefficients r values) were not significantly different at P < 0.05 between the two Maury treatments, the tilled blocks had a lower mean r value
than the sod blocks
. The higher r value for the sod blocks indicates a more stable flow pattern from start to finish than for the tilled blocks. For the Cecil soil, although the correlation coefficients for the two land-use treatments were not significantly different (P < 0.05), the Cecil tilled blocks had a lower mean r value
than the Cecil sod blocks
. The higher mean r value for the sod blocks shows a more steady and even flow rate from start to finish, while the lower mean r value for the tilled blocks indicates a more temporally variable flow pattern.
Tillage affects soil aggregation; in the absence of a protective cover like sod, aggregates are vulnerable to disruption by flowing water (Hillel, 1980). As aggregates collapse and flow lines change, the fluxes for some cells might change with time, while the overall flow rate for the block remains fairly constant from Tray 1 to Tray 15. The sod-covered soil had an accumulation of organic matter so that aggregates in this treatment should have been more stable than in the tilled soil. Greater aggregate stability would be expected to result in more stable flow patterns with time, as indicted by the higher r values for the sod vs. tilled blocks. The lower mean r values for the Cecil soil, relative to the Maury soil, may indicate a lower aggregate stability for this soil as compared with Maury soil.
Solute Breakthrough Curves
The breakthrough and the fitted curves of representative TDR and effluent data are presented in Fig. 1 to 4
. The CDE fitted the observed breakthrough curves very well, with R2 values ranging from 0.971 to 0.999. The average R2 was 0.988. Because of the goodness of fit obtained with the CDE, the MIM model was not employed in this study.
|
|
|
|
The 5-cm TDR probe was the first one to increase after the step increase in applied water conductivity (Fig. 1–4). Subsequent increases were registered at the 15- and 25-cm probe depths. The curves take on a more sigmoid shape as the solute moves farther from the source and deeper into the soil block. For most of the cases, the effluent breakthrough curve is increasing at later times compared with the TDR breakthrough curves. An exception was the Cecil tilled soil (Fig. 3), where the 25-cm depth TDR breakthrough curve was observed later than the effluent breakthrough curve. This may be due to the relative insensitivity of TDR resident concentrations to macropore flow. Previous field and laboratory experiments, with undisturbed Cecil soils, have shown that areas >50% of the total cross-sectional area within the Bt horizon have little or no flow at all, even after two or three pore volumes of solute have been eluted (Hatfield, 1988; Nelson, 1990). The TDR probes at the 25-cm depth in Fig. 3 may well have been located within such a zone.
Transport Parameters
The D and values for the Cecil soil and sod treatment were generally greater than the corresponding values for the Maury soil and tillage treatment, respectively (Table 4)
. However, the variation around these mean values was too high (see Table 4) to establish any statistically significant differences. The ANOVA results on the transport parameters for the effluent data showed no significant effects of soil type and land use on v, D, , and Peclet number.
|
|
|
under Maury sod. Thus, the depth effects discussed above only hold for the tilled soil in the case of . Results of ANOVA for the column Peclet number identified soil (P < 0.05), depth, and the soil x depth interaction (P < 0.01) as significant effects (Table 6). The Peclet number increased with increasing depth in the Maury soil, but showed no consistent trend with depth in the Cecil soil (Table 5). This result is due to the consistent increase in D with depth observed in the Cecil soil.
Relationships between Transport Parameters and Soil Properties
The regression analyses performed on all the data evaluated relationships between the soil water contents, structural characteristics, and solute transport parameters. For the steady-state, unsaturated flow conditions employed in this study, the increased as the Campbell b parameter increased (Fig. 5)
. This is because under unsaturated conditions more water is retained by soil with many small pores (and a large b value) than soil with few small pores (and small b value) at any given flow rate. For a constant water flux, the mean v should decrease with increasing water content. The expected result is apparent in Fig. 6
, although the R2 value was quite low.
|
|
, where m and n are fitting parameters and D0 is the effective diffusion coefficient. Biggar and Nielsen (1976) reported a strong relationship for this model
, with
. Bowman and Rice (1986) found weaker relationships, with n values ranging from 1.42 to 2.2. Other studies reported that correlations between D and v can be highly variable and in some cases apparently negative (Van Ommen et al., 1989).
In this study, the relationship between D and v was approximately linear, with D increasing with increasing v (Fig. 7)
. Because of the linear relationship between D and v, D also decreased with increasing water content (Fig. 8)
, indicating increased solute dispersion at low water contents as compared with high water contents, at the same flux. This result suggests an indirect influence of b on D through the relationship between and b. However, the relationship between b and D was not significant at P < 0.05. None of the other relationships between the solute transport parameters and soil properties were statistically significant.
|
|
|
Conclusions |
|---|
|
TOPABSTRACTINTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
Statistically significant linear relationships were found between and the Campbell b parameter, v and , D and v, and D and . However, no direct relationships were found between any of the soil structural properties and D or v. We conclude that structural controls on solute dispersion in unsaturated soils subjected to initial and boundary conditions similar to those in our experiments are likely to be indirect and related to differences in water content at a given flow rate produced by differences in pore-size distribution. More direct relationships between soil structural characteristics and transport parameters are to be expected under saturated conditions.
The effects of management on preferential flow, and hence solute dispersion, were not appreciable under the conditions of these experiments. The -2 kPa lower boundary condition ensured that unsaturated conditions were maintained throughout. Under saturated conditions, macropore flow would be expected to significantly increase solute dispersion. The rainfall rate of 1 cm h-1 was insufficient to generate any runoff. The occurrence of localized surface runoff could also greatly enhance solute dispersion by channeling water into macropores.
|
ACKNOWLEDGMENTS |
|---|
Received for publication March 1, 1999.
|
REFERENCES |
|---|
|
TOPABSTRACTINTRODUCTIONMaterials and methodsResults and discussionConclusionsREFERENCES |
|---|
This article has been cited by other articles:
|
|
 |
|
  J. Vanderborght and H. Vereecken Review of Dispersivities for Transport Modeling in Soils Vadose Zone J., January 24, 2007; 6(1): 29 - 52. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 L. Alletto, Y. Coquet, P. Vachier, and C. Labat Hydraulic Conductivity, Immobile Water Content, and Exchange Coefficient in Three Soil Profiles Soil Sci. Soc. Am. J., June 21, 2006; 70(4): 1272 - 1280. [Abstract] [Full Text] [PDF]
|
|
|
|
 |
|
 I. M. Young and J. W. Crawford Interactions and Self-Organization in the Soil-Microbe Complex Science, June 11, 2004; 304(5677): 1634 - 1637. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
M. Javaux and M. Vanclooster Robust Estimation of the Generalized Solute Transfer Function Parameters Soil Sci. Soc. Am. J., January 1, 2003; 67(1): 81 - 91. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
N. E. Derby, R. E. Knighton, and B. R. Montgomery Construction and Performance of Large Soil Core Lysimeters Soil Sci. Soc. Am. J., September 1, 2002; 66(5): 1446 - 1453. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. Vanderborght, P. Gahwiller, and H. Fluhler Identification of Transport Processes in Soil Cores Using Fluorescent Tracers Soil Sci. Soc. Am. J., May 1, 2002; 66(3): 774 - 787. [Abstract] [Full Text] [PDF]
|
|
|
|
|
|
J. S. Strock, D. K. Cassel, and M. L. Gumpertz Spatial Variability of Water and Bromide Transport Through Variably Saturated Soil Blocks Soil Sci. Soc. Am. J., November 1, 2001; 65(6): 1607 - 1617. [Abstract] [Full Text] [PDF]
|
|
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |
| The SCI Journals | Agronomy Journal | Crop Science | |||
| Vadose Zone Journal | Journal of Plant Registrations | ||||
| Journal of Natural Resources and Life Sciences Education |
Journal of Environmental Quality |
||||
