A two‐dimensional analytical solution for groundwater flow in a leaky confined aquifer system near open tidal water

Groundwater in coastal areas is commonly disturbed by tidal fluctuations. A two‐dimensional analytical solution is derived to describe the groundwater fluctuation in a leaky confined aquifer system near open tidal water under the assumption that the groundwater head in the confined aquifer fluctuates in response to sea tide whereas that of the overlying unconfined aquifer remains constant. The analytical solution presented here is an extension of the solution by Sun for two‐dimensional groundwater flow in a confined aquifer and the solution by Jiao and Tang for one‐dimensional groundwater flow in a leaky confined aquifer. The analytical solution is compared with a two‐dimensional finite difference solution. On the basis of the analytical solution, the groundwater head distribution in a leaky confined aquifer in response to tidal boundaries is examined and the influence of leakage on groundwater fluctuation is discussed. Copyright © 2001 John Wiley & Sons, Ltd.     

A two-dimensional analytical solution of groundwater responses to tidal loading in an estuary and ocean

Previous studies on tidal dynamics of coastal aquifers have focussed on the inland propagation of oceanic tides in the cross-shore direction, a configuration that is essentially one-dimensional. Aquifers at natural coasts can also be influenced by tidal waves in nearby estuaries, resulting in a more complex behaviour of head fluctuations in the aquifers. We present an analytical solution to the two-dimensional depth-averaged groundwater flow equation for a semi-infinite aquifer subject to oscillating head conditions at the boundaries. The solution describes the tidal dynamics of a coastal aquifer that is adjacent to a cross-shore estuary. Both the effects of oceanic and estuarine tides on the aquifer are included in the solution. The analytical prediction of the head fluctuations is verified by comparison with numerical solutions computed using a standard finite-difference method. An essential feature of the present analytical solution is the interaction between the cross- and along-shore tidal waves in the aquifer area near the estuary’s entry. As the distance from the estuary or coastline increases, the wave interaction is weakened and the aquifer response is reduced, respectively, to the one-dimensional solution for oceanic tides or the solution of Sun (Sun H. A two-dimensional analytical solution of groundwater response to tidal loading in an estuary, Water Resour Res 1997;33:1429–35) for two-dimensional non-interacting tidal waves.     

Series solutions for flow in stratified aquifers with natural geometry

An analytical series solution method is presented for modeling regional steady-state groundwater flow in a two-dimensional stratified aquifer cross-section where the water table is well-characterized. The aquifer system may have any number of contiguous or non-contiguous layers and the geometry of each layer is restricted only by the requirement that the elevation of the stratigraphic unconformities between layers is a function of the x-coordinate alone. Various techniques may be used to handle pinching layers, faults, and other discontinuities. The solutions are obtained by minimizing head and flow continuity errors between layers and errors in the Dirichlet surface at a set of control points along these unconformities; the governing equation is met exactly. The solutions are derived and demonstrated on multiple test cases. The errors for some specific, geometrically challenging cases are assessed and discussed.     

A semi‐analytical solution for groundwater responses to stream‐stage variations and tidal fluctuations in a coastal aquifer

A two‐dimensional semi‐analytical solution to analyse stream–aquifer interactions in a coastal aquifer where groundwater level responds to tidal effects is presented. The conceptual model considered is a two‐dimensional subsurface system with stream and coastline boundaries at right angles. The dimensional and non‐dimensional boundary value problems were solved for water level in the aquifer by successive application of Laplace and Fourier transform techniques, and the results were obtained by numerical inversion of the transformed solution. The solution was then verified by reducing the solutions to one‐dimensional known problems and comparing the results with those from previous studies. Hypothetical examples were used to examine the characteristics of water‐level variations due to the variations in stream stage and the fluctuations in tide level. Sensitivity analysis indicated that streambed leakance has no influence over the amplitude of groundwater fluctuations, but that the effect of stream stage increases with increasing leakance. Little difference was observed in the water level for different aquifer penetration ratios with narrow stream width. Increases in streambed leakance caused increases in the effect of aquifer penetration by the stream on the water level. An increased specific yield value resulted in decreased amplitude of water fluctuations and mean water level, and showed that water‐level variations due to stream and tidal boundaries are sensitive to specific yield. Copyright © 2006 John Wiley & Sons, Ltd.     

A conceptual–numerical model to simulate hydraulic head in aquifers that are hydraulically connected to surface water bodies

In this paper, we present a conceptual‐numerical model that can be deduced from a calibrated finite difference groundwater‐flow model, which provides a parsimonious approach to simulate and analyze hydraulic heads and surface water body–aquifer interaction for linear aquifers (linear response of head to stresses). The solution of linear groundwater‐flow problems using eigenvalue techniques can be formulated with a simple explicit state equation whose structure shows that the surface water body–aquifer interaction phenomenon can be approached as the drainage of a number of independent linear reservoirs. The hydraulic head field could be also approached by the summation of the head fields, estimated for those reservoirs, defined over the same domain set by the aquifer limits, where the hydraulic head field in each reservoir is proportional to a specific surface (an eigenfunction of an eigenproblem, or an eigenvector in discrete cases). All the parameters and initial conditions of each linear reservoir can be mathematically defined in a univocal way from the calibrated finite difference model, preserving its characteristics (geometry, boundary conditions, hydrodynamic parameters (heterogeneity), and spatial distribution of the stresses). We also demonstrated that, in practical cases, an accurate solution can be obtained with a reduced number of linear reservoirs. The reduced computational cost of these solutions can help to integrate the groundwater component within conjunctive use management models. Conceptual approximation also facilitates understanding of the physical phenomenon and analysis of the factors that influence it. A simple synthetic aquifer has been employed to show how the conceptual model can be built for different spatial discretizations, the parameters required, and their influence on the simulation of hydraulic head fields and stream–aquifer flow exchange variables. A real‐world case was also solved to test the accuracy of the proposed approaches, by comparing its solution with that obtained using finite‐difference MODFLOW code. Copyright © 2011 John Wiley & Sons, Ltd.     

Boundary depletion rate and drawdown in leaky wedge‐shaped aquifers

This study presents analytical solutions of the three‐dimensional groundwater flow to a well in leaky confined and leaky water table wedge‐shaped aquifers. Leaky wedge‐shaped aquifers with and without storage in the aquitard are considered, and both transient and steady‐state drawdown solutions are derived. Unlike the previous solutions of the wedge‐shaped aquifers, the leakages from aquitard are considered in these solutions and unlike similar previous work for leaky aquifers, leakage from aquitards and from the water table are treated as the lower and upper boundary conditions. A special form of finite Fourier transforms is used to transform the z‐coordinate in deriving the solutions. The leakage induced by a partially penetrating pumping well in a wedge‐shaped aquifer depends on aquitard hydraulic parameters, the wedge‐shaped aquifer parameters, as well as the pumping well parameters. We calculate lateral boundary dimensionless flux at a representative line and investigate its sensitivity to the aquitard hydraulic parameters. We also investigate the effects of wedge angle, partial penetration, screen location and piezometer location on the steady‐state dimensionless drawdown for different leakage parameters. Results of our study are presented in the form of dimensionless flux‐dimensionless time and dimensionless drawdown‐leakage parameter type curves. The results are useful for evaluating the relative role of lateral wedge boundaries and leakage source on flow in wedge‐shaped aquifers. This is very useful for water management problems and for assessing groundwater pollution. The presented analytical solutions can also be used in parameter identification and in calculating stream depletion rate and volume. Copyright © 2011 John Wiley & Sons, Ltd.     

Efficient conceptual model for simulating the effect of aquifer heterogeneity on natural groundwater discharge to rivers

When linearity can be assumed (linear response of heads to stresses), stream–aquifer flow exchange can be simulated as the drainage of a number of independent linear reservoirs. This conceptual model, which can be mathematically deduced in a univocal way from an eigenvalue solution of the linear groundwater flow problem, facilitates the understanding of the physical phenomenon and the analysis of influencing factors. The number of reservoirs required to simulate stream depletion in some ideal homogeneous cases of stream–aquifer connection was analyzed in detail in a previous investigation using analytical eigenvalue solutions [16]. However, most aquifers are heterogeneous in nature and numerical solutions must be employed to analyze whether they could also be simulated using few reservoirs. This paper presents a stochastic analysis of the influence of heterogeneity on the simulation of natural groundwater discharges in aquifers connected to rivers, as a series of linear reservoirs. A Monte-Carlo approach was employed to perform this study. The results show that, on a monthly time scale, many cases (even heterogeneous aquifers) can be simulated using just a few reservoirs with sufficient accuracy and at minimum computational cost. Therefore, this modeling technique can be useful to efficiently simulate the integrated management of complex water resources systems at the basin scale (with many aquifers, reservoirs, demands, etc.) that need to simultaneously consider surface and groundwater flow and stream–aquifer interaction.     

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.     

Effective hydraulic conductivity of nonstationary aquifers

Assuming that the ln hydraulic conductivity in an aquifer is mathematically approximated by a spatial deterministic surface, or trend, plus a stationary random noise, we treat the problem of finding what the effective hydraulic conductivity of that aquifer is. This problem is tackled by spectral methods applied to a type of diffusion equation of groundwater flow, together with suitable coordinate transformations. Analytical (exact) solutions in terms of elementary functions are presented for one- and three-dimensional finite and infinite domains. Stability criteria are obtained for the solutions, in terms of a critical parameter, that turns out to involve the product of correlation scale and trend gradient. For the case of finite and symmetrical domains, additional provisions to insure the stability of numerical calculations of effective hydraulic conductivity are provided. Effective hydraulic conductivity is an important property, with potential applications in the calibrations of groundwater and transport numerical models.     

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.     

Analytical and semi-analytical solutions of horizontal well capture times under no-flow and constant-head boundaries

The closed-form analytical solutions and semi-analytical solutions of capture times to horizontal wells are derived for different recovery scenarios. The capture time is the time a fluid particle takes to flow to the well. The first scenario is recovery from a confined aquifer in which the influence of regional groundwater flow upon the capture time is included. The second scenario is recovery from underneath a water reservoir in which the top boundary of the aquifer is constant-head. The third scenario is recovery from a low-permeability layer bounded above and below by much higher permeability media. Closed-form solutions are provided for the cases with: (1) a center or a bottom well for the first scenario; (2) a bottom well for the second scenario; and (3) a center well for the third scenario. Semi-analytical solutions are provided for general well locations for those scenarios. Solutions for both isotropic and anisotropic media are studied. These solutions can be used as quick references to calculate the capture times, and as benchmarks to validate numerical solutions. The limitations of the analytical solutions are analyzed. Our results show that the top and bottom no-flow boundaries of an aquifer constrain the vertical flow, but enhance the horizontal flow, resulting in elongated iso-capture time curves. When constant-head boundaries are presented, water can infiltrate vertically across those boundaries to replenish the aquifers, resulting in less elongated iso-capture time curves.     

Boundary condition effects on maximum groundwater withdrawal in coastal aquifers

Prevention of sea water intrusion in coastal aquifers subject to groundwater withdrawal requires optimization of well pumping rates to maximize the water supply while avoiding sea water intrusion. Boundary conditions and the aquifer domain size have significant influences on simulating flow and concentration fields and estimating maximum pumping rates. In this study, an analytical solution is derived based on the potential-flow theory for evaluating maximum groundwater pumping rates in a domain with a constant hydraulic head landward boundary. An empirical correction factor, which was introduced by Pool and Carrera (2011) to account for mixing in the case with a constant recharge rate boundary condition, is found also applicable for the case with a constant hydraulic head boundary condition, and therefore greatly improves the usefulness of the sharp-interface analytical solution. Comparing with the solution for a constant recharge rate boundary, we find that a constant hydraulic head boundary often yields larger estimations of the maximum pumping rate and when the domain size is five times greater than the distance between the well and the coastline, the effect of setting different landward boundary conditions becomes insignificant with a relative difference between two solutions less than 2.5%. These findings can serve as a preliminary guidance for conducting numerical simulations and designing tank-scale laboratory experiments for studying groundwater withdrawal problems in coastal aquifers with minimized boundary condition effects.     

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.     

Recent advances in modeling of well hydraulics

Well hydraulics is a discipline to understand the process of flow to the well in an aquifer which is regarded as a source of groundwater. A variety of analytical and numerical models have been developed over the last few decades to provide a framework for understanding and quantifying the flow behavior in aquifer systems. In this review, we first briefly introduce the background of the theory of well hydraulics and the concepts, methodologies, and applications of analytical, semi-analytical, numerical and approximate methods in solving the well-hydraulic problems. We then address the subjects of current interests such as the incorporation of effects of finite well radius, wellbore storage, well partial penetration, and the presence of skin into various practical problems of groundwater flow. Furthermore, we also summarize recent developments of flow modeling such as the flow in aquifers with horizontal wells or collector wells, the capture zone delineation, and the non-Darcian flow in porous media and fractured formations. Finally, we present a comprehensive review on the numerical calculations for five well functions frequently appearing in well-hydraulic literature and suggest some topics in groundwater flow for future research.     

Understanding and modelling spatial drain–aquifer interactions in a low‐lying coastal aquifer—the Thurne catchment,Norfolk, UK

Seawater intrusion into fresh groundwater formations generally results inadvertently from human activities, such as over‐abstraction from coastal aquifers. This article describes the data analysis to quantify drain–aquifer interactions in a low‐lying pump‐drained coastal aquifer, which is subject to saline intrusion due to widespread land drainage, and the resulting development and application of a numerical groundwater model to understand the spatial groundwater system behaviour (including groundwater salinity fluxes). Without measured flow data in this pump‐drained catchment, a novel groundwater head‐dependent approach to hydrograph separation is described. Time‐variant and time‐invariant MODFLOW analyses are utilised to examine the flow processes. A new approach to calculate drain coefficients, which represent the extensive network of drainage ditches in the regional model, using field information, is described; the sum of the drainage coefficients are close to the values independently estimated from the head‐dependent hydrograph separation. Results show that (1) the groundwater flows into the drainage systems are well reproduced using the new drain coefficients, (2) particle tracking of fresh and saline water can explain observed spatial salinity distribution within drainage networks and (3) the modelled flow of seawater across the coast is approximately 25% greater than that discharged by the pumps, demonstrating the need for drainage management to be aware of the slow response of groundwater systems to past drainage system changes. The article demonstrates that numerical groundwater modelling can produce the improved understanding needed to inform management decisions in such complex environments. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of heterogeneity on spatiotemporal variations of groundwater level in a bounded unconfined aquifer

Spatiotemporal variations of groundwater level due to a white noise recharge time series and a random transmissivity field in a bounded unconfined aquifer was studied. The analytical solutions for the variance and covariance of groundwater level were derived with non-stationary spectral analyses and superposition principle. It was found that the fluctuations of groundwater level are spatially non-stationary due to a fixed head boundary condition and temporal non-stationary at early time but gradually became stationary as time progresses due to effect of the initial condition. The variation in groundwater level is mainly caused by the random source/sink in the case of temporally random recharge and spatially random transmissivity. The effect of heterogeneity is to increase the variation of groundwater level and the maximum effect occurs close to the constant head boundary because of the linear mean hydraulic gradient. The heterogeneity also enhances the correlation of groundwater level, especially at large time intervals and small spatial distances.     

Conformal mapping of groundwater flow fields with internal boundaries

An analytical approach is presented for solving problems of steady, two-dimensional groundwater flow with inhomogeneity boundaries. A common approach for such problems is to separate the problem domain into two homogeneous domains, search for solutions in each domain, and then attempt to match conditions, either exactly or approximately, along the inhomogeneity boundary. Here, we use classical solutions to problems with inhomogeneity boundaries with simple geometries, and map conformally the entire domain onto a new one. In this way, existing solutions are used to solve problems with more complex, and more practical, boundary geometries. The approach is general, but subject to some restrictions on the mapping functions that may be used.Using this approach, we develop explicit analytical solutions for two problems of practical interest. The first problem addresses aquifer interaction across a gap in an impermeable separating layer; flow regimes are defined and the interaction is quantified. The second solution represents flow in the vertical plane to a partially clogged stream bed that is partially penetrating the aquifer; the stream bed is modeled as a thin layer of low-permeability silt. Flow regimes for groundwater surface–water interaction are quantified analytically.     

Analytical solutions for 2D topography-driven flow in stratified and syncline aquifers

Two analytical solution methods are presented for regional steady-state groundwater flow in a two-dimensional stratified aquifer cross section where the water table is approximated by the topographic surface. For the first solution, the surficial aquifer is represented as a set of dipping parallel layers with different, but piecewise constant, anisotropic hydraulic conductivities, where the anisotropy is aligned with the dip of the layered formation. The model may be viewed as a generalization of the solutions developed by [Tóth JA. A theoretical analysis of groundwater flows in small drainage basins. J Geophys Res 1963;68(16):4795–812; Freeze R, Witherspoon P. Theoretical analysis of regional groundwater flow 1) analytical and numerical solution to the mathematical model, water resources research. Water Resour Res 1966;2(4):641–56; Selim HM. Water flow through multilayered stratified hillside. Water Resour Res 1975;11:949–57] to an multi-layer aquifer with general anisotropy, layer orientation, and a topographic surface that may intersect multiple layers. The second solution presumes curved (syncline) layer stratification with layer-dependent anisotropy aligned with the polar coordinate system. Both solutions are exact everywhere in the domain except at the topographic surface, where a Dirichlet condition is met in a least-squared sense at a set of control points; the governing equation and no-flow/continuity conditions are met exactly. The solutions are derived and demonstrated on multiple test cases. The error incurred at the location where the layer boundaries intersect the surface is assessed.     

Modelling of the groundwater flow and of tracer movement in the porous and fissured media: Chalk Aquifer (Northern part of Paris Basin,France)

A groundwater flow model has been developed in order to study the chalk aquifer of Paris Basin, based on most of the geological and hydrological available data. The numerical processes are intended to modelling the groundwater flow in the Senonian (Late Cretaceous) formations and to visualize the tracer movement in groundwater resources in the experimental site of LaSalle Beauvais (northern part Paris Basin). Both objectives were achieved as follows: (i) the comprehension of the spatial distribution of the hydraulic conductivity in the chalk aquifer taking into account the characteristics of the hydrogeological system and (ii) the use of the analytical solution for describing one‐dimensional to two‐dimensional solute transport in a unidirectional steady‐state flow tracer with scale‐dependent dispersion. Advection and diffusion mechanisms are taken into account. Comparison between the breakthrough curves of the analytical and the numerical solutions provided an excellent agreement for various ranges of scale‐related transport parameters of interest. The developed power series solution facilitates fast prediction of the breakthrough curves at each observation point. Thus, the derived new solutions are widely applicable and are very useful for the validation of numerical transport. The numerical approach is carried out by MT3DMS, a Modular 3‐D Multi‐Species Transport Model for Simulation of Advection, Dispersion, and Chemical Reactions of Contaminants in Groundwater Systems, and based on total variation‐diminishing method using the ULTIMATE algorithm. The estimation of the infected surface could constitute an approach in water management and allows to prevent the risks of pollution and to manage the groundwater resource from a durable development perspective. Copyright © 2016 John Wiley & Sons, Ltd.     

Closed-form analytical solutions incorporating pumping and tidal effects in various coastal aquifer systems

Pumping wells are common in coastal aquifers affected by tides. Here we present analytical solutions of groundwater table or head variations during a constant rate pumping from a single, fully-penetrating well in coastal aquifer systems comprising an unconfined aquifer, a confined aquifer and semi-permeable layer between them. The unconfined aquifer terminates at the coastline (or river bank) and the other two layers extend under tidal water (sea or tidal river) for a certain distance L. Analytical solutions are derived for 11 reasonable combinations of different situations of the L-value (zero, finite, and infinite), of the middle layer’s permeability (semi-permeable and impermeable), of the boundary condition at the aquifer’s submarine terminal (Dirichlet describing direct connection with seawater and no-flow describing the existence of an impermeable capping), and of the tidal water body (sea and tidal river). Solutions are discussed with application examples in fitting field observations and parameter estimations.