Validity limits for the van Genuchten–Mualem model and implications for parameter estimation and numerical simulation

The calculation of the relative hydraulic conductivity function based on water retention data is an attractive and widely used approach, since direct measurements of unsaturated conductivities are difficult. We show theoretically under which conditions an air-entry value for water retention data is definitely required when using the statistical approach of Mualem. Moreover we rigorously specify the conditions for which the classical van Genuchten–Mualem model leads to wrong predictions of relative hydraulic conductivity and, hence, an alternative formulation including an air-entry value should be used. Significant consequences are demonstrated for the inverse parameter estimation based on multistep outflow experiments. Furthermore it is shown that the use of a physically correct formulation of the water retention curve including an air-entry value and the derived hydraulic conductivity function influences not only the stability of numerical simulations but also their final results. This is especially grave as simulations with van Genuchten–Mualem parameters are frequently used to compare experiments and simulations and to draw conclusions on the correctness of Richards’ equation.     

Inverse upscaling of hydraulic parameters during constant flux infiltration using borehole radar

Modeling unsaturated flow in porous media requires constitutive relations that describe the soil water retention and soil hydraulic conductivity as a function of either potential or water content. Often, the hydraulic parameters that describe these relations are directly measured on small soil cores, and many cores are needed to upscale to the entire heterogeneous flow field. An alternative to the forward upscaling method using small samples are inverse upscaling methods that incorporate soft data from geophysical measurements observed directly on the larger flow field. In this paper, we demonstrate that the hydraulic parameters can be obtained from cross borehole ground penetrating radar by measuring the first arrival travel time of electromagnetic waves (represented by raypaths) from stationary antennae during a constant flux infiltration experiment. The formulation and coupling of the hydrological and geophysical models rely on a constant velocity wetting front that causes critical refraction at the edge of the front as it passes by the antennae. During this critical refraction period, the slope of the first arrival data can be used to calculate (1) the wetting velocity and (2) the hydraulic conductivity of the wet (or saturated) soil. If the soil is undersaturated during infiltration, then an estimate of the saturated water content is needed before calculating the saturated hydraulic conductivity. The hydraulic conductivity value is then used in a nonlinear global optimization scheme to estimate the remaining two parameters of a Broadbridge and White soil.     

Estimation of Leakage Rates Through Flexible Membrane Liners

Leakage rate calculations for both low- permeability soil liners and composite liners using flexible membrane liners (FMLs) overlying low-permeability soil are developed. Latin-Hypercube simulations with uncertainly assigned to the soil liner hydraulic conductivity value and the spatial frequency of FML holes are used to examine the variability in the liner leakage rates. The low-permeability soil hydraulic conductivity is the parameter with the greatest effect on landfill liner leakages rates. Composite liners have a significant impact on reducing leakage rates through the landfill liner.     

Free surface and saturated-unsaturated analyses of borehole infiltration tests above the water table

Constant head borehole infiltration tests are widely used for the in situ evaluation of saturated hydraulic conductivities of unsaturated soils above the water table. The formulae employed in analysing the results of such tests disregard the fact that some of the infiltrating water may flow under unsaturated conditions. Instead, these formulae are based on various approximations of the classical free surface theory which treats the flow region as if it were fully saturated and enclosed within a distinct envelope, the so-called ‘free surface’. A finite element model capable of solving free surface problems is used to examine the mathematical accuracy of the borehole infiltration formulae. The results show that in the hypothetical case where unsaturated flow does not exist, the approximate formulae are reasonably accurate within·a practical range of borehole conditions. To see what happens under conditions closer to those actually encountered in the field, the effect of unsaturated flow on borehole infiltration is investigated by means of two different numerical models: a mixed explicit-implicit finite element model, and a mixed explicit-implicit integrated finite difference model. Both of these models give nearly identical results; however, the integrated finite difference model is considerably faster than the finite element model. The relatively low computational efficiency of the finite element scheme is attributed to the large number of operations required in order to re-evaluate the conductivity (stiffness) matrix at each iteration in this highly non-linear saturated-unsaturated flow problem. The saturated-unsaturated analysis demonstrates that the classical free surface approach provides a distorted picture of the flow pattern in the soil. Contrary to what one would expect on the basis of this theory, only a finite region of the soil in the immediate vicinity of the borehole is saturated, whereas a significant percentage of the flow takes place under unsaturated conditions. As a consequence of disregarding unsaturated flow, the available formulae may underestimate the saturated hydraulic conductivity of fine grained soils by a factor of two, three, or more. Our saturated-unsaturated analysis leads to an improved design of borehole infiltration tests and a more accurate method for interpreting the results of such tests. The analysis also shows how one can predict the steady state rate of infiltration from data collected during the early transient period of the test.     

Effects of experimental uncertainty on the calculation of hillslope flow paths

Measurement uncertainty is a key hindrance to the quantification of water fluxes at all scales of investigation. Predictions of soil‐water flux rely on accurate or representative measurements of hydraulic gradients and field‐state hydraulic conductivity. We quantified the potential magnitude of errors associated with the parameters and variables used directly and indirectly within the Darcy – Buckingham soil‐water‐flux equation. These potential errors were applied to a field hydrometric data set collected from a forested hillslope in central Singapore, and their effect on flow pathway predictions was assessed. Potential errors in the hydraulic gradient calculations were small, approximately one order of magnitude less than the absolute magnitude of the hydraulic gradients. However, errors associated with field‐state hydraulic conductivity derivation were very large. Borehole (Guelph permeameter) and core‐based (Talsma ring permeameter) techniques were used to measure field‐saturated hydraulic conductivity. Measurements using these two approaches differed by up to 3\9 orders of magnitude, with the difference becoming increasingly marked within the B horizon. The sensitivity of the shape of the predicted unsaturated hydraulic conductivity curve to ±5% moisture content error on the moisture release curve was also assessed. Applied moisture release curve error resulted in hydraulic conductivity predictions of less than ±0\2 orders of magnitude deviation from the apparent conductivity. The flow pathways derived from the borehole saturated hydraulic conductivity approach suggested a dominant near‐surface flow pathway, whereas pathways calculated from the core‐based measurements indicated vertical percolation to depth. Direct tracer evidence supported the latter flow pathway, although tracer velocities were approximately two orders of magnitude smaller than the Darcy predictions. We conclude that saturated hydraulic conductivity is the critical hillslope hydrological parameter, and there is an urgent need to address the issues regarding its measurement further. Copyright © 2000 John Wiley & Sons, Ltd.     

Calibration of high-resolution geophysical data with tracer test measurements to improve hydrological predictions

Relationships between porosity and hydraulic conductivity tend to be strongly scale- and site-dependent and are thus very difficult to establish. As a result, hydraulic conductivity distributions inferred from geophysically derived porosity models must be calibrated using some measurement of aquifer response. This type of calibration is potentially very valuable as it may allow for transport predictions within the considered hydrological unit at locations where only geophysical measurements are available, thus reducing the number of well tests required and thereby the costs of management and remediation. Here, we explore this concept through a series of numerical experiments. Considering the case of porosity characterization in saturated heterogeneous aquifers using crosshole ground-penetrating radar and borehole porosity log data, we use tracer test measurements to calibrate a relationship between porosity and hydraulic conductivity that allows the best prediction of the observed hydrological behavior. To examine the validity and effectiveness of the obtained relationship, we examine its performance at alternate locations not used in the calibration procedure. Our results indicate that this methodology allows us to obtain remarkably reliable hydrological predictions throughout the considered hydrological unit based on the geophysical data only. This was also found to be the case when significant uncertainty was considered in the underlying relationship between porosity and hydraulic conductivity.     

Quantifying the effects of three-dimensional subsurface heterogeneity on Hortonian runoff processes using a coupled numerical,stochastic approach

The impact of three-dimensional subsurface heterogeneity in the saturated hydraulic conductivity on hillslope runoff generated by excess infiltration (so-called Hortonian runoff) is examined. A fully coupled, parallel subsurface–overland flow model is used to simulate runoff from an idealized hillslope. Ensembles of correlated, Gaussian random fields of saturated hydraulic conductivity are used to create uncertainty in spatial structure. A large number of cases are simulated in a parametric manner with the variance of the hydraulic conductivity varied over orders of magnitude. These cases include rainfall rates above, equal and below the geometric mean of the hydraulic conductivity distribution. These cases are also compared to theoretical representations of runoff production based on simple assumptions regarding (1) the rainfall rate and the value of hydraulic conductivity in the surface cell using a spatially-indiscriminant approach; and (2) a percolation-theory type approach to incorporate so-called runon. Simulations to test the ergodicity of hydraulic conductivity on hillslope runoff are also performed. Results show that three-dimensional stochastic representations of the subsurface hydraulic conductivity can create shallow perching, which has an important effect on runoff behavior that is different than previous two-dimensional analyses. The simple theories are shown to be very poor predictors of the fraction of saturated area that might runoff due to excess infiltration. It is also shown that ergodicity is reached only for a large number of integral scales (∼30) and not achieved for cases where the rainfall rate is less than the geometric mean of the saturated hydraulic conductivity.     

Analysis of subsurface contaminant migration and remediation using high performance computing

Highly resolved simulations of groundwater flow, chemical migration and contaminant recovery processes are used to test the applicability of stochastic models of flow and transport in a typical field setting. A simulation domain encompassing a portion of the upper saturated aquifer materials beneath the Lawrence Livermore National Laboratory was developed to hierarchically represent known hydrostratigraphic units and more detailed stochastic representations of geologic heterogeneity within them. Within each unit, Gaussian random field models were used to represent hydraulic conductivity variation, as parameterized from well test data and geologic interpretation of spatial variability. Groundwater flow, transport and remedial extraction of two hypothetical contaminants were made in six different statistical realizations of the system. The effective flow and transport behavior observed in the simulations compared reasonably with the predictions of stochastic theories based upon the Gaussian models, even though more exacting comparisons were prevented by inherent nonidealities of the geologic model and flow system. More importantly, however, biases and limitations in the hydraulic data appear to have reduced the applicability of the Gaussian representations and clouded the utility of the simulations and effective behavior based upon them. This suggests a need for better and unbiased methods for delineating the spatial distribution and structure of geologic materials and hydraulic properties in field systems. High performance computing can be of critical importance in these endeavors, especially with respect to resolving transport processes within highly variable media.©1998 Elsevier Science Limited. All rights reserved     

On the Relation Between Saturated Conductivity and Capillary Retention Characteristics

Abstract. A simple closed-form expression relating saturated hydraulic conductivity to the van Genuchten capillary retention model parameters is derived. Application of this equation to an experimental data set shows reasonable agreement between measured and predicted saturated conductivity values. The proposed equation provides a consistent theoretical basis for estimating both saturated and unsaturated hydraulic conductivity from statistical pore structure models.     

Effects of crab burrows on pore water flows in salt marshes

Macro-pores such as crab burrows are found commonly distributed in salt marsh sediments. Their disturbance on the soil structure is likely to influence both pore water flows and solute transport in salt marshes; however, the effects of crab burrows are not well understood. Here, a three-dimensional model simulated tidally driven pore water flows subject to the influence of crab burrows in a marsh system. The model, based on Richards’ equation, considered variably saturated flow in the marsh with a two-layer soil configuration, as observed at the Chongming Dongtan wetland (Shanghai, China). The simulation results showed that crab burrows distributed in the upper low-permeability soil layer, acting as preferential flow paths, affected pore water flows in the marsh particularly when the contrast of hydraulic conductivity between the lower high-permeability soil layer and the overlying low-permeability soils was high. The burrows were found to increase the volume of tidally driven water exchange between the marsh soil and the tidal creek. The simulations also showed improvement of soil aeration conditions in the presence of crab burrows. These effects may lead to increased productivity of the marsh ecosystem and enhancement of its material exchange with coastal waters.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

A New Rapid Method for Measuring the Vertical Head Profile

This study describes a new technique for measuring the head profile in a geologic formation. The technique provides rapid, low cost information on the depth of water‐producing zones and aquitards in heterogeneous aquifers, yielding estimates of hydraulic heads in each zone while identifying any potential for cross contamination between zones. The measurements can be performed in a typical borehole in just a few hours. The procedure uses both the continuous transmissivity profile obtained by the installation (eversion) of a flexible borehole liner into an open borehole and the subsequent removal (inversion) of the same liner from the borehole. The method is possible because of the continuous transmissivity profile (T profile described by Keller et al. 2014) obtained by measuring the rate of liner eversion under a constant driving head. The hydraulic heads of producing zones are measured using the reverse head profile (RHP) method (patent no. 9,008,971) based on a stepwise inversion of the borehole liner. As each interval of the borehole is uncovered by inversion of the liner, the head beneath the liner is allowed to equilibrate to a steady‐state value. The individual hydraulic heads contributing to each measurement are calculated using the measured transmissivity for each zone. Application of the RHP method to a sedimentary bedrock borehole in New Jersey verified that it reproduced the head distribution obtained the same day in the same borehole instrumented with a multilevel sampling system.     

Usefulness of flux measurements during hydraulic tomographic survey for mapping hydraulic conductivity distribution in a fractured medium

Using the first-order analysis, we investigate the spatial cross-correlation between hydraulic conductivity variation and specific discharge (flux) as well as its components measured in a borehole under steady-state flow conditions during cross-hole pumping tests in heterogeneous aquifers. These spatial correlation patterns are found to be quite different from that between the hydraulic conductivity variation and the hydraulic head measurement in the same borehole. This finding suggests that a specific discharge measurement carries non-redundant information about the spatial distribution of heterogeneity, even this measurement is collected from the same location where the head measurement is taken. As such, specific discharge observations should be included in the analysis of hydraulic tomography to increase the resolution of estimated aquifer heterogeneity. Using numerical experiments, we demonstrate the effectiveness of the joint interpretation of both hydraulic heads and fluxes for mapping fracture distributions in a hypothetic geologic medium.     

Comparison of upscaling methods to estimate hydraulic conductivity

A hydraulic field test program was performed at a hard rock laboratory (Asp? HRL) on the Swedish east coast to test upscaling theories. The test program investigated the rock volume around a borehole located at a depth of approximately 340 m below sea level. Hydraulic packer tests were performed at various scales, from 2 m to the entire borehole length of 296 m. From this set of data the predictive ability of different upscaling methods could be evaluated. The comparison of the evaluated "true" field scale hydraulic conductivity with the upscaled hydraulic conductivity yielded that the majority of the upscaling methods tested in this paper predict the large scale values with significant accuracy. However, the ability to predict rapidly decreases when the variance of the natural logarithm of hydraulic conductivity of the subsamples is larger than one. Such a variance is consistently found in the crystalline rocks at the tested site at the 2 m scale. However, at scales of 10 m and larger, a variance larger than one is uncommon. Therefore, it is concluded that there exists a smallest possible scale for use of hydraulic pumping test results for estimating the effective hydraulic conductivity at scales typical for regional flow.     

Macroporosity and infiltration in blanket peat: the implications of tension disc infiltrometer measurements

Little is known about the processes of infiltration and water movement in the upper layers of blanket peat. A tension infiltrometer was used to measure hydraulic conductivity in a blanket peat in the North Pennines, England. Measurements were taken from the surface down to 20 cm in depth for peat under four different vegetation covers. It was found that macropore flow is a significant pathway for water in the upper layers of this soil type. It was also found that peat depth and surface vegetation cover were associated with macroporosity and saturated hydraulic conductivity. The proportion of macropore flow was found to be greater at 5 cm depth than at 0, 10 and 20 cm depth. Peat beneath a Sphagnum cover tends to be more permeable and a greater proportion of macropore flow can occur beneath this vegetation type. Functional macroporosity and matrix flow in the near‐surface layers of bare peat appear to have been affected by weathering processes. Comparision of results with rainfall records demonstrates that infiltration‐excess overland flow is unlikely to be a common runoff‐generating mechanism on blanket peat; rather, a saturation‐excess mechanism combined with percolation‐excess above much less permeable layers dominates the runoff response. Copyright © 2001 John Wiley & Sons, Ltd.     

A geostatistically based inverse model for three-dimensional variably saturated flow

We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.     

A geostatistically based inverse model for three-dimensional variably saturated flow

We present a geostatistically based inverse model for characterizing heterogeneity in parameters of unsaturated hydraulic conductivity for three-dimensional flow. Pressure and moisture content are related to perturbations in hydraulic parameters through cross-covariances, which are calculated to first-order. Sensitivities needed for covariance calculations are derived using the adjoint state sensitivity method. Approximations of the conditional mean parameter fields are then obtained from the cokriging estimator. Correlation between parameters and pressure – moisture content perturbations is seen to be strongly dependent on mean pressure or moisture content. High correlation between parameters and pressure data was obtained under saturated or near saturated flow conditions, providing accurate estimation of saturated hydraulic conductivity, while moisture content measurements provided accurate estimation of the pore size distribution parameter under unsaturated flow conditions.     

Effects of Improper Characterization of Aquifer Thickness on Estimates of Hydraulic Conductivity from Pumping Tests

Pumping test data for surficial aquifers are commonly analyzed under the assumption that the base of the aquifer corresponds to the bottom of the test wells (i.e., the aquifer is truncated). This practice can lead to inaccurate hydraulic conductivity estimates, resulting from the use of low saturated thickness values with transmissivity estimates, and not accounting for the effects of partially penetrating wells. Theoretical time-drawdown data were generated at an observation well in a hypothetical unconfined aquifer for various values of saturated thickness and were analyzed by standard curve-matching techniques. The base of the aquifer was assumed to be the bottom of the pumping and observation wells. The overestimation of horizontal hydraulic conductivity was found to be directly proportional to the error in assumed saturated thickness, and to the (actual) ratio of vertical to horizontal hydraulic conductivity (K_v/K_h). Inaccurately high estimates of hydraulic conductivity obtained by aquifer truncation can lead to overestimates of ground water velocity and contaminant plume spreading, narrow capture zone configuration estimates, and overestimates of available ground water resources.     

In situ Mixing of Organic Matter Decreases Hydraulic Conductivity of Denitrification Walls in Sand Aquifers

In a previous study, a denitrification wall was constructed in a sand aquifer using sawdust as the carbon substrate. Ground water bypassed around this sawdust wall due to reduced hydraulic conductivity. We investigated potential reasons for this by testing two new walls and conducting laboratory studies. The first wall was constructed by mixing aquifer material in situ without substrate addition to investigate the effects of the construction technique (mixed wall). A second, biochip wall, was constructed using coarse wood chips to determine the effect of size of the particles in the amendment on hydraulic conductivity. The aquifer hydraulic conductivity was 35.4 m/d, while in the mixed wall it was 2.8 m/d and in the biochip wall 3.4 m/d. This indicated that the mixing of the aquifer sands below the water table allowed the particles to re-sort themselves into a matrix with a significantly lower hydraulic conductivity than the process that originally formed the aquifer. The addition of a coarser substrate in the biochip wall significantly increased total porosity and decreased bulk density, but hydraulic conductivity remained low compared to the aquifer. Laboratory cores of aquifer sand mixed under dry and wet conditions mimicked the reduction in hydraulic conductivity observed in the field within the mixed wall. The addition of sawdust to the laboratory cores resulted in a significantly higher hydraulic conductivity when mixed dry compared to cores mixed wet. This reduction in the hydraulic conductivity of the sand/sawdust cores mixed under saturated conditions repeated what occurred in the field in the original sawdust wall. This indicated that laboratory investigations can be a useful tool to highlight potential reductions in field hydraulic conductivities that may occur when differing materials are mixed under field conditions.