Statistical estimation of the relative efficiency of natural attenuation mechanisms in contaminated aquifers

Using monitored natural attenuation is an increasingly popular strategy for dealing with contaminated aquifers. This paper provides a statistical methodology for the estimation of the relative efficiency of natural attenuation mechanisms. The methodology provides estimates, with associated measures of uncertainty, of the relative efficiency of four types of bio-degradation (oxidation using oxygen as the electron-acceptor, denitrification, iron reduction and sulfate reduction). A data set from Trecate, Italy, is analysed using the methodology. The analysis shows that sulfate is the main cause of hydrocarbon removal on this site. It is also seen that oxidation using oxygen seems to be more preferential than the other reactions, in the sense that this reaction is relatively more efficient than other reactions at locations where the hydrocarbon concentration is low.     

Isotope methods for management of shared aquifers in northern Africa

Access to fresh water is one of the major issues of northern and sub-Saharan Africa. The majority of the fresh water used for drinking and irrigation is obtained from large ground water basins where there is minor contemporary recharge and the aquifers cross national borders. These aquifers include the Nubian Aquifer System shared by Chad, Egypt, Libya, and Sudan; the Iullemeden Aquifer System, extending over Niger, Nigeria, Mali, Benin, and Algeria; and the Northwest Sahara Aquifer System shared by Algeria, Libya, and Tunisia. These resources are subject to increased exploitation and may be severely stressed if not managed properly as witnessed already by declining water levels. In order to make appropriate decisions for the sustainable management of these shared water resources, planners and managers in different countries need an improved knowledge base of hydrological information. Three technical cooperation projects related to aquifer systems will be implemented by the International Atomic Energy Agency, in collaboration with the United Nations Educational, Scientific and Cultural Organization and United Nations Development Programme/Global Environmental Facility. These projects focus on isotope hydrology studies to better quantify ground water recharge and dynamics. The multiple isotope approach combining commonly used isotopes 18O and 2H together with more recently developed techniques (chlorofluorocarbons, 36Cl, noble gases) will be applied to improve the conceptual model to study stratification and ground water flows. Moreover, the isotopes will be an important indicator of changes in the aquifer due to water abstraction, and therefore they will assist in the effort to establish a sustainable ground water management.     

Relationships between microbial communities and groundwater chemistry in two pristine confined groundwater aquifers in central China

This study explores linkages between the microbial composition and hydrochemical variables of pristine groundwater to identify active redox conditions and processes. Two confined aquifers underlying the city of Qianjiang in the Jianghan Plain in China were selected for this study, having different recharge sources and strong hydrochemical gradients. Typical methods for establishing redox processes according to threshold concentration criteria for geochemical parameters suggest iron or sulphate reduction processes. High‐throughput 16S rRNA sequencing was used to obtain diversity and taxonomic information on microbial communities. Instead of revealing iron‐ and sulphate‐reducing bacteria, salt‐ and alkali‐tolerant bacteria, such as the phylum Firmicutes and the class Gammaproteobacteria, and in particular, the family Bacillaceae, were dominant in the downstream groundwater of the first aquifer that had high ion concentrations caused by the dissolution of calcite and dolomite; meanwhile, the heterotrophic microaerophilic families Comamonadaceae and Rhodocyclaceae prevailed in the upstream groundwater of the first aquifer. Sulphate‐reducing bacteria were extremely abundant in the upstream groundwater of the second aquifer, as the SO₄^2? concentration was especially high. Methanogens and methanotrophs were predominant in the downstream groundwater of the second aquifer even though the concentration of SO₄^2? was much higher than 0.5 mg L^?1. The microbial communities, together with the geochemical parameters, indicated that the upstream region of the first aquifer was suboxic, that Fe(III) and Mn(IV) reductions were not the main redox processes in the downstream groundwater of the first aquifer with high Fe and Mn concentrations, and that the redox processes in the upstream and downstream regions of the second confined aquifer were SO₄^2? reduction and methanogenesis, respectively. This study expands understanding of the linkages between microbial communities and hydrogeochemistry in pristine groundwaters and provides more evidence for identifying active redox conditions and processes.     

Contaminant transport in heterogeneous aquifers: A critical review of mechanisms and numerical methods of non-Fickian dispersion

Natural aquifers usually exhibit complex physical and chemical heterogeneities, which are key factors complicating kinetic processes, such as contaminant transport and transformation, posing a great challenge in the remediation of contaminated groundwater. Aquifer heterogeneity usually leads to a distinct feature, the so-called “anomalous transport” in groundwater,which deviates from the phenomenon described by the classical advection-dispersion equation(ADE) based on Fick's Law.Anomalous transport, also known as non-Fickian dispersion or “anomalous dispersion” in a broad sense, can explain the hydrogeological mechanism that leads to the temporally continuous deterioration of water quality and rapid spatial expansion of pollutant plumes. Contaminants enter and then are retained in the low-permeability matrix from the high-permeability zone via molecular diffusion, chemical adsorption, and other mass exchange effects. This process can be reversed when the concentration of pollutants in high-permeability zones is relatively low. The contaminants slowly return to the high-permeability zones through reverse molecular diffusion, resulting in sub-dispersive anomalous transport leading to the chronic gradual deterioration of water quality. Meanwhile, some contaminants are rapidly transported along the interconnected preferential flow paths, resulting in super-dispersive anomalous transport, which leads to the rapid spread of contaminants. Aquifer heterogeneity is also an important factor that constrains the efficacy of groundwater remediation, while the development, application, and evaluation of groundwater remediation technologies are usually based on the Fickian dispersion process predicted by the ADE equation.Comprehensive studies of the impacts of non-Fickian dispersion on contaminant transport and remediation are still needed. This article reviews the non-Fickian dispersion phenomenon caused by the heterogeneity of geological media, summarizes the processes and current understanding of contaminant migration and transformation in highly heterogeneous aquifers, and evaluates mathematical methods describing the main non-Fickian dispersion features. This critical review also discusses the limitations of existing research and outlines potential future research areas to advance the understanding of mechanisms and modeling of non-Fickian dispersion in heterogeneous media.     

Biodegradation of crude oil using an efficient microbial consortium in a simulated marine environment

Ochrobactrum sp. N1, Brevibacillus parabrevis N2, B. parabrevis N3 and B. parabrevis N4 were selected when preparing a mixed bacterial consortium based on the efficiency of crude oil utilization. A crude oil degradation rate of the N-series microbial consortium reached upwards of 79% at a temperature of 25 °C in a 3.0% NaCl solution in the shake flask trial. In the mesocosm experiment, a specially designed device was used to simulate the marine environment. The internal tank size was 1.5 m (L)×0.8 m (W)×0.7 m (H). The microbial growth conditions, nutrient utilization and environmental factors were thoroughly investigated. Over 51.1% of the crude oil was effectively removed from the simulated water body. The escalation process (from flask trials to the mesocosm experiment), which sought to represent removal under conditions more similar to the field, proved the high efficiency of using N-series bacteria in crude oil degradation.     

Stochastic analysis of bounded unsaturated flow in heterogeneous aquifers: Spectral/perturbation approach

This paper describes a stochastic analysis of steady state flow in a bounded, partially saturated heterogeneous porous medium subject to distributed infiltration. The presence of boundary conditions leads to non-uniformity in the mean unsaturated flow, which in turn causes non-stationarity in the statistics of velocity fields. Motivated by this, our aim is to investigate the impact of boundary conditions on the behavior of field-scale unsaturated flow. Within the framework of spectral theory based on Fourier–Stieltjes representations for the perturbed quantities, the general expressions for the pressure head variance, variance of log unsaturated hydraulic conductivity and variance of the specific discharge are presented in the wave number domain. Closed-form expressions are developed for the simplified case of statistical isotropy of the log hydraulic conductivity field with a constant soil pore-size distribution parameter. These expressions allow us to investigate the impact of the boundary conditions, namely the vertical infiltration from the soil surface and a prescribed pressure head at a certain depth below the soil surface. It is found that the boundary conditions are critical in predicting uncertainty in bounded unsaturated flow. Our analytical expression for the pressure head variance in a one-dimensional, heterogeneous flow domain, developed using a nonstationary spectral representation approach [Li S-G, McLaughlin D. A nonstationary spectral method for solving stochastic groundwater problems: unconditional analysis. Water Resour Res 1991;27(7):1589–605; Li S-G, McLaughlin D. Using the nonstationary spectral method to analyze flow through heterogeneous trending media. Water Resour Res 1995; 31(3):541–51], is precisely equivalent to the published result of Lu et al. [Lu Z, Zhang D. Analytical solutions to steady state unsaturated flow in layered, randomly heterogeneous soils via Kirchhoff transformation. Adv Water Resour 2004;27:775–84].     

Analytical stochastic solutions of saltwater/freshwater interface in coastal aquifers

Three analyitcal solutions of saltwater intrusion under uncertain hydrologic, hydrogeologic, and pumping conditions are presented. These solutions extends the existing deterministic, sharp interface solutions to stochastic ones. The randomness is represented in terms of statistical measures of mean, standard deviation and covariance. The analysis is based on perturbation using Taylor series expansion. Simulations based on probability distribution are conducted for verification.     

Analytical stochastic solutions of saltwater/freshwater interface in coastal aquifers

Three analyitcal solutions of saltwater intrusion under uncertain hydrologic, hydrogeologic, and pumping conditions are presented. These solutions extends the existing deterministic, sharp interface solutions to stochastic ones. The randomness is represented in terms of statistical measures of mean, standard deviation and covariance. The analysis is based on perturbation using Taylor series expansion. Simulations based on probability distribution are conducted for verification.     

Comparison of three methods of hydrogeological parameter estimation in leaky aquifers using transient flow pumping tests

Accurate knowledge of hydrogeological parameters is essential for groundwater modeling, protection and remediation. Three methods (type curve fitting method, inflection point method and global curve‐fitting method (GCFM)) which are frequently applied in the estimation of leaky aquifer parameters were compared using synthetic pumping tests. The results revealed GCFM could provide best parameter estimation among the three methods with fewer uncertainties associated with the processes of parameter estimation. GCFM was also found to be both time saving and of low cost and is thus more preferable for hydrogeological parameter estimation than the other two methods. Copyright © 2013 John Wiley & Sons, Ltd.     

Arsenic in glacial aquifers: sources and geochemical controls

A total of 176 wells in sand-and-gravel glacial aquifers in central Illinois were sampled for arsenic (As) and other chemical parameters. The results were combined with archived and published data from several hundred well samples to determine potential sources of As and the potential geochemical controls on its solubility and mobility. There was considerable spatial variability in the As concentrations. High concentrations were confined to areas smaller than 1 km in diameter. Arsenic and well depth were uncorrelated. Arsenic solubility appeared to be controlled by oxidation-reduction (redox) conditions, especially the presence of organic matter. Geochemical conditions in the aquifers are typically reducing, but only in the most reducing water does As accumulate in solution. In wells in which total organic carbon (TOC) was below 2 mg/L and sulfate (SO4(2-)) was present, As concentrations were low or below the detection limit (0.5 microg/L). Arsenic concentrations >10 microg/L were almost always found in wells where TOC was >2 mg/L and SO4(2-) was absent or at low concentrations, indicating post-SO4 (2-)reducing conditions. Iron (Fe) is common in the aquifer sediments, and Fe oxide reduction appears to be occurring throughout the aquifers. Arsenic is likely released from the solid phase as Fe oxide is reduced.     

Describing near surface, transient flow processes in unconfined aquifers below irrigated lands: A deforming finite element model for heterogeneous aquifers

Hydrologic models of irrigated lands generally adopt either a basin-scale or a root-zone perspective. While basin-wide macro-scale models rely on the aggregation of important spatial and temporal data across large areas, micro-scale root-zone models depend on the definition of rigid boundaries around the zone of plant–soil–water interaction. In reality, irrigation management decisions are made on a field by field basis and can interact across field boundaries. This paper first describes a shallow water table model, based on deforming finite element (DFE) framework, to characterize the near-surface field-to-field hydrologic response to various irrigation and drainage management regimes along a gently sloping alluvial fan. The model is then enhanced through changing geometry of a fluctuating water table below a series of irrigated fields. Such an enhancement also offers computational flexibility relative to the saturated–unsaturated models commonly used in micro-scale studies. The model is designed with the alluvial fan aquifers of California’s western San Joaquin Valley as reference systems.     

An exploration of incorporating site response into PSHA—Part I: Issues related to site response analysis methods

This paper presents a series of analyses for the evaluation of the ground response of two NEHRP class D sites, subjected to shaking by a large number of strong ground-motion records. The two investigated sites have very distinct profiles, but they are characterised by almost identical V_s30 values. The site response analyses are performed using various methods of analysis and input parameters in order to explore the sensitivity of the ground response estimates and to identify the dominating parameters. Equivalent linear analysis is performed using different sets of dynamic soil properties curves, while nonlinear analysis is performed using different target dynamic soil curves, viscous damping formulations and fitting procedures for the constitutive model parameters. Particular focus is given to the sensitivity of the response when soil sites are subjected to high-intensity shaking, a subject of particular interest when the prediction of surface ground motions with low annual probabilities of exceedance is the target of probabilistic seismic hazard analyses (PSHA). The site response analysis results of this paper are incorporated into the probabilistic framework of Bazzurro and Cornell [1] in our companion paper in order to assess their impact on the final soil surface hazard calculation.     

Evidence of methylmercury production and modification of the microbial community structure in estuary sediments contaminated with wastewater treatment plant effluents

The Seine’s estuary (France) waters are the receptacle of effluents originating from wastewater treatment plants (WWTP). In this estuary, mudflats are deposition zones for sediments and their associated contaminants, and play an essential role in the mercury (Hg) biogeochemical cycle mainly due to indigenous microorganisms. Microcosms were used to assess the impact of WWTP-effluents on mercury methylation by monitoring Hg species (total dissolved Hg in porewater, methylmercury and total mercury) and on microbial communities in sediments. After effluent amendment, methylmercury (MeHg) concentrations increased in relation with the total Hg and organic matter content of the WWTP-effluents. A correlation was observed between MeHg and acid-volatile-sulfides concentrations. Quantification of sulfate-reducing microorganisms involved in Hg methylation showed no increase of their abundance but their activity was probably enhanced by the organic matter supplied with the effluents. WWTP-effluent spiking modified the bacterial community fingerprint, mainly influenced by Hg contamination and the organic matter amendment.     

High-resolution simulation and characterization of density-driven flow in CO2 storage in saline aquifers

Simulations are routinely used to study the process of carbon dioxide (CO₂) sequestration in saline aquifers. In this paper, we describe the modeling and simulation of the dissolution–diffusion–convection process based on a total velocity splitting formulation for a variable-density incompressible single-phase model. A second-order accurate sequential algorithm, implemented within a block-structured adaptive mesh refinement (AMR) framework, is used to perform high-resolution studies of the process. We study both the short-term and long-term behaviors of the process. It is found that the onset time of convection follows closely the prediction of linear stability analysis. In addition, the CO₂ flux at the top boundary, which gives the rate at which CO₂ gas dissolves into a negatively buoyant aqueous phase, will reach a stabilized state at the space and time scales we are interested in. This flux is found to be proportional to permeability, and independent of porosity and effective diffusivity, indicative of a convection-dominated flow. A 3D simulation further shows that the added degrees of freedom shorten the onset time and increase the magnitude of the stabilized CO₂ flux by about 25%. Finally, our results are found to be comparable to results obtained from TOUGH2-MP.     

Slug test in confined aquifers, the over-damped case: quasi-steady flow analysis

In the great majority of slug tests performed in wells fully penetrating confined geologic formations, and for over-damped conditions, the response data are evaluated with the transient-flow model of Cooper et al. (1967) when the radial hydraulic conductivity K(r) and the coefficient of specific storage S(s) are to be estimated. That particular analytical solution, however, is computationally involved and awkward to use. Thus, groundwater professionals often use a few pre-prepared type-curves to fit the data by a rough matching procedure, visually or computationally. On the other hand, the method of Hvorslev (1951), which assumes the flow to be quasi-steady, is much simpler but yields only K(r) estimates. In this work, we develop a complete quasi-steady flow model that includes a storage balance inside the aquifer and allows estimating K(r) and S(s) simultaneously, through matching of the well response data to a type-curve. The new model approximates the model of Cooper et al. closely and has the practical advantage that its solution type-curves are generated easily using an electronic spreadsheet, so that the optimal fit of data by a type-curve can be readily automated.     

Performance evaluation of interpolation methods for incorporating rain gauge measurements into NEXRAD precipitation data: a case study in the Upper Guadalupe River Basin

High spatial and temporal resolution of precipitation data is critical input for hydrological budget estimation and flash flood modelling. This study evaluated four methods [Bias Adjustment (BA), Simple Kriging with varying Local Means (SKlm), Kriging with External Drift (KED), and Regression Kriging (RK)] for their performances in incorporating gauge rainfall measurements into Next Generation Weather Radar (NEXRAD) multi‐sensor precipitation estimator (MPE; hourly and 4 × 4 km²). Measurements from a network of 50 gauges at the Upper Guadalupe River Basin, central Texas and MPE data for the year 2004 were used in the study. We used three evaluation coefficients percentage bias (PB), coefficient of determination (R²), and Nash–Sutcliffe efficiency (NSE) to examine the performance of the four methods for preserving regional‐ and local‐scale characteristics of observed precipitation data. The results show that the two Kriging‐based methods (SKlm and RK) are in general better than BA and KED and that the PB and NSE criteria are better than the R² criterion in assessing the performance of the four methods. It is also worth noting that the performance of one method at regional scale may be different from its performance at local scale. Critical evaluation of the performance of different methods at local or regional scale should be conducted according to the different purposes. The results obtained in this study are expected to contribute to the development of more accurate spatial rainfall products for hydrologic budget and flash flood modelling. Copyright © 2011 John Wiley & Sons, Ltd.     

Accuracy of numerical simulations of contaminant transport in heterogeneous aquifers: A comparative study

This work deals with a comparison of different numerical schemes for the simulation of contaminant transport in heterogeneous porous media. The numerical methods under consideration are Galerkin finite element (GFE), finite volume (FV), and mixed hybrid finite element (MHFE). Concerning the GFE we use linear and quadratic finite elements with and without upwind stabilization. Besides the classical MHFE a new and an upwind scheme are tested. We consider higher order finite volume schemes as well as two time discretization methods: backward Euler (BE) and the second order backward differentiation formula BDF (2). It is well known that numerical (or artificial) diffusion may cause large errors. Moreover, when the Péclet number is large, a numerical code without some stabilising techniques produces oscillating solutions. Upwind schemes increase the stability but show more numerical diffusion. In this paper we quantify the numerical diffusion for the different discretization schemes and its dependency on the Péclet number. We consider an academic example and a realistic simulation of solute transport in heterogeneous aquifer. In the latter case, the stochastic estimates used as reference were obtained with global random walk (GRW) simulations, free of numerical diffusion. The results presented can be used by researchers to test their numerical schemes and stabilization techniques for simulation of contaminant transport in groundwater.     

土壤—地下水系统非水相液体污染区的电性特征及电阻率法探测

从土壤-地下水系统被非水相液体污染后引起的导电性变化、污染区地电模型演化、电阻率异常认识及量化解释等方面分析了国内外利用电阻率法探测这一问题的研究现状.已有研究表明,污染区电性随污染程度、污染时间变化显示一定的动态性;污染区地电模型建立是个复杂过程,不仅反映污染物的扩散过程,还要考虑地质、地下水、土层含水量等影响因素,并随外部条件和时间变化显示一定的动态性;这也决定了实测电剖面异常特征的复杂性.由于基于电测数据的真电阻率获取及污染土电性基本模型建立都尚未解决,利用电性探测结果进行污染量化分析的研究还没有形成统一的认识和结论,这也是未来很长时间内需要攻克的难题.