Pollutant’s Horizontal Dispersion Along and Against Sinusoidally Varying Velocity from a Pulse Type Point Source

An analytical solution of a two-dimensional advection diffusion equation with time dependent coefficients is obtained by using Laplace Integral Transformation Technique. The horizontal medium of solute transport is considered of semi-infinite extent along both the longitudinal and lateral directions. The input concentration is assumed at an intermediate position of the domain. It helps to evaluate concentration level along the flow as well as against the flow through one model only. The source of the input concentration is considered to be of pulse type. In the presence of the source, it is assumed to be decreasing very slowly with time, and just after the elimination of the source it is assumed to be zero. The dispersion coefficient and the advection parameter are considered directly proportional to each other. The analytical solution may be used to predict the solute concentration level with position and time in an open medium as well as in a porous medium. The effect of heterogeneity on the solute transport may also be predicted.     

Two‐dimensional solute transport for periodic flow in isotropic porous media: an analytical solution

In this article, a mathematical model is presented for the dispersion problem in finite porous media in which the flow is two‐dimensional, the seepage flow velocity is periodic, and dispersion parameter is proportional to the flow velocity. In addition to these, first‐order decay and zero‐order production parameters have also been considered directly proportional to the velocity. Retardation factor is taken into account in the present problem. First‐type boundary condition of periodic nature is considered at the extreme end of the boundary. Mixed‐type boundary condition is assumed at the origin of the domain. A classical mathematical substitution transforms the original advection–dispersion equation into diffusion equation in terms of other dependent and independent variables, with constant coefficients. Laplace transform technique is used to obtain the analytical solution. Copyright © 2011 John Wiley & Sons, Ltd.     

Two-dimensional concentration distribution for mixing-controlled bioreactive transport in steady state

Under steady-state conditions, the degradation of contaminant plumes introduced continuously into an aquifer is controlled by transverse dispersion when the other reacting compound is provided from ambient groundwater. Given that the reaction is instantaneous and longitudinal dispersion can be neglected, the length of the plume is inversely proportional to the transverse dispersion coefficient. In typical scenarios of natural attenuation, however, the considered reaction is biotic and kinetic. The standard model of bioreactive transport relies on double-Monod kinetics and pseudo first-order biomass decay. Under these conditions, a fraction of the injected mass flux remains beyond the length of the plume determined for the instantaneous reaction. We present an analytical framework to derive the steady-state concentration distributions of the dissolved compounds and the biomass from the concentration distribution of a conservative compound, assuming double-Monod kinetics and two different models describing biomass decay. The first biomass-decay model assumes a constant first-order decay coefficient, while the second assumes that the decay coefficient depends upon the electron-acceptor concentration. We apply the method to the case of a line-injection in two-dimensional uniform flow. In general, the bioreactive concentration distributions are similar to the distributions computed for an instantaneous reaction. The similarity is greater when the biomass decay coefficient is assumed to depend on the electron-acceptor concentration rather than being constant.     

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.     

Upstream Dispersion in Solute Transport Models: A Simple Evaluation and Reduction Methodology

This article outlines analytical solutions to quantify the length scale associated with “upstream dispersion,” the artificial movement of solutes in the opposite direction to groundwater flow, in solute transport models. Upstream dispersion is an unwanted artifact in common applications of the advection-dispersion equation (ADE) in problems involving groundwater flow in the direction of increasing solute concentrations. Simple formulae for estimating the one-dimensional distance of upstream dispersion are provided. These show that under idealized conditions (i.e., steady-state flow and transport, and a homogeneous aquifer), upstream dispersion may be a function of only longitudinal dispersivity. The scale of upstream dispersion in a selection of previously presented situations is approximated to highlight the utility of the presented formulae and the relevance of this ADE anomaly in common transport problems. Additionally, the analytical solution is applied in a hypothetical scenario to guide the modification of dispersion parameters to minimize upstream dispersion.     

Analytical solution for drainage and recession from an unconfined aquifer

One-dimensional transient groundwater flow from a divide to a river in an unconfined aquifer described by the Boussinesq equation was studied. We derived the analytical solution for the water table recession and drainage change process described with a linearized Boussinesq equation with a physically based initial condition. A method for determining the average water table in the solutions was proposed. It is shown that the solution derived in the form of infinite series can be well approximated with the simplified solution which contains only the leading term of the original solution. The solution and their simplification can be easily evaluated and used by others to study the groundwater flow problems, such as drainage and base flow estimation, in an unconfined aquifer.     

One-dimensional solute dispersion along unsteady flow through a heterogeneous medium,dispersion being proportional to the square of velocity

Abstract

One-dimensional solute transport, originating from a continuous uniform point source, is studied along unsteady longitudinal flow through a heterogeneous medium of semi-infinite extent. Velocity is considered as directly proportional to the linear spatially-dependent function that defines the heterogeneity. It is also assumed temporally dependent. It is expressed in both the independent variables in degenerate form. The dispersion parameter is considered to be proportional to square of the velocity. Certain new independent variables are introduced through separate transformations to reduce the variable coefficients of the advection–diffusion equation to constant coefficients. The Laplace Transformation Technique (LTT) is used to obtain the desired solution. The effects of heterogeneity and unsteadiness on the solute transport are investigated.

Editor D. Koutsoyiannis; Associate editor F.F. Hattermann

Citation Kumar, A., Jaiswal, D.K., and Kumar, N., 2012. One-dimensional solute dispersion along unsteady flow through a heterogeneous medium, dispersion being proportional to the square of velocity. Hydrological Sciences Journal, 57 (6), 1223–1230.     

Vectorized simulation of groundwater flow and contaminant transport using analytic element method and random walk particle tracking

Groundwater contaminant transport processes are usually simulated by the finite difference (FDM) or finite element methods (FEM). However, they are susceptible to numerical dispersion for advection‐dominated transport. In this study, a numerical dispersion‐free coupled flow and transport model is developed by combining the analytic element method (AEM) with random walk particle tracking (RWPT). As AEM produces continuous velocity distribution over the entire aquifer domain, it is more suitable for RWPT than FDM/finite element methods. Using the AEM solutions, RWPT tracks all the particles in a vectorized manner, thereby improving the computational efficiency. The present model performs a convolution integral of the response of an impulse contaminant injection to generate concentration distributions due to a permanent contaminant source. The RWPT model is validated with an available analytical solution and compared to an FDM solution, the RWPT model more accurately replicates the analytical solution. Further, the coupled AEM‐RWPT model has been applied to simulate the flow and transport in hypothetical and field aquifer problems. The results are compared with the FDM solutions and found to be satisfactory. The results demonstrate the efficacy of the proposed method.     

THE USE OF VARIATIONAL METHODS AND OF ERROR DISTRIBUTION PRINCIPLES IN GROUNDWATER HYDRAULICS

Abstract

The aim of the present paper is to present some mathematical techniques for the solution of problems connected with three-dimensional steady-state groundwater flow with a free surface. The validity of Darcy's law is assumed. As no use is made of the Dupuit-Forschheimer approximation, the shape of the free surface and the velocity potential must be determined simultaneously from a non-linear boundary value problem. In order to demonstrate the use of a variational method and of error distribution principles for the solution of those problems by an example as simple as possible, we investigate the gravity flow of incompressible, homogeneous groundwater towards a circular well completely penetrating an isotropic, homogeneous, inelastic aquifer resting on a horizontal, impermeable substratum.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.     

A general solution for groundwater flow in an estuarine's leaky aquifer system after considering aquifer anisotropy

This paper presents an analytical model for describing the tidal effects in a two‐dimensional leaky confined aquifer system in an estuarine delta where ocean and river meet. This system has an unconfined aquifer on top and a confined aquifer on the bottom with an aquitard in between the two. The unconfined and confined aquifers interact with each other through leakage. It was assumed that the aquitard storage was negligible and that the leakage was linearly proportional to the head difference between the unconfined and confined aquifers. This model's solution was based on the separation of variables method. Two existing solutions that deal with the head fluctuation in one‐dimensional or two‐dimensional leaky confined aquifers are shown as special cases in the present solution. Based on this new solution, the dynamic effect of the water table's fluctuations can be clearly explored, as well as the influence of leakage on the behaviour of fluctuations in groundwater levels in the leaky aquifer system. Copyright © 2015 John Wiley & Sons, Ltd.     

Monte Carlo evaluation of microbial-mediated contaminant reactions in heterogeneous aquifers

Monte Carlo simulations are conducted to evaluate microbial-mediated contaminant reactions in an aquifer comprised of spatially variable microbial biomass concentrations, aquifer hydraulic conductivities, and initial electron donor/acceptor concentrations. A finite element simulation model is used that incorporates advection, dispersion, and Monod kinetic expressions to describe biological processes. Comparisons between Monte Carlo simulations of heterogeneous systems and simulations using homogeneous formulation of the same two-dimensional transport problem are presented. For the assumed set of parameters, physical aquifer heterogeneity is found to have a minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system; however, it noticeably changes the dispersion, skewness, and peakness of contaminant concentration distributions. Similarly, for low microbial growth rate, given favorable microbial growth characteristics, biological heterogeneity has minor effect on the mass of contaminant biodegraded/transformed when compared to a homogeneous system. On the other hand, when higher effective growth rates are assumed, biological heterogeneity and spatial heterogeneities in essential electron donor/acceptors reduce the efficiency of biotic contaminant reactions; consequently, model simulations derived from heterogeneous biomass distributions predict remediation time scales that are longer than those simulated for homogeneous systems. When correlations between physical aquifer and biological heterogeneities are considered, the assumed correlation affects predicted mean and variance of contaminant concentration and biomass distributions. For example, an assumed negative correlation between hydraulic conductivity and the initial biomass distribution produces a plume where less efficient biotic contaminant reactions occur at the leading edge of the plume; this is consistent with less degradation/transformation occurring over regions of higher groundwater velocities. However, the presence and absence of these correlations do not appear to affect the efficiency of microbial-mediated contaminant attenuation.     

Point dilution tests to calculate groundwater velocity: an example in a porous aquifer in northeast Italy

ABSTRACT

The point dilution test is a single-well technique for estimating horizontal flow velocity in the aquifer surrounding a well. The test is conducted by introducing a tracer into a well section and monitoring its decreasing concentration over time. When using a salt tracer, the method is easy and inexpensive. Traditionally, the horizontal Darcy velocity is calculated as a function of the rate of dilution and is based on the simple assumption that the decreasing tracer concentration is proportional both to the apparent velocity into the test section and to the Darcy velocity in the aquifer. In this article, an alternative approach to analyse the results of point dilution tests is proposed and verified using data acquired at a test site in the middle Venetian plain, northeast Italy. In this approach, the one-dimensional equilibrium advection–dispersion equation is inverted using the CXTFIT model to estimate the apparent velocity inside the test section. Analysis of the field data obtained by the two approaches shows good agreement between the methods and suggests that it is possible to use the equilibrium advection–dispersion equation to estimate apparent velocity over a wide range of velocities.

Editor D. Koutsoyiannis; Associate editor K. Heal     

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.     

Dissolution of nonaqueous phase liquid pools in anisotropic aquifers

A two-dimensional numerical transport model is developed to determine the effect of aquifer anisotropy and heterogeneity on mass transfer from a dense nonaqueous phase liquid (DNAPL) pool. The appropriate steady state groundwater flow equation is solved implicitly whereas the equation describing the transport of a sorbing contaminant in a confined aquifer is solved by the alternating direction implicit method. Statistical anisotropy in the aquifer is introduced by two-dimensional, random log-normal hydraulic conductivity field realizations with different directional correlation lengths. Model simulations indicate that DNAPL pool dissolution is enhanced by increasing the mean log-transformed hydraulic conductivity, groundwater flow velocity, and/or anisotropy ratio. The variance of the log-transformed hydraulic conductivity distribution is shown to be inversely proportional to the average mass transfer coefficient.     

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.     

Tracing groundwater flow in the Borden aquifer using krypton-85

Krypton-85 was measured in air, soil gas, and ground water at the Borden aquifer in Ontario in October 1989. The measured specific activities in air and soil gas were 52.0 ± 2.0 and 53.6 ± 1.8 disintegrations per min (dpm) cm⁻³ krypton. These measurements are in excellent agreement with the global atmospheric trend and demonstrate that krypton-85 enters the water table at the Borden site without a lag in the soil gas reservoir. The krypton-85 specific activity in five groundwater samples ranged from 44.9 to 9.5 dpm cm⁻³ corresponding to groundwater ages of 2–17 years with a monotonic decrease in specific activity (increase in age) along the groundwater flow path. Travel times calculated from a two-dimensional steady-state model of groundwater flow agree well with the krypton-85 ages in the main recharge region of the aquifer where flow is predominantly vertical, but were 30–40% older than the krypton-85 age downstream of the main recharge area where the flow is mainly horizontal. The effect of dispersion on the distribution of krypton-85 was determined by modelling the transport of krypton-85 in the Borden aquifer with a two-dimensional time-dependent advection dispersion model using the steady-state flow field. Agreement between model specific activity and observed specific activity was excellent for samples in the main recharge region, but the model specific activities were 30–50% lower than observed specific activities in the region of horizontal flow. The differences in travel times and krypton-85 ages and in model krypton-85 and observed krypton-85 specific activities are considered to be small given the heterogeneities that exist in the hydraulic conductivity and aquifer geometry and hence in the groundwater flow field. The model simulated krypton-85 distribution was not sensitive to changes in longitudinal dispersivity and was only weakly sensitive to changes in transverse dispersivity. The geochemical inertness, well-defined source function, and insensitivity to dispersion of krypton-85 allow estimates of groundwater age to be made in a straightforward manner and measurement of krypton-85 can significantly enhance the characterization of groundwater flow in many shallow subsurface systems.     

Direct simulation of solute recycling in irrigated areas

Solute recycling from irrigation can be described as the process that occurs when the salt load that is extracted from irrigation wells and distributed on the fields is returned to the groundwater below irrigated surfaces by deep percolation. Unless the salt load leaves the system by means of drains or surface runoff, transfer to the groundwater will take place, sooner or later. This can lead to solute accumulation and thus to groundwater degradation, particularly in areas where extraction rates exceed infiltration rates (semi-arid and arid regions). Thus, considerable errors can occur in a predictive solute mass budget if the recycling process is not accounted for in the calculation. A method is proposed which allows direct simulation of solute recycling. The transient solute response at an extraction well is shown to be a superposition of solute mass flux contributions from n recycling cycles and is described as a function of the travel time distribution between a recycling point and a well. This leads to an expression for a transient ‘recycling source’ term in the advection–dispersion equation, which generates the effect of solute recycling. At long times, the ‘recycling source’ is a function of the local capture probability of the irrigation well and the solute mass flux captured by the well from the boundaries. The predicted concentration distribution at steady state reflects the maximum spatial concentration distribution in response to solute recycling and can thus be considered as the solute recycling potential or vulnerability of the entire domain for a given hydraulic setting and exploitation scheme. Simulation of the solute recycling potential is computationally undemanding and can therefore, for instance, be used for optimisation purposes. Also, the proposed method allows transient simulation of solute recycling with any standard flow and transport code.     

Time–space fractional governing equations of transient groundwater flow in confined aquifers: Numerical investigation

In order to model non‐Fickian transport behaviour in groundwater aquifers, various forms of the time–space fractional advection–dispersion equation have been developed and used by several researchers in the last decade. The solute transport in groundwater aquifers in fractional time–space takes place by means of an underlying groundwater flow field. However, the governing equations for such groundwater flow in fractional time–space are yet to be developed in a comprehensive framework. In this study, a finite difference numerical scheme based on Caputo fractional derivative is proposed to investigate the properties of a newly developed time–space fractional governing equations of transient groundwater flow in confined aquifers in terms of the time–space fractional mass conservation equation and the time–space fractional water flux equation. Here, we apply these time–space fractional governing equations numerically to transient groundwater flow in a confined aquifer for different boundary conditions to explore their behaviour in modelling groundwater flow in fractional time–space. The numerical results demonstrate that the proposed time–space fractional governing equation for groundwater flow in confined aquifers may provide a new perspective on modelling groundwater flow and on interpreting the dynamics of groundwater level fluctuations. Additionally, the numerical results may imply that the newly derived fractional groundwater governing equation may help explain the observed heavy‐tailed solute transport behaviour in groundwater flow by incorporating nonlocal or long‐range dependence of the underlying groundwater flow field.     

A New Capture Fraction Method to Map How Pumpage Affects Surface Water Flow

All groundwater pumped is balanced by removal of water somewhere, initially from storage in the aquifer and later from capture in the form of increase in recharge and decrease in discharge. Capture that results in a loss of water in streams, rivers, and wetlands now is a concern in many parts of the United States. Hydrologists commonly use analytical and numerical approaches to study temporal variations in sources of water to wells for select points of interest. Much can be learned about coupled surface/groundwater systems, however, by looking at the spatial distribution of theoretical capture for select times of interest. Development of maps of capture requires (1) a reasonably well-constructed transient or steady state model of an aquifer with head-dependent flow boundaries representing surface water features or evapotranspiration and (2) an automated procedure to run the model repeatedly and extract results, each time with a well in a different location. This paper presents new methods for simulating and mapping capture using three-dimensional groundwater flow models and presents examples from Arizona, Oregon, and Michigan.