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, δ13C-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 δ13C 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.

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.

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 37Cl 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 1H 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).