Identification of Degradation Pathways of Chlorohydrocarbons in Saturated Low-Permeability Sediments Using Compound-Specific Isotope Analysis

Risk classification of contaminated sites - Comparison of the Swedish and the German method

To classify contaminated sites into different risk classes, many different methods exist in Europe and worldwide. However, no systematic comparison of European risk classification methods has been carried out so far to carve out the advantages and disadvantages of the methods and to homogenize them. To address this research gap, this study aims at comparing the Swedish Method for Inventories of Contaminated Sites (MIFO) with the German Individual Assessment of Contaminated Sites Method (EB) from the Hessian Agency for Nature Conservation, Environment and Geology (HLNUG) regarding the risk class categorization of 51 contaminated sites. The results revealed that with the MIFO 39% fewer contaminated sites are assigned to risk classes 1 and 2 and thus, subject to remediation compared to the EB. Moreover, in comparison to the EB, the MIFO showed a lower comparability, traceability, and a larger room for interpretation, which could be related to the lack of a quantitative approach such as a point or ranking system in the MIFO. Hence, we recommend providing the MIFO and other methods that lack a quantitative approach with a point and/or ranking system, similar to the EB, to increase their objectivity for the risk class categorization of contaminated sites.

Determining effective diffusion coefficients of chlorohydrocarbons in natural clays: Unique results from highly resolved controlled release field experiments

This study aims to precisely determine the effective diffusion coefficients of chlorohydrocarbons in low permeable units under in-situ field conditions. To this end, two controlled release field experiments using TCE and PCE as dense non-aqueous phase liquids (DNAPLs) were conducted in two natural clayey deposits. Several months to years after the controlled DNAPL release, highly resolved concentration profiles were determined for the chlorohydrocarbons that had diffused into the clayey deposits. Effective diffusion coefficients for TCE and PCE were then determined by calibrating a 3D numerical and 1D analytical model, respectively, to the measured high-resolution concentration profiles. The simulations revealed that the effective diffusion coefficients vary by as much as a factor of four within the same low permeability unit being consistent with observed small-scale heterogeneities. The determined chlorohydrocarbon effective diffusion coefficients were further used to determine the equivalent thickness of DNAPL that would completely dissolve in an idealized, parallel-plate fracture by diffusion transport into clayey deposits for the time periods of the controlled release field experiments. The equivalent TCE and PCE DNAPL film thicknesses ranged between 36 and 581 μm, respectively, comparable and exceeding fracture apertures measured in naturally fractured clay rich deposits. Hence, films of DNAPL initially contained within fractures in clay-rich deposits can completely dissolve away within a few months to a few years due to diffusion. This stored contaminant mass poses a risk to adjacent aquifers if it is re-released due to diffusion out of the matrix after source depletion or remediation.

Investigating the roles of advection and degradation in chlorinated solvent back-diffusion from multi-layer aquitards: A novel analytical approach

Aquitards contaminated by chlorinated solvents may act as a secondary source slowly releasing contaminants into adjacent aquifers, thus severely hampering the remediation of groundwater systems. Accurate predicting the long-term exchange of solvents between aquifers and aquitards can more effectively guide site management and remediation. This study presented a general analytical model for the back-diffusion of chlorinated solvents through multilayer aquitards. This model considers the slow advection and local degradation of dissolved constituents in natural aquitards and the dynamic depletion of dense nonaqueous phase liquid (DNAPL) source zone in aquifers. Transient solutions for the proposed multilayer model were derived using Duhamel's Theorem, the separation of variables method, and the transfer matrix method, verified against experimental and numerical concentration data. Results reveal that advection in aquitards can significantly shorten the trailing time of chlorinated solvent plumes, and highly adsorptive soils may reduce this effect in layered aquitards. The previous no-degradation model is no longer applicable to predict the back-diffusion behavior of chlorinated solvents when the extent and rate of solvent degradation are large, giving a "strong-effect zone". Based on numerous example simulations and data fitting, the forecast functions for the back-diffusion onset time and plume trailing time were proposed, greatly facilitating remediation decisions and risk assessment of chlorinated-solvent contaminated sites.

Identification of chlorohydrocarbon degradation pathways in aquitards using dual element compound-specific isotope measurements in aquifers

Compound-specific isotope analysis (CSIA) has been increasingly used to understand and quantify the (bio)degradation processes affecting chlorohydrocarbons in aquifer-aquitard systems. In this study, we aimed at investigating through reactive transport simulations if dual element (C, Cl) CSIA in aquifer samples can provide information about the occurring (bio)degradation pathways in the underlying aquitard. To that end, we modeled the continous dissolution of a 1,1,2,2-tetrachloroethane (TeCA) dense nonaqueous phase liquid (DNAPL) source in an aquifer as well as the resulting TeCA groundwater plume formation and diffusion into the underlying aquitard. The (bio)degradation of TeCA in the aquifer-aquitard system was simulated in four scenarios: TeCA biodegradation via dehydrohalogenation to trichloroethene (TCE) and TeCA dichloroelimination to dichloroethene (DCE) in the aquifer as well as in the aquitard. The simulations revealed that dual element (C, Cl) CSIA in the aquifer allows the disentanglement of whether TeCA degradation occurs in the aquifer or the aquitard and which (bio)degradation pathways occur in the aquitard. This demonstrates that chlorohydrocarbon (bio)degradation pathways in aquitards can be identified based on CSIA aquifer measurements only, which is an advantage as aquifers are easier to monitor than aquitards.

Forward and back diffusion of reactive contaminants through multi-layer low permeability sediments

Contaminants stored in the low permeability sediments will continue to threaten the adjacent shallow groundwater system after the aquifer is remediated. Understanding the storage and discharge behavior of contaminants in the aquitards is essential for the efficient remediation of contaminated sites, but most of the previous analytical studies focused on nonreactive solutes in a single homogenous aquitard. This study presents novel analytical solutions for the forward and back diffusion of contaminants through multi-layer low permeability sediments considering abiotic and biotic environmental degradation. Three representative source depletion patterns (i.e., instantaneous, linear, and exponential patterns) were selected to describe the dissolution of dense non-aqueous phase liquids (DNAPL) in the aquifer more realistically. At the forward diffusion stage, the mass storage of contaminants in the aquitards with the instantaneous pattern is the largest, nearly twice that with the exponential pattern. A simple equivalent homogeneous model is generally adopted in the risk assessment. However, relative to the proposed multi-layer model, it will significantly underestimate the onset of the back-diffusion of heterogeneous aquitards and overestimate the persistence of aquifer plumes. The previously-reported semi-infinite boundary assumption is also not applicable, with a maximum error of over 200% in the long-term prediction of back diffusion behavior of a thin aquitard. Moreover, when the degradation half-life is less than 16 years, less than 10% of the contaminants stored in the aquitards will diffuse into the overlying aquifer, suggesting that biostimulation or bioaugmentation can effectively mitigate back-diffusion risk. Overall, the proposed diffusion-reaction coupled model with multi-layer media is of great value and high demand in predicting the back-diffusion behavior of heterogeneous aquitards and guiding the soil bioremediation.

Evaluating the impact of back diffusion on groundwater cleanup time

Back diffusion of groundwater contaminants from low permeability (K) zones can be a major factor controlling the time to reach cleanup goals in downgradient monitor wells. We identify the aquifer and contaminant characteristics that have the greatest influence on the time (T_OoM) after complete source removal for contaminant concentrations to decline by 1, 2 and 3 Orders-of-Magnitude (T₁, T₂ and T₃). Two aquifer configurations are evaluated: (a) layered geometry (LG) with finite thickness low K layers; and (b) boundary geometry (BG) with thick semi-infinite low K boundaries. A semi-analytical modeling approach (Muskus and Falta, 2018) is used to simulate the concentration decline following source removal for a range of conditions and generate ≈21,000 independent values of T₁, T₂ and T₃. Linear regression is applied to interpret this large dataset and develop simple relationships to estimate T_OoM from three characteristic parameters - the mass residence time (T_M), diffusion time (T_D), and ratio of low K to high K mass storage (γ). T_M is most important predictor of T₁, T₂ and T₃ for both geometries and is equal to the combined high and low K contaminant mass divided by the mass flux, at the end of the loading period (T_L). For LG, T₃ is strongly influenced by T_D = R_LL_D²/(4D*), where R_L is the low K retardation factor, L_D is the half-thickness of the embedded low K layers, and D* is the effective diffusion coefficient. For BG, T₃ is strongly influenced by γ. Contaminant decay in low K zones can significantly reduce cleanup times when λ_LT_D > 0.01, where λ_L is the effective first order decay rate in the low K zone. The 1st Damköhler (Da), equal to T_M/T_D, provides a useful indicator of the relative importance of back diffusion on T_OoM. Back diffusion impacts are greatest on T₃ when 0.01 > Da > 0.1, then decrease with increasing Da. Back diffusion has less impacts on T₂, with limited influence on T₁. The results are summarized in a simple conceptual model to aid in evaluating the impact of back diffusion on the time for concentrations to decline by 1-3 OoM.

Compound-specific carbon isotope analysis of volatile organic compounds in complex soil extracts using purge and trap concentration coupled to heart-cutting two-dimensional gas chromatography-isotope ratio mass spectrometry

Compound-specific carbon isotope analysis (CSIA) is a powerful tool to track the origin and fate of organic subsurface contaminants including petroleum and chlorinated hydrocarbons and is typically applied to water samples. However, soil can form a significant contaminant reservoir. In soil samples, it can be challenging to recover sufficient amounts of volatile organic compounds (VOC) to perform CSIA. Soil samples often contain complex contaminant mixtures and gas chromatography combustion isotope ratio mass spectrometry (GC-C-IRMS) is highly dependent on good chromatographic separation due to the conversion to a single analyte. To extend the applicability of CSIA to complex volatile organic compound mixtures in soil samples, and to recover sufficient amounts of target compounds for carbon CSIA, we compared two soil extraction solvents, tetraglyme (TGDE) and methanol, and developed a heart-cutting two-dimensional GC-GC-C-IRMS method. We used purge & trap concentration of solvent-water mixtures to increase the amount of analyte delivered to the column and thus lower method detection limits. We optimized purge & trap and chromatographic parameters for twelve target compounds, including one suffering from poor purge efficiency. By using a 30 m thick-film non-polar column in the first and a 15 m polar column in the second dimension, we achieved good chromatographic separation for the target compounds in difficult matrices and high accuracy (trueness and precision) for carbon isotopic analysis. Tetraglyme extraction was shown to offer advantages over methanol for purge & trap concentration, leading to lower target compound method detection limits for CSIA of soil samples. The applicability of the developed method was demonstrated for a case study on soil extracts from a former manufacturing facility. Our approach extends the applicability of CSIA to an important matrix that often controls the long-term fate of contaminants in the subsurface.

Assessment of chlorinated ethenes degradation after field scale injection of activated carbon and bioamendments: Application of isotopic and microbial analyses

Over the last decade, activated carbon amendments have successfully been applied to retain chlorinated ethene subsurface contamination. The concept of this remediation technology is that activated carbon and bioamendments are injected into aquifer systems to enhance biodegradation. While the scientific basis of the technology is established, there is a need for methods to characterise and quantify the biodegradation at field scale. In this study, an integrated approach was applied to assess in situ biodegradation after the establishment of a cross sectional treatment zone in a TCE plume. The amendments were liquid activated carbon, hydrogen release donors and a Dehalococcoides containing culture. The integrated approach included spatial and temporal evaluations on flow and transport, redox conditions, contaminant concentrations, biomarker abundance and compound-specific stable isotopes. This is the first study applying isotopic and microbial techniques to assess field scale biodegradation enhanced by liquid activated carbon and bioamendments. The injection enhanced biodegradation from TCE to primarily cis-DCE. The Dehalococcoides abundances facilitated characterisation of critical zones with insufficient degradation and possible explanations. A conceptual model of isotopic data together with distribution and transport information improved process understanding; the degradation of TCE was insufficient to counteract the contaminant input by inflow into the treatment zone and desorption from the sediment. The integrated approach could be used to document and characterise the in situ degradation, and the isotopic and microbial data provided process understanding that could not have been gathered from conventional monitoring tools. However, quantification of degradation through isotope data was restricted for TCE due to isotope masking effects. The combination of various monitoring tools, applied frequently at high-resolution, with system understanding, was essential for the assessment of biodegradation in the complex, non-stationary system. Furthermore, the investigations revealed prospects for future research, which should focus on monitoring contaminant fate and microbial distribution on the sediment and the activated carbon.

Contaminant occurrence and migration between high- and low-permeability zones in groundwater systems: A review

In recent decades, water quality problems that impact human health, especially groundwater pollution, have been intensely studied, and this has contributed to new ideas and policies around the world such as Low Impact Development (LID) and Superfund legislation. The fundamental to many of these problems is pollutant occurrence and migration in saturated porous media, especially in groundwater. Such environments often contain contrasting zones of high and low permeability with significant differences in hydraulic conductivity (~10^-4 and 10^-8 m/s, respectively). High-permeability zones (HPZs) represent the primary pathways for pollutant transport in groundwater, while low-permeability zones (LPZs) are often diffusion dominated and serve as both sinks and sources (i.e., via back-diffusion) of pollutants over many decades. In this review, concepts and mechanisms of solute source depletion, contaminant accumulation, and back-diffusion in high- and low-permeability systems are presented, and new insights gained from both experimental and numerical studies are analyzed and summarized. We find that effluent monitoring and novel image analysis techniques have been adroitly used to investigate temporal and spatial evolutions of contaminant concentration; simultaneously, mathematical models are constantly upscaled to verify, optimize and extend the experimental data. However, the spatial concentration data during back-diffusion lacks diversity due to the limitations of pollutant species in studies, the microscopic mechanisms controlling pollutant transformation are poorly understood, and the impacts of these reactions on contaminant back-diffusion are rarely considered. Hence, most simulation models have not been adequately validated and are not capable of accurately predicting pollutant fate and cleanup in realistic heterogeneous aquifers. Based on these, some hypotheses and perspectives are mentioned to promote the investigation of contaminant migration in high- and low-permeability systems in groundwater.

The importance of transects for characterizing aged organic contaminant plumes in groundwater

A complex mixture of dissolved organic contaminants, emanating from a many decades-old, residual, dense non-aqueous phase liquid (DNAPL) source, migrates through unconfined, moderately heterogeneous, glacial-derived sediments and sedimentary rock in a residential area of Dane County, Wisconsin, USA. A portion of this contaminant plume intersects a large man-made pond, roughly 400 m downgradient of the source zone. Depth-discrete, multilevel groundwater sampling, detailed sedimentological logs, and hydraulic head profiles were used to delineate the spatial distribution of hydraulic, geologic, organic contaminant, and redox hydrochemical conditions within the established plume along two transects immediately upgradient of the pond. Twenty-one contaminants were detected and classified into four major contaminant groups: chlorinated ethenes, chlorinated ethanes, aromatics (BTEX: benzene, toluene, ethylbenzene, xylene), and aliphatic ketones. Within the glacial sediments and shallow bedrock, zones of reductive dechlorination of chlorinated ethenes and ethanes were juxtaposed with zones of BTEX and ketone degradation. Spatial heterogeneity in the concentration and distribution of contaminant groups and redox conditions was observed over lateral distances of tens of meters and vertical distances of tens of centimeters along the two transects. Although the site was situated in a complex glacial depositional environment, lithologic and hydraulic heterogeneity surprisingly only had a modest influence on the spatial distribution of plume contaminants. Depth-discrete sampling along paired, closely spaced transects (~20 m apart) was essential to assess internal plume composition/concentration evolution along flow paths with strong attenuation over short migration distances. This study shows how paired, highly resolved transects can enhance understanding of transverse and longitudinal variability in areas where contaminant-induced redox conditions control reaction zones and strong plume attenuation.

Controls on the persistence of aqueous-phase groundwater contaminants in the presence of reactive back-diffusion

The persistence of groundwater contaminants is influenced by several interacting processes. Physical, physico-chemical, and (bio-)chemical processes all influence the transport of contaminants in the subsurface. However, for a given hydrogeological system, it is generally unclear to which degree each of these phenomena acts as a control on plume behaviour. Here, we present a comprehensive investigation of these processes and their influences on plume behaviour and persistence in layered sedimentary systems. We investigate different scenarios that represent fundamental configurations of common contaminant situations. A confined aquifer over- and underlain by aquitard layers is investigated in a source-removal scenario and a constant-source equilibrium scenario. Additionally, an aquitard overlain and underlain by high permeability units is investigated in a source-removal scenario. In these investigations, we vary layer thickness, as well as parameters governing advection, (back-)diffusion, sorption, and degradation. Extensive analysis of these results enables quantification of the influence of these parameters on maximum down-gradient concentration, plume persistence duration, and plume spatial extent. Finally, parameter space dimensionality reduction is used to establish trends and regimes in which certain processes dominate as controls. A lower limit to plume extent as a function of a novel constructed parameter is also determined. These results provide valuable quantitative information for the assessment of the fate of groundwater contaminants and are applicable to a wide range of aqueous-phase solutes.

Virtual experiments to assess opportunities and pitfalls of CSIA in physical-chemical heterogeneous aquifers

Degradation of chlorinated ethenes (CEs) in low conductivity layers of aquifers reduces pollution plume tailing and accelerates natural attenuation timeframes. The degradation pathways involved are often different from those in the higher conductive layers and might go undetected when only highly conductive layers are targeted in site assessments. Reactive transport model simulations (PHT3D in FloPy) were executed to assess the performance of dual carbon and chlorine compound specific stable isotope analysis (CSIA) in degradation pathway identification and quantification in a coupled physical-chemical heterogeneous virtual aquifer. Degradation rate constants were assumed correlated to the hydraulic conductivity: positively for oxidative transformation (higher oxygen availability in coarser sands) and negatively for chemical reduction (higher content of reducing solids in finer sediments). Predicted carbon isotope ratios were highly heterogeneous. They generally increased downgradient of the pollution source but the large variation across depth illustrates that monotonously increasing isotope ratios downgradient, as were associated with the oxidative component, are not necessarily a common situation when degradation is favorable in low conductivity layers. Dual carbon-chlorine CSIA performed well in assessing the occurrence of the spatially separated degradation pathways and the overall degradation, provided appropriate enrichment factors were known and sufficiently different. However, pumping to obtain groundwater samples especially from longer well screens causes a bias towards overestimation of the contribution of oxidative transformation associated with the higher conductive zones. As degradation was less intense in these highly conductive zones under the simulated conditions, overall degradation was underestimated. In contrast, in the usual case of limited CSIA data, dual CSIA plots may rather indicate dominance of chemical reduction, while oxidative transformation could go unnoticed, despite being an equally important degradation pathway.

Assessing microbial degradation degree and bioavailability of BDE-153 in natural wetland soils: Implication by compound-specific stable isotope analysis

Microbial degradation is an important pathway for the attenuation of polybrominated diphenyl ethers (PBDEs) in natural soils. In this study, the compound-specific stable isotope analysis (CSIA) was applied to characterize microbial degradation of BDE-153, one of the prevailing and toxic PBDE congeners, in natural wetland soils. During the 45-day incubation, the residual percentages of BDE-153 decreased to 67.9% and 73.6% in non-sterilized soils spiked with 1.0 and 5.0 μg/g, respectively, which were both much lower than those in sterilized soils (96.0% and 97.2%). This result indicated that microbial degradation could accelerate BDE-153 elimination in wetland soils. Meanwhile, the significant carbon isotope fractionation was observed in non-sterilized soils, with δ¹³C of BDE-153 shifting from -29.4‰ to -26.7‰ for 1.0 μg/g and to -27.2‰ for 5.0 μg/g, respectively, whilst not in sterilized soils. This phenomenon indicated microbial degradation could induce stable carbon isotope fractionation of BDE-153. The carbon isotope enrichment factor (ε_c) for BDE-153 microbial degradation was first determined as -7.58‰, which could be used to assess the microbial degradation and bioavailability of BDE-153 in wetland soils. Based on δ¹³C and ε_c, the new methods were developed to dynamically and quantitatively estimate degradation degree and bioavailability of BDE-153 during degradation process, respectively, which could exclude interference of physical processes. This work revealed that CSIA was a promising method to investigate in situ microbial degradation of PBDEs in field studies.

Acceleration and centralization of a back-diffusion process: Effects of EDTA-2Na on cadmium migration in high- and low-permeability systems

Pollutant accumulation in the low-permeability zones (LPZs) in groundwater systems is regarded as a secondary source, and its consequent back-diffusion can extend the timeframe of pump-and-treat remediation. However, the bioavailability and mobility of heavy metals and the medium characteristics can be changed during the process. This study investigated the accumulation and back-diffusion law of toxic metals and the effects of ethylenediaminetetraacetic acid disodium salt (EDTA-2Na) on them by implementing a series of tank experiments. In these experiments, a cadmium solution was injected first, and deionized water or EDTA-2Na constantly washed the system consisting of different medium layers. The experimental results showed that the cadmium breakthrough curves had some concentration gradient reverse points where the curves fluctuated with elution by deionized water, which did not exist when EDTA-2Na was the eluent. In these scenarios, the mass of accumulated cadmium in the media before elution was large, with a value of 931 mg (153 mg/kg), when the low-permeability medium was clay. However, when EDTA-2Na was injected together with cadmium, the value dropped to 319 mg (52.3 mg/kg), greatly reducing the cadmium accumulation. Additionally, the use of EDTA-2Na as an eluent resulted in the appearance of a secondary peak in the breakthrough curve, showing that EDTA-2Na accelerated and centralized the back-diffusion. Notably, the reduced cadmium accumulation in LPZs with the elution by EDTA-2Na was partly due to a reduced adsorption capacity of the clay minerals. The above results can advance the technology related to pump-and-treat remediation.

Assessing toluene biodegradation under temporally varying redox conditions in a fractured bedrock aquifer using stable isotope methods

In complex hydrogeological settings little is known about the extent of temporally varying redox conditions and their effect on aromatic hydrocarbon biodegradation. This study aims to assess the impact of changing redox conditions over time on aromatic hydrocarbon biodegradation in a fractured bedrock aquifer using stable isotope methods. To that end, four snapshots of highly spatio-temporally resolved contaminant and redox sensitive species concentrations, as well as stable isotope ratio profiles, were determined over a two-years time period in summer 2016, spring 2017, fall 2017 and summer 2018 in a toluene contaminated fractured bedrock aquifer. The concentration profiles of redox sensitive species and stable isotope ratio profiles for dissolved inorganic carbon (DIC) and sulfate (δ¹³C_DIC, δ³⁴S_SO4, δ¹⁸O_SO4) revealed that the aquifer alternates between oxidising (spring 2017/summer 2018) and reducing conditions (summer 2016/fall 2017). This alternation was attributed to a stronger aquifer recharge with oxygen-rich meltwater in spring 2017/summer 2018 compared to summer 2016/fall 2017. The temporally varying redox conditions coincided with various extents of toluene biodegradation revealed by the different magnitude of heavy carbon (¹³C) and hydrogen (²H) isotope enrichment in toluene. This indicated that the extent of toluene biodegradation and its contribution to plume attenuation was controlled by the temporally changing redox conditions. The highest toluene biodegradation was observed in summer 2016, followed by spring 2017 and fall 2017, whereby these temporal changes in biodegradation occurred throughout the whole plume. Thus, under temporally varying recharge conditions both the core and the fringe of a contaminant plume can be replenished with terminal electron acceptors causing biodegradation in the whole plume and not only at its distal end as previously suggested by the plume fringe concept. Overall, this study highlights the importance of highly temporally resolved groundwater monitoring to capture temporally varying biodegradation rates and to accurately predict biodegradation-induced contaminant attenuation in fractured bedrock aquifers.

Do CSIA data from aquifers inform on natural degradation of chlorinated ethenes in aquitards?

Back-diffusion of chlorinated ethenes (CEs) from low-permeability layers (LPLs) causes contaminant persistence long after the primary spill zones have disappeared. Naturally occurring degradation in LPLs lowers remediation time frames, but its assessment through sediment sampling is prohibitive in conventional remediation projects. Scenario simulations were performed with a reactive transport model (PHT3D in FloPy) accounting for isotope effects associated with degradation, sorption, and diffusion, to evaluate the potential of CSIA data from aquifers in assessing degradation in aquitards. The model simulated a trichloroethylene (TCE) DNAPL and its pollution plume within an aquifer-aquitard-aquifer system. Sequential reductive dechlorination to ethene and sorption were uniform in the aquitard and did not occur in the aquifer. After 10 years of loading the aquitard through diffusion from the plume, subsequent source removal triggered release of TCE by back-diffusion. In the upper aquifer, during the loading phase, δ¹³C-TCE was slightly enriched (up to 2‰) due to diffusion effects stimulated by degradation in the aquitard. In the upper aquifer, during the release phase, (i) source removal triggered a huge δ¹³C increase especially for higher CEs, (ii) moreover, downstream decreasing isotope ratios (caused by downgradient later onset of the release phase) with temporal increasing isotope ratios reflect aquitard degradation (as opposed to downstream increasing and temporally constant isotope ratios in reactive aquifers), and (iii) the carbon isotope mass balance (CIMB) enriched up to 4‰ as lower CEs (more depleted, less sorbing) have been transported deeper into the aquitard. Thus, enriched CIMB does not indicate oxidative transformation in this system. The CIMB enrichment enhanced with more sorption and lower aquitard thickness. Thin aquitards are quicker flushed from lower CEs leading to faster CIMB enrichment over time. CIMB enrichment is smaller or nearly absent when daughter products accumulate. Aquifer CSIA patterns indicative of aquitard degradation were similar in case of linear decreasing rate constants but contrasted with previous simulations assuming a thin bioactive zone. The Rayleigh equation systematically underestimates the extent of TCE degradation in aquifer samples especially during the loading phase and for conditions leading to long remediation time frames (low groundwater flow velocity, thicker aquitards, strong sorption in the aquitard). The Rayleigh equation provides a good and useful picture on aquitard degradation during the release phase throughout the sensitivity analysis. This modelling study provides a framework on how aquifer CSIA data can inform on the occurrence of aquitard degradation and its pitfalls.

Isotope fractionation due to aqueous phase diffusion - What do diffusion models and experiments tell? - A review

For the interpretation of stable isotope ratio trends in saturated geochemical systems, the magnitude of aqueous phase diffusion-induced isotope fractionation needs to be known. This study reviews how five diffusion models (Fick, Maxwell-Stefan, Einstein, Langevin, Mode-Coupling Theory Analysis (MCTA) of diffusion) predict isotope fractionation due to aqueous phase diffusion and compares them with experimental results. The reviewed diffusion models were not consistent regarding the prediction of the mass (m) dependency of the aqueous phase diffusion coefficient (D). The predictions range from a square root power law (D ∝ m^-0.5) to an opposite mass dependency of D (D ∝ m^β). Experimental studies exhibited consistently a weak power law mass dependency of the diffusion coefficient (D ∝ m^-β with β < 0.5) for the vast majority of dissolved species and a larger diffusion-induced isotope effect for low weight noble gases (D ∝ m^-0.5). The weak power law mass dependency of D for the species other than low weight noble gases is consistent with the MCTA of diffusion. The MCTA suggests that the weak power law mass dependency of D originates from interplays between strongly mass dependent short-term and mass independent long-term solute-solvent interactions. The larger isotope fractionation for low weight noble gases could be attributed to quantum isotope effects significantly magnifying the aqueous phase diffusion-induced isotope fractionation. Our review shows, that except for low weight noble gases a weak power law mass dependency of D is likely the most adequate assumption for aqueous phase diffusion-induced isotope fractionation in geochemical systems.