Soil Science Society of America Journal 66:696-701 (2002)
© 2002 Soil Science Society of America
DIVISION S-1SOIL PHYSICS
Prediction of Dispersivity for Undisturbed Soil Columns from Water Retention Parameters
E. Perfect*,a,
M. C. Sukopb and
G. R. Haszlerc
a Dep. of Geological Sciences, Univ. of Tennessee, Knoxville, TN 37996
b Dep. of Plants, Soils and Biometeorology, Utah State Univ., Logan, UT 84322
c Dep. of Agronomy, Univ. of Kentucky, Lexington, KY 40546
* Corresponding author (eperfect{at}utk.edu)
 |
ABSTRACT
|
|---|
Dispersivity (
) is a required input parameter in solute-transport models based on the advection-dispersion equation (ADE). Normally
is obtained from miscible-displacement experiments. This dependency on inverse procedures imposes a severe limitation on our predictive capability. If solute breakthrough curves and soil hydraulic properties were measured simultaneously, pedotransfer functions could be developed to predict
from independent measurements. In this study, short (6 cm long) undisturbed columns were employed to investigate the relationship between
and the water-retention curve as parameterized by the air-entry value (
a) and Campbell exponent (b). We worked with 69 columns from six soil types ranging in texture from loamy sand to silty clay, conventional-till and no-till management practices, steady-state saturated flow conditions, and a step decrease in CaCl2 concentration from 0.009 to 0.001 M. Breakthrough curves were measured by monitoring changes in effluent electrical conductivity using a computerized data acquisition system. Estimates of
(calculated using the method of moments) ranged from 1 to 192 mm for the six soil types. Stepwise multiple-regression analysis explained
50% of the total variation in
, and indicated that dispersion increased as
a and b increased. Since both
a and b increase with increasing clay content,
also increases moving from coarse- to fine-textured soils. Our regression equation can be used as a pedotransfer function to predict
from existing databases of soil hydraulic properties. Further research is needed to independently validate its predictive capability, and to develop strategies for upscaling the model predictions.
Abbreviations: ADE, advection-dispersion equation b, Campbell exponent D, dispersive coefficient v, mean pore-water velocity
, dispersivity
, total porosity
a, air-entry value **, significant at the 0.01 probability level
 |
INTRODUCTION
|
|---|
NUMEROUS MATHEMATICAL MODELS are available for describing solute transport in soil. Of these, the ADE is the most widely used. For steady state, one-dimensional water flow, the ADE for a nonreactive solute is given by (Fried and Combarnous, 1971):
 | [1] |
where C is concentration of solute in the soil water (M L-3), t is time (T), x is distance (L), D is the dispersion coefficient (L2 T-1), and v is the mean pore-water velocity (L T-1). Equation [1] is normally applied to large volumes of homogeneous soil. Parker and van Genuchten (1984) have shown that, with appropriate boundary conditions, it is able to reproduce the highly asymmetric breakthrough curves obtained from short laboratory columns, even when continuous macropores are present. However, the use of Eq. [1] under such conditions remains a contentious issue (Germann, 1991; Feyen et al., 1998).
Because D in Eq. [1] depends upon v, it is desirable to define an alternative mixing parameter, the dispersivity
D/v (L), that can be related solely to characteristics of the porous medium (Fried and Combarnous, 1971). Dispersivity is a required input parameter in contaminant transport models based on the ADE (Zheng and Bennett, 1995). Except for a few simple systems such as packed beds of uniformly sized particles,
cannot be obtained from independent measurements (e.g., Koch and Flühler, 1993). Since most natural porous media are heterogeneous, soil physicists are currently unable to predict solute dispersion in undisturbed soil without first conducting a miscible-displacement experiment to measure it. This dependency on inverse procedures imposes a severe limitation on our predictive capability.
If solute breakthrough curves and static physical properties were determined on the same sample, empirical relations could be developed to predict
from independent measurements. However, such studies are remarkably rare, and most of them have used packed beds of disturbed media (Passioura and Rose, 1971; Han et al., 1985; Xu and Eckstein, 1997). Under saturated conditions, the magnitude of solute dispersion at any given flow rate is controlled by the pore-space geometry (Perfect and Sukop, 2001). Since it is the geometrical characteristics of solids or aggregates rather than the pore space that are measured in the packed bed approach, the resulting relationships are not directly applicable to undisturbed soil. While studies involving artificial macropores (Kanchanasut et al., 1978; Li and Ghodrati, 1997) may provide more information on the relationship between solute spreading and pore-space geometry, they are subject to similar criticisms regarding their applicability to natural systems.
In terms of heterogeneous systems, Anderson and Bouma (1977) observed greater Cl dispersion in undisturbed soil samples with subangular blocky structure as compared with those with 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. Network model simulations by Bruderer and Bernabé (2001) indicate that for a given flow rate, solute dispersion increases logarithmically as the normalized geometric standard deviation for a log-normal distribution of pores increases.
Several authors have investigated the relationship between pore structures revealed by dye-staining patterns and solute-transport parameters. Seyfried and Rao (1987) and Vervoort et al. (1999) report increasing solute dispersivity with decreasing percentage of dyed area. The percentage of dyed area can be thought of as a flow weighted measure of pore size and connectivity, with small values corresponding to high macropore connectivity and vice versa. Hatano et al. (1992) investigated the relationship between solute dispersion and the fractal geometry of dye-staining patterns. These authors conducted Cl-miscible displacement experiments on undisturbed soil columns that were later destructively sampled, and their mass and surface fractal dimensions determined by image analysis of methylene blue dye stains. An empirical equation was obtained by regression analysis relating the dimensionless Brenner number (v multiplied by column length divided by D) for Cl breakthrough to both the mass and surface fractal dimensions of the dye-stained pore space.
Pore-space geometry also determines the retention and transport of water in soil. Since hydraulic properties such as the saturated hydraulic conductivity and parameters describing the water-retention curve are more frequently measured than pore characteristics, it seems logical to seek empirical relations between such properties and
. However, very few studies have explored this line of reasoning. Vervoort et al. (1999) investigated the relationship between
and soil hydraulic properties experimentally. They reported that
increased logarithmically as the reciprocal of the slope parameter from the water-retention curve increased. Computer simulations by Vogel (2000) suggest that different pore-size distributions can result in similar water-retention curves, and that solute dispersion is more sensitive to the water-retention curve than the pore-size distribution.
Databases of soil hydraulic properties are widely available (e.g., Rawls et al., 1991; Leij et al., 1996). Thus, if solute breakthrough curves and soil hydraulic properties were measured simultaneously on samples from a wide range of soil types, pedotransfer functions (Bouma, 1989) could be developed to predict
from information stored in existing databases. Pedotransfer functions are regression equations used to predict difficult-to-obtain parameters from more easily measured soil properties. They have been widely used to predict input parameters for soil hydrological models from basic soil properties such as particle-size distribution, bulk density, and organic C content (e.g., Rawls et al., 1991). To our knowledge however, this approach has not been previously applied to predict solute-transport parameters.
The major objectives of this study were: (i) to simultaneously measure solute-transport parameters and soil hydraulic properties on undisturbed columns from a wide range of soil types and management practices, and (ii) to establish an empirical relationship (pedotransfer function) between
and soil hydraulic properties. To facilitate these goals, we performed miscible-displacement experiments with a step change in electrolyte solution concentration under saturated, steady-state flow conditions, combined with measurements of total porosity, saturated hydraulic conductivity, and the soil water retention curve.
 |
MATERIALS AND METHODS
|
|---|
Undisturbed columns were obtained from the 2- to 8-cm depth of six different soil types with an Uhland sampler (Table 1). Sampling locations were selected on the basis of textural class and management history to give a wide range of soil hydraulic properties (Table 2). The Allegheny (fine-loamy, mixed, semiactive, mesic Typic Hapludults), Bruno (sandy, mixed, mesic Typic Udifluvents), Jefferson (fine-loamy, siliceous, semiactive, mesic Typic Hapludults), and Pope (coarse-loamy, mixed, active, mesic Fluventic Dystrudepts) soil columns were all obtained from the Robinson Agricultural Experiment Station, Quicksand, KY in May 1999. The Karnak (fine, smectitic, nonacid, mesic Vertic Endoaquepts) soil columns were collected from an agricultural field near Utica, KY in July 1998.
The Maury (fine, mixed, semiactive, mesic Typic Paleudalfs) soil columns came from long-term research plots at the University of Kentucky's Spindletop Agricultural Experiment Station, Lexington, KY. The treatments sampled were conventional-till and no-till corn (Zea mays L.), with ryegrass (Secale cereale L.) as a winter-cover crop, under two N fertilizer inputs (0 and 336 kg N ha-1) (see Frye and Blevins, 1997 for details). Three subsamples were obtained from each of the four field replications. Sampling took place in May 1998, before tillage in the conventional-till treatment.
The 5.37-cm diam. by 6-cm long aluminum sleeves containing the soil columns were sealed in plastic bags, transported to the laboratory and stored at 4°C until required for analysis. Upon removal from storage, the columns were trimmed to remove any excess soil and were saturated from below with a 0.009 M CaCl2 solution. The ends of each column were covered with Nylon screens (Nitex 53-µm mesh, Sefar America Inc., Briarcliff Manor, NY) and fitted into a plexiglass Tempe cell (Soilmoisture Equipment Corp., Santa Barbara, CA), with its porous ceramic base plate removed. The top inlet on the Tempe cell was connected to a precision peristaltic pump (Technicon AutoAnalyzer Proportioning Pump III, Technicon Corp., Tarrytown, NY). The bottom outlet on the Tempe cell was connected to a plexiglass conductivity cell, with electrodes formed by the faces of two stainless steel bolts. The conductivity cell was specially designed to minimize apparatus-induced dispersion relative to that produced by the soil columns.
The experimental setup for the miscible-displacement experiments is illustrated in Fig. 1
. A steady flow of
1.4 mL min-1 of 0.009 M CaCl2 solution was pumped through the soil columns from top to bottom. A standpipe was connected to the inflow side of the Tempe cell to measure the head (
h) developed across the column (Fig. 1). If the head exceeded that available in the laboratory, a pressure transducer (Tensimeter, Soil Measurement Systems, Tucson, AZ) was attached in place of the standpipe. Saturated conditions and a constant outflow head were maintained by piping the conductivity cell discharge back up to the top of the column (Fig. 1). The head was controlled at the outlet by a drip point. Assuming no head loss over the conductivity cell, the head loss over the sample is known and was used to determine the saturated hydraulic conductivity during the solute transport experiments.

View larger version (26K):
[in this window]
[in a new window]
|
Fig. 1. Schematic of the experimental set up for the miscible-displacement experiments; h is the hydraulic head.
|
|
After
20 pore volumes, the concentration of the influent solution was changed to 0.001 M CaCl2. When breakthrough was complete, the influent solution concentration was switched back to 0.009 M CaCl2. These concentrations were selected to produce a mean ionic strength similar to field conditions (Seaman et al., 1995), and to maintain the structural integrity of the samples. Because of possible density-driven unstable mixing during a step increase in solute concentration with downward flow, only the step decrease was used to compute the breakthrough curve.
Four columns were run at a time, with the conductivity cells connected to a YSI Model 35 conductivity bridge (YSI Scientific, Yellow Springs, OH) via a multiplexer. The multiplexer was controlled by a computer, and switched to each column once every minute. The voltage output of the conductivity bridge (proportional to electrical conductivity) was similarly switched between one of four analog-digital conversion channels on a computer data acquisition board (CIO-DAS08/Jr-AO, Computer Boards, Inc., Mansfield, MA). The analog input voltages were converted to digital values by the board. Data collection was managed using a program written with the LABTECH NOTEBOOK software package (LABTECH Corporation, Andover, MA). The outputs of this program consisted of ASCII files containing the voltage over time for each cell.
The voltage data were normalized to the voltage prior to the solution step change (taken as C/C0 = 0) and to that after complete breakthrough was attained (taken as C/C0 = 1). Typical breakthrough curves are shown in Fig. 2
. The breakthrough curves were analyzed using the method of moments (Skopp, 1984; Jury and Sposito, 1985; Leij and Dane, 1991). The first, M1 [T], and second, M2 [T2], central time moments for a step increase in solute concentration were computed using the trapezoidal rule (Yu et al., 1999). The ADE parameters, v and D, were estimated from M1 and M2 using the following expressions (Yu et al., 1999):
 | [2] |
 | [3] |
where
is the length of the column [L]. By substituting Eq. [2] and [3] into our previous definition of dispersivity, an estimate of
is obtained.
Upon completion of the miscible-displacement experiments, a water-retention curve was determined for each column. For the Maury and Karnak soils, the retention measurements were performed on the undisturbed columns using hanging water columns (Klute, 1986) over the tension range 0 to 5.0 x 100 kPa, and pressure plates (Klute, 1986) from 5.0 x 100 to 1.7 x 102 kPa. Pressure plates were used in conjunction with disturbed subsamples over the range 1.7 x 102 to 1.5 x 103 kPa (Klute, 1986). For the other soils, the undisturbed columns were employed down to a tension of 1.0 x 102 kPa using the same techniques as before. A dewpoint water activity meter (Gee et al., 1992) was then used in conjunction with disturbed subsamples to obtain the water-retention curve from 1.0 x 102 to 4.8 x 105 kPa.
The water-retention curves were parameterized with the Campbell (1974) model:
 | [4] |
where S is relative saturation [L3 L-3] and
is the soil water tension [M L-1 T-2]. This model has two fitting parameters: the air-entry value,
a [M L-1 T-2], which is inversely related to the size of the largest pores present, and the dimensionless constant, b, which is directly related to the width of the pore-size distribution. Equation [4] was fitted to the experimental data using segmented nonlinear regression analysis. The mean coefficient of determination (R2) associated with these fits was >0.95. The resulting estimates of
a and b, along with the total porosity (
) and saturated hydraulic conductivity (Ksat), were used to characterize the hydraulic properties of each column.
 |
RESULTS AND DISCUSSION
|
|---|
The textural classes of the six sampled soil types ranged from loamy sand to silty clay (Table 1). This relatively wide range in soil texture, coupled with the different management practices, produced a variety of pore structures and associated hydraulic properties (Table 2). Total porosity was the least variable property measured; there was no obvious relationship between the mean
and soil texture (Table 2). In contrast, both
a and b generally increased with increasing clay content (Table 2). The saturated hydraulic conductivity varied over three orders of magnitude, with values generally increasing with increasing porosity (r = 0.41**).
The estimates of v and D for the different soils and management practices are summarized in Table 3. Because of the constant flow rate produced by the peristaltic pump, the v parameter exhibited relatively little variability, both between and within soil types. The smallest and largest values of D were obtained for the loamy sand and silty clay soils, respectively (Table 3). Within soil variability was extremely high for the D parameter, with coefficients of variation ranging from 30 to 700%.
The estimates of D generally increased as v increased (Fig. 3) . Numerous authors have reported positive correlations between D and v. The most frequently encountered models are linear and power law relationships (Perfect and Sukop, 2001). Network simulations suggest a transition from Taylor dispersion (D
v2) in homogenous porous media to mechanical dispersion (D
v) in heterogeneous porous media (Bruderer and Bernabé, 2001). In our study involving undisturbed soils, a linear model gave the best overall fit to the data (Fig. 3).
Following Elprince and Day (1977), we treated both v and D in Eq. [1] as adjustable parameters. This enabled us to evaluate any solute adsorption or exclusion effects without introducing an additional retardation factor. An effective transport porosity,
eff [L3 L-3], can be defined as (Elprince and Day, 1977):
 | [5] |
where q is the measured Darcy flux. Values of
eff <
denote solute exclusion, while values of
eff >
signify solute adsorption.
The mean values of
eff are given in Table 3. Comparison of the calculated
eff values in Table 3 with the corresponding measured values of
in Table 2 suggests that effluent electrical conductivity was a conservative indicator of the water flux in most cases. For the Allegheny, Jefferson, and Pope soils however,
eff >>
indicating retardation of ions relative to the water flux. Cation adsorption in transport experiments is well documented (e.g., Schweich et al., 1983; Kool et al., 1989). Less well known is anion adsorption, which can occur in highly weathered acid soils such as Ultisols (Bellini et al., 1996; Korom, 2000). The Allegheny and Jefferson soils are both classified as Ultisols (see Table 1), and we speculate that they had appreciable anion as well as cation-exchange capacities which, combined with the low ionic strength of the solutions employed, served to retard the movement of both cations and anions. Seaman et al. (1995) observed a similar phenomenon with highly weathered, sandy aquifer material. Use of higher ionic strength solutions, much greater than would normally be encountered in the field, is recommended to eliminate this possibility in future laboratory studies.
The mean dispersivities ranged from <0.5 cm for Bruno loamy sand to >20 cm for Karnak silty clay (Table 3). Figure 4
shows the mean
values ranked from low to high, along with their corresponding textural classes. It can be seen that dispersivity generally increased sequentially moving from the coarser to finer textural classes.
Simple linear correlations between the solute-transport parameters and soil hydraulic properties are given in Table 4. The
was the most sensitive solute-transport parameter, producing significant correlations with three out of the four soil hydraulic properties investigated. The dispersivity increased as
decreased, as
a increased, and as b increased (Table 4). The inverse relationship between
and
is to be expected, since the tortuosity of connected flow paths increases as the porosity approaches the percolation threshold (Eggleston and Pierce, 1995). The positive relationship between
and
a implies that fine-textured soils are more dispersive than coarse-textured soils. Similarly, the positive relationship between
and b means that dispersivity increases as the width of the pore-size distribution increases. A similar conclusion was reached by Vervoort et al. (1999).
The results of stepwise multiple-regression analyses relating the solute-transport parameters to soil hydraulic properties are summarized in Table 5. Once again dispersivity was the most sensitive parameter, with
a and b combining to explain
50% of the total variation in
(Table 5). Dispersivity increased as both
a and b increased. Neither
nor Ksat contributed significantly to the prediction of
once these two water-retention parameters were accounted for. Comparison of the observed and predicted dispersivities (Fig. 5)
indicated the regression model tends to over-predict small values of
, and under-predict large values of
.
The fact that multiple-regression analysis could explain only half of the total variation in
may be related to the relatively narrow range of soil textures investigated. Inclusion of additional data for sands and clays might improve the overall goodness of fit. Another source of unexplained variability is the different methods used to measure the water-retention curves for the different soils. It is also possible that the suite of predictor variables investigated (i.e.,
,
a, b, and Ksat) did not provide a sufficiently complete or sensitive description of the actual pore-space geometry.
In an example application we used the regression coefficients in Table 5 as a pedotransfer function to predict
(
pred) from the soil-water retention database of Cosby et al. (1984). These authors presented mean values of
a and b for 11 soil textural classes ranging from sand to clay. Their data are reproduced in Table 6 along with our predicted dispersivities. All of the predictions were physically reasonable (i.e., no negative
values). Furthermore, there was a clear trend towards increasing dispersivity with increasing clay content. The range in
pred for the 11 textural classes was 12 cm (Table 6). Given Fig. 5, this value is probably somewhat lower than the actual range in dispersivities that might be encountered.
A potential restriction on the use of our pedotransfer function is the dependence of
upon the measurement scale. Estimates of
tend to increase as the size of the solute transport experiment increases (Neuman, 1990; Gehlar et al., 1992). In this context, we note that all of the predictions in Table 6 are for 6-cm long columns. Further research is needed to develop appropriate strategies for upscaling these values to larger soil volumes.
 |
CONCLUSIONS
|
|---|
We have developed a pedotransfer function for predicting the dispersivity of 6-cm long soil columns from the Campbell (1974) water-retention parameters:
a and b. The model explained
50% of the observed variability in
. Our analysis suggests that
increases as
a and b increase. Since both
a and b increase with increasing clay content,
also increases moving from coarse- to fine-textured soils. In an example application, we used the pedotransfer function to predict
for the 11 textural classes in the soil water-retention database of Cosby et al. (1984). The predictions ranged from 0.8 cm for sands to 12.8 cm for clays. The actual range in
at this scale is likely to be somewhat wider. Research is currently underway on refining and independently validating this pedotransfer function. We are also developing strategies for upscaling the model predictions.
Received for publication April 24, 2001.
 |
REFERENCES
|
|---|
- Anderson, J.L., and J. Bouma. 1977. Water movement through pedal soils. I saturated flow. Soil Sci. Soc. Am. J. 41:413418.[Abstract/Free Full Text]
- Bellini, G., M.E. Sumner, D.E. Radcliffe, and N.P. Qafoku. 1996. Anion transport through columns of highly weathered acid soil: Adsorption and retardation. Soil Sci. Soc. Am. J. 60:132137.[Abstract/Free Full Text]
- Bouma, J. 1989. Using soil survey data for quantitative land evaluation. Adv. Soil Sci. 9:177213.
- Bruderer, C., and Y. Bernabé. 2001. Network modeling of dispersion: Transition from Taylor dispersion in homogeneous networks to mechanical dispersion in very heterogeneous ones. Water Resour. Res. 37:897908.
- Campbell, G.S. 1974. A simple method for determining unsaturated conductivity from moisture retention data. Soil Sci. 117:311314.
- Cosby, B.J., G.M. Hornberger, R.B. Clapp, and T.R. Ginn. 1984. A statistical exploration of the relationships of soil moisture characteristics to the physical properties of soils. Water Resour. Res. 20:682690.
- Eggleston, J.R., and J.J. Peirce. 1995. Dynamic programming analysis of pore space. Eur. J. Soil Sci. 46:581590.
- Elprince, A.M., and P.R. Day. 1977. Fitting solute breakthrough equations to data using two adjustable parameters. Soil Sci. Soc. Am. J. 41:3941.
- Feyen, J., D. Jacques, A. Timmerman, and J. Vanderborght. 1998. Modelling water flow and solute transport in heterogeneous soils: A review of recent approaches. J. Agric. Eng. Res. 70:231256.
- Fried, J.J., and M.A. Combarnous. 1971. Dispersion in porous media. Adv. Hydroscience 7:169282.
- Frye, W.W., and R.L. Blevins. 1997. Soil organic matter under long-term no-tillage and conventional tillage corn production in Kentucky. p. 227234. In E.A. Paul et al. (ed.) Soil organic matter in temperate agroecosystems: Long-term experiments in North America. CRC Press, Boca Raton, FL.
- Gee, G.W., M.D. Campbell, G.S. Campbell, and J.H. Campbell. 1992. Rapid measurement of low water potentials using a water activity meter. Soil Sci. 56:10681070.
- Gelhar, L.W., C. Welty, and K.R. Rehfeldt. 1992. A critical review of data on field-scale dispersion in aquifers. Water Resour. Res. 28:19551974.
- Germann, P.F. 1991. Length scales of convection-dispersion approaches to flow and transport in porous media. J. Contam. Hydrol. 7:3949.
- Gist, G.A., A.H. Thompson, A.J. Katz, and R.L. Higgins. 1990. Hydrodynamic dispersion and pore geometry in consolidated rock. Phys. Fluids A 2:15331544.
- Han, N., J. Bhakta, and R.G. Carbonell. 1985. Longitudinal and lateral dispersion in packed beds: Effect of column length and particle size distribution. AIChE. J. 31:277288.
- Hatano, R., N. Kawamura, J. Ikeda, and T. Sakuma. 1992. Evaluation of the effect of morphological features of flow paths on solute transport by using fractal dimensions of methylene blue staining pattern. Geoderma 53:3144.
- Jury, W.A., and G. Sposito. 1985. Field calibration and validation of solute transport models for the unsaturated zone. Soil Sci. Soc. Am. J. 49:13311341.[Abstract/Free Full Text]
- Kanchanasut, P., D.R. Scotter, and R.W. Tillman. 1978. Preferential solute movement through larger soil voids. II. Experiments with saturated soil. Aust. J. Soil Res. 16:269276.
- Klute, A. 1986. Water retention: Laboratory methods. p. 635662. In A. Klute (ed.) Methods of soil analysis. Part 1. 2nd ed. Agron. Monogr. 9. ASA and SSSA, Madison, WI.
- Koch, S., and H. Flühler. 1993. Solute transport in aggregated porous media: comparing model independent and dependent parameter estimation. Water Air Soil Pollut. 68:275289.
- Kool, J.B., J.C. Parker, and L.W. Zelazny. 1989. On the estimation of cation exchange parameters from column displacement experiments. Soil Sci. Soc. Am. J. 53:13471355.[Abstract/Free Full Text]
- Korom, S.F. 2000. An adsorption isotherm for bromide. Water Resour. Res. 36:19691974.
- Leij, F.J., W.J. Alves, M.Th. van Genuchten, and J.R. Williams. 1996. Unsaturated soil hydraulic database, UNSODA 1.0 User's manual. Rep. EPA/600/R-96/095, US EPA, Ada, OK.
- Leij, F.J., and J.H. Dane. 1991. Solute transport in a two-layer medium investigated with time moments. Soil Sci. Soc. Am. J. 55:15291535.[Abstract/Free Full Text]
- Li, Y., and M. Ghodrati. 1997. Preferential transport of solute through soil columns containing constructed macropores. Soil Sci. Soc. Am. J. 61:13081317.[Abstract/Free Full Text]
- Neuman, S.P. 1990. Universal scaling of hydraulic conductivities and dispersivities in geologic media. Water Resour. Res. 26:17491758.
- Parker, J.C., and M.Th. van Genuchten. 1984. Determining transport parameters from laboratory and field tracer experiments. Virginia Agric. Exp. Stn. Bull. 84-3, Virginia Polytechnic Inst. and State Univ., Blacksburg, VA.
- Passioura, J.B., and D.A. Rose. 1971. Hydrodynamic dispersion in aggregated media 2. Effects of velocity and aggregate size. Soil Sci. 111:345351.
- Perfect, E., and M.C. Sukop. 2001. Models relating solute dispersion to pore space geometry in saturated media: A review. p. 77146. In H.M. Selim and D.L. Sparks (ed.) Physical and chemical processes of water and solute transport/retention in soils. SSSA Special Publ. 56. SSSA, Madison WI.
- Rawls, W.J., T.J. Gish, and D.L. Brakensiek. 1991. Estimating soil water retention from soil physical properties and characteristics. Adv. Agron. 16:213234.
- Schweich, D., M. Sardin, and J.-P. Gaudet. 1983. Measurement of a cation exchange isotherm from elution curves obtained in a soil column: Preliminary results. Soil Sci. Soc. Am. J. 47:3237.[Abstract/Free Full Text]
- Seaman, J.C., P.M. Bertsch, and W.P. Miller. 1995. Ionic tracer movement through highly weathered sediments. J. Contam. Hydrol. 20:127143.
- Seyfried, M.S., and P.S.C. Rao. 1987. Solute transport in undisturbed columns of an aggregated tropical soil: Preferential flow effects. Soil Sci. Soc. Am. J. 51:14341444.[Abstract/Free Full Text]
- Skopp, J. 1984. Analysis of solute movement in structured soils. p. 220227. In J. Bouma and P.A.C. Raats (ed.) Proc. ISSS Symp. Water and solute movement in heavy clay soils. ILRI, Wageningen, the Netherlands.
- Vervoort, R.W., D.E. Radcliffe, and L.T. West. 1999. Soil structure development and preferential solute flow. Water Resour. Res. 35:913928.
- Vogel, H.J. 2000. A numerical experiment on pore size, pore connectivity, water retention, permeability, and solute transport using network models. Eur. J. Soil Sci. 51:99105.
- Walker, P.J.C., and S.T. Trudgill. 1983. Quantimet image analysis of soil pore geometry: comparison with tracer breakthrough curves. Earth Surf. Processes Landforms 8:465472.
- Xu, M., and Y. Eckstein. 1997. Statistical analysis of the relationships between dispersivity and other physical properties of porous media. Hydrogeology J. 5:420.
- Yu, C., A.W. Warrick, and M.H. Conklin. 1999. A moment method for analyzing breakthrough curves of step inputs. Water Resour. Res. 35:35673572.
- Zheng, C., and G.D. Bennett. 1995. Applied contaminant transport modeling. Van Nostrand Reinhold, New York.
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]
|
 |
|

|
 |

|
 |
 
J. L. Starr, A. M. Sadeghi, and Y. A. Pachepsky
Monitoring and Modeling Lateral Transport through a Large In Situ Chamber
Soil Sci. Soc. Am. J.,
October 27, 2005;
69(6):
1871 - 1880.
[Abstract]
[Full Text]
[PDF]
|
 |
|

|
 |

|
 |
 
E. Perfect, A. B. Kenst, M. Diaz-Zorita, and J. H. Grove
Fractal Analysis of Soil Water Desorption Data Collected on Disturbed Samples with Water Activity Meters
Soil Sci. Soc. Am. J.,
July 1, 2004;
68(4):
1177 - 1184.
[Abstract]
[Full Text]
[PDF]
|
 |
|