Comparison of Field Data and Water-Balance Predictions for a Capillary Barrier Cover

Predictions of surface runoff (R), evapotranspiration (ET), soil–water storage (S), and percolation obtained using three numerical codes (LEACHM, HYDRUS, and UNSAT-H) employed to simulate the hydrology of water-balance covers are compared to measured water-balance data from a lysimeter used to monitor a capillary barrier cover profile in a subhumid climate. All of the codes captured the seasonal variations in water-balance quantities observed in the field. LEACHM and HYDRUS predicted total R during the monitoring period with reasonable accuracy (within 18?mm using general mean parameters), but the timing of predicted and observed R events was different. In contrast, UNSAT-H consistently overpredicted R by at least 239?mm. Evapotranspiration was predicted reliably (within 60?mm) with all three codes when data from the first year were excluded. However, all three codes overpredicted ET in late winter and early spring, when snowmelt was occurring and S was accumulating in the field. Consequently, S generally was underpredicted by all three codes. Predicted and measured percolation were in good agreement (within 1?mm/year), except during the first year. Results of the comparison indicate that cover modelers should scrutinize runoff predictions for reasonableness and carefully account for snow accumulation, snow melt, and ET during snow cover.     

Capillary Barriers: Design Variables and Water Balance

Water balance simulations were conducted with the unsaturated flow model UNSAT-H to assess how layer thicknesses, unsaturated hydraulic properties, and climate affect the performance of capillary barriers. Simulations were conducted for four locations in semiarid and arid climates. Hydraulic properties of four finer-grained and two coarser-grained soils were selected to study how saturated and unsaturated hydraulic properties affect the water balance. Results of the simulations indicate that thickness and hydraulic properties of the surface layer significantly affect the water balance of capillary barriers. As expected, increasing the thickness or reducing the saturated hydraulic conductivity of the finer-grained surface layer reduces percolation. Unsaturated hydraulic properties of the coarser layer also affect the water balance, including the storage capacity of the surface layer as well as the onset and amount of percolation from the cover. Thickness of the coarser layer has a much smaller impact on the water balance. Climate also affects the water balance. Greater soil water storage capacity is required at sites where the season with more frequent and less intense precipitation does not coincide with the season having highest evapotranspiration.     

Field Data from a Capillary Barrier and Model Predictions with UNSAT-H

Water balance data are presented from a capillary barrier test section located on the final cover of a municipal solid waste landfill in a semiarid region (E. Wenatchee, Washington, U.S.). Water balance and meteorological data were collected from November 1992 to August 1995. Estimates of the water balance were made using the program UNSAT-H, with input consisting of meteorological data, soil properties, and vegetative information. Estimates of evapotranspiration and soil-water storage by UNSAT-H agree reasonably well with the field data. Peak soil-water storage was underestimated during the winter and evapotranspiration was overestimated in late winter. Water contents were estimated reasonably, although the changes in water content of the sand obtained from UNSAT-H were not as large as, and occurred less quickly than, that in the field. Percolation was generally overestimated, with the greatest overestimation occurring during Winter 1993, which had substantial snowfall. Surface runoff was underestimated; no runoff was obtained from UNSAT-H, whereas 7.4 cm of runoff was measured in the field. The overestimates in percolation appear to be closely related to underestimates in runoff and extra storage in the sand layer caused by the geocomposite drain used in the test section. Snowmelt, freezing of the soil surface, and hysteresis in soil hydraulic properties also appear to have had an effect on the differences between estimated and measured water balances.     

Evolution of the Hydraulic Conductivity of Reclamation Covers over Sodic/Saline Mining Overburden

The evolution of the field saturated hydraulic conductivity of four covers located on a reclaimed saline-sodic shale overburden from oil sands mining is presented. Three covers consisted of a surface layer of peat/glacial topsoil over a mineral, soil. and one cover was a single layer of mixed peat and mineral soil. Measurements of the field saturated hydraulic conductivity of the cover and shale materials were made with a Guelph permeameter between 2000 and 2004. The hydraulic conductivity of the cover materials in the multilayered covers increased by one to two orders of magnitude over the first few monitoring seasons. The hydraulic conductivity of the single-layer cover system, which was placed three years before the multilayered covers, marginally increased from 2000 to 2002 and then remained relatively unchanged. The hydraulic conductivity of the shale underlying all four covers increased approximately one order of magnitude. Soil temperature measurements indicated that one freeze/thaw cycle occurred each year within all cover soils and the surficial overburden. This suggests that freeze/thaw effects were the cause of the observed increases in hydraulic conductivity, as previously observed by other researchers working on compacted clays.     

Comparison of Soil Hydraulic Property Measurement Methods

Unsaturated and saturated soil hydraulic properties were determined and compared for three sandy soils at adjacent field sites. Drying soil–water retention curves were measured on soil specimens using a pressure plate apparatus. Saturated hydraulic conductivities (Ks) were measured with a Guelph permeameter and falling head tests. Parameter optimization was used to simultaneously estimate the drying and wetting soil–water retention and hydraulic conductivity curves from cone permeameter and multistep inflow/outflow data. Ks values from all test methods were within an order of magnitude of each other at each site and, as expected, trended with bulk density. The Guelph permeameter generally yielded the highest Ks values. The soil–water retention curves were similar in shape, except for the cone permeameter curves, which had steeper slopes due to rapid flow of water into the soil. Relative hydraulic conductivity curves were similar in character to the soil–water retention curves. Each method provided important information about the soil hydraulic properties. No one method provided the entire range of information provided by all of the tests combined, and no one method was found to be superior to the others.     

Shear Strength and Shear-Induced Volume Change Behavior of Unsaturated Soils from a Triaxial Test Program

A series of unsaturated soil triaxial tests were performed on four soils including sand, silt, and a low plasticity clay. Attempts were made to correlate unsaturated soil properties from these tests and data from the literature with soil-water characteristics curve (SWCC), soil gradation, and saturated soil properties. The feasibility of estimating unsaturated soil property functions from saturated soil properties, SWCCs and gradation data, is demonstrated. A hyperbolic model for estimation of the unsaturated soil parameter, ?b, versus matric suction is presented. Shear induced volume change behavior was also studied, and results are included in this paper. Although not correlated with soil index properties, these shear-induced volume change data are important to complete stress-deformation analyses, and represent a significant addition to the existing data base of unsaturated soil properties.     

Postconstruction Changes in the Hydraulic Properties of Water Balance Cover Soils

Hydraulic properties of soils used for water balance covers measured at the time of construction and one to four years after construction are compared to assess how the hydraulic properties of cover soils change over time as a result of exposure to field conditions. Data are evaluated from ten field sites in the United States that represent a broad range of environmental conditions. The comparison shows that the saturated hydraulic conductivity (Ks) can increase by a factor of 10,000, saturated volumetric water content (θs) by a factor of 2.0, van Genuchten’s α parameter by a factor of 100, and van Genuchten’s n parameter can decrease by a factor of 1.4. Larger changes occur for denser or more plastic fine-textured soils that have lower as-built Ks, α, and θs and higher as-built n, resulting in a reduction in the variation in hydraulic properties that can be attributed to compaction. After two to four years, many water balance cover soils can be assumed to have Ks between 10?5 and 10?3?cm/s, θs between 0.36 and 0.40, α between 0.002 and 0.2?kPa?1, and n between 1.2 and 1.5. The data may be used to estimate changes in hydraulic properties for applications such as waste containment, where long-term maintenance of hydraulic properties in shallow engineered soil layers is important.     

Hydraulic Performance of Geosynthetic Clay Liners in a Landfill Final Cover

Percolation from a landfill final cover containing a geosynthetic clay liner (GCL) as the hydraulic barrier is described. The GCL was covered with 760?mm of vegetated silty sand and underlain with two gravel-filled lysimeters to monitor percolation from the base of the cover. Higher than anticipated percolation rates were recorded in both lysimeters within 4–15?months after installation of the GCL. The GCL was subsequently replaced with a GCL laminated with a polyethylene geofilm on one surface (a “composite” GCL). The composite GCL was installed in two ways, with the geofilm oriented upwards or downwards. Low percolation rates (2.6–4.1?mm/year) have been transmitted from the composite GCL for more than 5?years regardless of the orientation of the geofilm. Samples of the conventional GCL that were exhumed from the cover ultimately had hydraulic conductivities on the order of 5×10?5?cm/s. These high hydraulic conductivities apparently were caused by exchange of Ca and Mg for Na on the bentonite combined with dehydration. The overlying and underlying soils likely were the source of the Ca and Mg involved in the exchange. Column experiments and numerical modeling indicated that plant roots and hydraulic anomalies caused by the lysimeters were not responsible for the high hydraulic conductivity of the GCL. Despite reports by others, the findings of this study indicate that a surface layer 760?mm thick is unlikely to protect conventional GCLs from damage caused by cation exchange and dehydration. Accordingly, GCLs should be used in final covers with caution unless if cation exchange and dehydration can be prevented or another barrier layer is present (geomembrane or geofilm).     

Hydraulic Conductivity of Soils from Grain-Size Distribution: New Models

This paper presents new developments of regression-based models to predict the saturated hydraulic conductivity of compacted soils from grain-size distribution. The models incorporate parameter values that adequately represent the distribution of grain sizes. Alternative representations of the grain-size distribution, the fractal dimension and entropy of the distributions, as well as porosity, soil density, and fines content are used in the models to estimate the hydraulic conductivity. These parameters that characterize the textural and hydraulic properties of the soil are combined and used in a multidimensional analysis to estimate the hydraulic conductivity. The predictions of the developed models are compared with those of existing models and laboratory measurements of hydraulic conductivity. The results suggest that the newly developed models outperform the existing models in predicting hydraulic conductivity using information from grain-size distribution. The presented models are suggested as alternatives to, for example, laboratory measurements of the hydraulic conductivity of certain soils that may be difficult to prepare or that may take several days or perhaps weeks to perform. In certain circumstances it may also be used to give first-hand information about the hydraulic properties in a field environment.     

Rate of Capillary Rise in Soil

A rigorous closed-form analytical solution is developed for analyzing the rate of capillary rise in soils. The new solution can be reduced to Terzaghi’s classical solution if the nonlinearity in the hydraulic conductivity with changing soil suction is ignored. Results obtained using the new solution are compared with Terzaghi’s classical solution and a series of previously documented experimental data from open-tube capillary rise tests. The new solution is a significant improvement over the previous solution, thus providing more realistic and practical predictions for the rate of capillary rise in unsaturated soils.     

Process Modeling of Storm-Water Flow in a Bioretention Cell

A two-dimensional variable saturated flow model was developed to simulate subsurface flow in bioretention facilities employing the Richards’ equation. Variable hydrologic performances of bioretention are evaluated using the underdrain outflow hydrographs, outflow volumes for 10 storms with various duration and depth, and flow duration curves for 25 different storms. The effects of some important design parameters and elements are tested, including media type, surrounding soils, initial water content, ratio of drainage area to bioretention surface area, and ratio of cell length to width. Model results indicate that the outflow volume via underdrain is less than the inflow; the flow peak is significantly reduced and delayed. Underdrain outflow volume from loamy sand media (with larger Ks) is larger than that from sandy clay loam media. The saturated hydraulic conductivity, storage capacity, and exfiltration into surrounding soils contribute to the hydrologic performance of a bioretention cell. Initial media storage capacity is affected by the hydraulic properties of media soils, initial water content, and bioretention surface area. The exfiltration volume is determined by the surrounding soil type and exfiltration area, dominated by flow through the bottom of the media.     

Centrifuge Permeameter for Unsaturated Soils. I: Theoretical Basis and Experimental Developments

A new centrifuge permeameter was developed with the specific objective of expediting the measurement of the hydraulic characteristics of unsaturated soils. The development, theoretical basis, and typical results associated with using the centrifuge permeameter for concurrent determination of the soil-water retention curve (SWRC) and hydraulic conductivity function (K function) of unsaturated soils are presented in this paper. Components developed for the centrifuge permeameter are described, including the centrifuge, permeameter, water flow control system, and instrumentation used to concurrently and nondestructively measure the infiltration rate (flow pump and outflow transducer), volumetric water content (time domain reflectometry), and matric suction (tensiometers) in flight during steady-state infiltration. A companion paper focuses on definition of the SWRC and K function for a clay soil using the procedures described in this paper. While conventional geotechnical centrifuges are used to reproduce the response of earth structure prototypes, the centrifuge developed in this study is used to accelerate flow processes. Accordingly, it required a comparatively small radius (0.7 m) but high angular velocity (up to 875 rpm or 600 g’s) to impart a wide range of hydraulic gradients to an unsaturated soil specimen. Analytical solutions to Richards’ equation in the centrifuge indicate that steady-state infiltration allows direct determination of the relationships between suction, volumetric water content, and hydraulic conductivity from the instrumentation results. Typical instrumentation results during a drying stage are presented to illustrate determination of data points on the SWRC and K function at steady state. These results were found to be consistent with analytical flow solutions.     

Case Study of a Full-Scale Evapotranspiration Cover

The design, construction, and performance analyses of a 6.1?ha evapotranspiration (ET) landfill cover at the semiarid U.S. Army Fort Carson site, near Colorado Springs, Colo. are presented. Initial water-balance model simulations, using literature reported soil hydraulic data, aided selection of borrow-source soil type(s) that resulted in predictions of negligible annual drainage ( ? 1?mm/year). Final construction design was based on refined water-balance simulations using laboratory determined soil hydraulic values from borrow area natural soil horizons that were described with USDA soil classification methods. Cover design components included a 122?cm thick clay loam (USDA), compaction ? 80% of the standard Proctor maximum dry density (dry bulk density ～ 1.3?Mg/m3), erosion control measures, top soil amended with biosolids, and seeding with native grasses. Favorable hydrologic performance for a 5?year period was documented by lysimeter-measured and Richards’-based calculations of annual drainage that were all <0.4?mm/year. Water potential data suggest that ET removed water that infiltrated the cover and contributed to a persistent driving force for upward flow and removal of water from below the base of the cover.     

Combined Brooks-Corey/Burdine and van Genuchten/Mualem Closed-Form Model for Improving Prediction of Unsaturated Conductivity

The performance of the widely used conventional closed-form models of Brooks and Corey and van Genuchten is restricted to soils characterized by a specific shape form of the water retention curve. Hence, the van Genuchten model more accurately describes “S”-shaped retention curves characterizing finer-textured soils, whereas the Brooks-Corey model is much better adapted for “J”-shaped retention curves characterizing relatively coarse-textured soils. In this work, a new closed-form soil hydraulic model is proposed. The suggested continuous-form function accurately describes soil retention curves irrespective of their specific shape form. New algebraic expressions based on Mualem’s statistical model and another new model that is a combination of the Mualem and Burdine theories were derived for the prediction of the unsaturated hydraulic conductivity function. Comparisons of ten soils known from international bibliography were performed. It is concluded that the proposed water retention curve expression as well as the hydraulic conductivity predictions showed significant improvement over the conventional van Genuchten and Brooks-Corey closed-form models, particularly for conductivity values near the residual water content and saturation.     

Centrifuge Permeameter for Unsaturated Soils. II: Measurement of the Hydraulic Characteristics of an Unsaturated Clay

This paper presents the hydraulic characteristics of an unsaturated, compacted clay, including its soil-water retention curve (SWRC) and hydraulic conductivity function (K function), determined using a new centrifuge permeameter developed at the University of Texas at Austin. A companion paper describes the apparatus, its instrumentation layout, and data reduction procedures. Three approaches are evaluated in this study to define the SWRC and K function of the compacted clay under both drying and wetting paths, by varying the inflow rate, the g level, or both. For imposed inflow rates ranging from 20 to 0.1 mL/h and g levels ranging from 10 to 100 g, the measured matric suction ranged from 5 to 70 kPa, the average volumetric water content ranged from 23 to 33%, and the hydraulic conductivity ranged from 2×10?7 to 8×10?11?m/s. The SWRCs and K functions obtained using the three different testing approaches were very consistent, and yielded suitable information for direct determination of the hydraulic characteristics. The approaches differed in the time required to complete a testing stage and in the range of measured hydraulic conductivity values. The g level had a negligible effect on the measured hydraulic characteristics of the compacted clay. The SWRCs and K functions defined using the centrifuge permeameter are consistent with those obtained using pressure chamber and column infiltration tests. The K functions defined using the centrifuge permeameter follow the same shape as those obtained from predictive relationships, although the measured and predicted K functions differ by two orders of magnitude at the lower end of the volumetric water content range.     

Transient Rainfall-Runoff Loadings to a Partial Exfiltration System: Implications for Urban Water Quantity and Quality

Modification of rainfall-runoff processes by urban infrastructure and anthropogenic activities impacts receiving waters and the surrounding terrestrial environment. Infiltration–exfiltration systems such as a partial exfiltration reactor (PER) when loaded by transient sheet flow have the potential to attenuate the impact of both the quantity and quality of urban runoff. These in situ systems are subject to highly variable water quality and quantity while functioning under variably saturated flow conditions. To improve the understanding of field-scale PER performance as a rainfall-runoff unit operation and process, a two-dimensional (2D) numerical model was used to simulate the effluent hydrograph and water content profiles under transient hydraulic loadings. Richard’s equation was applied in the 2D model using parameters estimated from laboratory experiments and hydrographs measured for an in situ PER. The temporal dynamics of the water content illustrated the ability of the PER to lower peak flow, redistribute volume, and attenuate temporal aspects of the inflow hydrograph. Results demonstrated the role of the PER to attenuate runoff water quantity, while also providing water quality improvements, as illustrated for suspended solids and dissolved Cu. Simulation of historical events for different surrounding soils illustrated the critical role of surrounding soil conditions on PER performance. While the PER demonstrated water quantity attenuation benefits for design storms (1, 2, 5?year return periods), results also illustrate how a given PER design for clayey soils conditions can be limiting for intense events. Evaporation was a dominant mechanism for the drying process in the PER upper layer; with a residual moisture content in the porous pavement layer achieved in less than 2?days in summer for Cincinnati, Ohio.     

Modulus-Suction-Moisture Relationship for Compacted Soils in Postcompaction State

Despite clear evidence, changes in mechanical properties (i.e., stiffness or modulus) of compacted subgrades in response to subgrade moisture regime changes after construction have rarely been investigated in the geotechnical profession. In particular, when in-service assessment of pavement subgrade is made, the modulus-moisture variation should be addressed on the basis of unsaturated soil mechanics. This study presents the unsaturated small-strain modulus behavior of five predominately fine-grained compacted subgrade soils. The small-strain shear modulus (Go) of saturated compacted specimens subjected to a desorption soil-water characteristic curve (SWCC) was evaluated using bender elements. A test apparatus was designed to apply two stress state variables, the net confining pressure and matric suction, during the Go measurements. The relationship between Go and the SWCC under a constant mean net stress was developed. Additionally, the effect of compaction moisture content, compaction energy, and soil type on the Go-SWCC relationship was investigated. Finally, a relationship describing the small-strain modulus behavior of unsaturated compacted soils is proposed.     

Hydraulic Property Estimation Using Piezocone Results

Governing underground water flow, hydraulic properties such as hydraulic conductivity or coefficient of consolidation are major geotechnical parameters. Determination of hydraulic properties, however, is traditionally time consuming and expensive. This research proposes an easy and economical way of determining the hydraulic properties of soils through piezocone penetration tests. Pore pressure responses of soils from piezocone penetration tests are numerically analyzed herein by the coupled theory of mixtures, which is based on the large strain elastoplasticity. Using the numerical results, the effects of input parameters are evaluated. Simple equations are also derived for a faster estimation of the hydraulic conductivity or the coefficient of consolidation of soils. The hydraulic properties predicted by these derived equations agree reasonably with the measured results.     

Relative Abundance of Monovalent and Divalent Cations and the Impact of Desiccation on Geosynthetic Clay Liners

Laboratory experiments were conducted on a geosynthetic clay liner (GCL) containing Na–bentonite to determine how the swell index and hydraulic conductivity of GCLs are affected by wet-dry cycling with solutions having different relative abundance of monovalent and multivalent cations. Relative abundance of monovalent and multivalent cations was characterized by the RMD of the test solution, which is defined as the ratio of the total molarity of monovalent cations to the square root of the total molarity of multivalent cations at a given ionic strength. RMD was found to control the final swell index, relative abundance of monovalent and divalent cations in the final exchange complex, and the final hydraulic conductivity of bentonite exposed to wet-dry cycling. Ionic strength affects the number of wet-dry cycles required for a change in hydraulic conductivity to occur and the rate of change in swell index. Large increases in hydraulic conductivity and loss of swelling capacity occurred for solutions having RMD ? 0.07?M1/2. Modest or small changes in hydraulic conductivity and swell index were obtained when the RMD was ≥ 0.14?M1/2. These findings suggest that chemical analysis of the pore water in cover soils may prove useful in evaluating the compatibility of GCLs and cover soils used in applications where wet-dry cycling may occur.