Aerobic biodegradation of chlorinated ethenes in a fractured bedrock aquifer: quantitative assessment by compound-specific isotope analysis (CSIA) and reactive transport modeling

Sorption- and diffusion-induced isotopic fractionation in chloroethenes

Isotopic fractionation of groundwater contaminants can occur due to degradation, diffusion and sorption. Of these, only degradation has been extensively explored, yet diffusive isotopic fractionation (DIF) and sorptive isotopic fractionation (SIF) can have significant effects on the isotopic enrichment of groundwater contaminants. Understanding how to mathematically describe and model these processes is vital to the correct interpretation of compound-specific isotope analysis (CSIA) data in the field. Here, models for these physical fractionation processes are developed and described, including the definition of a sorption enrichment factor. These models are then implemented numerically using inverse and finite-element methods to investigate two scenarios, diffusion-sorption and diffusion-sorption-advection, that have been measured in the laboratory. Concentration, δ³⁷Cl, and δ²H data from cis-dichloroethene (cDCE) and trichloroethene (TCE) are used as inputs to the models. Unknown transport parameters including diffusive fractionation exponents are determined from an inverse modelling approach. DIF is shown to have a stronger influence on chlorine isotopologues than on hydrogen isotopologues. For both cDCE and TCE, the sorption enrichment factor of chlorine is found to be negative while that of hydrogen is positive. The presented approach and results provide novel tools and insight into DIF and SIF and underline that these processes should be taken into account when using CSIA to assess contaminant fate.

Dynamics of chlorinated aliphatic hydrocarbons in the Chalk aquifer of northern France

The Chalk aquifer used for drinking-water production in the southwest of the Lille European Metropolis is threatened by the presence of chlorinated aliphatic hydrocarbons (CHCs), their concentrations in groundwater regularly exceeding the regulatory limits for drinking water in France. This hinders its use for drinking-water production. Understanding the dynamics and spatial distribution of CHC in the aquifer is a key factor for resource sustainability. For that purpose, an intensive monitoring was undertaken in several well fields and at different depths over eight years. To assess a possible migration and/or degradation of the compounds, the water column in several wells was sampled at various depths with passive samplers. Furthermore, CHC degradation mechanisms were investigated with compound-specific carbon-isotope analysis. The CHC concentrations and their distributions in the area depend on past and current industrial activity, causing plumes emphasized by pumping in the wells, such plumes being multi-source with no identified origin in most wells. In the south area of Les Ansereuilles, reductive dechlorination of tetrachloroethylene from a former industrial laundry highly impacted the surrounding area with its main degradation product cis-1,2-dichloroethylene. The same area is also affected by tetrachlroroethylene from several industrial laundries, textile factories and dyeing industries with also an anaerobic degradation. In the northern part of Les Ansereuilles, tetrachloroethylene, trichloroethane, trichloroethylene and 1,1-dichloroethylene were found as primary products, whereas cis-1,2-dichloroethylene appears to be an anaerobic degradation product of TCE. The other well fields (Houplin-Ancoisne, Seclin and Emmerin) are less impacted by CHC pollution, and it was shown that no CHC degradation occurred in the wells. However, the stratification of CHCs in the well-water columns, their constant concentration values over time caused by the large amount of available CHCs, and the minor degradation occurring in wells are of concern for water operators in the future.

Modeling 3D-CSIA data: Carbon, chlorine, and hydrogen isotope fractionation during reductive dechlorination of TCE to ethene

Reactive transport modeling of multi-element, compound-specific isotope analysis (CSIA) data has great potential to quantify sequential microbial reductive dechlorination (SRD) and alternative pathways such as oxidation, in support of remediation of chlorinated solvents in groundwater. As a key step towards this goal, a model was developed that simulates simultaneous carbon, chlorine, and hydrogen isotope fractionation during SRD of trichloroethene, via cis-1,2-dichloroethene (and trans-DCE as minor pathway), and vinyl chloride to ethene, following Monod kinetics. A simple correction term for individual isotope/isotopologue rates avoided multi-element isotopologue modeling. The model was successfully validated with data from a mixed culture Dehalococcoides microcosm. Simulation of Cl-CSIA required incorporation of secondary kinetic isotope effects (SKIEs). Assuming a limited degree of intramolecular heterogeneity of δ³⁷Cl in TCE decreased the magnitudes of SKIEs required at the non-reacting Cl positions, without compromising the goodness of model fit, whereas a good fit of a model involving intramolecular CCl bond competition required an unlikely degree of intramolecular heterogeneity. Simulation of H-CSIA required SKIEs in H atoms originally present in the reacting compounds, especially for TCE, together with imprints of strongly depleted δ²H during protonation in the products. Scenario modeling illustrates the potential of H-CSIA for source apportionment.

Validation of GC-IRMS techniques for δ13C and δ2H CSIA of organophosphorus compounds and their potential for studying the mode of hydrolysis in the environment

Compound-specific stable isotope analysis (CSIA) is among the most promising tools for studying the fate of organic pollutants in the environment. However, the feasibility of multidimensional CSIA was limited by the availability of a robust method for precise isotope analysis of heteroatom-bearing organic compounds. We developed a method for δ ¹³C and δ ²H analysis of eight organophosphorus compounds (OPs) with different chemical properties. In particular, we aimed to compare high-temperature conversion (HTC) and chromium-based HTC (Cr/HTC) units to explore the limitations of hydrogen isotope analysis of heteroatom-bearing compounds. Analysis of the amount dependency of the isotope values (linearity analysis) of OPs indicated that the formation of HCl was a significant isotope fractionation process leading to inaccurate δ ²H analysis in HTC. In the case of nonchlorinated OPs, by-product formation of HCN, H₂S, or PH₃ in HTC was observed but did not affect the dynamic range of reproducible isotope values above the limit of detection. No hydrogen-containing by-products were found in the Cr/HTC process by use of ion trap mass spectrometry analysis. The accuracy of gas chromatography - isotope ratio mass spectrometry was validated in comparison with elemental analyzer - isotope ratio mass spectrometry. Dual-isotope fractionation yielded Λ values of 0 ± 0 at pH 7, 7 ± 1 at pH 9, and 30 ± 6 at pH 12, indicating the potential of 2D CSIA to characterize the hydrolysis mechanisms of OPs. This is the first report on the combination of δ ²H and δ ¹³C isotope analysis of OPs, and this is the first study providing a systematic evaluation of HTC and Cr/HTC for hydrogen isotope analysis using OPs as target compounds. Graphical Abstract Comparison of δ²H measurement of non-chlorinated and chlorinated OPs via GC-Cr/HTC-IRMS and GC-HTC-IRMS system.

Dual Mechanism Conceptual Model for Cr Isotope Fractionation during Reduction by Zerovalent Iron under Saturated Flow Conditions

Chromium isotope analysis is rapidly becoming a valuable complementary tool for tracking Cr(VI) treatment in groundwater. Evaluation of various treatment materials has demonstrated that the degree of isotope fractionation is a function of the reaction mechanism, where reduction of Cr(VI) to Cr(III) induces the largest fractionation. However, it has also been observed that uniform flow conditions can contribute complexity to isotope measurements. Here, laboratory batch and column experiments were conducted to assess Cr isotope fractionation during Cr(VI) reduction by zerovalent iron under both static and saturated flow conditions. Isotope measurements were accompanied by traditional aqueous geochemical measurements (pH, Eh, concentrations) and solid-phase analysis by scanning electron microscopy and X-ray absorption spectroscopy. Increasing δ(53)Cr values were associated with decreasing Cr(VI) concentrations, which indicates reduction; solid-phase analysis showed an accumulation of Cr(III) on the iron. Reactive transport modeling implemented a dual mechanism approach to simulate the fractionation observed in the experiments. The faster heterogeneous reaction pathway was associated with minimal fractionation (ε=-0.2‰), while the slower homogeneous pathway exhibited a greater degree of fractionation (ε=-0.9‰ for the batch experiment, and ε=-1.5‰ for the column experiment).

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.

Aerobic biodegradation of trichloroethene without auxiliary substrates

Trichloroethene (TCE) represents a priority pollutant and is among the most frequently detected contaminants in groundwater. The current bioremediation measures have certain drawbacks like e.g. the need for auxiliary substrates. Here, the aerobic biodegradation of TCE as the sole growth substrate is demonstrated. This new process of metabolic TCE degradation was first detected in groundwater samples. TCE degradation was stable in an enriched mixed bacterial culture in mineral salts medium for over five years and repeated transfers of the culture resulting in a 10(10) times dilution of the original groundwater. Aerobic TCE degradation resulted in stoichiometric chloride formation. Stable carbon isotope fractionation was observed providing a reliable analytical tool to assess this new biodegradation process at field sites. The results suggest that aerobic biodegradation of TCE without auxiliary substrate could be considered as an option for natural attenuation or engineered bioremediation of contaminated sites.

Multi-isotope (carbon and chlorine) analysis for fingerprinting and site characterization at a fractured bedrock aquifer contaminated by chlorinated ethenes

The use of compound specific multi-isotope approach (C and Cl) in the characterization of a chlorinated ethenes contaminated fractured aquifer allows the identification of several sources and contaminant plumes, as well as the occurrence of biodegradation and mixing processes. The study site is located in Spain with contamination resulting in groundwater concentrations of up to 50mg/L of trichloroethene (TCE), the most abundant chlorinated ethene, and 7 mg/L of tetrachloroethene (PCE). The potential sources of contamination including abandoned barrels, an underground tank, and a disposal lagoon, showed a wide range in δ(13)C values from -15.6 to -40.5‰ for TCE and from -18.5 to -32.4‰ for PCE, allowing the use of isotope fingerprinting for tracing of the origin and migration of these contaminants in the aquifer. In contrast, there is no difference between the δ(37)Cl values for TCE in the contaminant sources, ranging from +0.53 to +0.66‰. Variations of δ(37)Cl and δ(13)C in the different contaminant plumes were used to investigate the role of biodegradation in groundwater. Moreover, the isotopic data were incorporated into a reactive transport model for determination of whether the isotope pattern observed downstream from the tank's source could be explained by the simultaneous effect of mixing and biodegradation. The results demonstrate that a multi-isotope approach is a valuable tool for characterization of complex sites such as fractured bedrock aquifer contaminated by multiple sources, providing important information which can be used by consultants and site managers to prioritize and design more successful remediation strategies.

Carbon and nitrogen isotope analysis of atrazine and desethylatrazine at sub-microgram per liter concentrations in groundwater

Environmental degradation of organic micropollutants is difficult to monitor due to their diffuse and ubiquitous input. Current approaches-concentration measurements over time, or daughter-to-parent compound ratios-may fall short, because they do not consider dilution, compound-specific sorption characteristics or alternative degradation pathways. Compound-specific isotope analysis (CSIA) offers an alternative approach based on evidence from isotope values. Until now, however, the relatively high limits for precise isotope analysis by gas chromatography-isotope ratio mass spectrometry (GC-IRMS) have impeded CSIA of sub-microgram-per-liter scale micropollutant concentrations in field samples. This study presents the first measurements of C and N isotope ratios of the herbicide atrazine and its metabolite desethylatrazine at concentrations of 100 to 1,000 ng/L in natural groundwater samples. Solid-phase extraction and preparative HPLC were tested and validated for preconcentration and cleanup of groundwater samples of up to 10 L without bias by isotope effects. Matrix interferences after solid-phase extraction could be greatly reduced by a preparative HPLC cleanup step prior to GC-IRMS analysis. Sensitivity was increased by a factor of 6 to 8 by changing the injection method from large-volume to cold-on-column injection on the GC-IRMS system. Carbon and nitrogen isotope values of field samples showed no obvious correlation with concentrations or desethylatrazine-to-atrazine ratios. Contrary to expectations, however, δ (13) C values of desethylatrazine were consistently less negative than those of atrazine from the same sites. Potentially, this line of evidence may contain information about further desethylatrazine degradation. In such a case, the common practice of using desethylatrazine-to-atrazine ratios would underestimate natural atrazine degradation.

Reactive transport modeling of chromium isotope fractionation during Cr(VI) reduction

Chromium isotope fractionation is indicative of mass-transfer processes, such as reduction of Cr(VI) to Cr(III) during groundwater remediation. Laboratory experiments comparing batch and column treatment of Cr(VI) using organic carbon suggest that the associated isotope fractionation may be influenced by solute-transport mechanisms. These batch and column experiments were simulated using the reactive transport model MIN3P to further evaluate the effects of Cr reduction and transport on isotope fractionation under saturated flow conditions. Simulation of the batch experiment provided a good fit to the experimental data, where a fractionation factor (α₅₃) of 0.9965 was attributed to a single, dominant Cr(VI) removal mechanism. Calibration of the column simulations to the experimental results suggested the presence of a second, more rapid Cr(VI) removal mechanism with α₅₃ = 0.9992. Results from this study demonstrate that the interpretation of Cr isotope fractionation during reduction can be complex, particularly where multiple removal mechanisms are evident. Reactive transport modeling of Cr isotope fractionation can provide a quantitative assessment of the contaminant removal mechanisms, thus improving the application of Cr isotope measurements as a tool to track Cr(VI) migration and attenuation in groundwater.

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.

Origin and fate of organic compounds in water: characterization by compound-specific stable isotope analysis

Within the past 15 years, compound-specific stable isotope analysis has continued to increase in popularity in the area of contaminant hydrology of organic molecules. In particular, in cases where concentration data alone are insufficient to elucidate environmental processes unequivocally, the isotope signature can provide additional unique information. Specifically, it can help answer questions about contaminant source apportionment, quantification of biotic and abiotic processes, and identification of transformation reactions on a mechanistic level. We review advances in laboratory and field investigations and exemplary applications in contaminant hydrology via stable isotope analysis. We also highlight future directions in the field.

Effectiveness of nanoscale zero-valent iron for treatment of a PCE-DNAPL source zone

Nanoscale zero-valent iron (nZVI) has received considerable attention as a potential in situ remediation technology for treating chlorinated solvent source zones. Experimental and mathematical modeling studies were conducted to investigate the performance of nZVI in the transformation of tetrachloroethene (PCE) entrapped as a dense nonaqueous phase liquid (DNAPL). Injection of a 60 g/L suspension of nZVI into a column containing 20-30 mesh Ottawa sand and PCE-DNAPL at a residual saturation of 5.5% resulted in a uniform distribution of nZVI and minimal displacement of PCE. Subsequent flushing with 267 pore volumes of water containing 3mM CaCl(2) at a Darcy velocity of 0.75 m/day resulted in steady-state effluent concentrations of PCE near the solubility limit (ca. 200mg/L) and production of dissolved-phase ethene (10-30 mg/L). Over the duration of the experiment, approximately 30% of the initial PCE-DNAPL mass reacted to form ethene, 50% was eluted as dissolved-phase PCE, and 20% remained in the column as PCE-DNAPL. To further explore the implications of the nZVI column results, a multiphase transport model was developed that incorporated rate-limited PCE-DNAPL dissolution and reactions with nZVI. Using a fitted pseudo first-order transformation rate coefficient of 1.421/h, the model accurately captured observed trends in effluent concentrations of PCE and ethene and overall mass balance. A model sensitivity study reveals a strong dependence of treatment effectiveness on system characteristics. The sensitivity analysis suggests that an increase in the extent of PCE transformation is facilitated by decreasing flow rate, emplacement of nZVI down-gradient of the DNAPL source zone, and decreasing length of the DNAPL source zone. These findings indicate that, although emplacement of high concentrations of nZVI within a PCE-DNAPL source zone can result in substantial transformation of the parent compound, careful attention to design parameters (e.g. flow rate, location and amount nZVI delivered) will be required to achieve complete conversion to benign reaction products.

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.