Modification of the Bouwer and Rice Method to a Cutoff Wall with a Filter Cake

The Bouwer and Rice method is a line-fitting method used to estimate the hydraulic conductivity of an aquifer by means of a slug test. When considering a relatively impermeable layer, called a filter cake, which may form at the interface between a cutoff wall and the natural soil formation, the assumptions of the Bouwer and Rice method are violated. A modification of the Bouwer and Rice method is proposed that incorporates the concept of a flow net, whereby the geometry of the cutoff wall and filter cake is effectively considered in estimating the hydraulic conductivity of a vertical cutoff wall.     

Tensor Hydraulic Conductivity Estimation for Heterogeneous Aquifers under Unknown Boundary Conditions

A physically based inverse method is developed using hybrid formulation and coordinate transform to simultaneously estimate hydraulic conductivity tensors, steady‐state flow field, and boundary conditions for a confined aquifer under ambient flow or pumping condition. Unlike existing indirect inversion techniques, the physically based method does not require forward simulations to assess model‐data misfits. It imposes continuity of hydraulic head and Darcy fluxes in the model domain while incorporating observations (hydraulic heads, Darcy fluxes, or well rates) at measurement locations. Given sufficient measurements, it yields a well‐posed inverse system of equations that can be solved efficiently with coarse grids and nonlinear optimization. When pumping and injection are active, well rates are used as measurements and flux sampling is not needed. The method is successfully tested on synthetic aquifer problems with regular and irregular geometries, different hydrofacies and flow patterns, and increasing conductivity anisotropy ratios. All problems yield stable inverse solutions under increasing head measurement errors. For a given set of observations, inversion accuracy is strongly affected by the conductivity anisotropy ratio. Conductivity estimation is also affected by flow pattern: within a hydrofacies, when Darcy flux component is very small, the corresponding directional conductivity perpendicular to streamlines becomes less identifiable. Finally, inversion is successful even if the location of aquifer boundaries is unknown. In this case, the inversion domain is defined by the location of the measurements.     

General Steady-State Shape Factor for a Partially Penetrating Well

We present the closed form of a general steady-state shape factor for a partially penetrating well in a uniform anisotropic aquifer. Our simple analytical expression for the shape factor has a uniform representation for full range of parameters and meets or exceeds the accuracy of previous results obtained through semiempirical methods (e.g., Bouwer and Rice [1976] equations). This general shape factor pertains to the flow of fluids (water or air) in subsurface formations when the upper formation boundary has constant potential and the lower boundary is impermeable. The results of our investigation are directly applicable to analyses of (1) slug tests with falling or rising head and (2) injection/extraction tests with constant head, essential techniques for the characterization of hydraulic conductivity of aquifers, streambeds, or lakebeds as well as the design of aquifer and soil remediation systems.     

Semi‐analytical solution for a slug test in partially penetrating wells including the effect of finite‐thickness skin

This paper presents a new semi‐analytical solution for a slug test in a well partially penetrating a confined aquifer, accounting for the skin effect. This solution is developed based on the solution for a constant‐flux pumping test and a formula given by Peres and co‐workers in 1989. The solution agrees with that of Cooper and co‐workers and the KGS model when the well is fully penetrating. The present solution can be applied to simulate the temporal and spatial head distributions in both the skin and formation zones. It can also be used to demonstrate the influences of skin type or skin thickness on the well water level and to estimate the hydraulic parameters of the skin and formation zones using a least‐squares approach. The results of this study indicate that the determination of hydraulic conductivity using a conventional slug‐test data analysis that neglects the presence of a skin zone will give an incorrect result if the aquifer has a skin zone. Copyright © 2008 John Wiley & Sons, Ltd.     

Groundwater response to tidal fluctuations in a leaky confined coastal aquifer with a finite length

The study on the hydraulic properties of coastal aquifers has significant implications both in hydrological sciences and environmental engineering. Although many analytical solutions are available, most of them are based on the same basic assumption that assumes aquifers extend landward semi‐infinitely, which does not necessarily reflect the reality. In this study, the general solutions for a leaky confined coastal aquifer have been developed that consider both finitely landward constant‐head and no‐flow boundaries. The newly developed solutions were then used to examine theoretically the joint effects of leakage and aquifer length on hydraulic head fluctuations within the leaky confined aquifer, and the validity of using the simplified solution, which assumes the aquifer is semi‐infinite. The results illustrated that the use of the simplified solution may cause significant errors, depending on joint effects of leakage and aquifer length. A dimensionless characteristic parameter was then proposed as an index for judging the applicability of the simplified solution. In addition, practical application of the general solution for the constant‐head inland boundary was used to characterize the hydraulic properties of a leaky confined aquifer using the data collected from a field site at the Seine River estuary, France, and the versatility of the general solution was further justified.     

Stochastic analysis of flux and head moments in a heterogeneous aquifer system

This study investigates the behavior of flux and head in a strongly heterogeneous three-dimensional aquifer system. The analyses relied on data from 520 slug tests together with 38,000 one-foot core intervals lithological data from the site of the General Separations Area in central Savannah River Site, South Carolina, USA. The skewness in the hydraulic conductivity histograms supported the geologic information for the top two aquifers, but revealed stronger clay content, than was reported for the bottom aquifer. The log-normal distribution model described adequately the hydraulic conductivity measurements for all three aquifers although, other distributions described equally well the bottom aquifer measurements. No apparent anisotropy on the horizontal plane was found for the three aquifers, but ratios of horizontal to vertical correlation lengths between 33 and 75 indicated a strong stratification at the site. Three-dimensional Monte Carlo stochastic simulations utilized a grid with larger elements than the support volume of measurements, but of sub-REV (representative elementary volume) dimensions. This necessitated, on one hand, the use of upscaled hydraulic conductivity expressions, but on the other hand did not allow for the use of anisotropic effective hydraulic conductivity expressions (Sarris and Paleologos in J Stoch Environ Res Risk Assess 18: 188–197, 2004). Flux mean and standard deviations components were evaluated on three vertical cross-sections. The mean and variance of the horizontal flux component normal to a no-flow boundary tended to zero at approximately two to three integral scales from that boundary. Close to a prescribed head boundary both the mean and variance of the horizontal flux component normal to the boundary increased from a stable value attained at a distance of about five integral scales from that boundary. The velocity field 〈q_x〉 was found to be mildly anisotropic in the top two aquifers, becoming highly anisotropic in the bottom aquifer; 〈q_y〉 was anisotropic in all three aquifers with directions of high continuity normal to those of the 〈q_x〉 field; finally, 〈q_z〉 was highly anisotropic in all three aquifers, with higher continuity along the east–west direction. The mean head field was found to be continuous, despite the high heterogeneity of the underlying hydraulic conductivity field. Directions of high continuity were in alignment with field boundaries and mean flow direction. Conditioning did not influence significantly the expected value of the flux terms, with more pronounced being the effect on the standard deviation of the flux vector components. Conditioning reduced the standard deviations of the horizontal flux components by as much as 50% in the bottom aquifer. Variability in the head cross-sections was affected only marginally, with an average 10% reduction in the respective standard deviation. Finally, the location of the conditioning data did not appear to have a significant effect on the surrounding area, with uniform reduction in standard deviations.     

Obtaining the steady‐state drawdown solutions of constant‐head and constant‐flux tests

The solutions of constant‐head and constant‐flux tests are commonly used to predict the temporal or spatial drawdown distribution or to determine aquifer parameters. Theis and Thiem equations, for instance, are well‐known transient and steady‐state drawdown solutions, respectively, of the constant‐flux test. It is known that the Theis equation is not applicable to the case where the aquifer has a finite boundary or the pumping time tends to infinity. On the other hand, the Thiem equation does not apply to the case where the aquifer boundary is infinite. However, the issue of obtaining the Thiem equation from the transient drawdown solution has not previously been addressed. In this paper, the drawdown solutions for constant‐head and constant‐flux tests conducted in finite or infinite confined aquifers with or without consideration of the effect of the well radius are examined comprehensively. Mathematical verification and physical interpretation of the solutions to these two tests converging or not converging to the Thiem equation are presented. The result shows that there are some finite‐domain solutions for these two tests that can converge to the Thiem equation when the time becomes infinitely large. In addition, the time criteria to give a good approximation to the finite‐domain solution by the infinite‐domain solution and the Thiem equation are investigated and presented. Copyright © 2008 John Wiley & Sons, Ltd.     

New analytical solutions for groundwater flow in wedge-shaped aquifers with various topographic boundary conditions

The hydraulic head distribution in a wedge-shaped aquifer depends on the wedge angle and the topographic and hydrogeological boundary conditions. In addition, an equation in terms of the radial distance with trigonometric functions along the boundary may be suitable to describe the water level configuration for a valley flank with a gentle sloping and rolling topography. This paper develops a general mathematical model including the governing equation and a variety of boundary conditions for the groundwater flow within a wedge-shaped aquifer. Based on the model, a new closed-form solution for transient flow in the wedge-shaped aquifer is derived via the finite sine transform and Hankel transform. In addition, a numerical approach, including the roots search scheme, the Gaussian quadrature, and Shanks’ method, is proposed for efficiently evaluating the infinite series and the infinite integral presented in the solution. This solution may be used to describe the head distribution for wedges that image theory is inapplicable, and to explore the effects of the recharge from various topographic boundaries on the groundwater flow system within a wedge-shaped aquifer.     

A Leaky Aquifer below Champlain Sea Clay: Closed‐Form Solutions for Natural Seepage

Closed‐form solutions are proposed for natural seepage in semiconfined (leaky) aquifers such as those existing below the massive Champlain Sea clay layers in the Saint‐Lawrence River Valley. The solutions are for an ideal horizontal leaky aquifer below an ideal aquitard that may have either a constant thickness and a constant hydraulic head at its surface, or a variable thickness and a variable hydraulic head at its surface. A few simplifying assumptions were needed to obtain the closed‐form solutions. These have been verified using a finite element method, which did not make any of the assumptions but gave an excellent agreement for hydraulic heads and groundwater velocities. For example, the difference between the two solutions was smaller than 1 mm for variations in the 5 to 8 m range for the hydraulic head in the semiconfined aquifer. Note that fitting the hydraulic head data of monitoring wells to the theoretical solutions gives only the ratio of the aquifer and aquitard hydraulic conductivities, a clear case of multiple solutions for an inverse problem. Consequently, field permeability tests in the aquitard and the aquifer, and pumping tests in the aquifer, are still needed to determine the hydraulic conductivity values.     

Modeling Steady Sea Water Intrusion with Single‐Density Groundwater Codes

Steady interface flow in heterogeneous aquifer systems is simulated with single‐density groundwater codes by using transformed values for the hydraulic conductivity and thickness of the aquifers and aquitards. For example, unconfined interface flow may be simulated with a transformed model by setting the base of the aquifer to sea level and by multiplying the hydraulic conductivity with 41 (for sea water density of 1025 kg/m³). Similar transformations are derived for unconfined interface flow with a finite aquifer base and for confined multi‐aquifer interface flow. The head and flow distribution are identical in the transformed and original model domains. The location of the interface is obtained through application of the Ghyben‐Herzberg formula. The transformed problem may be solved with a single‐density code that is able to simulate unconfined flow where the saturated thickness is a linear function of the head and, depending on the boundary conditions, the code needs to be able to simulate dry cells where the saturated thickness is zero. For multi‐aquifer interface flow, an additional requirement is that the code must be able to handle vertical leakage in situations where flow in an aquifer is unconfined while there is also flow in the aquifer directly above it. Specific examples and limitations are discussed for the application of the approach with MODFLOW. Comparisons between exact interface flow solutions and MODFLOW solutions of the transformed model domain show good agreement. The presented approach is an efficient alternative to running transient sea water intrusion models until steady state is reached.     

Hydraulics of a partially penetrating well with skin zone in a confined aquifer

A steady/quasi-steady model is developed for predicting flow into a partially penetrating well with skin zone in a confined aquifer overlying an impervious layer. The model takes into account flow through the bottom of the wellbore, finite skin thickness and finite horizontal and vertical extent of the aquifer. Moreover, the solution can be easily extended to include the mixed-type boundary condition at the well face, where a Dirichlet in the form of a specified hydraulic head and a Neumann in the form of zero flux coexist at the same time at different portions of the well face. The validity of the proposed solution is tested by comparing a few results obtained from the developed model with corresponding results obtained by analytical and numerical means. The study shows that, among other factors remaining constant, both the horizontal and vertical extent of an artesian aquifer, thickness of the skin zone, bottom flow and conductivity contrast of the skin and formation zones, play an important part in deciding flow to a well dug in the aquifer, and hence these factors must be considered while analyzing the problem. The model proposed here can be used to estimate skin thickness as well as hydraulic conductivities of the skin and formation zones of a well with skin zone in an artesian aquifer underlain by an impervious layer by utilizing pumping test data falling in the steady or quasi-steady state of a typical pumping test. As the proposed solution is of a general nature in the sense that it can handle, apart from partial penetration and bottom flow, the finite size skin zone and finite horizontal and vertical extent of an artesian aquifer together with the mixed-type boundary condition at the well face, it is hoped that the predictions coming out of the model will be more realistic than those obtained using solutions developed with more stringent assumptions.     

Estimating finite difference block equivalent hydraulic conductivity for numerically solving the Richards' equation

A method is proposed for calculating the equivalent hydraulic conductivity (EHC) within a finite difference block (FDB). Application of the constant‐flux assumption of Darcy's Law, the EHC equals to the integration of effective hydraulic conductivity (K_w) as a function of pressure head (h_w) divided by the head difference at the ends of the FDB. Error analysis show that the constant‐flux (CF) EHC estimates are better than those computed by the commonly used arithmetic‐mean (AM), geometric‐mean (GM), and harmonic‐mean (HM) techniques. CF EHC results are even more superior at larger interblock head difference situations. Simulations of water infiltration experiments show that simulations using the CF EHC or AM or GM weighting technique have only slight difference while applying the Neumann type boundary condition at the ground surface. In case of the Dirichlet type boundary condition, however, the CF EHC is superior to the other two in correctly estimating the depth of infiltration while enlarging the grid size. Therefore, it is recommended to adopt the CF EHC with a larger grid size to the more stable and more efficient results. Copyright © 2007 John Wiley & Sons, Ltd.     

Wellbore flow‐rate solution for a constant‐head test in two‐zone finite confined aquifers

The solution describing the wellbore flow rate in a constant‐head test integrated with an optimization approach is commonly used to analyze observed wellbore flow‐rate data for estimating the hydrogeological parameters of low‐permeability aquifers. To our knowledge, the wellbore flow‐rate solution for the constant‐head test in a two‐zone finite‐extent confined aquifer has never been reported so far in the literature. This article is first to develop a mathematical model for describing the head distribution in the two‐zone aquifer. The Laplace domain solutions for the head distributions and wellbore flow rate in a two‐zone finite confined aquifer are derived using the Laplace transform, and their corresponding time domain solutions are then obtained using the Bromwich integral method and residue theorem. These new solutions are expressed in terms of an infinite series with Bessel functions and not straightforward to calculate numerically. A large‐time solution for the wellbore flow rate is therefore developed by employing the relationship of small Laplace variable versus large time variable and L'Hospital's rule. The result shows that the large‐time solution is identical to the steady‐state solution obtained after applying the Tauberian theorem into the Laplace domain solution. This large‐time solution can reduce to the Thiem equation in the case of no skin. Finally, the newly developed solution is used to investigate the effects of outer boundary distance and conductivity ratio on the wellbore flow rate. Copyright © 2011 John Wiley & Sons, Ltd.     

A Correction on Coastal Heads for Groundwater Flow Models

We introduce a simple correction to coastal heads for constant‐density groundwater flow models that contain a coastal boundary, based on previous analytical solutions for interface flow. The results demonstrate that accurate discharge to the sea in confined aquifers can be obtained by direct application of Darcy's law (for constant‐density flow) if the coastal heads are corrected to ((α + 1)/α)h_s ? B/2α, in which h_s is the mean sea level above the aquifer base, B is the aquifer thickness, and α is the density factor. For unconfined aquifers, the coastal head should be assigned the value . The accuracy of using these corrections is demonstrated by consistency between constant‐density Darcy's solution and variable‐density flow numerical simulations. The errors introduced by adopting two previous approaches (i.e., no correction and using the equivalent fresh water head at the middle position of the aquifer to represent the hydraulic head at the coastal boundary) are evaluated. Sensitivity analysis shows that errors in discharge to the sea could be larger than 100% for typical coastal aquifer parameter ranges. The location of observation wells relative to the toe is a key factor controlling the estimation error, as it determines the relative aquifer length of constant‐density flow relative to variable‐density flow. The coastal head correction method introduced in this study facilitates the rapid and accurate estimation of the fresh water flux from a given hydraulic head measurement and allows for an improved representation of the coastal boundary condition in regional constant‐density groundwater flow models.     

Numerical method of moments for solute transport in a nonstationary flow field conditioned on hydraulic conductivity and head measurements

The unconditional stochastic studies on groundwater flow and solute transport in a nonstationary conductivity field show that the standard deviations of the hydraulic head and solute flux are very large in comparison with their mean values (Zhang et al. in Water Resour Res 36:2107–2120, 2000; Wu et al. in J Hydrol 275:208–228, 2003; Hu et al. in Adv Water Resour 26:513–531, 2003). In this study, we develop a numerical method of moments conditioning on measurements of hydraulic conductivity and head to reduce the variances of the head and the solute flux. A Lagrangian perturbation method is applied to develop the framework for solute transport in a nonstationary flow field. Since analytically derived moments equations are too complicated to solve analytically, a numerical finite difference method is implemented to obtain the solutions. Instead of using an unconditional conductivity field as an input to calculate groundwater velocity, we combine a geostatistical method and a method of moment for flow to conditionally simulate the distributions of head and velocity based on the measurements of hydraulic conductivity and head at some points. The developed theory is applied in several case studies to investigate the influences of the measurements of hydraulic conductivity and/or the hydraulic head on the variances of the predictive head and the solute flux in nonstationary flow fields. The study results show that the conditional calculation will significantly reduce the head variance. Since the hydraulic head measurement points are treated as the interior boundary (Dirichlet boundary) conditions, conditioning on both the hydraulic conductivity and the head measurements is much better than conditioning only on conductivity measurements for reduction of head variance. However, for solute flux, variance reduction by the conditional study is not so significant.     

Constraints on hydraulic parameters and implications for groundwater flux across the upland–estuary interface

Analyses of independent laboratory- and field-scale measurements from two sites on Sapelo Island, Georgia reveal heterogeneity in hydraulic parameters across the upland–estuary interface. Regardless of the method used (short-duration pumping tests, amplitude attenuation of tidal pumping data, sediment grain size distributions, and falling head permeameter tests), we obtain hydraulic conductivity of 10⁻⁴ m s⁻¹ for the fine-grained, well-sorted, clean sands that make up the upland areas. Proximal to the upland–estuary boundary, the tidal pumping analyses and permeameter tests suggest that hydraulic conductivities decrease by more than two orders of magnitude, a result consistent with the presence of a clogging layer. Such a clogging layer may arise due to a variety of physical, chemical, or biological processes. The extent and orientation of the layers of reduced hydraulic conductivity near the upland–estuary boundary influence the nature of the aquifer's response to tidal forcing. Where the lower conductivity layer forms a relatively flat creek bank, tidal pumping produces a primarily mechanical response in the adjacent aquifer. Where the creek bank is nearly vertical, there is a more direct hydraulic connection between the tidal creek and the adjacent aquifer. The clogging layer likely contributes to the development of complicated flow pathways across the upland–estuary boundary. Effective flow paths calculated from tidal pumping data terminate within the marsh, beyond the boundary of the upland aquifer, suggesting a diffuse regime of groundwater discharge in the marsh. We postulate that, in many settings, submarsh flow may be as important as seepage faces for groundwater discharge into the marsh–estuary complex.     

Transient Solutions to Groundwater Mounding in Bounded and Unbounded Aquifers

In this study, the well‐known Hantush solution procedure for groundwater mounding under infinitely long infiltration strips is extended to finite and semi‐infinite aquifer cases. Initially, the solution for infinite aquifers is presented and compared to those available in literature and to the numerical results of MODFLOW. For the finite aquifer case, the method of images, which is commonly used in well hydraulics, is used to be able to represent the constant‐head boundaries at both sides. It is shown that a finite number of images is enough to obtain the results and sustain the steady state. The effect of parameters on the growth of the mound and on the time required to reach the steady state is investigated. The semi‐infinite aquifer case is emphasized because the growth of the mound is not symmetric. As the constant‐head boundary limits the growth, the unbounded side grows continuously. For this reason, the groundwater divide shifts toward the unbounded side. An iterative solution procedure is proposed. To perform the necessary computations a code was written in Visual Basic of which the algorithm is presented. The proposed methodology has a wide range of applicability and this is demonstrated using two practical examples. The first one is mounding under a stormwater dispersion trench in an infinite aquifer and the other is infiltration from a flood control channel into a semi‐infinite aquifer. Results fit very well with those of MODFLOW.     

Combining 3D Hydraulic Tomography with Tracer Tests for Improved Transport Characterization

Hydraulic tomography (HT) is a method for resolving the spatial distribution of hydraulic parameters to some extent, but many details important for solute transport usually remain unresolved. We present a methodology to improve solute transport predictions by combining data from HT with the breakthrough curve (BTC) of a single forced‐gradient tracer test. We estimated the three dimensional (3D) hydraulic‐conductivity field in an alluvial aquifer by inverting tomographic pumping tests performed at the Hydrogeological Research Site Lauswiesen close to Tübingen, Germany, using a regularized pilot‐point method. We compared the estimated parameter field to available profiles of hydraulic‐conductivity variations from direct‐push injection logging (DPIL), and validated the hydraulic‐conductivity field with hydraulic‐head measurements of tests not used in the inversion. After validation, spatially uniform parameters for dual‐domain transport were estimated by fitting tracer data collected during a forced‐gradient tracer test. The dual‐domain assumption was used to parameterize effects of the unresolved heterogeneity of the aquifer and deemed necessary to fit the shape of the BTC using reasonable parameter values. The estimated hydraulic‐conductivity field and transport parameters were subsequently used to successfully predict a second independent tracer test. Our work provides an efficient and practical approach to predict solute transport in heterogeneous aquifers without performing elaborate field tracer tests with a tomographic layout.     

Flow to an eccentric well in a leaky circular aquifer with varied lateral replenishment

Summary An analytical solution is obtained for the flow to an eccentric well in a leaky circular aquifer with lateral replenishment, both for steady and unsteady cases. The flows for external boundary conditions of constant head and zero flux, which were treated previously, follow in the limit from a more general boundary condition. Graphs are developed to show the influence of vertical leakage and lateral replenishment on the relationship between drawdown at the well and eccentricity.Other symbols are defined in the text as they occur.     

Applicability of spectral analysis to determine hydraulic diffusivity

This study is to evaluate the applicability of estimating the one-dimensional horizontal hydraulic diffusivity of an unconfined aquifer with time-dependent fluctuation of lateral head and vertical recharge boundaries using observed water level spectra. Different models of boundary condition are imposed to evaluate the statistical significance between the calculated hydraulic diffusivity (ξ^′) with the given hydraulic diffusivity (ξ). The auto-spectra of the water level in observation wells tapping the same aquifer are closely related to those at the disturbed boundaries. For an aquifer with a constant hydraulic diffusivity, the water level fluctuation in the monitoring wells is linearly related to the water level spectra observed at the boundaries. The spectral density function of aquifer hydraulic head varies inversely with specific yield (S _y) and directly with recharge. Given small variation in water level spectra at the disturbed boundaries, the water level fluctuation in the aquifer is affected by the recharge condition and the aquifer spectral density function is sensitive to S _y. Using an iterative technique to estimate ξ from 1400 sets of given parameters, 99% of the ξ^′/ξ values deviated within only one order of magnitude with the model length (L) being equal to 1 km and 10 km. For L equal to 100 m, approximately 82% of the ξ^′/ξ population falls within two orders of magnitude. Therefore, spectral analysis of aquifer hydraulic head response can be used to estimate the hydraulic diffusivity of an unconfined aquifer which is affected by periodic variations in recharge and head at boundaries.