One-dimensional uniform and time varying solute dispersion along transient groundwater flow in a semi-infinite aquifer

An analytical solution for the space-time variation of contaminant concentration in one-dimensional transient groundwater flow in a homogenous semi-infinite aquifer, subjected to time-dependent source contamination, is derived. The uniform and time varying dispersion along transient groundwater flow is investigated under two conditions. First, the flow velocity distribution in the aquifer is considered as a sinusoidally varying function, and second, the flow velocity distribution is treated as an exponentially increasing function of time. It is assumed that initially the aquifer is not solute free, so the initial background concentration is considered as an exponentially decreasing function of the space variable which is tending to zero at infinity. It is assumed that dispersion is directly proportional to the square of the velocity, noting that experimental observations indicate that dispersion is directly proportional to the velocity with a power ranging from 1 to 2. The analytical solution is illustrated using an example and may help benchmark numerical codes and solutions.     

Matrix Diffusion Effects on Contaminant Migration from an Injection Well in Fractured Sandstone

Deep-well injection into fractured sandstone is an option for the disposal of contaminated mine dewatering discharge from an open pit uranium mine. As part of the assessment of potential contaminant migration from deep-well injection, the effect of matrix diffusion was evaluated. An analytical mathematical model was developed for the simulation of the radial movement of a contaminant front away from an injection point under steady flow conditions in a planar fracture with uniform properties. The model includes the effects of advection in the fracture, diffusion of contaminants from the fracture into the rock matrix, and equilibrium adsorption on the fracture surface as well as in the rock matrix. Effective diffusion coefficients obtained from laboratory experiments on 11 intact core samples varied from 3.4 × 10⁻⁸ to 3.2 × 10⁻⁷ cm²/s. Model simulations were made with diffusion coefficient values in this range and with single-fracture injection rates estimated from fracture frequencies in boreholes, and from bulk hydraulic conductivity values obtained from field tests. Because of matrix diffusion, the rate of outward movement of the front of the nonreactive contaminants from the injection well is much slower than the rate of water flow in the fractures. Simulations of the movement of contaminants that undergo adsorption indicate that even a small distribution coefficient for the rock matrix causes the contaminants to remain very close to the injection well during the one-year period. The results of the simplified model demonstrate that matrix diffusion is an important process that cannot be neglected in the assessment of a waste disposal scheme located in fractured porous rock. However, in order to make a definitive assessment of the capability of matrix diffusion and associated matrix adsorption to significantly limit the extent of contaminant migration around injection wells, it would be necessary to conduct field tests such as a preliminary or experimental injection.     

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.     

Numerical modelling of stream–aquifer interaction: Quantifying the impact of transient streambed permeability and aquifer heterogeneity

Stream–aquifer interaction plays a vital role in the water cycle, and a proper study of this interaction is needed for understanding groundwater recharge, contaminants migration, and for managing surface water and groundwater resources. A model‐based investigation of a field experiment in a riparian zone of the Schwarzbach river, a tributary of the Rhine River in Germany, was conducted to understand stream–aquifer interaction under alternative gaining and losing streamflow conditions. An equivalent streambed permeability, estimated by inverting aquifer responses to flood waves, shows that streambed permeability increased during infiltration of stream water to aquifer and decreased during exfiltration. Aquifer permeability realizations generated by multiple‐point geostatistics exhibit a high degree of heterogeneity and anisotropy. A coupled surface water groundwater flow model was developed incorporating the time‐varying streambed permeability and heterogeneous aquifer permeability realizations. The model was able to reproduce varying pressure heads at two observation wells near the stream over a period of 55 days. A Monte Carlo analysis was also carried out to simulate groundwater flow, its age distribution, and the release of a hypothetical wastewater plume into the aquifer from the stream. Results of this uncertainty analysis suggest (a) stream–aquifer exchange flux during the infiltration periods was constrained by aquifer permeability; (b) during exfiltration, this flux was constrained by the reduced streambed permeability; (c) the effect of temporally variable streambed permeability and aquifer heterogeneity were found important to improve the accurate capture of the uncertainty; and (d) probabilistic infiltration paths in the aquifer reveal that such pathways and the associated prediction of the extent of the contaminant plume are highly dependent on aquifer heterogeneity.     

Groundwater monitoring strategies for variable contaminant migration boundaries

Abstract

Detection efficiencies of alternative groundwater monitoring networks were evaluated in relation to distance to a buffer zone (contaminant migration) boundary. This boundary establishes a distance limit within which contaminant plumes should pass through monitoring wells, located on curvilinear segments (monitoring loci) near a waste storage facility. Alternative strategies allocated monitoring wells to loci at specified distances, measured parallel to groundwater flow, from the downgradient boundaries of a landfill. One approach constrained wells to equal spacing, measured perpendicular to groundwater flow. Compressing well locations 10% closer to the downgradient corner of the landfill rendered alternative monitoring configurations. Computations by a monitoring efficiency model indicated: (a) networks largely maintained detection efficiency for different contaminant migration boundaries; (b) one network most efficiently attained a target detection capability for all contaminant migration boundaries; and (c) compressed networks slightly outperformed equal-spaced counterparts. Compressed networks with more wells along closer monitoring loci best maintained the detection efficiency when shifting the contaminant migration boundary closer to the landfill. Procedures described in this paper may be useful for examining trade-offs between monitoring efficiency and distance limits of contaminant travel at landfills posing potential hazards to underlying groundwater.     

A Remedial Investigation of an Organically Polluted Outwash Aquifer

A glacial outwash aquifer underlying the Gloucester Landfill near Ottawa, Canada, has become polluted with various organic chemicals following the disposal of laboratory solvents in shallow trenches immediately above the aquifer. Several remedial alternatives have been considered by the government of Canada. Impermeable barrier walls were rejected as being unsuitable given the permeable nature of the underlying bedrock. It appears improbable that pools of liquid organic chemicals (DNAPLs) exist within the aquifer, although ganglia are likely present. Therefore, much of the contaminant plume can be removed hydraulically over a period of five years by the operation of four purge wells discharging to an on-site treatment plant from which the purified water is returned to the aquifer by recharge wells. The residual contamination is anticipated to be cleaned up by in situ biorestoration techniques currently under development.     

Multiple Working Transport Hypotheses in a Heterogeneous Glacial Aquifer System

Multiple working hypotheses can be used to evaluate permissible alternative hydrogeological interpretations at sites with limited subsurface control. This approach was applied to test the viability of three conceptual aquifer system architecture models coupled with three hypothesized source locations for a 1,4-dioxane plume in a heterogeneous glacial aquifer system in Washtenaw County, Michigan. The three alternative conceptual models characterized the site hydrogeology with increasingly complex distributions of hydrostratigraphic units: (A) an effective aquifer, (B) a layered confined aquifer, and (C) a discretely heterogeneous aquifer model. Each was incorporated into an independently calibrated numerical ground water flow (MODFLOW) model. Steady-state and transient flow simulations of the alternative models were evaluated using both hydraulic flow field characteristics observed under natural conditions and the perturbed response after local remedial pumping activity began. Three plausible locations where 1,4-dioxane could have entered the aquifer system were identified using historical information at the site: (1) manufacturing waste water disposal lagoons, (2) a 60 foot (18 m) deep kettle lake, and (3) a shallow impoundment on a local stream. Advective transport modeling (MODPATH) was used to assess the consistency of the hypothesized source locations with observed contaminant migration pathways inferred from the mapped location of the plume. Evaluation of the nine combinations of hydrogeologic conceptualizations and 1,4-dioxane source locations led to elimination of four working hypotheses and discounting of two others, leading to reduced overall uncertainty and the development of new insights into the system behavior.     

Analytical models of contaminant transport in coastal aquifers

The Henry formulation, which couples subsurface flow and salt transport via a variable-density flow formulation, can be used to evaluate the extent of sea water intrusion into coastal aquifers. The coupling gives rise to nontrivial flow patterns that are very different from those observed in inland aquifers. We investigate the influence of these flow patterns on the transport of conservative contaminants in a coastal aquifer. The flow is characterized by two dimensionless parameters: the Péclet number, which compares the relative effects of advective and dispersive transport mechanisms, and a coupling parameter, which describes the importance of the salt water boundary on the flow. We focus our attention on two regimes – low and intermediate Péclet number flows. Two transport scenarios are solved analytically by means of a perturbation analysis. The first, a natural attenuation scenario, describes the flushing of a contaminant from a coastal aquifer by clean fresh water, while the second, a contaminant spill scenario, considers an isolated point source.     

Simulation of liquid waste buoyancy in a deep heterogeneous aquifer

The results of 3D numerical simulation of coupled flow and transport processes of light-waste migration in heterogeneous aquifer with brine are described. The study is based on the geological data and monitoring of deep-well injection of liquid waste at Chepetsk Mechanical Plant (ChMZ, Glazov). The buoyancy of light waste is simulated for the post-injection period of 300 years using different conceptual models of heterogeneity. The procedure of 3D geostatistical simulation based on transition probability analysis with TSIM code of the reservoir heterogeneity is presented. The results of the ascending waste migration simulation using SEAWAT2000 code indicate that both the variable fluid density and the hydraulic heterogeneity must be considered in the forecast models for determining the possible waste distribution during the injection and the post injection periods.     

Effects of shield brine on the safe disposal of waste in deep geologic environments

The salinity of groundwater increases with depth in the Canadian Shield (up to 1.3 kg/L of density). The existence of brine can be critically important for the safe geologic disposal of radioactive wastes, as dense brine can significantly retard the upward migration of radionuclides released from repositories. Static and flushing conditions of the deep brine are analyzed using a U-tube analogy model. Velocity reduction due to the presence of dense brine is derived under flushing conditions. A set of illustrative numerical simulations in a two-dimensional cross section is presented to demonstrate that dense brine can significantly influence regional groundwater flow patterns in a shield environment. It is implied from the results that (1) the existence of Shield brine can be an indicator of a hydrogeologically stable environment, (2) activities near ground surface may not perturb the stable groundwater environment in the deep brine region, and thus, (3) the deep brine region can be considered as a candidate geologic site for the safe disposal of waste. In addition to brine, other issues associated with long-term waste disposal, such as geological, glacial and seismic events, may need to be considered for the safe storage of spent nuclear fuel in a shield environment.     

Assessing the Transport of Pharmaceutical Compounds in a Layered Aquifer Discharging to a Stream

A groundwater plume containing high concentrations of pharmaceutical compounds, mainly sulfonamides, barbiturates, and ethyl urethane, in addition to chlorinated ethenes and benzene was investigated. The contamination originating from a former pharmaceutical industry discharges into a multilayered aquifer system and a downgradient stream. In this study, geological and hydrogeological data were integrated into a numerical flow model to examine identified trends using statistical approaches, including principal component analysis and hierarchal cluster analysis. A joint interpretation of the groundwater flow paths and contaminant concentrations in the different compartments (i.e., groundwater and hyporheic zone) provided insight on the transport processes of the different contaminant plumes to the stream. The analysis of historical groundwater concentrations of pharmaceutical compounds at the site suggested these compounds are slowly degrading. The pharmaceutical compounds migrate in both a deep semiconfined aquifer, as well as in the shallow unconfined aquifer, and enter the stream along a 2-km stretch. This contrasted with the chlorinated ethenes, which mainly discharge to the stream as a focused plume from the unconfined aquifer. The integrated approach developed here, combining groundwater flow modeling and statistical analyses of the contaminant concentration data collected in groundwater and the hyporheic zone, lead to an improved understanding of the observed distribution of contaminants in the unconfined and semiconfined aquifers, and thus to their discharge to the stream. This approach is particularly relevant for large and long-lasting contaminant sources and plumes, such as abandoned landfills and industrial production sites, where field investigations may be very expensive.     

Vegetation-groundwater dynamics at a former uranium mill site following invasion of a biocontrol agent: A time series analysis of Landsat normalized difference vegetation index data

Because groundwater recharge in dry regions is generally low, arid and semiarid environments have been considered well-suited for long-term isolation of hazardous materials (e.g., radioactive waste). In these dry regions, water lost (transpired) by plants and evaporated from the soil surface, collectively termed evapotranspiration (ET), is usually the primary discharge component in the water balance. Therefore, vegetation can potentially affect groundwater flow and contaminant transport at waste disposal sites. We studied vegetation health and ET dynamics at a Uranium Mill Tailings Radiation Control Act (UMTRCA) disposal site in Shiprock, New Mexico, where a floodplain alluvial aquifer was contaminated by mill effluent. Vegetation on the floodplain was predominantly deep-rooted, non-native tamarisk shrubs (Tamarix sp.). After the introduction of the tamarisk beetle (Diorhabda sp.) as a biocontrol agent, the health of the invasive tamarisk on the Shiprock floodplain declined. We used Landsat normalized difference vegetation index (NDVI) data to measure greenness and a remote sensing algorithm to estimate landscape-scale ET along the floodplain of the UMTRCA site in Shiprock prior to (2000–2009) and after (2010–2018) beetle establishment. Using groundwater level data collected from 2011 to 2014, we also assessed the role of ET in explaining seasonal variations in depth to water of the floodplain. Growing season scaled NDVI decreased 30% (p < .001), while ET decreased 26% from the pre- to post-beetle period and seasonal ET estimates were significantly correlated with groundwater levels from 2011 to 2014 (r² = .71; p = .009). Tamarisk greenness (a proxy for health) was significantly affected by Diorhabda but has partially recovered since 2012. Despite this, increased ET demand in the summer/fall period might reduce contaminant transport to the San Juan River during this period.     

Bioremediation in Fractured Rock: 1. Modeling to Inform Design,Monitoring, and Expectations

Field characterization of a trichloroethene (TCE) source area in fractured mudstones produced a detailed understanding of the geology, contaminant distribution in fractures and the rock matrix, and hydraulic and transport properties. Groundwater flow and chemical transport modeling that synthesized the field characterization information proved critical for designing bioremediation of the source area. The planned bioremediation involved injecting emulsified vegetable oil and bacteria to enhance the naturally occurring biodegradation of TCE. The flow and transport modeling showed that injection will spread amendments widely over a zone of lower‐permeability fractures, with long residence times expected because of small velocities after injection and sorption of emulsified vegetable oil onto solids. Amendments transported out of this zone will be diluted by groundwater flux from other areas, limiting bioremediation effectiveness downgradient. At nearby pumping wells, further dilution is expected to make bioremediation effects undetectable in the pumped water. The results emphasize that in fracture‐dominated flow regimes, the extent of injected amendments cannot be conceptualized using simple homogeneous models of groundwater flow commonly adopted to design injections in unconsolidated porous media (e.g., radial diverging or dipole flow regimes). Instead, it is important to synthesize site characterization information using a groundwater flow model that includes discrete features representing high‐ and low‐permeability fractures. This type of model accounts for the highly heterogeneous hydraulic conductivity and groundwater fluxes in fractured‐rock aquifers, and facilitates designing injection strategies that target specific volumes of the aquifer and maximize the distribution of amendments over these volumes.     

Effect of variable‐density groundwater flow on nitrate flux to coastal waters

A two‐dimensional variable‐density groundwater flow and transport model was developed to provide a conceptual understanding of past and future conditions of nitrate (NO₃) transport and estimate groundwater nitrate flux to the Gulf of Mexico. Simulation results show that contaminant discharge to the coast decreases as the extent of saltwater intrusion increases. Other natural and/or artificial surface waters such as navigation channels may serve as major sinks for contaminant loading and act to alter expected transport pathways discharging contaminants to other areas. Concentrations of NO₃ in the saturated zone were estimated to range between 30 and 160 mg?L^?1 as NO₃. Relatively high hydraulic vertical gradients and mixing likely play a significant role in the transport processes, enhancing dilution and contaminant migration to depth. Residence times of NO₃ in the deeper aquifers vary from 100 (locally) to about 300 years through the investigated aquifer system. NO₃ mass fluxes from the shallow aquifers (0 to 5.7 × 10⁴ mg?m^?2?day^?1) were primarily directed towards the navigation channel, which intersects and captures a portion of the shallow groundwater flow/discharge. Direct NO₃ discharge to the sea (i.e. Gulf of Mexico) from the shallow aquifer was very low (0 to 9.0 × 10¹ mg · m^?2?day^?1) compared with discharge from the deeper aquifer system (0 to 8.2 × 10³ mg?m^?2?day^?1). Both model‐calibrated and radiocarbon tracer‐determined contaminant flux estimates reveal similar discharge trends, validating the use of the model for density‐dependent flow conditions. The modelling approach shows promise to evaluate contaminant and nutrient loading for similar coastal regions worldwide. Copyright © 2015 John Wiley & Sons, Ltd.     

An Idealized Ground-Water Flow and Chemical Transport Model (S-PATHS)

The number of studies on the actual and potential environmental consequences of contaminated ground water is growing. One means of studying these consequences is through an idealized flow and transport model, S-PATHS, which allows the hydrologist to determine the salient features of contaminant migration with a minimum of data. The transport of contaminants by ground water from many waste disposal sites can be geometrically idealized as flow between a line and a circle. The flow system adjacent to the disposal site can be represented as a contaminant line source, and a downgradient pumping well as a circular sink. To study waste disposal sites on a larger scale the model geometry is reversed and the disposal site is represented as a circular source, and a river or other convenient line of evaluation is represented as a line sink. This idealization allows S-PATHS to describe the flow and transport process directly by a single partial differential expression. S-PATHS considers transmissivity, effective porosity, sorption, source strength, source concentration, decay, potentiometric gradient, circle size, and distance to the line. Coding for the model is not lengthy and can be run on a large-capacity, hand-held calculator.     

Groundwater Flow Field Distortion by Monitoring Wells and Passive Flux Meters

Due to differences in hydraulic conductivity and effects of well construction geometry, groundwater lateral flow through a monitoring well typically differs from groundwater flow in the surrounding aquifer. These differences must be well understood in order to apply passive measuring techniques, such as passive flux meters (PFMs) used for the measurement of groundwater and contaminant mass fluxes. To understand these differences, lab flow tank experiments were performed to evaluate the influences of the well screen, the surrounding filter pack and the presence of a PFM on the natural groundwater flux through a monitoring well. The results were compared with analytical calculations of flow field distortion based on the potential theory of Drost et al. (1968). Measured well flow field distortion factors were found to be lower than calculated flow field distortion factors, while measured PFM flow field distortion factors were comparable to the calculated ones. However, this latter is not the case for all conditions. The slotted geometry of the well screen seems to make a correct analytical calculation challenging for conditions where flow field deviation occurs, because the potential theory assumes a uniform flow field. Finally, plots of the functional relationships of the distortion of the flow field with the hydraulic conductivities of the filter screen, surrounding filter pack and corresponding radii make it possible to design well construction to optimally function during PFM applications.     

The effects of the 2004 tsunami on a coastal aquifer in Sri Lanka

On December 26, 2004, the earthquake off the southern coast of Sumatra in the Indian Ocean generated far-reaching tsunami waves, resulting in severe disruption of the coastal aquifers in many countries of the region. The objective of this study was to examine the impact of the tsunami on groundwater in coastal areas. Field investigations on the east coast of Sri Lanka were carried out along a transect located perpendicular to the coastline on a 2.4 km wide sand stretch bounded by the sea and a lagoon. Measurements of groundwater table elevation and electrical conductivity (EC) of the groundwater were carried out monthly from October 2005 to August 2007. The aquifer system and tsunami saltwater intrusion were modeled using the variable-density flow and solute transport code HST3D to understand the tsunami plume behavior and estimate the aquifer recovery time. EC values reduced as a result of the monsoonal rainfall following the tsunami with a decline in reduction rate during the dry season. The upper part of the saturated zone (down to 2.5 m) returned to freshwater conditions (EC < 1000 μS/cm) 1 to 1.5 years after the tsunami, according to field observations. On the basis of model simulations, it may take more than 15 years for the entire aquifer (down to 28 m) to recover completely, although the top 6 m of the aquifer may become fresh in about 5 years.     

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 field study of the coupled effects of aquifer stratification,fluid density,and groundwater fluctuations on dispersivity assessments

A number of experimental studies have tackled the issue of solute transport parameter assessments either in the laboratory or in the field. But yet, the behavior of a plume in the field under density driven forces, is not well known due to possible development of instabilities. Some field tracer tests on the fate of plumes denser than native groundwater such as those encountered under waste disposal facilities, have pointed out the processes of sinking and splitting at the early stage of migration. The process of dispersion was widely investigated, but the range of dispersivity values obtained from either experimental tests, or numerical and theoretical calculations is still very large, even for the same type of aquifers. These discrepancies were considered to be essentially caused by soil heterogeneities and scale effects. In the meantime, studies on the influence of sinking and fingering have remained more scarce. The objective of the work is to analyze how transport parameters such as dispersivities can be affected by unstable conditions, which lead to plume sinking and fingering. A series of tracer tests were carried out to study under natural conditions, the transport of a dense chloride solution injected in a shallow two-layered aquifer. Two types of experiments were performed: in the first type, source injection was such that the plume could travel downward from one layer to the other of higher pore velocity, and in the second one, the migration took place only in the faster layer. The results suggest some new insights in the processes occurring at the early stages of a dense plume migration moving in a stratified aquifer under groundwater fluctuations, which can be summarized through the following points: (i) Above a stability criterion threshold, a fingering process and a multi modal plume transport take place, but local dispersivities can be cautiously derived, using breakthrough curves matching. (ii) When water table is subject to some cycling or rising, the plume can be significantly distorted in the transverse direction, leading to unusual values of the ratio between longitudinal and transverse dispersivities. (iii) Under stable conditions, for example in the case of straightforward injection in the faster aquifer layer, longitudinal dispersivity is greater than the transverse component as usually encountered, and the obtained transport parameters are closed to macro dispersivity values, which reach their asymptotic limit at very short distances. (iv) The classical scale effect about the varying dispersivity at short distances could be a process mainly due to the distance required for a plume stabilization.     

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.