Hydraulic Conductivity Measurement from On-the-Fly uCPT Sounding and from VisCPT

Detailed profiles of hydraulic conductivity are recovered from the deployment of direct-push permeameters at the Geohydrologic Experimental and Monitoring Site, Kansas. Measurements with thin tapered tips, and with standard cone penetration test (uCPT) tips, show only minor differences, suggesting that tip-local disturbance effects are small, and that routine uCPT measurements are therefore representative of pristine conditions. Permeameter measurements are correlated against closely deployed uCPT measurements, estimates of hydraulic conductivity from uCPT sounding correlations, and from grain size correlations derived from both vision CPT (VisCPT) and from cone metrics. On-the-fly evaluations of hydraulic conductivity require that the tip-local pressure field is both steady and partially drained. Continuous penetration is shown to yield pore pressures sufficiently close to steady to enable conductivities to be directly determined. Cone metrics of cone resistance, sleeve friction, and pore pressure ratio are shown to be sufficient to discriminate between partially drained and undrained behavior, and therefore to define the permissible regime where conductivities may be determined from uCPT sounding data. Estimates of hydraulic conductivities from uCPT sounding data are shown to correlate with independently measured magnitudes of hydraulic conductivity recovered using the permeameter tests. However, most of hydraulic conductivities from the permeameter tests (4.5?cm length screen) are underpredicted, suggesting that storage effects, the inability to reach a steady state, or the effects of dilation may influence the response. Profiles of hydraulic conductivities evaluated from the on-the-fly method also correlate well with the permeameter measurements. Predictions from soil classification and from VisCPT methods are also capable of estimating conductivities, with soil classifications giving the closest correlations of these two for this particular suite of data.     

Comparison of Models for Computing Drainage Discharge

The WAVE model describes the transport and transformations of matter and energy in the soil, crop, and vadose environment. A lateral field drainage subprogram was added to the WAVE model to simulate lateral subsurface drainage flow. The subsurface drainage is considered as the drainage provided by evenly spaced parallel drains with a free outlet: drain tubing or ditch. The rate of subsurface water movement into drain tubes or ditches depends on the hydraulic conductivity of the soil, drain or ditch spacing, hydraulic head in the drains, profile depth, and water table elevation. Hooghoudt's steady-state equation was selected for incorporation in the WAVE model. The subsurface drainage subprogram was calibrated and validated by comparison with the SWAP model (The Netherlands) and DRAINMOD (the United States) and partially by using 7 years of drain outflow data from an experimental field under fallow and cropped conditions. The comparative study revealed that the three models performed equally well and that the models were reliable and accurate tools for predicting the drainage flux as a function of rainfall-evapotranspiration and local conditions. The WAVE model, in comparison to the SWAP and DRAINMOD model, provided as good a prediction of the lateral subsurface drainage flow to drains. The statistical analysis between each model and observed data revealed that the three models were able to predict with sufficient accuracy the observed drainage discharge. The DRAINMOD model, however, has the advantage of giving a more accurate estimate of the discharge, resulting in a more precise modeling. The models were consistent in predicting water table levels, but they could not be verified against field data because of a lack of suitable measurements.     

Laboratory Evaluation of Sand Underdrains

Constant-head hydraulic conductivity tests are performed on layered heterogeneous porous media to evaluate the use of underdrains to calculate the hydraulic conductivity of an overlying, less permeable medium. The layered profiles consist of a barrier layer comprising sand mixed with 10% kaolin, overlying a foundation layer comprising sand mixed with only 5% kaolin. Underdrains are evaluated by replacing excavated portions of the foundation layer with only sand. The results indicate that preferential flow of water occurs around, rather than through, the sand underdrains resulting in an underestimate of the measured hydraulic conductivity of the barrier layer assuming 1D, saturated flow in accordance with standard practice. The observed preferential flow effect is consistent with previously published numerical simulations of unsaturated flow through similarly layered heterogeneous soil profiles that indicate lateral flow around underdrains due to the contrast in unsaturated properties of the soils. The results of this study have important ramifications with respect to the use of underdrains to measure in situ hydraulic conductivity of compacted clay liners for waste containment.     

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.     

Contaminant Transport Model for Unsaturated Soil Using Fuzzy Approach

Contaminant transport in the unsaturated zone is important for managing water resources and assessing the damage due to contamination in the field of irrigation, water management, wastewater management, and urban and agricultural drainage systems. Deterministic modeling which is widely used for contaminant transport is not adequate because it considers model input parameters as well-defined crisp values and hence does not account for uncertainties and imprecision. This paper presents a contaminant transport model based on fuzzy set theory to simulate water flow and contaminant transport in the unsaturated soil zone under surface ponding condition. Among all soil hydraulic parameters that have uncertainty associated with them, saturated hydraulic conductivity was found to be the most sensitive to model outputs. Trapezoidal fuzzy numbers were used to express the uncertainties associated with saturated hydraulic conductivity. The incorporation of uncertainties into contaminant transport model is useful in decision making, as it yields scientifically and practically based estimates of contaminant concentration.     

Sorbent Wicking Device for Sampling Hydrophobic Organic Compounds in Unsaturated Soil Pore Water. I: Design and Hydraulic Characteristics

The hydraulic characteristics of horizontally installed sorbent wick sampling devices were evaluated through wick tracer studies and laboratory soil column experiments to assess the influence of horizontal wick length and sampler interface design on sampling pore water in unsaturated soils. The nominal sampler design consisted of a cylindrical porous metal interface packed with granular-activated carbon encapsulating the end of a fiberglass wick that extended 100 cm horizontally from the interface before dropping 100 cm vertically to a collection vessel. The maximum sampling rate of horizontally installed wick systems declines exponentially with increasing horizontal wick length, while the vertical length influences the range of soil–water pressures that may be sampled. The nominal design sampled pore water from clay loam laboratory columns at 8 to 14 mL?h?1 under steady-state infiltration conditions and 2 to 5 mL?h?1 under draining conditions across a ?10 to ?45 cm H2O soil–water pressure range. Sampling rates in medium-grained sand under similar flow conditions were less than that of the clay loam due to reduced water content and reduced interface/soil contact area. An analysis of observed sampling velocities versus calculated soil water contents and hydraulic conductivities indicated that the design performs best when the soil water content is greater than 0.15 and unsaturated hydraulic conductivity is greater than 0.2 cm?h?1. A hydraulic model was developed that estimates the sampling velocity of the nominal design based on sampler interface pressure, which was linearly correlated with soil pressure.     

Foundry Green Sands as Hydraulic Barriers: Field Study

A field study was conducted to determine if the field hydraulic conductivity of barrier layers constructed with foundry green sand is comparable to the hydraulic conductivity measured in the laboratory on laboratory-compacted specimens normally used for testing during design. Three test pads were constructed with foundry green sand. Their field hydraulic conductivity was measured using sealed double ring infiltrometers, two-stage borehole permeameters, and on large block specimens. Additional field hydraulic conductivity tests were conducted on the test pads after exposure to winter weather causing freeze-thaw cycling and summer weather causing desiccation. The field hydraulic conductivity data followed the same trends with bentonite content and liquid limit observed in the laboratory. When the bentonite content is greater than 6% (by weight), the plasticity index is greater than 3, or the liquid limit is greater than 20, the hydraulic conductivity is less than 10?7?cm/s. Testing after winter exposure showed that the field hydraulic conductivity was unaffected by winter weather, even though the test pads underwent up to six freeze-thaw cycles (depending on depth). Similarly, exposing the test pads to summer weather had no measurable effect on the field hydraulic conductivity. The field study validated that foundry sand is a useful industrial by-product that can be beneficially used as a hydraulic barrier material.     

Foundry Green Sands as Hydraulic Barriers: Laboratory Study

A laboratory testing program was conducted to assess the use of foundry sands from gray iron foundries, typically called green sands, as hydraulic barrier materials. Foundry green sands are mixtures of fine uniform sand, bentonite, and other additives. Specimens of foundry sand were compacted in the laboratory at a variety of water contents and compactive efforts and then permeated in rigid-wall and flexible-wall permeameters to define relationships between hydraulic conductivity, compaction water content, and dry unit weight. Additional tests were conducted to assess how hydraulic conductivity of compacted foundry sand is affected by environmental stresses such as desiccation, freeze-thaw, and chemical permeation. Results of the tests show that the hydraulic conductivity of foundry sand is sensitive to the same variables that affect hydraulic conductivity of compacted clays (i.e., compaction water content, and compactive effort). However, hydraulic conductivities <10?7 cm∕s can be obtained for many foundry sands using a broad range of water contents and compactive efforts, including water contents dry of optimum and at lower compactive effort. The hydraulic conductivity of foundry sand was generally unaffected by freeze-thaw, desiccation, or permeation with 0.1 N salt solution or municipal solid waste leachate, but was incompatible with acetic acid (pH = 3.5). Hydraulic conductivity of foundry sands correlates well with bentonite content and liquid limit, with hydraulic conductivity less than 10?7 cm∕s being achieved for bentonite content ≥6% and∕or liquid limit >20.     

Effects of Soil Water Salinity on Field Soil Hydraulic Functions

Unsaturated soil hydraulic parameters and functions used in numerical models to simulate water flow and solute transport in the unsaturated zone are generally considered invariant of soil water salinity levels. This study uses 5 years of field soil water salinity levels at three observation sites from the Land Retirement Demonstration Project (LRDP) (20069) located in western Fresno County, California, to test the hypothesis that field unsaturated soil hydraulic properties are also a function of soil water salinity level. The HYDRUS-1D software package for simulating one-dimensional (1D) movement of water, heat, and multiple solutes in variably saturated media, and Parameter Estimation (PEST), a model-independent parameter optimizer, is used to optimize the soil hydraulic parameters and downward bottom flux corresponding to three different average soil salinity levels at each site. The results show that at the same pressure head, soil water content is less with higher soil water salinity as compared with lower soil water salinity. It is thus concluded that the use of soil water salinity invariant soil water hydraulic parameters in numerical modeling can seriously compromise predictions, especially for a variable soil water salinity environment.     

Field Performance of a Compacted Clay Landfill Final Cover at a Humid Site

A study was conducted in southern Georgia, USA, to evaluate how the hydraulic properties of the compacted clay barrier layer in a final landfill cover changed over a 4-year service life. The cover was part of a test section constructed in a large drainage lysimeter that allowed continuous monitoring of the water balance. Patterns in the drainage (i.e., flow from the bottom of the cover) record suggest that preferential flow paths developed in the clay barrier soon after construction, apparently in response to desiccation cracking. After four years, the clay barrier was excavated and examined for changes in soil structure and hydraulic conductivity. Tests were conducted in situ with a sealed double-ring infiltrometer and two-stage borehole permeameters and in the laboratory on hand-carved blocks taken during construction and after four years of service. The in situ and laboratory tests indicated that the hydraulic conductivity increased approximately three orders of magnitude (from ≈ 10?7?to? ≈ 10?4?cm?s?1) during the service life. A dye tracer test and soil structure analysis showed that extensive cracking and root development occurred throughout the entire depth of the barrier layer. Laboratory tests on undisturbed specimens of the clay barrier indicated that the hydraulic conductivity of damaged clay barriers can be underestimated significantly if small specimens (e.g., tube samples) are used for hydraulic conductivity assessment. The findings also indicate that clay barriers must be protected from desiccation and root intrusion if they are expected to function as intended, even at sites in warm, humid locations.     

Mathematical Formulation and Validation of a Mixed Finite Element–Finite Difference Model for Simulating Phreatic Surfaces

The phreatic surface in an unconfined aquifer exists as a movable interface between the saturated and unsaturated zones. The movement of the phreatic surface depends on recharge, hydraulic conductivity, porosity, and horizontal and vertical flows. The location of the phreatic surface helps define the variably saturated flow domain in the subsurface. The variably saturated flow process in the subsurface is described by a parabolic partial differential equation. In this equation, the hydraulic conductivity and soil moisture capacity are used as the subsurface characteristics. The location of the phreatic surface is governed by a first-order partial differential equation. The governing parabolic partial differential equation is solved using a variational finite element formulation. The first order phreatic surface equation is then solved by loosely coupling with the governing parabolic partial differential equation describing the variably saturated flow. In the present study, a two-dimensional space is used to investigate the movement of the phreatic surface in a variably saturated unconfined flow domain. Based on the time-varying solutions of hydraulic heads, the location of the phreatic surface is simulated in a finite two-dimensional space.     

Infiltration Considerations for Ground-Water Recharge with Waste Effluent

Water reuse and ground-water recharge can be used to meet the growing demands for water, particularly in arid regions. Ground-water recharge using fresh water or treated wastewater is most often accomplished by infiltration from surface basins. The water percolates through the unsaturated soil region to an underlying aquifer for storage and future use. In the case of wastewater, additional treatment occurs as the effluent flows through the soil. The system hydraulics of recharge basins have been examined through a combination of field and laboratory investigations. These studies indicate that infiltration rates and soil aquifer treatment of wastewater are influenced by soil type and soil profile characteristics, surface clogging material, pond depth, and wetting∕drying cycle times. The surface-clogging layer was found to be susceptible to consolidation and to associated reduction in hydraulic conductivity under seepage forces.     

Study of Variation of Saturated Hydraulic Conductivity with Time

The variation of the saturated hydraulic conductivity with time, as a function of temperature, has been studied involving field measurement at the College of Aboureyhan Research farm. The College of Aboureyhan is a part of the University of Tehran and the above-mentioned farm is located in the lowland of the southeast of Tehran, Iran. For the purposes of this research study it was planned to measure and record the field data in a plot of 18?m2. The hydraulic conductivity data were measured in 18 test points inside the study area using the inverse auger-hole method. Experiments were carried out from August 16, 2005 to June 14, 2006. The recorded filed data were then used to calculate the saturated hydraulic conductivity data using Excel software. The relationship between soil temperature, water temperature, and also water’s viscosity with hydraulic conductivity, respectively, were determined. The results of the statistical analysis involving SAS software demonstrated that the variation of temperature can considerably affect the saturated hydraulic conductivity values. The results showed that the lowest values of Ks were obtained in the winter when the water, soil, and air temperature were minimal and these values increased when the temperature increased. The effect of the means of daily K values (values measured in each sample date Kd) and their corresponding soil temperature adjusted values on drain spacing was determined compared with the average of total K values.?The maximum over- and underestimation of drain spacing was 19.1 and 23.3%, respectively, for measured values. These estimations were 9 and 16% for adjusted ones. Also, using the average values of K measured in the soil temperature range of 17–23°C resulted in a lower over- or underestimation of drain spacing.     

Design of Two-Layered Porous Landscaping Detention Basin

Under the mandate of the Federal Clean Water Act, porous landscaping detention (PLD) has been widely used to increase on-site infiltration. A PLD system consists of a surface storage basin and subsurface filtering layers. The major design parameters for a PLD system are the infiltration rate on the land surface and the seepage rate through the subsurface medium. A low infiltration rate leads to a sizable storage basin while a high infiltration rate results in standing water if the subsurface seepage does not sustain the surface loading. In this study, the design procedure of a PLD basin is revised to take both detention flow hydrology and seepage flow hydraulics into consideration. The design procedure begins with the basin sizing according to the on-site water quality control volume. The ratio of design infiltration rate to sand-mix hydraulic conductivity is the key factor to select the thickness of sand-mix layer underneath a porous bed. The total filtering thickness for both sand-mix and gravel layers is found to be related to the drain time and infiltration rate. The recommended sand-mix and granite gravel layers underneath a PLD basin are reproduced in the laboratory for infiltration tests. The empirical decay curve for sand-mix infiltration rate was derived from the laboratory data and then used to maximize the hydraulic efficiency through the subsurface filtering layers. In this study, it is recommended that a PLD system be designed with the optimal performance to consume the hydraulic head available and then evaluated using the prolonged drain time for potential standing water problems under various clogging conditions.     

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.     

Simulation of Ground-Water Flow in Steep Basin with Shallow Surface Soil

A coupled ground-water∕channel flow distributed model has been developed for continuous simulation in a 123-km2 basin. The aim was to analyze the streamflow generation processes in natural vegetated environments. Finite-difference schemes have been used to solve conservation equations of the 2D saturated subsurface flow and the 1D kinematic surface flow. Because of the high hydraulic conductivity of the surface soil, only the saturation excess mechanism of runoff production has been considered. Parameter sensitivity analysis showed the overriding influence of soil storage capacity and conductivity. A grid discretization >100 m produces a hydraulic conductivity greater than physically meaningful, which considerably increases as the space-grid step increases. Results indicate that the model can satisfactorily simulate the water-flow behavior of the catchment after fitting the three parameters of surface hydraulic conductivity, effective porosity, and evapotranspiration losses. These are done after calculating the conductivity as a function of the height of the water table. The simulation efficiency has varied from 87% in the first 5-year calibration period to 85.8% in the subsequent 5-year validation period.     

Ohmic Conductivity of a Compacted Silty Clay

In its natural state, loess can be considered as an unstable soil, which develops large deformations when moistened. In Argentina, loess is used in most Geotechnical constructions, including embankments and liners. The interest of this work to evaluate the potential application of electrical conductivity measurements for monitoring the effects introduced by remolding and compaction in the soil. Samples of loess were compacted at varied densities and mixed with electrolytes of different concentrations. Electrical conductivity was measured with a two electrode cell. The effects introduced on the measured conductivity by frequency, degree of saturation, soil density, temperature, and electrolyte type and concentration are addressed. Additionally, hydraulic permeability tests were performed on compacted specimens of loess and the relationship between electrical and hydraulic conductivity was determined. It is concluded here that the ohmic conductivity of compacted specimens depends mainly on the salt concentration in the pore fluid, and volumetric water content. The effect of compaction density was observed to be less significant. The whole behavior of electric conductivity of loess is well described by the Archie’s law.     

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.     

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 Modeling of Light Nonaqueous Phase Liquids Transport in Unsaturated Soils

Centrifuge modeling appears useful for studying geo-environmental problems such as pollutant migration in subsurface systems. In this study, centrifuge tests were conducted to simulate a gasoline spill from a leaking underground storage tank (UST) and the subsequent subsurface migration of the gasoline. When the centrifugal acceleration reached the desired g level, the gasoline was released from the UST and then it migrated in the unsaturated soil for a prototype time equivalent to 1 year. After the centrifuge tests, soil samples were collected using sampling tubes and the concentrations of individual constituents in the light nonaqueous phase liquids (LNAPLs) were directly measured by means of gas chromatograph analysis. Two types of unsaturated soils were used to study the migration patterns of LNAPLs in unsaturated porous media. Centrifuge test data show that the migration pattern of LNAPLs is related to the soil type and the physical properties of individual constituents in the LNAPLs.