Biogeochemical and isotopic gradients in a BTEX/PAH contaminant plume: model-based interpretation of a high-resolution field data set

Multicomponent and Surface Charge Effects on PFOS Sorption and Transport in Goethite-Coated Porous Media under Variable Hydrochemical Conditions

Perfluorooctanesulfonate (PFOS), a toxic anionic perfluorinated surfactant, exhibits variable electrostatic adsorption mechanisms on charge-regulated minerals depending on solution hydrochemistry. This work explores the interplay of multicomponent interactions and surface charge effects on PFOS adsorption to goethite surfaces under flow-through conditions. We conducted a series of column experiments in saturated goethite-coated porous media subjected to dynamic hydrochemical conditions triggered by step changes in the electrolyte concentration of the injected solutions. Measurements of pH and PFOS breakthrough curves at the outlet allowed tracking the propagation of multicomponent reactive fronts. We performed process-based reactive transport simulations incorporating a mechanistic network of surface complexation reactions to quantitatively interpret the geochemical processes. The experimental and modeling outcomes reveal that the coupled spatio-temporal evolution of pH and electrolyte fronts, driven by the electrostatic properties of the mineral, exerts a key control on PFOS mobility by determining its adsorption and speciation reactions on goethite surfaces. These results illuminate the important influence of multicomponent transport processes and surface charge effects on PFOS mobility, emphasizing the need for mechanistic adsorption models in reactive transport simulations of ionizable PFAS compounds to determine their environmental fate and to perform accurate risk assessment.

Microscale Chaotic Mixing as a Driver for Chemical Reactions in Porous Media

Mixing-induced reactions play a key role in a large range of biogeochemical and contaminant transport processes in the subsurface. Fluid flow through porous media was recently shown to exhibit chaotic mixing dynamics at the pore scale, enhancing microscale concentration gradients and controlling mixing rates. While this phenomenon is likely ubiquitous in environmental systems, it is not known how it affects chemical reactions. Here, we use refractive index matching and laser-induced fluorescence imaging of a bimolecular redox reaction to investigate the consequence of pore scale chaotic mixing on the reaction rates. The overestimation of measured reaction rates by the classical macrodispersion model highlights the persistence of incomplete mixing on the pore scale. We show that the reaction product formation is controlled by microscale chaotic mixing, which induces an exponential increase of the mixing interface and of the reaction rates. We derive a reactive transport model that captures experimental results and predicts that chaotic mixing has a first order control on reaction rates across a large range of time scales and Péclet and Damköhler numbers. These findings provide a new framework for understanding, assessing, and predicting mixing-induced reactions and their role on the fate and mobility of environmental compounds in natural porous media.

Experimental research on the transport-transformation of organic contaminants under the influence of multi-field coupling at a site scale

Organic-contaminated shallow aquifers have become a global concern of groundwater contamination, yet little is known about the coupled effects of hydrodynamic-thermal-chemical-microbial (HTCM) multi-field on organic contaminant transport and transformation over a short time in aquifers. Therefore, this study proposed a quick and efficient field experimental method for the transport-transformation of contaminants under multi-field coupling to explore the relationship between organic contaminants (total petroleum hydrocarbon (TPH), polycyclic aromatic hydrocarbons (PAHs), benzene-toluene-ethylbenzene-xylene (BTEX) and phthalates acid esters (PAEs)) and multi-field factors. The results showed that hydrodynamics (affecting pH, p < 0.001) and temperature (affecting dissolved oxygen, pH and HCO3-, p < 0.05) mainly affected the organic contaminants indirectly by influencing the hydrochemistry to regulate redox conditions in the aquifer. The main degradation reactions of the petroleum hydrocarbons (TPH, PAHs and BTEX) and PAEs in the aquifer were sulfate reduction and nitrate reduction, respectively. Furthermore, the organic contamination was directly influenced by microbial communities, whose spatial patterns were shaped by the combined effects of the spatial pattern of hydrochemistry (induced by the organic contamination pressure) and other multi-field factors. Overall, our findings imply that the spatiotemporal patterns of organic contaminants are synergistically regulated by HTCM, with distinct mechanisms for petroleum hydrocarbons and PAEs.

Enhancing in situ biodegradation in groundwater using pump and treat remediation: a proof of concept and modelling analysis of controlling variables

A remediation approach which uses pump and treatment (PAT) to enhance the biodegradation of organic contaminants by increasing dispersive mixing between plumes and groundwater was evaluated for a phenol-contaminated aquifer, using a reactive transport model which simulates kinetic reactions between an electron donor (ED) in the plume and electron acceptor (EA) in the groundwater. The influence of system design and operation was examined in six modelling scenarios. Injection or extraction of groundwater increases biodegradation above no action and the location, pumping rate, and distance between well(s) are important variables which influence biodegradation. An increase in pumping rate, distance of the wells from the plume centreline, and changing the flow direction increase dispersive mixing between the plume and groundwater. This increases plume spreading and the plume fringe interface, providing a greater flux of dissolved EAs for biodegradation. In general, injection of groundwater containing natural EAs enhances biodegradation more than extraction. The enhancement of biodegradation is sensitive to the relative fluxes of ED and EA, as controlled by the arrangement of the wells. In the best performing scenario, biodegradation was enhanced by 128%, compared with no action.

Quantification of contaminant mass discharge from point sources in aquitard/aquifer systems based on vertical concentration profiles and 3D modeling

Point sources with contaminants, such as chlorinated solvents, per- and polyfluoroalkyl substances (PFAS), or pesticides, are often located in low-permeability aquitards, where they can act as long-term sources and threaten underlying groundwater resources. We demonstrate the use of a 3D numerical model integrating comprehensive hydrogeological and contamination data to determine the contaminant mass discharge (CMD) from an aquitard into the underlying aquifer. A mature point source with a dissolved chlorinated solvent in a clayey till is used as an example. The quantitative determination is facilitated by model calibration to high-resolution vertical concentration profiles obtained by direct-push sampling techniques in the aquifer downgradient of the contaminant source zone. The concentration profiles showed a plume sinking with distance from the source characteristic for such aquitard/aquifer settings. The sinking is caused by the interplay between infiltrating water and horizontal groundwater flow. The application of 3D solute transport modeling on high-resolution profiles allowed for determining the infiltration rate, the hydraulic conductivity in the aquitard, and, ultimately, the CMD. Different source zone conceptualizations demonstrate the potential effects of fractures and sorption in source zones in aquitards on CMD development. Fractures in the aquitard had a minor influence on the current CMD determined with the presented approach. Still, fractures with hydraulic apertures larger than 10 μm were crucial for the temporal development of the CMD and plume. A thorough characterization of the source zone conditions combined with high-resolution concentration profiles and detailed modeling is valuable for shedding light on the probable future development of groundwater contamination arising from sources in aquitard/aquifer settings and evaluating remedial actions.

Bioremediation of chlorinated ethenes contaminated groundwater and the reactive transport modeling - A review

Improper disposal of chlorinated ethenes (CEs), a class of widely used solvents in chemical manufacturing and cleaning industries, often leads to severe groundwater contamination. In situ bioremediation of CE-contaminated groundwater has received continuous attention in recent years. The reactive transport simulation is a valuable tool for planning and designing in situ bioremediation systems. This paper presents a detailed and comprehensive review on the main biotransformation pathways of CEs in aquifers, the mathematical modeling of bioremediation processes, and the available software developed for the simulation of reactive transport of CEs over past three decades. The aim of this research is to provide guidance on the selection of appropriate models and software suitable for systems of varying scales, and to discern prevailing research trends while identifying areas worthy of further study. This paper provides a detailed summary of the equations, parameters, and applications of existing biotransformation models from literature studies, highlighting the operation, benefits, and limitations of software available for CEs reactive transport simulations. Lastly, the support of reactive transport simulation programs for the design of full-scale in situ bioremediation systems was elucidated. Further research is needed for incorporating the effects of key subsurface environmental factors on biodegradation processes into models and balancing model complexity with computer data processing power to better support the development and application of reactive transport modeling software.

Isotope forensics of polycyclic aromatic compounds (PACs) in a contaminated shallow aquifer

Diesel was accidently released into the shallow subsurface at an industrial site in the province of Québec, Canada, in the late 1980s. Subsequent remediation efforts removed much of the contamination; however, traces of petroleum hydrocarbons continue to impact the local aquifer. In addition to the historical diesel spill, more recent yet unconfirmed accidental releases from ongoing on-site and neighbouring industrial activities may have potentially contributed to elevated levels of polycyclic aromatic compounds (PACs) in groundwater. To identify the main source(s) of contamination, compound-specific stable carbon isotope ratios (δ13C) of PACs in groundwater monitoring wells were compared to those in asphalt produced from a nearby plant and in fuel oil #6 oil being used by local industry. The δ13C values of five individual compounds (biphenyl, C2-naphthalene, C1-fluorene, dibenzothiophene and phenanthrene) and two groups of combined C1-phenanthrenes/anthracenes in all groundwater samples were within analytical uncertainty (±0.5‰). Moreover, the δ13CPAC values in groundwater samples were distinct from those in asphalt and fuel oil #6, indicating negligible contributions from these sources. The similarity in δ13CPAC values across monitoring wells, including one situated in the former source zone containing a floating hydrocarbon phase, pointed to a common source of subsurface contamination that was attributed to the historical diesel spill. These results thus demonstrate that δ13CPAC values can be used for source apportionment in shallow aquifers decades after the original spill event.

Evolution of plume geometry, dilution and reactive mixing in porous media under highly transient flow fields at the surface water-groundwater interface

Highly transient boundary conditions affect mixing of dissolved solutes in groundwater. An example of these transient boundary conditions occurs at the surface water-groundwater interface, where the water level in rivers can change rapidly due to the operation of hydropower plants, leading to a regime known as hydropeaking. Inspired by this phenomenon, this work studies at laboratory scale the effects of fluctuating surface water bodies on solute transport in aquifers. We performed flow-through experiments at two different flow velocities and under steady and transient flow conditions where a conservative tracer was injected in the system and its concentration measured with optical imaging methods. The experimental results were quantitatively interpreted with numerical simulations implementing a non-linear velocity-dependent dispersive transport model. We estimated plume dilution by computing the dilution index and its evolution as well as two key geometrical metrics of the transient plumes: the perimeter and the area. We further investigated reactive mixing and mixing enhancement considering mixing-controlled bimolecular reactions using different critical mixing ratios. In general, highly transient boundary conditions lead to a larger area, perimeter and plume dilution and the results show greater relative enhancement for the scenarios with low groundwater flow velocity. A linear relationship was observed between the evolution of the area and the dilution index of the plumes for the transient flow scenarios investigated. Considering reactive transport and mixing-limited reactions at the surface water-groundwater interface, we identified different dilution and reaction dominated regimes, characterized, respectively, by increasing and decreasing plume entropies at different mixing ratios of the reactants. Furthermore, reactive mixing was enhanced by transient flows leading to a faster degradation of contaminant plumes compared to corresponding steady flow conditions.

Sulfate reduction accelerates groundwater arsenic contamination even in aquifers with abundant iron oxides

Groundwater contamination by geogenic arsenic is a global problem affecting nearly 200 million people. In South and Southeast Asia, a cost-effective mitigation strategy is to use oxidized low-arsenic aquifers rather than reduced high-arsenic aquifers. Aquifers with abundant oxidized iron minerals are presumably safeguarded against immediate arsenic contamination, due to strong sorption of arsenic onto iron minerals. However, preferential pumping of low-arsenic aquifers can destabilize the boundaries between these aquifers, pulling high-arsenic water into low-arsenic aquifers. We investigate this scenario in a hybrid field-column experiment in Bangladesh where naturally high-arsenic groundwater is pumped through sediment cores from a low-arsenic aquifer, and detailed aqueous and solid-phase measurements are used to constrain reactive transport modelling. Here we show that elevated groundwater arsenic concentrations are induced by sulfate reduction and the predicted formation of highly mobile, poorly sorbing thioarsenic species. This process suggests that contamination of currently pristine aquifers with arsenic can occur up to over 1.5 times faster than previously thought, leading to a deterioration of urgently needed water resources.

Pore-Scale Heterogeneities Improve the Degradation of a Self-Inhibiting Substrate: Insights from Reactive Transport Modeling

In situ bioremediation is a common remediation strategy for many groundwater contaminants. It was traditionally believed that (in the absence of mixing-limitations) a better in situ bioremediation is obtained in a more homogeneous medium where the even distribution of both substrate and bacteria facilitates the access of a larger portion of the bacterial community to a higher amount of substrate. Such conclusions were driven with the typical assumption of disregarding substrate inhibitory effects on the metabolic activity of enzymes at high concentration levels. To investigate the influence of pore matrix heterogeneities on substrate inhibition, we use a numerical approach to solve reactive transport processes in the presence of pore-scale heterogeneities. To this end, a rigorous reactive pore network model is developed and used to model the reactive transport of a self-inhibiting substrate under both transient and steady-state conditions through media with various, spatially correlated, pore-size distributions. For the first time, we explore on the basis of a pore-scale model approach the link between pore-size heterogeneities and substrate inhibition. Our results show that for a self-inhibiting substrate, (1) pore-scale heterogeneities can consistently promote degradation rates at toxic levels, (2) the effect reverses when the concentrations fall to levels essential for microbial growth, and (3) an engineered combination of homogeneous and heterogeneous media can increase the overall efficiency of bioremediation.

Reclaimed Water Reuse for Groundwater Recharge: A Review of Hot Spots and Hot Moments in the Hyporheic Zone

As an alternative resource, reclaimed water is rich in the various nutrients and organic matter that may irreparably endanger groundwater quality through the recharging process. During groundwater recharge with reclaimed water, hot spots and hot moments (HSHMs) in the hyporheic zones, located at the groundwater–reclaimed water interface, play vital roles in cycling and processing energy, carbon, and nutrients, drawing increasing concern in the fields of biogeochemistry, environmental chemistry, and pollution treatment and prevention engineering. This paper aims to review these recent advances and the current state of knowledge of HSHMs in the hyporheic zone with regard to groundwater recharge using reclaimed water, including the generation mechanisms, temporal and spatial characteristics, influencing factors, and identification indicators and methods of HSHMs in the materials cycle. Finally, the development prospects of HSHMs are discussed. It is hoped that this review will lead to a clearer understanding of the processes controlling water flow and pollutant flux, and that further management and control of HSHMs can be achieved, resulting in the development of a more accurate and safer approach to groundwater recharge with reclaimed water.

Impact of solute charge and diffusion coefficient on electromigration and mixing in porous media

The application of electrokinetic techniques in porous media has great potential to enhance mass transfer rates and, thus, to mobilize contaminants and effectively deliver reactants and amendments. However, the transport mechanisms induced by the application of an external electric field are complex and entail the coupling of physical, chemical and electrostatic processes. In this study we focus on electromigration and we provide experimental evidence of the impact of compound-specific properties, such as the aqueous diffusivity and the valence of charged species, on the macroscopic electrokinetic transport. We performed a series of multidimensional experiments considering the displacement of three different tracer plumes (i.e., permanganate, allura red and new coccine) in different background electrolyte solutions. The outcomes of the experiments clearly show that both the compound-specific diffusivity and the charge of the injected and resident ions impact the transport of the selected color tracer plumes, whose evolution was monitored with image analysis. The investigated experimental scenarios led to distinct plume behavior characterized by different mass distribution, average displacement velocities, longitudinal and lateral plume spreading, shape of the invading and receding fronts, as well as dilution of the injected solutes. A numerical simulator, based on the Nernst-Planck-Poisson equations and on aqueous speciation reactions in the pore water, allowed us to quantitatively interpret the experimental results, to capture the observed patterns of plume evolution, and to illuminate the coupling between the governing physico-chemical mechanisms and the controlling role of small scale compound-specific and electrostatic properties. Finally, the model was also extended to a typical configuration of in situ electrokinetic remediation of contaminated groundwater to show the impact of such mechanisms at larger scale.

Mass-Transfer-Limited Biodegradation at Low Concentrations-Evidence from Reactive Transport Modeling of Isotope Profiles in a Bench-Scale Aquifer

Organic contaminant degradation by suspended bacteria in chemostats has shown that isotope fractionation decreases dramatically when pollutant concentrations fall below the (half-saturation) Monod constant. This masked isotope fractionation implies that membrane transfer is slow relative to the enzyme turnover at μg L^-1 substrate levels. Analogous evidence of mass transfer as a bottleneck for biodegradation in aquifer settings, where microbes are attached to the sediment, is lacking. A quasi-two-dimensional flow-through sediment microcosm/tank system enabled us to study the aerobic degradation of 2,6-dichlorobenzamide (BAM), while collecting sufficient samples at the outlet for compound-specific isotope analysis. By feeding an anoxic BAM solution through the center inlet port and dissolved oxygen (DO) above and below, strong transverse concentration cross-gradients of BAM and DO yielded zones of low (μg L^-1) steady-state concentrations. We were able to simulate the profiles of concentrations and isotope ratios of the contaminant plume using a reactive transport model that accounted for a mass-transfer limitation into bacterial cells, where apparent isotope enrichment factors *ε decreased strongly below concentrations around 600 μg/L BAM. For the biodegradation of organic micropollutants, mass transfer into the cell emerges as a bottleneck, specifically at low (μg L^-1) concentrations. Neglecting this effect when interpreting isotope ratios at field sites may lead to a significant underestimation of biodegradation.

Process-based modeling of electrokinetic-enhanced bioremediation of chlorinated ethenes

This study presents a process-based modeling analysis of electrokinetic-enhanced bioremediation (EK-Bio) to illuminate the complex interactions between physical, electrostatic and biogeochemical processes occurring during the application of this remediation technique. The features of the proposed model include: (i) multidimensional electrokinetic transport in saturated porous media by electromigration and electroosmosis, (ii) charge interactions, (iii) degradation kinetics, (iv) microbial populations dynamics of indigenous and specialized exogenous degraders, (v) mass transfer limitations, and (vi) geochemical reactions. A scenario modeling investigation is presented, which was inspired by an EK-Bio pilot application conducted in a clayey aquitard at the Skuldelev site (Denmark) contaminated by chlorinated ethenes. Lactate and specialized degraders are delivered under conservative and reactive transport conditions. In the considered setup, transport of lactate using electrokinetics results in more than fourfold increase in the distribution efficiency with respect to a diffusion-only scenario. Moreover, EK transport by electromigration and electroosmosis yields fluxes at least two orders of magnitude larger than diffusive fluxes. Quantitative metrics are also defined and used to assess the amendment distribution and the enhanced contaminant biodegradation in the different conservative and reactive transport scenarios.

Effects of local transverse dispersion on macro-scale coefficients of oxygen-limited biodegradation in a stratified formation

The correct characterization of macro-scale contaminant transport and transformation rates is an important issue for modeling reactive transport in heterogeneous aquifers. While previous studies have investigated field-scale heterogeneity of transport and biochemical properties, the effects of local transverse dispersion on macro-scale transport and transformation rates have not been well understood. In this paper, the process of oxygen-limited biodegradation in a stratified aquifer is analysed by spectral perturbation approach, and longitudinal macrodispersivity, effective biodegradation rate, effective retardation factor and effective velocity are derived for the coupled transport equations of a system consisting of a contaminant and an oxidizing agent (oxygen). The effects of local transverse dispersion on these macro-scale coefficients are studied. It is shown that local transverse dispersion can smooth the heterogeneity in biodegradation and sorption processes and enlarge effective biodegradation rate and retardation factor. The local transverse dispersion can also limit the effects of heterogeneity in biodegradation process on longitudinal macrodispersivities and effective velocities for the contaminant and dissolved oxygen. But the effects of heterogeneity in sorption process on the contaminant macrodispersivity is likely to be magnified by local transverse dispersion.

A modeling approach integrating microbial activity, mass transfer, and geochemical processes to interpret biological assays: An example for PCE degradation in a multi-phase batch setup

The rate at which organic contaminants can be degraded in aquatic environments is not only dependent upon specific degrading bacteria, but also upon the composition of the microbial community, mass transfer of the contaminant, and abiotic processes that occur in the environment. In this study, we present three-phase batch experiments of tetrachloroethene (PCE) degradation by a consortium of organohalide-respiring bacteria, cultivated alone or in communities with iron- and/or sulfate-reducers. We developed a modeling approach to quantitatively evaluate the experimental results, comprised of chemical and biomolecular time series data. The model utilizes the IPhreeqc module to couple multi-phase mass transfer between gaseous, organic and aqueous phases with microbial and aquatic geochemical processes described using the geochemical code PHREEQC. The proposed approach is able to capture the contaminant degradation, the microbial population dynamics, the effects of multi-phase kinetic mass transfer and sample removal, and the geochemical reactions occurring in the aqueous phase. The model demonstrates the importance of aqueous speciation and abiotic reactions on the bioavailability of the substrates. The model-based interpretation allowed us to quantify the reaction kinetics of the different bacterial guilds. The model further revealed that the inclusion of sulfate-reducing bacteria lowers the rate of PCE degradation and that this effect is moderated in the presence of iron-reducing bacteria.

The fate of DNAPL contaminants in non-consolidated subsurface systems - Discussion on the relevance of effective source zone geometries for plume propagation

Dense non-aqueous phase liquids, i.e., DNAPLs and the evolving contaminant plumes in aquifers provide significant potential to pose hazards affecting both environment and human health. Therefore, a proper assessment of contaminant spreading within the subsurface is critical. This includes a sufficient characterization of governing parameters describing both the subsurface and the contaminant itself. Thereby, knowledge on the contaminant source zone and especially the source zone geometry, i.e., SZG is critically required, yet very uncertain. This study identifies current limitations and open research questions in the formation and shape determination of source zone geometry, as well as its relevance for contaminant plumes. Our literature review reveals that existing characterization methods are subject to data interpretation uncertainties, while the application of these methods on field scale is limited by technical demands and accompanied efforts. In a next step, methods to implement increased source zone information into calculation methods are discussed. By means of an exemplary application of selected assessment tools, i.e., plume response models, results clearly proof the relevance of SZGs for site assessment. However, existing plume response models consider over-simplified geometries that may compromise their suitability. Our findings identify the demand for improved characterization of complex SZGs and the need to better evaluate the dependency of DNAPL migration on system properties and external influences. With emphasized knowledge on the most relevant SZG features, the delineation of "effective" SZGs allowing for straightforward implementation into plume response models and an adaption of the latter to incorporate more information on SZGs should be possible.

Recent Advances in Experimental Studies of Steady-State Dilution and Reactive Mixing in Saturated Porous Media

Transverse dispersive mixing plays an important role in controlling natural attenuation of contaminant plumes and the performance of engineered remediation strategies. The extent of transverse mixing can be significantly affected by porous media heterogeneity and anisotropy. For instance, flow focusing in the high-permeability inclusions leads to an enhancement of dilution and reactive mixing in steady-state solute transport. Numerous modeling studies have been performed to understand the mechanism of conservative and reactive transport in homogeneous and complex heterogeneous porous media. However, experimental investigations are necessary to show an intuitive phenomenon and to validate the modeling results. This paper briefly reviews recent laboratory experimental studies on dilution and reactive mixing of steady-state transport in saturated homogeneous and heterogeneous porous media. In this context, setups and measuring techniques are described in pore-scale and Darcy-scale experiments. Parameters quantifying dilution and reactive mixing in the experiments are also introduced. Finally, we discuss the further experimental works necessary to deepen our understanding of dilution and reactive mixing in natural aquifers.

Applying the Rayleigh Approach for Stable Isotope-Based Analysis of VOC Biodegradation in Diffusion-Dominated Systems

Compound-specific stable isotope analysis (CSIA) has become an established tool for assessing biodegradation in the subsurface. Diffusion-dominated vapor phase transport thereby is often excluded from quantitative assessments due to the problem of diffusive mixing of concentrations with different isotopic signatures for CSIA interpretation. In soils and other unsaturated porous media volatile organic compounds (VOCs) however, are mainly transported via gas-phase diffusion and may thus prohibit a CSIA-based quantitative assessment of the fate of VOCs. The present study presents and verifies a concept for the assessment of biodegradation-induced stable isotope fractionation along a diffusive transport path of VOCs in unsaturated porous media. For this purpose data from batch and column toluene biodegradation experiments in unsaturated porous media were combined with numerical reactive transport simulations; both addressing changes of concentration and stable isotope fractionation of toluene. The numerical simulations are in good agreement with the experiment data, and our results show that the presented analytically derived assessment concept allows using the slope of the Rayleigh plot to obtain reasonable estimates of effective in situ fractionation factors in spite of diffusion-dominated transport. This enlarges the application range of CSIA and provides a mean for a better understanding of VOC fate in the unsaturated subsurface.

Contaminant concentration versus flow velocity: drivers of biodegradation and microbial growth in groundwater model systems

Aromatic hydrocarbons belong to the most abundant contaminants in groundwater systems. They can serve as carbon and energy source for a multitude of indigenous microorganisms. Predictions of contaminant biodegradation and microbial growth in contaminated aquifers are often vague because the parameters of microbial activity in the mathematical models used for predictions are typically derived from batch experiments, which don’t represent conditions in the field. In order to improve our understanding of key drivers of natural attenuation and the accuracy of predictive models, we conducted comparative experiments in batch and sediment flow-through systems with varying concentrations of contaminant in the inflow and flow velocities applying the aerobic Pseudomonas putida strain F1 and the denitrifying Aromatoleum aromaticum strain EbN1. We followed toluene degradation and bacterial growth by measuring toluene and oxygen concentrations and by direct cell counts. In the sediment columns, the total amount of toluene degraded by P. putida F1 increased with increasing source concentration and flow velocity, while toluene removal efficiency gradually decreased. Results point at mass transfer limitation being an important process controlling toluene biodegradation that cannot be assessed with batch experiments. We also observed a decrease in the maximum specific growth rate with increasing source concentration and flow velocity. At low toluene concentrations, the efficiencies in carbon assimilation within the flow-through systems exceeded those in the batch systems. In all column experiments the number of attached cells plateaued after an initial growth phase indicating a specific “carrying capacity” depending on contaminant concentration and flow velocity. Moreover, in all cases, cells attached to the sediment dominated over those in suspension, and toluene degradation was performed practically by attached cells only. The observed effects of varying contaminant inflow concentration and flow velocity on biodegradation could be captured by a reactive-transport model. By monitoring both attached and suspended cells we could quantify the release of new-grown cells from the sediments to the mobile aqueous phase. Studying flow velocity and contaminant concentrations as key drivers of contaminant transformation in sediment flow-through microcosms improves our system understanding and eventually the prediction of microbial biodegradation at contaminated sites.

High Resolution Monitoring Above and Below the Groundwater Table Uncovers Small-Scale Hydrochemical Gradients

Hydrochemical solute concentrations in the shallow subsurface can be spatially highly variable within small scales, particularly at interfaces. However, most monitoring systems fail to capture these small scale variations. Within this study, we developed a high resolution multilevel well (HR-MLW) with which we monitored water across the interface of the unsaturated and saturated zone with a vertical resolution of 0.05-0.5 m. We installed three of these 4 m deep HR-MLWs in the riparian zone of a third-order river and analyzed for hydrochemical parameters and stable water isotopes. The results showed three distinct vertical zones (unsaturated zone, upper saturated zone, lower saturated zone) within the alluvial aquifer. A 2 m thick layer influenced by river water (upper saturated zone) was not captured by existing monitoring wells with higher screen length. Hydrochemical data (isotopes, total ions) were consistent in all HR-MLWs and showed similar variation over time emphasizing the reliability of the installed monitoring system. Further, the depths zones were also reflected in the NO₃-N concentrations; with high spatial variabilities between the three wells. The zonation was constant over time, with seasonal variability in the upper saturated zone due to the influence of river water. This study highlights the use of high resolution monitoring for identifying the spatial and temporal variability of hydrochemical parameters present in many aquifer systems. Possible applications range from riparian zones, agricultural field sites to contaminated site studies, wherever an improved understanding of biogeochemical turnover processes is necessary.

Joint interpretation of enantiomer and stable isotope fractionation for chiral pesticides degradation

Chiral pesticides are important contaminants affecting the health and functioning of aquatic systems. The combination of stable isotope and enantiomer analysis techniques has been recently proposed to better characterize the fate of these contaminants in natural and engineered settings. We introduce a modeling approach with the aim of unifying and integrating the interpretation of isotopic and enantiomeric fractionation. The model is based on the definition of enantiomer-specific isotopologues and jointly predicts the evolution of concentration, enantiomer fractionation, as well as changes in stable isotope ratios of different elements. The method allows evaluating different transformation pathways and was applied to investigate enzymatic degradation of dichlorprop (DCPP), enzymatic degradation of mecoprop methyl ester (MCPPM), and microbial degradation of α-hexachlorocyclohexane (α-HCH) by different bacterial strains and under different redox conditions. The model accurately reproduces the isotopic and enantiomeric data observed in previous experimental studies and precisely captures the dual-dimensional trends characterizing different reaction pathways. Furthermore, the model allows testing possible combinations of enantiomer analysis (EA), compound specific isotope analysis (CSIA), and enantiomer specific isotope analysis (ESIA) to identify and assess isotope and enantiomer selective reaction mechanisms.

Hydrogen Isotope Fractionation As a Tool to Identify Aerobic and Anaerobic PAH Biodegradation

Aerobic and anaerobic polycyclic aromatic hydrocarbon (PAH) biodegradation was characterized by compound specific stable isotope analysis (CSIA) of the carbon and hydrogen isotope effects of the enzymatic reactions initiating specific degradation pathways, using naphthalene and 2-methylnaphtalene as model compounds. Aerobic activation of naphthalene and 2-methylnaphthalene by Pseudomonas putida NCIB 9816 and Pseudomonas fluorescens ATCC 17483 containing naphthalene dioxygenases was associated with moderate carbon isotope fractionation (εC = -0.8 ± 0.1‰ to -1.6 ± 0.2‰). In contrast, anaerobic activation of naphthalene by a carboxylation-like mechanism by strain NaphS6 was linked to negligible carbon isotope fractionation (εC = -0.2 ± 0.2‰ to -0.4 ± 0.3‰). Notably, anaerobic activation of naphthalene by strain NaphS6 exhibited a normal hydrogen isotope fractionation (εH = -11 ± 2‰ to -47 ± 4‰), whereas an inverse hydrogen isotope fractionation was observed for the aerobic strains (εH = +15 ± 2‰ to +71 ± 6‰). Additionally, isotope fractionation of NaphS6 was determined in an overlaying hydrophobic carrier phase, resulting in more reliable enrichment factors compared to immobilizing the PAHs on the bottle walls without carrier phase. The observed differences especially in hydrogen fractionation might be used to differentiate between aerobic and anaerobic naphthalene and 2-methylnaphthalene biodegradation pathways at PAH-contaminated field sites.

PHT3D-UZF: A Reactive Transport Model for Variably-Saturated Porous Media

A modified version of the MODFLOW/MT3DMS-based reactive transport model PHT3D was developed to extend current reactive transport capabilities to the variably-saturated component of the subsurface system and incorporate diffusive reactive transport of gaseous species. Referred to as PHT3D-UZF, this code incorporates flux terms calculated by MODFLOW's unsaturated-zone flow (UZF1) package. A volume-averaged approach similar to the method used in UZF-MT3DMS was adopted. The PHREEQC-based computation of chemical processes within PHT3D-UZF in combination with the analytical solution method of UZF1 allows for comprehensive reactive transport investigations (i.e., biogeochemical transformations) that jointly involve saturated and unsaturated zone processes. Intended for regional-scale applications, UZF1 simulates downward-only flux within the unsaturated zone. The model was tested by comparing simulation results with those of existing numerical models. The comparison was performed for several benchmark problems that cover a range of important hydrological and reactive transport processes. A 2D simulation scenario was defined to illustrate the geochemical evolution following dewatering in a sandy acid sulfate soil environment. Other potential applications include the simulation of biogeochemical processes in variably-saturated systems that track the transport and fate of agricultural pollutants, nutrients, natural and xenobiotic organic compounds and micropollutants such as pharmaceuticals, as well as the evolution of isotope patterns.

A Post-Genomic View of the Ecophysiology, Catabolism and Biotechnological Relevance of Sulphate-Reducing Prokaryotes

Dissimilatory sulphate reduction is the unifying and defining trait of sulphate-reducing prokaryotes (SRP). In their predominant habitats, sulphate-rich marine sediments, SRP have long been recognized to be major players in the carbon and sulphur cycles. Other, more recently appreciated, ecophysiological roles include activity in the deep biosphere, symbiotic relations, syntrophic associations, human microbiome/health and long-distance electron transfer. SRP include a high diversity of organisms, with large nutritional versatility and broad metabolic capacities, including anaerobic degradation of aromatic compounds and hydrocarbons. Elucidation of novel catabolic capacities as well as progress in the understanding of metabolic and regulatory networks, energy metabolism, evolutionary processes and adaptation to changing environmental conditions has greatly benefited from genomics, functional OMICS approaches and advances in genetic accessibility and biochemical studies. Important biotechnological roles of SRP range from (i) wastewater and off gas treatment, (ii) bioremediation of metals and hydrocarbons and (iii) bioelectrochemistry, to undesired impacts such as (iv) souring in oil reservoirs and other environments, and (v) corrosion of iron and concrete. Here we review recent advances in our understanding of SRPs focusing mainly on works published after 2000. The wealth of publications in this period, covering many diverse areas, is a testimony to the large environmental, biogeochemical and technological relevance of these organisms and how much the field has progressed in these years, although many important questions and applications remain to be explored.

Dynamics of suspended and attached aerobic toluene degraders in small-scale flow-through sediment systems under growth and starvation conditions

The microbially mediated reactions, that are responsible for field-scale natural attenuation of organic pollutants, are governed by the concurrent presence of a degrading microbial community, suitable energy and carbon sources, electron acceptors, as well as nutrients. The temporal lack of one of these essential components for microbial activity, arising from transient environmental conditions, might potentially impair in situ biodegradation. This study presents results of small scale flow-through experiments aimed at ascertaining the effects of substrate-starvation periods on the aerobic degradation of toluene by Pseudomonas putida F1. During the course of the experiments, concentrations of attached and mobile bacteria, as well as toluene and oxygen were monitored. Results from a fitted reactive-transport model, along with the observed profiles, show the ability of attached cells to survive substrate-starvation periods of up to four months and suggest a highly dynamic exchange between attached and mobile cells under growth conditions and negligible cell detachment under substrate-starvation conditions. Upon reinstatement of toluene, it was readily degraded without a significant lag period, even after a starvation period of 130 days. Our experimental and modeling results strongly suggest that aerobic biodegradation of BTEX-hydrocarbons at contaminated field sites is not hampered by intermittent starvation periods of up to four months.

Evolution of carbon isotope signatures during reactive transport of hydrocarbons in heterogeneous aquifers

Compound-specific isotope analysis (CSIA) of organic pollutants has become a well-established tool for assessing the occurrence and extent of biodegradation processes in contaminated aquifers. However, the precision of CSIA is influenced by the degree to which assumptions underlying CSIA data interpretation hold under realistic field-scale conditions. For the first time this study demonstrates how aquifer analogs combined with reactive transport models offer an underexplored way to develop generic process understanding, evaluate monitoring and quantification strategies in highly heterogeneous subsurface settings. Data from high-resolution aquifer analogs were used in numerical experiments to track the propagation of a representative oxidizable organic compound (toluene) within a variety of realistic heterogeneous aquifers and to investigate its detailed fate. The simulations were used to analyze (1) the effects of physical aquifer heterogeneities on spatiotemporal patterns of contaminant concentrations and isotope signatures, (2) the performance of the commonly applied Rayleigh equation and (3) the applicability of an extension of the Rayleigh equation for complex hydrogeological conditions. The results indicate that if field-derived enrichment factors are applied without corrections for dilution, the conventional Rayleigh equation is inaccurate and estimates for biodegradation are typically overestimated and unreliable in heterogeneous aquifers. Underestimations can occur due to the partial source zone depletion. In contrast, if dilution can be accurately accounted for, field-derived enrichment factors comprise a suitable alternative to laboratory-derived and redox-specific enrichment factors. The study also examines to what extent variations in monitoring/sampling strategies influence the obtained results. Especially measurements from long-screened wells (>1 m) reveal to be inappropriate for the application of the Rayleigh equation in the investigated aquifer analogs, as low resolution data sampled from the simulated scenarios only enable a qualitative assessment of biodegradation. Measurements from both long- and short-screened wells employing the Rayleigh equation streamline approach are only partly viable for in situ biodegradation measurements in heterogeneous systems.

Highly parameterized inversion of groundwater reactive transport for a complex field site

In this study a numerical groundwater reactive transport model of a shallow groundwater aquifer contaminated with volatile organic compounds is developed. In addition to advective-dispersive transport, the model includes contaminant release from source areas, natural attenuation, abiotic degradation by a permeable reactive barrier at the site, and dilution by infiltrating rain. Aquifer heterogeneity is parameterized using pilot points for hydraulic conductivity, specific yield and groundwater recharge. A methodology is developed and applied to estimate the large number of parameters from the limited data at the field site (groundwater levels, groundwater concentrations of multiple chemical species, point-scale measurements of soil hydraulic conductivity, and lab-scale derived information on chemical and biochemical reactions). The proposed methodology relies on pilot point parameterization of hydraulic parameters and groundwater recharge, a regularization procedure to reconcile the large number of spatially distributed model parameters with the limited field data, a step-wise approach for integrating the different data sets into the model, and high performance computing. The methodology was proven to be effective in reproducing multiple contaminant plumes and in reducing the prior parameter uncertainty of hydraulic conductivity and groundwater recharge. Our results further indicate that contaminant transport predictions are strongly affected by the choice of the groundwater recharge model and flow parameters should be identified using both head and concentration measurements.

Fringe-controlled biodegradation under dynamic conditions: quasi 2-D flow-through experiments and reactive-transport modeling

Biodegradation in contaminated aquifers has been shown to be most pronounced at the fringe of contaminant plumes, where mixing of contaminated water and ambient groundwater, containing dissolved electron acceptors, stimulates microbial activity. While physical mixing of contaminant and electron acceptor by transverse dispersion has been shown to be the major bottleneck for biodegradation in steady-state plumes, so far little is known on the effect of flow and transport dynamics (caused, e.g., by a seasonally fluctuating groundwater table) on biodegradation in these systems. Towards this end we performed experiments in quasi-two-dimensional flow-through microcosms on aerobic toluene degradation by Pseudomonas putida F1. Plume dynamics were simulated by vertical alteration of the toluene plume position and experimental results were analyzed by reactive-transport modeling. We found that, even after disappearance of the toluene plume for two weeks, the majority of microorganisms stayed attached to the sediment and regained their full biodegradation potential within two days after reappearance of the toluene plume. Our results underline that besides microbial growth, also maintenance and dormancy are important processes that affect biodegradation performance under transient environmental conditions and therefore deserve increased consideration in future reactive-transport modeling.

Experimental investigation of compound-specific dilution of solute plumes in saturated porous media: 2-D vs. 3-D flow-through systems

Dilution of solute plumes in groundwater strongly depends on transverse mixing. Thus, the correct parameterization of transverse dispersion is of critical importance for the quantitative description of solute transport. In this study we perform flow-through laboratory experiments to investigate the influence of transport dimensionality on transverse mixing. We present a high-resolution experimental setup to study solute dilution and transverse dispersion in three-dimensional porous media. We conduct multi-tracer experiments in the new 3-D setup and compare the results with the outcomes of analogous tracer experiments performed in a quasi 2-D system. We work under steady-state flow and transport conditions and consider a range of velocities relevant for groundwater flow (0.5-8 m/day). Transverse dispersion coefficients are determined from high-resolution concentration profiles at the outlet of the flow-through chambers (7×7 ports in the 3-D setup and 7 ports in the quasi 2-D system), considering conservative tracers with significantly different aqueous diffusion coefficients, namely fluorescein and dissolved oxygen. To quantify dilution in the 2-D and 3-D systems, we experimentally determine the flux-related dilution index using the flow rates and the concentrations measured at the inlet and outlet ports, and we propose semi-analytical expressions to predict its evolution with travel distance in uniform groundwater flow. The experimental results in the quasi 2-D and 3-D flow-through systems are consistent and show a compound-specific behavior of the transverse dispersion coefficient and its non-linear dependence on the seepage velocity in both setups. The degree of dilution and the compound-specific effects of transverse dispersion are considerably more pronounced in 3-D than in quasi 2-D transport systems.

Evaluation of the role of heterogeneities on transverse mixing in bench-scale tank experiments by numerical modeling

In this work, numerical modeling is used to evaluate and interpret a series of detailed and well-controlled two-dimensional bench-scale conservative tracer tank experiments performed to investigate transverse mixing in porous media. The porous medium used consists of a fine matrix and a more permeable lens vertically aligned with the tracer source and the flow direction. A sensitivity analysis shows that the tracer distribution after passing the lens is only slightly sensitive to variations in transverse dispersivity, but strongly sensitive to the contrast of hydraulic conductivities. A unique parameter set could be calibrated to closely fit the experimental observations. On the basis of calibrated and validated model, synthetic experiments with different contrasts in hydraulic conductivity and more complex setups were performed and the efficiency of mixing evaluated. Flux-related dilution indices derived from these simulations show that the contrasts in hydraulic conductivity between matrix and high-permeable lenses as well as the spatial configuration of tracer plumes and lenses dominate mixing, rather than the actual pore scale dispersivities. These results indicate that local material distributions, the magnitude of permeability contrasts, and their spatial and scale relation to solute plumes are more important for macro-scale transverse dispersion than the micro-scale dispersivities of individual materials. Local material characterization by thorough site investigation hence is of utmost importance for the evaluation of mixing-influenced or -governed problems in groundwater, such as tracer test evaluation or an assessment of contaminant natural attenuation.

On the importance of diffusion and compound-specific mixing for groundwater transport: an investigation from pore to field scale

Mixing processes significantly affect and limit contaminant transport and transformation rates in the subsurface. The correct quantification of mixing in groundwater systems must account for diffusion, local-scale dispersion and the flow variability in heterogeneous flow fields (e.g., flow-focusing in high-conductivity and de-focusing in low-conductivity zones). Recent results of multitracer laboratory experiments revealed the significant effect of compound-specific diffusive properties on the physical displacement of dissolved species across a representative range of groundwater flow velocities. The goal of this study is to investigate the role of diffusion and compound-specific mixing for solute transport across a range of scales including: (i) pore-scale (~10⁻² m), (ii) laboratory bench-scale (~10⁰ m) and (iii) field-scale (~10² m). We investigate both conservative and mixing-controlled reactive transport using pore-scale modeling, flow-through laboratory experiments and simulations, and field-scale numerical modeling of complex heterogeneous hydraulic conductivity fields with statistical properties similar to the ones reported for the extensively investigated Borden aquifer (Ontario, Canada) and Columbus aquifer (Mississippi, USA, also known as MADE site). We consider different steady-state and transient transport scenarios. For the conservative cases we use as a metric of mixing the exponential of the Shannon entropy to quantify solute dilution either in a given volume (dilution index) or in a given solute flux (flux-related dilution index). The decrease in the mass and the mass-flux of the contaminant plumes is evaluated to quantify reactive mixing. The results show that diffusive processes, occurring at the small-scale of a pore channel, strongly affect conservative and reactive solute transport at larger macroscopic scales. The outcomes of our study illustrate the need to consider and properly account for compound-specific diffusion and mixing limitations in order to accurately describe and predict conservative and reactive transport in porous media.

Model-based integration and analysis of biogeochemical and isotopic dynamics in a nitrate-polluted pyritic aquifer

Leaching of nitrate from agricultural land to groundwater and the resulting nitrate pollution are a major environmental problem worldwide. Its impact is often mitigated in aquifers hosting sufficiently reactive reductants that can promote autotrophic denitrification. In the case of pyrite acting as reductant, however, denitrification is associated with the release of sulfate and often also with the mobilization of trace metals (e.g., arsenic). In this study, reactive transport modeling was used to reconstruct, quantify and analyze the dynamics of the dominant biogeochemical processes in a nitrate-polluted pyrite-containing aquifer and its evolution over the last 50 years in response to changing agricultural practices. Model simulations were constrained by measured concentration depth profiles. Measured (3)H/(3)He profiles were used to support the calibration of flow and conservative transport processes, while the comparison of simulated and measured sulfur isotope signatures acted as additional calibration constraint for the reactive processes affecting sulfur cycling. The model illustrates that denitrification largely prevented an elevated discharge of nitrate to surface waters, while sulfate discharges were significantly increased, peaking around 15 years after the maximum nitrogen inputs.

3D-CSIA: carbon, chlorine, and hydrogen isotope fractionation in transformation of TCE to ethene by a Dehalococcoides culture

Carbon (C), chlorine (Cl), and hydrogen (H) isotope effects were determined during dechlorination of TCE to ethene by a mixed Dehalococcoides (Dhc) culture. The C isotope effects for the dechlorination steps were consistent with data published in the past for reductive dechlorination (RD) by Dhc. The Cl effects (combined with an inverse H effect in TCE) suggested that dechlorination proceeded through nucleophilic reactions with cobalamin rather than by an electron transfer mechanism. Depletions of (37)Cl in daughter compounds, resulting from fractionation at positions away from the dechlorination center (secondary isotope effects), further support the nucleophilic dechlorination mechanism. Determination of C and Cl isotope ratios of the reactants and products in the reductive dechlorination chain offers a potential tool for differentiation of Dhc activity from alternative transformation mechanisms (e.g., aerobic degradation and reductive dechlorination proceeding via outer sphere mechanisms), in studies of in situ attenuation of chlorinated ethenes. Hydrogenation of the reaction products (DCE, VC, and ethene) showed a major preference for the (1)H isotope. Detection of depleted dechlorination products could provide a line of evidence in discrimination between alternative sources of TCE (e.g., evolution from DNAPL sources or from conversion of PCE).

Model complexity needed for quantitative analysis of high resolution isotope and concentration data from a toluene-pulse experiment

Separating microbial- and physical-induced effects on the isotope signals of contaminants has been identified as a challenge in interpreting compound-specific isotope data. In contrast to simple analytical tools, such as the Rayleigh equation, reactive-transport models can account for complex interactions of different fractionating processes. The question arises how complex such models must be to reproduce the data while the model parameters remain identifiable. In this study, we reanalyze the high-resolution data set of toluene concentration and toluene-specific δ(13)C from the toluene-pulse experiment performed by Qiu et al. (this issue). We apply five reactive-transport models, differing in their degree of complexity. We uniquely quantify degradation and sorption properties of the system for each model, estimate the contributions of biodegradation-induced, sorption-induced, and transverse-dispersion-induced isotope fractionation to the overall isotope signal, and investigate the error introduced in the interpretation of the data when individual processes are neglected. Our results show that highly resolved data of both concentration and isotope ratios are needed for unique process identification facilitating reliable model calibration. Combined analysis of these highly resolved data demands reactive transport models accounting for nonlinear degradation kinetics and isotope fractionation by both reactive and physical processes such as sorption and transverse dispersion.

Numerical simulation of isotope fractionation in steady-state bioreactive transport controlled by transverse mixing

Compound-specific stable isotope analysis (CSIA) has increasingly been used as a tool to assess intrinsic biodegradation at contaminated field sites. Typically, the Rayleigh equation is used to estimate the extent of biodegradation from measured isotope ratios of the contaminant. However, if the rate-limiting step in overall degradation is not the microbial reaction itself, the Rayleigh equation may no more be applicable. In this study we simulate biodegradation of continuously emitted petroleum hydrocarbons in groundwater systems. These contaminants are effectively degraded at the plume fringe where transverse dispersion makes them mix with dissolved electron acceptors present in the ambient groundwater. We simulate reactive transport to study the coupled effects of transverse mixing and biodegradation on the spatial patterns of carbon isotope signatures and their interpretation based on depth-integrated sampling which represents the most common setup in the assessment of contaminated sites. We present scenarios mimicking a hydraulically uniform laboratory experiment and a field-scale application considering heterogeneous conductivity fields. We compare cases in which isotopologue-specific transverse dispersion is considered to those where this additional fractionation process is neglected. We show that these effects cause significant shifts in the isotopic signals and may lead to overestimation of biodegradation. Moreover, our results provide evidence that the rate-limiting effect of transverse mixing on the overall degradation spatially varies along the length of a steady-state contaminant plume. The control of biodegradation by transverse dispersion and the fractionating effect of dispersion modulate the fractionation caused by the microbial reaction alone. As a consequence, significantly nonlinear isotopic patterns are observed in a Rayleigh plot. Simulations in heterogeneous flow domains show that these effects persist at larger field scales and are sensitive to the degree of mixing enhancement, determined by the heterogeneity of the hydraulic conductivity fields, and to the groundwater flow velocity.

Stochastic evaluation of mixing-controlled steady-state plume lengths in two-dimensional heterogeneous domains

We study plumes originating from continuous sources that require a dissolved reaction partner for their degradation. The length of such plumes is typically controlled by transverse mixing. While analytical expressions have been derived for homogeneous flow fields, incomplete characterization of the hydraulic conductivity field causes uncertainty in predicting plume lengths in heterogeneous domains. In this context, we analyze the effects of three sources of uncertainty: (i) The uncertainty of the effective mixing rate along the plume fringes due to spatially varying flow focusing, (ii) the uncertainty of the volumetric discharge through (and thus total mass flux leaving) the source area, and (iii) different parameterizations of the Darcy-scale transverse dispersion coefficient. The first two are directly related to heterogeneity of hydraulic conductivity. In this paper, we derive semi-analytical expressions for the probability distribution of plume lengths at different levels of complexity. The results are compared to numerical Monte Carlo simulations. Uncertainties in mixing and in the source strength result in a statistical distribution of possible plume lengths. For unconditional random hydraulic conductivity fields, plume lengths may vary by more than one order of magnitude even for moderate degrees of heterogeneity. Our results show that the uncertainty of volumetric flux through the source is the most relevant contribution to the variance of the plume length. The choice of different parameterizations for the local dispersion coefficient leads to differences in the mean estimated plume length.

Transverse hydrodynamic dispersion effects on isotope signals in groundwater chlorinated solvents' plumes

The effects of transverse hydrodynamic dispersion on altering transformation-induced compound-specific isotope analysis (CSIA) signals within groundwater pollution plumes have been assessed with reactive transport modeling accommodating diffusion-induced isotope fractionation (DIF) and implementing different parameterizations of local transverse dispersion. The model reproduced previously published field data showing a negative carbon isotope pattern (-2 ‰) at the fringes of a nondegrading PCE plume. We extended the study to reactive transport scenarios considering vinyl chloride as a model compound and assessing, through a detailed sensitivity analysis, the coupled effects of transverse hydrodynamic dispersion (with and without DIF) and aerobic fringe degradation on the evolution of carbon and chloride isotope ratios. Transformation-induced positive isotope signals were increasingly attenuated with distance from the source and higher degradation rate. The effect of DIF on the overall isotope signal attenuation was greatest near the source and for low values of groundwater flow velocity, transverse dispersion coefficient, molecular weight, rate constant, and isotope fractionation factor, α, of the degradation reaction. Models disregarding DIF underestimate the actual α. The approximately twice larger DIF effect for chlorine than for carbon together with the low α for oxidation resulted in strong chlorine CSIA depletions for VC at the plume fringe.

Field applicability of Compound-Specific Isotope Analysis (CSIA) for characterization and quantification of in situ contaminant degradation in aquifers

Microbial processes govern the fate of organic contaminants in aquifers to a major extent. Therefore, the evaluation of in situ biodegradation is essential for the implementation of Natural Attenuation (NA) concepts in groundwater management. Laboratory degradation experiments and biogeochemical approaches are often biased and provide only indirect evidence of in situ degradation potential. Compound-Specific Isotope Analysis (CSIA) is at present among the most promising tools for assessment of the in situ contaminant degradation within aquifers. One- and two-dimensional (2D) CSIA provides qualitative and quantitative information on in situ contaminant transformation; it is applicable for proving in situ degradation and characterizing degradation conditions and reaction mechanisms. However, field application of CSIA is challenging due to a number of influencing factors, namely those affecting the observed isotope fractionation during biodegradation (e.g., non-isotope-fractionating rate-limiting steps, limited bioavailability), potential isotope effects caused by processes other than biodegradation (e.g., sorption, volatilization, diffusion), as well as non-isotope-fractionating physical processes such as dispersion and dilution. This mini-review aims at guiding practical users towards the sound interpretation of CSIA field data for the characterization of in situ contaminant degradation. It focuses on the relevance of various constraints and influencing factors in CSIA field applications and provides advice on when and how to account for these constraints. We first evaluate factors that can influence isotope fractionation during biodegradation, as well as potential isotope-fractionating and non-isotope-fractionating physical processes governing observed isotope fractionation in the field. Finally, the potentials of the CSIA approach for site characterization and the proper ways to account for various constraints are illustrated by means of a comprehensive CSIA field study at the benzene, toluene, ethylbenzene, and xylene (BTEX)-contaminated site Zeitz.

Methane oxidation in a crude oil contaminated aquifer: Delineation of aerobic reactions at the plume fringes

High resolution direct-push profiling over short vertical distances was used to investigate CH(4) attenuation in a petroleum contaminated aquifer near Bemidji, Minnesota. The contaminant plume was delineated using dissolved gases, redox sensitive components, major ions, carbon isotope ratios in CH(4) and CO(2), and the presence of methanotrophic bacteria. Sharp redox gradients were observed near the water table. Shifts in δ(13)C(CH4) from an average of -57.6‰ (±1.7‰) in the methanogenic zone to -39.6‰ (±8.7‰) at 105m downgradient, strongly suggest CH(4) attenuation through microbially mediated degradation. In the downgradient zone the aerobic/anaerobic transition is up to 0.5m below the water table suggesting that transport of O(2) across the water table is leading to aerobic degradation of CH(4) at this interface. Dissolved N(2) concentrations that exceeded those expected for water in equilibrium with the atmosphere indicated bubble entrapment followed by preferential stripping of O(2) through aerobic degradation of CH(4) or other hydrocarbons. Multivariate and cluster analysis were used to distinguish between areas of significant bubble entrapment and areas where other processes such as the infiltration of O(2) rich recharge water were important O(2) transport mechanisms.

The application of a simplified method to map the aerobic acetate mineralization rates at the groundwater table of the Netherlands

A simplified method is used to assess the microbial activity of subsoils and soils across a broad geographic scale. Acetate was selected because it is a major intermediate in catabolic biochemical pathways. In order to get minimal disturbance, only a small amount of tritium labelled acetate and water is added to the subsoil material. After an incubation time, the subsoil material is separated from the water by centrifugation and the formed tritium labelled water is separated from the remaining acetate by evaporation. The data of 128 locations in the Netherlands were plotted in a soil map and were also compared with the depth, dry weight, electric conductivity, pH and nitrate concentration. The peat areas consisted of limed meadows with a high groundwater level whereas the sand areas often showed deeper groundwater levels and a lower pH. The subsoils at the groundwater table of the peat areas, which are in contact with soil air, showed a higher mineralization rate compared with the surface soils in our study. In contrast, the mineralization rate of the subsoil at the groundwater table of sandy soils showed on average a factor 30 lower rate. Nevertheless, the self purification capacity of the subsoil can be vital under weather conditions where the surface soil becomes less active.

Metabolites indicate hot spots of biodegradation and biogeochemical gradients in a high-resolution monitoring well

Anaerobic degradation processes play an important role in contaminated aquifers. To indicate active biodegradation processes signature metabolites can be used. In this study field samples from a high-resolution multilevel well in a tar oil-contaminated, anoxic aquifer were analyzed for metabolites by liquid chromatography-tandem mass spectrometry and time-of-flight mass spectrometry. In addition to already known specific degradation products of toluene, xylenes, and naphthalenes, the seldom reported degradation products benzothiophenemethylsuccinic acid (BTMS), benzofuranmethylsuccinic acid (BFMS), methylnaphthyl-2-methylsuccinic acid (MNMS), and acenaphthene-5-carboxylic acid (AC) could be identified (BFMS, AC) and tentatively identified (BTMS, MNMS). The occurrence of BTMS and BFMS clearly show that the fumarate addition pathway, known for toluene and methylnaphthalene, is also important for the anaerobic degradation of heterocyclic contaminants in aquifers. The molar concentration ratios of metabolites and their related parent compounds differ over a wide range which shows that there is no simple and consistent quantitative relation. However, generally higher ratios were found for the more recalcitrant compounds, which are putatively cometabolically degraded (e.g., 2-carboxybenzothiophene and acenaphthene-5-carboxylic acid), indicating an accumulation of these metabolites. Vertical concentration profiles of benzylsuccinic acid (BS) and methyl-benzylsuccinic acid (MBS) showed distinct peaks at the fringes of the toluene and xylene plume indicating hot spots of biodegradation activity and supporting the plume fringe concept. However, there are some compounds which show a different vertical distribution with the most prominent concentrations where also the precursor compounds peaked.

Isotopic fractionation by transverse dispersion: flow-through microcosms and reactive transport modeling study

Flow-through experiments were carried out to investigate the role of transverse dispersion on the isotopic behavior of an organic compound during conservative and bioreactive transport in a homogeneous porous medium. Ethylbenzene was selected as model contaminant and a mixture of labeled (perdeuterated) and light isotopologues was continuously injected in a quasi two-dimensional flow-through system. We observed a significant fractionation of ethylbenzene isotopologues during conservative transport at steady state. This effect was particularly pronounced at the plume fringe and contrasted with the common assumption that physical processes only provide a negligible contribution to isotope fractionation. Under the experimental steady state conditions, transverse hydrodynamic dispersion was the only process that could have caused the observed fractionation. Therefore, the measured isotope ratios at the outlet ports were interpreted with different parameterizations of the transverse dispersion coefficient. A nonlinear compound-specific parameterization showed the best agreement with the experimental data. Successively, bioreactive experiments were performed in two subsequent stages: a first oxic phase, involving a single strain of ethylbenzene degraders and a second phase with aerobic and anaerobic (i.e., ethylbenzene oxidation coupled to nitrate reduction) degradation. Significant fractionation through biodegradation occurred exclusively due to the metabolic activity of the anaerobic degraders. We performed analytical and numerical reactive transport simulations of the different experimental phases which confirmed that both the effects of physical processes (diffusion and dispersion) and microbially mediated reactions have to be considered to match the observed isotopic fractionation behavior.

Evidence of compound-dependent hydrodynamic and mechanical transverse dispersion by multitracer laboratory experiments

Mass transfer, mixing, and therefore reaction rates during transport of solutes in porous media strongly depend on dispersion and diffusion. In particular, transverse mixing is a significant mechanism controlling natural attenuation of contaminant plumes in groundwater. The aim of the present study is to gain a deeper understanding of vertical transverse dispersive mixing of reaction partners in saturated porous media. Multitracer laboratory experiments in a quasi two-dimensional tank filled with glass beads were conducted and transverse dispersion coefficients were determined from high-resolution vertical concentration profiles. We investigated the behavior of conservative tracers (i.e., fluorescein, dissolved oxygen, and bromide), with different aqueous diffusion coefficients, in a range of grain-related Peclet numbers between 1 and 562. The experimental results do not agree with the classical linear parametric model of hydrodynamic dispersion, in which the transverse component is approximated as the sum of pore diffusion and a compound-independent mechanical dispersion term. The outcome of the multitracer experiments clearly indicates a nonlinear relation between the dispersion coefficient and the average linear velocity. More importantly, we show that transverse mechanical dispersion depends on the diffusion coefficient of the compound, at least at the experimental bench-scale. This result has to be considered in reactive-transport models, because the typical assumption that two reactants with different aqueous diffusive properties are characterized by the same dispersive behavior does not hold anymore.

Enhancement of dilution and transverse reactive mixing in porous media: experiments and model-based interpretation

Transport and natural attenuation of contaminant plumes in groundwater are often controlled by transverse dispersion. The extent of mixing between dissolved reaction partners at the fringe of a plume determines its length and depends strongly on the groundwater flow field. Transient flow conditions as well as the focusing of the flow in high-permeability zones may enhance transverse mixing of dissolved species and, therefore, create favorable conditions for the natural attenuation of contaminant plumes. The aim of the present study is to experimentally test the influence of these processes on solute mixing and to directly compare the results with those under analogous homogeneous and steady-state conditions. We have performed conservative and reactive tracer experiments in a quasi two-dimensional tank filled with glass beads of different sizes. The experiments have been carried out in both homogeneous and heterogeneous porous media under steady-state and transient (i.e. oscillating) flow fields. We used fluorescein as conservative tracer; whereas an alkaline solution (NaOH) was injected into ambient acidic water (HCl) in the reactive experiments. A pH indicator was added to the reacting solutions in order to visualize the emerging plume. We simulated the laboratory experiments with a numerical model and compared the outcomes of the model with the measured concentrations at the outlet of the tank and with the observed tracer plumes. Spatial moments, a newly defined flux-related dilution index, the product mass fluxes and the reaction enhancement factors were calculated to quantify the differences in mixing and reaction extent under various experimental conditions. The results show that flow focusing in heterogeneous porous media significantly enhances transverse mixing and mixing-controlled reactions, whereas temporally changing flow fields appear to be of minor importance.