Natural Gradient Tracer Test to Evaluate Natural Attenuation of MTBE Under Anaerobic Conditions

A natural gradient tracer test using perdeuterated MTBE was conducted in an anaerobic aquifer to determine the relative importance of dispersion and degradation in reducing MTBE concentrations in ground water. Preliminary ground water chemistry and hydraulic conductivity data were used to place the tracer within an existing dissolved MTBE plume at Port Hueneme, California. Following one year of transport, the tracer plume was characterized in detail.
Longitudinal dispersion was identified as the dominant mechanism for lowering the perdeuterated MTBE concentrations. The method of moments was used to determine the longitudinal and lateral dispersion coefficients (0.85 m²/day and 0.08 m²/day, respectively). A mass-balance analysis, carried out after one year of transport, accounted for 110% of the injected mass and indicated that no significant mass loss occurred. The plume structure created by zones of higher and lower hydraulic conductivity at the site was complex, consisting of several localized areas of high tracer concentration in a lower concentration plume. This is important because the aquifer has generally been characterized as exhibiting fairly minor heterogeneity. In addition, the tracer plume followed a curved flowpath that deviated from the more macroscopic direction of ground water flow inferred from local ground water elevation measurements and the behavior of the existing plume. Understanding the mass balance, plume structure, curvature of the tracer plume, and consequently natural attenuation behavior required the detailed sampling approach employed in this study. These data imply that a detailed understanding of site hydrogeology and an extensive sampling network may be critical for the correct interpretation of monitored natural attenuation of MTBE.     

Tide‐induced groundwater head fluctuation in a coastal aquifer system with a submarine outcrop covered by a thin silt layer

We present an analytical solution of groundwater head response to tidal fluctuation in a coastal multilayered aquifer system consisting of an unconfined aquifer, a leaky confined aquifer and a semi‐permeable layer between them. The submarine outcrop of the confined aquifer is covered by a thin silt layer. A mathematical model and the analytical solution of this model are given. The silt layer reduces the amplitude of the hydraulic head fluctuation by a constant factor, and shifts the phase by a positive constant (time lag), both of which depend on the leakances of the silt layer and the semi‐permeable layer. The time lag is less than 1·5 h and 3·0 h for semi‐diurnal and diurnal sea tides respectively. When the leakance of the semi‐permeable layer or the silt layer assumes certain special values, the solution becomes the existing solutions derived by previous researchers. The amplitude of the hydraulic head fluctuation in the confined aquifer increases with the leakance of the silt layer and decreases with the leakance of the semi‐permeable layer, whereas the phase shift of the fluctuation decreases with both of them. A hypothetical example shows that neglecting the silt layer may result in significant parameter estimation discrepancy between the amplitude attenuation and the time‐lag fittings. Copyright © 2007 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.     

Atoll island hydrogeology: flow and freshwater occurrence in a tidally dominated system

A layered-aquifer model of groundwater occurrence in an atoll island was tested with a solute-transport numerical model. The computer model used, SUTRA, incorporates density-dependent flow. This can be significant in freshwater-saltwater interactions associated with the freshwater lens of an atoll island. Boundary conditions for the model included ocean and lagoon tidal variations. The model was calibrated to field data from Enjebi Island, Enewetak Atoll, and tested for sensitivity to a variety of parameters. This resulted in a hydraulic conductivity of 10 m day⁻¹ for the surficial aquifer and 1000 m day⁻¹ for the deeper aquifer; this combination of values gave an excellent reproduction of the tidal response data from test wells. The average salinity distribution was closely reproduced using a dispersivity of 0.02m. The computer simulation quantitatively supports the layered-aquifer model, including under conditions of density-dependent flow, and shows that tidal variations are the predominant driving force for flow beneath the island. The oscillating, vertical flow produced by the tidal variations creates an extensive mixing zone of brackish water. The layered-aquifer model with tidally driven flow is a significant improvement over the Ghyben-Herzberg-Dupuit model as it is conventionally applied to groundwater studies for many Pacific reef islands.     

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.     

Eolian Transport of Geogenic Hexavalent Chromium to Ground Water

A conceptual model of eolian transport is proposed to address the widely distributed, high concentrations of hexavalent chromium (Cr⁺⁶) observed in ground water in the Emirate of Abu Dhabi, United Arab Emirates. Concentrations (30 to more than 1000 μg/L Cr⁺⁶) extend over thousands of square kilometers of ground water systems. It is hypothesized that the Cr is derived from weathering of chromium-rich pyroxenes and olivines present in ophiolite sequence of the adjacent Oman (Hajar) Mountains. Cr⁺³ in the minerals is oxidized to Cr⁺⁶ by reduction of manganese and is subsequently sorbed on iron and manganese oxide coatings of particles. When the surfaces of these particles are abraded in this arid environment, they release fine, micrometer-sized, coated particles that are easily transported over large distances by wind and subsequently deposited on the surface. During ground water recharge events, the readily soluble Cr⁺⁶ is mobilized by rain water and transported by advective flow into the underlying aquifer. Chromium analyses of ground water, rain, dust, and surface (soil) deposits are consistent with this model, as are electron probe analyses of clasts derived from the eroding Oman ophiolite sequence. Ground water recharge flux is proposed to exercise some control over Cr⁺⁶ concentration in the aquifer.     

Periodic forcing in composite aquifers

Observations of periodic components of measured heads have long been used to estimate aquifer diffusivities. The estimations are often made using well-known solutions of linear differential equations for the propagation of sinusoidal boundary fluctuations through homogeneous one-dimensional aquifers. Recent field data has indicated several instances where the homogeneous aquifer solutions give inconsistent estimates of aquifer diffusivity from measurements of tidal lag and attenuation. This paper presents new algebraic solutions for tidal propagation in spatially heterogeneous one-dimensional aquifers. By building on existing solutions for homogeneous aquifers, comprehensive solutions are presented for composite aquifers comprising of arbitrary (finite) numbers of contiguous homogeneous sub-aquifers and subject to sinusoidal linear boundary conditions. Both Cartesian and radial coordinate systems are considered. Properties of the solutions, including rapid phase shifting and attenuation effects, are discussed and their practical relevance noted. Consequent modal dispersive effects on tidal waveforms are also examined via tidal constituent analysis. It is demonstrated that, for multi-constituent tidal forcings, measured peak heights of head oscillations can seem to increase, and phase lags seem to decrease, with distance from the forcing boundary unless constituents are separated and considered in isolation.     

A Field Test of Electromigration as a Method for Remediating Sulfate from Shallow Ground Water

Eloctromigraiion offers a potential tool for remediating ground water contaminated with highly soluble components, such as Na⁺, Cl, NO₃ and SO₄⁻. A field experiment was designed to lest the efficacy of electromigration for preconcontrating dissolved SO₄² in ground water associated with a fossil-fuel power plant. Two shallow wells, 25 feel apart (one 25 feel deep, the other 47 feet deep), were constructed in the upper portion of an unconfined alluvial aquifer. The wells were constructed with a double-wall design, with an outer casing of 4-inch PVC and an inner lube of 2-inch FVC; both were fully slotted (0.01 inch). Electrodes were constructed by wrapping the inner lulling with a 100-foot length of rare-earth metal oxide/copper wire. An electrical potential of 10.65 volts DC Was applied, and tests were run for periods of 12, 44, and 216 hours. Results showed large changes in the pH from the initial pH of ground water of about 7.5 to values of approximately 2 and 12 at the anode and cathode, respectively. Despite the fact that the test conditions were far from ideal, dissolved SO₄^2-; was significantly concentrated at the anode. Over a period of approximately nine days, the concentration of SO₄^2- at the anode reached what appeared to he a steady-state value of 2200 mg/L. compared lo the initial value in ground water of approximately 1150 mg/L. The results of this field lest should encourage further investigation of electromigration as a tool in the remediation of contaminated ground water.     

Cyclodextrin-Enhanced Vertical Flushing of a Trichloroethene Contaminated Aquifer

Pilot-scale testing of an innovative ground water remediation technology was conducted in a source zone of a trichloroethene-contaminated Superfund site in Tucson, Arizona. The technology is designed to enhance the removal of low-solubility organic contaminants from heterogeneous sedimentary aquifers by using a dual-screened vertical circulation well to inject and extract solutions containing a complexing sugar (hydroxypropyl-beta-cyclodextrin (HPCD]). Prior to initiating the pilot test, tracer tests were conducted to determine hydraulic characteristics of the vertical flow field and to evaluate trichloroethene-elution behavior during water flushing. The pilot test involved injecting approximately 4 m³ of a 20% HPCD solution into the upper screened interval of the well and extracting from the lower screened interval. The results of the pilot test indicate that the cyclodextrin solution increased the rate of trichloroethene removal from the aquifer. The concentrations of trichloroethene in the ground water extracted from the lower screened interval of the well increased by a factor of three (∼750 μg/L) in the presence of the cyclodextrin pulse, compared to concentrations obtained during previous water flushing (∼250 μg/L). Furthermore, the concentration of trichloroethene in water collected from the circulation well under static conditions was reduced to 6% of the levels measured prior to the test.     

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.     

Impact of coastal land reclamation on ground water level and the sea water interface

Land reclamation in coastal areas may have a significant effect on local ground water systems. Steady-state analytic solutions based on Dupuit and Ghyben-Herzberg assumptions are derived to evaluate this effect. Two situations are considered, both with ground water flow resulting from precipitation recharge: the coastal aquifer of an extensive landmass and an island. The results show that after reclamation, the water table rises and the salt water-fresh water interface moves seaward. The degree of these changes depends on the extent of reclamation and the hydraulic conductivity of the fill material. For the island situation, the reclamation displaces the ground water divide and changes the ground water conditions in the entire island. An unintended advantage of the reclamation is an increase of fresh ground water resource because the reclaimed land can be an additional aquifer and rain recharge takes place over a larger area.     

Urban Ground Water Pollution: A Case Study in Cuttack City, India

A large number of ground water samples (360) was collected from 60 stations over six consecutive seasons to study the influence of the main sewerage drain on shallow ground water table beneath the municipal area of Cuttack, India. A majority of the samples collected from stations close to the drain exceeded the maximum permissible limits set by the World Health Organization (WHO). Almost all the samples near the drain exceeded the WHO limit for NO₃^- and Na⁺. However, the concentrations decreased as the distance from the drain increased. The winter season registered the maximum concentrations of NH₄⁺, NO₃^-, and SO₄^2- ions whereas the minimum values always coincided with the rainy season. R-mode factor analysis was conducted to find relationships amongst the 16 chemical parameters studied. Fluoride showed a negative correlation with Cl^-, Na⁺, NO₃^-, SO₄^2-, and PO₄^3-. The concentration of F^- may be lower in raw waste water than naturally occurs in the ground water. Therefore, a decrease in the concentration of F^- near the drain may be attributed to dilution by contributions of waste water to the ground water. The rest of the parameters were found to be directly related to the distance of collection points to the sewerage. The distribution of nutrients is strongly affected by leaching of waste water into the ground water.     

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.     

Steady flow to a well near a stream with a leaky bed

We present an explicit analytic solution for steady, two-dimensional ground water flow to a well near a leaky streambed that penetrates the aquifer partially. Leakage from the stream is approximated as occurring along the centerline of the stream. The problem domain is infinite and pumping on one side of the stream induces flow on the other side. The solution includes the effects of uniform flow in the far field and a sloping hydraulic head in the stream. We use the solution to investigate the interaction between ground water and surface water in the stream, the effects of pumping on the opposite side of the stream, and the effects of the leaky streambed on the capture zone envelope of the well. We develop a relationship between parameters such that the pumping well will not capture water from the stream, or from the opposite side of the stream. When the discharge of the well is large enough to capture water from the stream, the shape of the capture zone envelope depends on flow conditions on the side of the stream opposite the well.     

Tide-induced head fluctuations in coastal aquifers of variable thickness

In this work, a new analytical solution to describe tide-induced head fluctuations in aquifers of variable thickness is presented. The proposed model assumes a finite and confined aquifer with a thickness that increases or decreases quadratically with the distance to the coast. A closed-form analytical solution is obtained by solving a boundary-value problem with both a separation of variables method and a change of variables method. This solution is a generalization of the solution obtained by Cuello et al., Hydrogeological Journal, 2017, 25, 1509–1515. The analytical solution is expressed in terms of the wedging parameter, a parameter that depends on the length and thicknesses at the coast and at the inland edge of the aquifer. Positive values of the wedging parameter describe aquifers with increasing thickness towards land and negative values describe aquifers with a decreasing thickness in the inland direction. The comparison of the new solution and the solution for a finite aquifer with constant thickness indicates that the sign of the wedging parameter enhances or decreases the amplitude of the tide-induced signal. However, the differences in time-lag between both solutions are negligible near the coast. The slope factor, which quantifies the inconsistencies between aquifer diffusivities estimated from attenuation and time-lag data, is computed and analysed. Near the coast, slope factor values greater than one are obtained for negative wedging parameters while slope factor values less than one are obtained for positive wedging parameters. The analysis of the new solution also indicates that more reliable estimates of the hydraulic diffusivity can be obtained from time-lag data.     

潮汐水位振幅和位相变化研究及其在地下水异常分析中的应用

将多孔介质中井-含水层-隔水层的潮汐水位振幅和位相的计算公式推广到裂隙饱水岩体潮汐分析中, 分析了裂隙含水层中井与裂隙, 裂隙与微裂隙的潮汐孔压响应原理和水流交换过程, 提取了影响裂隙含水层潮汐水位振幅和位相的主要因素, 应用井-裂隙排水产生的井水位-引潮高的振幅比A和位相差α₂主要随径向等效导水系数T同向变化, 裂隙和微裂隙(孔隙)排水产生的孔压-引潮高的振幅比D和位相差α₁主要随不排水条件下微裂隙与裂隙间振幅比E'/E反向变化的规律, 提出了潮汐井水位振幅和位相的8种不同变化类型, 分析了不同类型所反映的含水层形变, 并用于分析东川、 弥勒和西昌川03等3口井井水位振幅和位相变化的成因.      

Application of spectral analysis to determine hydraulic diffusivity of a sandy aquifer (Pingtung County,Taiwan)

The study demonstrates spectral relationships in the time–frequency domain for one‐dimensional groundwater flow in aquifers bounded by fluctuating boundaries. By nature, the solutions of spectral equations are non‐linear complex functions. To determine hydraulic diffusivity in the governing equations, it is required that the data are collected from the spectra of water levels at the fluctuating boundaries and observation wells. Hydraulic diffusivity thus can be obtained by an iterative inverse approach. This paper presents an application in Pingtung County of Taiwan to determine the hydraulic diffusivity of a sandy aquifer under confined conditions. Spectral density function of water level obtained from tidal boundaries and observation wells are used to approximate hydraulic diffusivity, which yields an averaged value of 1·26 × 10⁶ m²/h. Copyright © 2004 John Wiley & Sons, Ltd.     

Bacterial Production Stimulated Across the Zone of Influence of a Ground Water Circulation Well in a BTEX-Contaminated Aquifer

Benzene, toluene, ethylbenzene, and xylene (BTEX) hydrocarbons are typically the most abundant carbon source for bacteria in gasoline-contaminated ground water. In situ bioremediation strategies often involve stimulating bacterial heterotrophic production in an attempt to increase carbon demand of the assemblage. This may, in turn, stimulate biodegradation of contaminant hydrocarbons. In this study, ground water circulation wells (GCWs) were used as an in situ treatment for a fuel-contaminated aquifer to stimulate bacterial production, purportedly by increasing oxygen transfer to the subsurface, circulating limiting nutrients, enhancing bioavailability of hydrocarbons, or by removing metabolically inhibitory volatile organics. Bacterial production, as measured by rates of bacterial protein synthesis, was stimulated across the zone of influence (ZOI) of a series of GCWs. Productivity increased from ∼10² to >10⁵ ng C/L hour across the ZOI, suggesting that treatment stimulated overall biodegradation of carbon sources present in the ground water. However, even if BTEX carbon met all bacterial carbon demand, biodegradation would account for <4.3% of the total estimated BTEX removed from the ground water. Although bacterial productivity measurements alone cannot prove the effectiveness of in situ bioremediation, they can estimate the maximum amount of contaminant that may be biodegraded by a treatment system.     

Delayed impact of new volcanic ejecta on ground water quality

We observed long-term changes in the concentrations of dissolved ions in ground water caused by leachate from new volcanic ejecta deposited on the ground surface of the volcanic Miyakejima Island, Japan. Water samples were collected from nine wells and two rain collectors over a period of more than 10 years, and samples of runoff water were collected periodically. The samples were analyzed for temperature, pH, alkalinity, Cl⁻, and SO₄²⁻; some of the samples were also analyzed for δ¹³C. Because the leachate from the volcanic ejecta contained sulfate, we recorded an increase in SO₄²⁻ concentrations in the (unconfined) well water. The increase in SO₄²⁻ was initially detected between less than 1.4 and 5.2 years after the eruption, showing peak concentrations from 2.4 to 6.4 years after the eruption. This delayed response reflects the transit time of downward-moving SO₄²⁻ in the vadose zone, corresponding to an apparent movement rate of 0.4 to 7.2 cm/d. The rate relates to the mean recharge, represented as a fraction of local mean rainfall, and is calculated using the Cl⁻ balance method. The magnitude of the recorded increases reflects the volume of volcanic mudflow on the ground surface within the basin. For the management of ground water after an eruption, it is therefore important to know the chemical properties of the volcanic ejecta and the spatial distribution of mudflow to estimate the magnitude of the effect of ejecta on ground water quality.     

Fate of MTBE Relative to Benzene in a Gasoline-Contaminated Aquifer (1993–98)

Methyl tert -butyl ether (MTBE) and benzene have been measured since 1993 in a shallow, sandy aquifer contaminated by a mid-1980s release of gasoline containing fuel oxygenates. In wells downgradient of the release area, MTBK was detected before benzene, reflecting a chromatographic-like separation of these compounds in the direction of ground water flow. Higher concentrations of MTBE and benzene were measured in the deeper sampling ports of multilevel sampling wells located near the release area, and also up to 10 feet (3 m) below the water table surface in nested wells located farther from the release area. This distribution of higher concentrations at depth is caused by recharge events that deflect originally horizontal ground water flowlines. In the laboratory, microcosms containing aquifer material incubated with uniformly labeled ¹⁴C-MTBE under aerobic and anaerobic. Fe(III)-reducing conditions indicated a low but measurable biodegradation potential (<3%¹⁴C-MTBW as ¹⁴CO₂) after a seven-month incubation period, Tert -butyl alcohol (TBA), a proposed microbial-MTBE transformation intermediate, was detected in MTBE-contaminated wells, but TBA was also measured in unsaturated release area sediments. This suggests that TBA may have been present in the original fuel spilled and does not necessarily reflect microbial degradation of MTBE. Combined, these data suggest that milligram per liter to microgram per liter decreases in MTBE concentrations relative to benzene are caused by the natural attenuation processes of dilution and dispersion with less-contaminated ground water in the direction of flow rather than biodegradation at this point source gasoline release site.