Back diffusion from thin low permeability zones

Assessing the persistence of a contaminant plume generated by linear aquifer source depletion and back diffusion from an aquitard

This study presents analytical solutions for describing contaminant storage and release from an aquitard with linear source depletion (LSD) boundary conditions. We investigated three scenarios for trichloroethylene (TCE) mass exchange before and after the LSD period in an aquifer bounded by an adjacent aquitard based on the LSD dynamics, a resistance coefficient, and the aquitard thickness. The developed analytical solutions showed good agreement with measured profiles and breakthrough curves from a previous study. In three scenarios, the factors delaying the onset of TCE release into the aquifer were a decrease in the resistance coefficient, an increase in LSD period and aquitard thickness. The changes in the duration, mass, and rate of TCE storage in the aquitard during LSD loading process affected the equilibrium of the aquifer-aquitard concentration gradient. After TCE loading, the period maintained above the maximum contaminant level was directly related to the three variables; the longest plume persistence occurred when TCE penetration distance at transition point from storage to release coincided with the aquitard thickness. Overall, the developed analytical solution aids in evaluating the risk of plume persistence, enhancing site management efficiency, and reducing remediation costs.

Diffusion in Porous Rock Is Anomalous

Molecular diffusion of chemical species in subsurface environments─rock formations, soil sediments, marine, river, and lake sediments─plays a critical role in a variety of dynamic processes, many of which affect water chemistry. We investigate and demonstrate the occurrence of anomalous (non-Fickian) diffusion behavior, distinct from classically assumed Fickian diffusion. We measured molecular diffusion through a series of five chalk and dolomite rock samples over a period of about two months. We demonstrate that in all cases, diffusion behavior is significantly different than Fickian. We then analyze the results using a continuous time random walk framework that can describe anomalous diffusion in heterogeneous porous materials such as rock. This methodology shows extreme long-time tailing of tracer advance as compared to conventional Fickian diffusion processes. The finding that distinct anomalous diffusion occurs ubiquitously implies that diffusion-driven processes in subsurface zones should be analyzed using tools that account for non-Fickian diffusion.

Mechanistic insights into contaminant transport dynamics in the saturated porous system in the presence of low permeability region using numerical simulations and temporal moment analysis

The influence of low permeability porous media (LPPM) on contaminant transport dynamics in saturated porous systems was investigated using numerical simulations and temporal moments of contaminant concentrations. Two-dimensional flow and contaminant transport simulations were conducted, considering various parameters such as longitudinal dispersivity (ranging from 15 to 60 m), the ratio of transverse to longitudinal dispersivity (ranging from 0.05 to 0.2), retardation factor (ranging from 1 to 4), and hydraulic gradient (ranging from 0.005 to 0.02) for both homogeneous and heterogeneous porous systems. The findings revealed significant differences in the transport behavior of conservative and highly reactive contaminants between the porous systems without and with the LPPM region. The center of mass of contaminant plume and peak concentration zone were observed inside the LPPM region for the heterogeneous porous system, especially during the source off-loading period. Furthermore, asymmetric distributions of the zeroth temporal moment (ZTM), mean residence time (MRT), and variance of the breakthrough curve (BTC) were observed along the longitudinal distance within the LPPM region for heterogeneous porous system, highlighting the impact of heterogeneity on contaminant plume evolution dynamics. The moment analysis results provided insights into the influence of LPPM region on time-averaged contaminant transport dynamics in adjacent porous systems. These findings can help risk managers understand the complex fate and transport dynamics in heterogeneous porous systems. Future studies could explore the modelling of multispecies contaminants in heterogeneous saturated porous systems subjected to fluctuating water table.

Contaminant back-diffusion from layered aquitards subjected to barrier-controlled source zones

Low-permeability aquitards may serve as secondary sources of slow-releasing contaminants into the adjacent aquifer system, creating considerable obstacles to groundwater cleanup. Accurately capturing the exchange of contaminant mass between aquitards and aquifers can facilitate site management and remediation. Previous simulation studies were mainly limited to one-dimensional (1D) back diffusion from aquitards during the remediation of the source zone. In this study, a novel two-dimensional (2D) back-diffusion model is developed to investigate the storage and release of contaminants in aquitards after source isolation. This model coupled the dynamical decay of isolated sources and the diffusion-sorption process of contaminants in the layered aquitards. Exact analytical solutions for the present 2D multilayer model were derived using the finite cosine transform, Duhamel Theorem, separation of variables, and transfer matrix method. Results indicated that the previous 1D model would overestimate the contaminant concentration in the aquitard and the back-diffusion risk when the source zone was isolated. The proposed 2D back-diffusion model enables quantitative prediction of how source zone width, source concentration, and aquitard heterogeneity impact plume trailing time, thus aiding in understanding the mechanisms of aquifer contamination beyond barrier-controlled source zones.

Enhanced semi-volatile DNAPL accessibility at sub-boiling temperature during electrical resistance heating in heterogeneous porous media

Successful remediation of semi-volatile contaminants using electrical resistance heating (ERH) coupled technologies requires a deep understanding of contaminant migration and accessibility, especially with stratigraphic heterogeneity and dense nonaqueous phase liquid (DNAPL) occurrence. Here, we chose nitrobenzene (NB) as a model contaminant of semi-volatile DNAPL and uniquely demonstrated that temperature variation during ERH could induce NB DNAPL migration out of the low permeability zone (LPZ) even below water boiling temperature. When heating the system using alternating current (AC) of 140 V to a temperature range of 50-79 °C, obvious DNAPL migration was visually observed. The upward migration of DNAPL would considerably increase the mass of accessible contaminant by other remediation measures. The downstream cumulative NB mass of 1092 mg in 140 V system raised 56-folds compared to that of 19 mg in the control experiment with only groundwater flow. This migration was mainly attributed to a complex natural convection caused by temperature gradient. Comparing with traditional AC heating, ERH powered by pulsed direct current (PDC-ERH) showed a higher and more uneven heating pattern, resulting in a stronger convection at the same voltage that enhanced the DNAPL migration out of LPZ. These results revealed the importance of natural convection in the ERH process, which could be further optimized to improve the energy efficiency of remediation.

PFOS Mass Flux Reduction/Mass Removal: Impacts of a Lower-Permeability Sand Lens within Otherwise Homogeneous Systems

Perfluorooctane sulfonic acid (PFOS) is one of the most common per- and polyfluoroalkyl substances (PFAS) and is a significant risk driver for these emerging contaminants of concern. A series of two-dimensional flow cell experiments was conducted to investigate the impact of flow field heterogeneity on the transport, attenuation, and mass removal of PFOS. A simplified model heterogeneous system was employed consisting of a lower-permeability fine sand lens placed within a higher-permeability coarse sand matrix. Three nonreactive tracers with different aqueous diffusion coefficients, sodium chloride, pentafluorobenzoic acid, and β-cyclodextrin, were used to characterize the influence of diffusive mass transfer on transport and for comparison to PFOS results. The results confirm that the attenuation and subsequent mass removal of the nonreactive tracers and PFOS were influenced by mass transfer between the hydraulically less accessible zone and the coarser matrix (i.e., back diffusion). A mathematical model was used to simulate flow and transport, with the values for all input parameters determined independently. The model predictions provided good matches to the measured breakthrough curves, as well as to plots of reductions in mass flux as a function of mass removed. These results reveal the importance of molecular diffusion and pore water velocity variability even for systems with relatively minor hydraulic conductivity heterogeneity. The impacts of the diffusive mass transfer limitation were quantified using an empirical function relating reductions in contaminant mass flux (MFR) to mass removal (MR). Multi-step regression was used to quantify the nonlinear, multi-stage MFR/MR behavior observed for the heterogeneous experiments. The MFR/MR function adequately reproduced the measured data, which suggests that the MFR/MR approach can be used to evaluate PFOS removal from heterogeneous media.

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.

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.

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.

Enhancing DNAPL removal from low permeability zone using electrical resistance heating with pulsed direct current

Electrical resistance heating (ERH) has been widely applied for contaminant remediation in heterogeneous sites especially when low permeability zones exist, yet requires high energy input. To address the low energy-efficiency of ERH using conventional alternating current (AC), pulsed direct current (PDC) obtained by current rectification was introduced for heating to enhance dense nonaqueous phase liquid (DNAPL) migration in low permeability zones. Here we showed the proof-of-concept in a lab-scale two-dimensional heterogeneous sand system (40 cm × 30 cm) with trichloroethylene (TCE) DNAPL in the central low permeability zone. Applying PDC achieved faster temperature increase compared to that with the conventional AC of the same voltage gradient. The overall TCE removal efficiency from the cell increased from 79.0% to 89.6% with increasing PDC voltage gradient from 3 to 3.75 V cm^-1, compared to that of 9.4-91.1% with conventional AC. The lowest energy consumption of PDC was 390 kWh kg^-1 at a medium voltage gradient of 3.5 V cm^-1, which was 27.8% lower compared to that with AC at the same voltage gradient. These results suggest that remediation using pulsed direct current is a promising approach to improve the energy-efficiency and effectiveness of ERH.

Increasing electron donor concentration does not accelerate complete microbial reductive dechlorination in contaminated sediment with native organic carbon

Experiments with Fe(III)-rich, chloroethene-contaminated sediment demonstrated that trichloroethylene (TCE) and vinyl chloride (VC) were completely reduced to ethene regardless of whether electron donor(s) were added at 1 × stoichiometry or 10 × stoichiometry relative to all-electron acceptors. Unamended controls uniformly reduced TCE to ethene with a mean time to complete dechlorination (operationally defined as the presence of stoichiometric ethene production) of 79 days. Adding 1 × and 10 × acetate hindered the rate and extent of TCE and VC reduction relative to unamended controls, with several only partially reduced when the experiments were terminated. Adding high molecular mass (soybean oil derivative) substrates did not increase microbial reductive dechlorination relative to unamended incubations, and in many cases, hindered microbial dechlorination in favor of methanogenesis. The mean time to complete dechlorination was comparable between low (× 1) and high (× 10) electron donor concentration for all lipid-based electron donors tested. Those tested included Newman Zone® Standard without sodium lactate (96 vs. 75 days, respectively), CAP 18 ME (85 vs. 94 days, respectively), EOS 598B42 (68 vs. 72 days, respectively), and acetate (134 vs. 125 days, respectively). These data suggest that the addition of an electron donor does not always increase the rate and extent of reductive dechlorination but will increase costs. In particular, increasing the concentration of electron donors higher than the stoichiometric demand only decreased complete microbial reductive dechlorination, which is the opposite of most standard "more time and more electrons" approaches. These data argue that site-specific electron donor demands must be evaluated, and in some cases, a monitored natural attenuation (MNA) approach is most favorable.

Using MODFLOW and RT3D to simulate diffusion and reaction without discretizing low permeability zones

Low permeability zones (LPZs) are major sources of groundwater contamination after active remediation to remove pollutants in adjacent high permeability zones (HPZs). Slow back diffusion from LPZs to HPZs can extend management of polluted sites by decades. Numerical models are often used to simulate back diffusion, estimate cleanup times, and develop site management strategies. Sharp concentration gradients of pollutants are present at the interface between HPZs and LPZs, and hence accurate simulation requires fine grid sizes resulting in high computational burden. Since the MODFLOW family of codes is widely used in practice, we develop a new approach for modeling pollutant back diffusion using MODFLOW/RT3D that eliminates the need for fine discretization of the LPZ. Instead, the LPZ is treated as an impermeable region in MODFLOW, while in RT3D the LPZ is conceptualized as a series of immobile zones coupled with a mobile zone at the HPZ/LPZ interface. Finite volume discretization of diffusion and reaction within the LPZ is then modeled as mass transfer and reaction among several immobile species. This results in a simulation domain with significantly fewer grid cells compared to that required if all LPZs are discretized, providing potential for improved computational efficiency. Cases, including a layer of HPZ over an LPZ, a thin/thick lens of LPZ embedded in HPZ, and multiple lens of LPZs embedded in HPZ are tested by the new approach for tracer and reactive scenarios.

Impact of diffuse layer processes on contaminant forward and back diffusion in heterogeneous sandy-clayey domains

Low-permeability aquitards can significantly affect the transport, distribution, and persistence of contaminant plumes in subsurface systems. Although such low-permeability materials are often charged, the key role of charge-induced electrostatic processes during contaminant transport has not been extensively studied. This work presents a detailed investigation exploring the coupled effects of heterogeneous distribution of physical, chemical and electrostatic properties on reactive contaminant transport in field-scale groundwater systems including spatially distributed clay zones. We performed an extensive series of numerical experiments in three distinct heterogeneous sandy-clayey domains with different levels of complexity. The flow and reactive transport simulations were performed by explicitly resolving the complex velocity fields, the small-scale electrostatic processes, the compound-specific diffusive/dispersive fluxes and the chemical processes utilizing a multi-continua based reactive transport code (MMIT-Clay). In each particular domain, numerical experiments were performed focusing on both the forward and back diffusion through the sandy-clayey interfaces. The results illuminate the control of microscopic electrostatic mechanisms on macroscopic mass transfer. Coulombic interactions in the clay's diffuse layer can significantly accelerate or retard a particular species depending on its charge. Furthermore, the chemical heterogeneity plays a major role in mass storage and release during reactive transport. Neglecting such processes can lead to substantial over- or underestimation of the overall transport behavior, which underlines the need for integrated physical, chemical and electrostatic approaches to accurately describe mass transfer processes in systems including low-permeability inclusions.

Mass-Based, Field-Scale Demonstration of PFAS Retention within AFFF-Associated Source Areas

Transport of poly- and perfluoroalkyl substances (PFAS) at aqueous film-forming foam (AFFF)-impacted sites is limited by various processes that can retain PFAS mass within the source area. This study used concentration data obtained via a high-resolution sampling and analytical protocol to estimate the PFAS mass distribution in source and downgradient areas of a former firefighter training area. The total PFAS mass present at the site was approximately 222 kg, with 106 kg as perfluoroalkyl acids (PFAAs) and 116 kg as polyfluorinated precursors. Zwitterionic and cationic PFAS represented 83% of the total precursor mass and were found primarily in the source and up/side-gradient areas (75%), likely due to preferential hydrophobic partitioning, electrostatic interactions, and diffusion into lower-permeability soils. Based on the release history and the high percentage of total PFAS mass represented by precursors (primarily electrochemical fluorination-derived compounds), the estimated conversion rate of precursors to PFAAs was less than 2% annually. Eighty-two percent of the total PFAS mass was encountered in lower-permeability soils, which limited the potential for advection and transformation. This contributed to a 99% decrease in the mass discharge rate at the far-downgradient plume (0.048 kg/yr compared to the near-source area (3.6 kg/yr)). The results provide field-scale evidence of the importance of these PFAS retention processes at sites where AFFF has been released.

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.

Phase-transfer catalysis enhanced remediation of trichloroethylene polluted groundwater by potassium permanganate

As one of the remediation reagents, potassium permanganate (KMnO₄) is injected to the aquifer, degrading trichloroethylene (TCE) by chemical oxidation. This study investigated the kinetics of TCE degradation by series of batch experiments, as well as the influence of medium size. Moreover, phase-transfer catalyst (PTCs), such as pentyltriphenylphosphonium bromide (PTPP) and sodium hexametaphosphate (SHMP) were used for enhancing oxidation. The batch experimental results showed that in the absence of PTC, the removal efficiency of TCE was 36.14% and 86.79% within 4 and 30 min, respectively. However, the removal rate of TCE was up to 67.48% and 49.90% within 4 min for 15 mol% PTPP- and SHMP-added system, respectively. The results indicated that PTPP and SHMP promoted the depletion of M n O 4 - to oxidize DNAPL TCE, but its effectiveness varied with the addition ratio of PTPP or SHMP. Its promotion was more remarkable when PTC added with a higher proportion. The alleviation of MnO₂ by phosphates ( P O 4 3 - , H P O 4 2 - and H 2 P O 4 - ) or PTC in the presence of media was qualitatively investigated. Results showed that the content of MnO₂ in the dissolved phase during the same reaction period decreased by PTC. Moreover, H P O 4 2 - and SHMP have apparent beneficial effects of reducing MnO₂ formation. The presence of aquifer media has a pH buffer and a negative influence on the reaction between TCE and the oxidant; moreover, as particle size of media decreased, the negative effect increased.

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.

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.

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.

Contributions of biotic and abiotic pathways to anaerobic trichloroethene transformation in low permeability source zones

Low permeability source zones sustain long-term trichloroethene (TCE) groundwater contamination. In anaerobic environments, TCE is transformed by both biological reductive dechlorination and abiotic reactions with reactive minerals. Little is known about the relative contribution of these two pathways as TCE diffuses from low permeability zones (LPZs) into high permeability zones (HPZs). This study combines a flow cell experiment, batch experiments, and a diffusion-reaction model to evaluate the contributions of biotic and abiotic TCE transformation in LPZs. Natural clay (LPZ) and sand (HPZ) from a former Air Force base were used in all experiments. In batch, the LPZ material transformed TCE and cis-1,2-dichloroethene (cis-DCE) to acetylene with pseudo first-order rate constants of 8.57 × 10^-6 day^-1 and 1.02 × 10^-6 day^-1, respectively. Biotic and abiotic pathways were then evaluated together in a bench-scale flow cell (16.5 cm × 2 cm × 16.5 cm) that contained a LPZ layer, with a source of TCE at the base, overlain by a HPZ continuously purged with lactate-amended groundwater. Diffusion controlled mass transfer in the LPZ, while advection controlled migration in the HPZ. The mass discharge rate of TCE and its biotic (cis-DCE and vinyl chloride) and abiotic (acetylene) transformation products were measured over 180 days in the flow cell effluent. Depth profiles of these compounds through the LPZ were determined after terminating the experiment. A one-dimensional diffusion-reaction model was used to interpret the effluent and depth profile data and constrain reaction parameters. Abiotic transformation rate constants for TCE to acetylene, normalized to in situ solids loading, were approximately 13 times greater in batch than in the flow cell. Slower transformation rates in the flow cell indicate elevated TCE concentration and/or further degradation of acetylene to other reduced gas compounds in the flow cell LPZ (thereby partially masking TCE abiotic transformation). Biotic and abiotic parameters used to interpret the flow cell data were then used to simulate a field site with a 300 cm thick LPZ. Abiotic processes contributed to a 2% reduction in TCE flux after 730 days. When abiotic rate constants were changed to that observed in batch, or to rate constants previously reported for a pyrite rich mudstone, the TCE flux reduction was 21% and 53%, respectively, after 730 days. Though biotic processes dominated TCE transformation in the flow cell experiment, the simulations indicate that abiotic processes have potential to significantly contribute to TCE attenuation in electron donor limited environments provided suitable reactive minerals are present.

Experimental and numerical evaluation of Groundwater Circulation Wells as a remediation technology for persistent, low permeability contaminant source zones

Contaminants removal stoked inside low permeability zones of aquifers is one of the most important challenge of groundwater remediation process today. Low permeability layers can be considered persistent secondary sources of contamination because they release pollutants by molecular diffusion after primary source of contamination is reduced, causing long plum tails (Back-Diffusion). In this study, the Groundwater Circulation Well (GCW) system was investigated as an alternative remediation technology to the low efficient traditional pumping technologies to restore contaminated low permeability layers of aquifers. The GCW system creates vertical groundwater circulation cells by drawing groundwater through a screen of a multi-screen well and discharging it through another screen. The suitability of this technology to remediate contaminated low permeability zones was investigated by laboratory test and numerical simulations. The collected data were used to calibrate a model created to simulate the Back-Diffusion process and to evaluate the effect of different pumping technologies on the depletion time of that process. Results show that the efficiency of the GCW is dependent on the position and on the geometry of the low permeability zones, however the GCW system appears more suitable to restore contaminated low permeability layers of aquifers than the traditional pumping technology.

Flow field dynamics and high ethanol content in gasohol blends enhance BTEX migration and biodegradation in groundwater

Gasohol spills may easily descend through the soil column down and impact sensitive receptors as contaminants dissolve into the groundwater. Gasoline formulations are commonly blended with ethanol to alleviate environmental and economic issues associated with fossil fuels. However, the amount of ethanol added to gasoline and the groundwater hydraulic regime can significantly affect BTEX plume dynamics and lifespan. In this study, two long-term (5 and 10 years) field-scale gasohol releases with ethanol contents of 85% (E85) and 24% (E24), respectively, were assessed to discern the different dynamics undergone by gasohol blends. Statistical, geochemical, microbiological and trend approaches were employed to estimate the influence of groundwater flow variations on ethanol and dissolved BTEX transport, and the associated biodegradation rates of different gasohol blend spills. Ethanol and BTEX groundwater flow were quantified in terms of breakthrough curve characteristics, plume centroid positions and spreading, source depletion and mass degradation rates. In addition, bromide migration was evaluated to address the contribution of flow-driven dissolution. Results revealed that the high amount of ethanol along with a fast and dynamic flow exerted a flushing behavior that enhanced BTEX dissolution, migration (vertical and horizontal) and concentrations in groundwater. The higher amount of ethanol in E85 enhanced BTEX dissolution (and bioavailability) relative to E24 site and led to faster biodegradation rates, which can be explained by the cosolvency effect and metabolic flux dilution. Therefore, flow field dynamics and high ethanol content in gasohol blends enhance BTEX migration and biodegradation in gasohol-contaminated sites. The balance of these factors is crucial to determine fate and transport of contaminants in field sites. These findings suggest that hydraulic regime should be spatially and temporally characterized to support decisions on appropriate monitoring plan and remedial strategies for gasohol spills.

Diffusion of solutes from depleting sources into and out of finite low-permeability zones

Two important factors that affect groundwater contaminant persistence are the temporal pattern of contaminant source depletion and solute diffusion into and out of aquitards. This study provides a framework to evaluate the relative importance of these effects on contaminant persistence, with emphasis on the importance of thin aquitards. We developed one-dimensional (1D) analytical solutions for forward and back diffusion in a finite domain with a no flux boundary using the method of images and demonstrated their applicability to measured data from three well-controlled laboratory diffusion experiments with exponentially depleting sources. We used both in situ aquitard solute concentrations and aquifer breakthrough curves for sorbing and non-sorbing solutes. The finite-domain no flux boundary solutions showed better agreement with measured data than was available with semi-infinite approaches, with increasing discrepancy for dimensionless relative diffusion length scale beyond a critical threshold value (Z_d > 0.7). We also used a mass balance to demonstrate that the temporal pattern of contaminant source depletion controls the duration of solute mass accumulation in the aquitard, as well as the total solute mass release back into the aquifer. Lower rates of source depletion result in a longer period of mass accumulation in the aquitard and later back diffusion initiation time. The amount of solute mass stored in the aquitard increases with longer loading duration, thereby contributing to overall longer contaminant persistence in aquifers. This study entails widespread implications for anthropogenic waste and contamination sites, which are all dependent on efficient and cost-effective contaminant management strategies.

Semi-analytical method for matrix diffusion in heterogeneous and fractured systems with parent-daughter reactions

A semi-analytical/numerical method for modeling matrix diffusion in heterogeneous and fractured groundwater systems is developed. This is a significant extension of the Falta and Wang (2017) method that only applied to diffusion in an aquitard of infinite thickness. The current solution allows for the low permeability matrix to be embedded within a numerical gridblock, having finite average thickness, a specified volume fraction and a specified interfacial area with the high permeability domain. The new formulation also allows for coupled parent-daughter decay reactions with multiple species that each have independent retardation factors, decay rates, and yield coefficients in both the high and low permeability parts of the system. The method uses a fitting function to approximate the transient concentration profile in the low permeability part of each gridblock so that the matrix diffusion flux into the high permeability part of the gridblock can be computed as a concentration dependent source-sink term. This approach is efficient because the only unknowns in each gridblock are the concentrations in the high permeability domain, so there is practically no increase in computational effort compared to a conventional transport simulation. The method is shown to compare favorably with an analytical solution for matrix diffusion in fractured media with parallel fractures, with an analytical solution for matrix diffusion with parent-daughter decay reactions, with laboratory experiments of matrix diffusion in a layered system, with a laboratory experiment involving lens shaped inclusions, and with fine grid numerical simulations of transport in highly heterogeneous systems.

Natural attenuation of chlorinated ethenes in hyporheic zones: A review of key biogeochemical processes and in-situ transformation potential

Chlorinated ethenes (CEs) are legacy contaminants whose chemical footprint is expected to persist in aquifers around the world for many decades to come. These organohalides have been reported in river systems with concerning prevalence and are thought to be significant chemical stressors in urban water ecosystems. The aquifer-river interface (known as the hyporheic zone) is a critical pathway for CE discharge to surface water bodies in groundwater baseflow. This pore water system may represent a natural bioreactor where anoxic and oxic biotransformation process act in synergy to reduce or even eliminate contaminant fluxes to surface water. Here, we critically review current process understanding of anaerobic CE respiration in the competitive framework of hyporheic zone biogeochemical cycling fuelled by in-situ fermentation of natural organic matter. We conceptualise anoxic-oxic interface development for metabolic and co-metabolic mineralisation by a range of aerobic bacteria with a focus on vinyl chloride degradation pathways. The superimposition of microbial metabolic processes occurring in sediment biofilms and bulk solute transport delivering reactants produces a scale dependence in contaminant transformation rates. Process interpretation is often confounded by the natural geological heterogeneity typical of most riverbed environments. We discuss insights from recent field experience of CE plumes discharging to surface water and present a range of practical monitoring technologies which address this inherent complexity at different spatial scales. Future research must address key dynamics which link supply of limiting reactants, residence times and microbial ecophysiology to better understand the natural attenuation capacity of hyporheic systems.

Pore-scale simulations of concentration tails in heterogeneous porous media

The retention of contaminants in the finest and less-conductive regions of natural aquifer is known to strongly affect the decontamination of polluted aquifers. In fact, contaminant transfer from low to high mobility regions at the back end of a contaminant plume (i.e. back diffusion) is responsible for the long-term release of contaminants during remediation operation. In this paper, we perform pore-scale calculations for the transport of contaminant through heterogeneous porous media composed of low and high mobility regions with two objectives: (i) study the effect of permeability contrast and solute transport conditions on the exchange of solutes between mobile and immobile regions and (ii) estimate the mass of contaminants sequestered in low mobility regions based on concentration breakthrough curves.

Field-scale forward and back diffusion through low-permeability zones

Understanding the effects of back diffusion of groundwater contaminants from low-permeability zones to aquifers is critical to making site management decisions related to remedial actions. Here, we combine aquifer and aquitard data to develop recommended site characterization strategies using a three-stage classification of plume life cycle based on the solute origins: aquifer source zone dissolution, source zone dissolution combined with back diffusion from an aquitard, and only back diffusion. We use measured aquitard concentration profile data from three field sites to identify signature shapes that are characteristic of these three stages. We find good fits to the measured data with analytical solutions that include the effects of advection and forward and back diffusion through low-permeability zones, and linearly and exponentially decreasing flux resulting from source dissolution in the aquifer. Aquifer contaminant time series data at monitoring wells from a mature site were well described using analytical solutions representing the combined case of source zone and back diffusion, while data from a site where the source had been isolated were well described solely by back diffusion. The modeling approach presented in this study is designed to enable site managers to implement appropriate remediation technologies at a proper timing for high- and low-permeability zones, considering estimated plume life cycle.

A semi-analytical method for simulating matrix diffusion in numerical transport models

A semi-analytical approximation for transient matrix diffusion is developed for use in numerical contaminant transport simulators. This method is an adaptation and extension of the heat conduction method of Vinsome and Westerveld (1980) used to simulate heat losses during thermally enhanced oil recovery. The semi-analytical method is used in place of discretization of the low permeability materials, and it represents the concentration profile in the low permeability materials with a fitting function that is adjusted in each element at each time-step. The resulting matrix diffusion fluxes are added to the numerical model as linear concentration-dependent source/sink terms. Since only the high permeability zones need to be discretized, the numerical formulation is extremely efficient compared to traditional approaches that require discretization of both the high and low permeability zones. The semi-analytical method compares favorably with the analytical solution for transient one-dimensional diffusion with first order decay, with a two-layer aquifer/aquitard solution, with the solution for transport in a fracture with matrix diffusion and decay, and with a fully numerical solution for transport in a thin sand zone bounded by clay with variable decay rates.

Solute source depletion control of forward and back diffusion through low-permeability zones

Solute diffusive exchange between low-permeability aquitards and high-permeability aquifers acts as a significant mediator of long-term contaminant fate. Aquifer contaminants diffuse into aquitards, but as contaminant sources are depleted, aquifer concentrations decline, triggering back diffusion from aquitards. The dynamics of the contaminant source depletion, or the source strength function, controls the timing of the transition of aquitards from sinks to sources. Here, we experimentally evaluate three archetypical transient source depletion models (step-change, linear, and exponential), and we use novel analytical solutions to accurately account for dynamic aquitard-aquifer diffusive transfer. Laboratory diffusion experiments were conducted using a well-controlled flow chamber to assess solute exchange between sand aquifer and kaolinite aquitard layers. Solute concentration profiles in the aquitard were measured in situ using electrical conductivity. Back diffusion was shown to begin earlier and produce larger mass flux for rapidly depleting sources. The analytical models showed very good correspondence with measured aquifer breakthrough curves and aquitard concentration profiles. The modeling approach links source dissolution and back diffusion, enabling assessment of human exposure risk and calculation of the back diffusion initiation time, as well as the resulting plume persistence.

Implications of matrix diffusion on 1,4-dioxane persistence at contaminated groundwater sites

Management of groundwater sites impacted by 1,4-dioxane can be challenging due to its migration potential and perceived recalcitrance. This study examined the extent to which 1,4-dioxane's persistence was subject to diffusion of mass into and out of lower-permeability zones relative to co-released chlorinated solvents. Two different release scenarios were evaluated within a two-layer aquifer system using an analytical modeling approach. The first scenario simulated a 1,4-dioxane and 1,1,1-TCA source zone where spent solvent was released. The period when 1,4-dioxane was actively loading the low-permeability layer within the source zone was estimated to be <3years due to its high effective solubility. While this was approximately an order-of-magnitude shorter than the loading period for 1,1,1-TCA, the mass of 1,4-dioxane stored within the low-permeability zone at the end of the simulation period (26kg) was larger than that predicted for 1,1,1-TCA (17kg). Even 80years after release, the aqueous 1,4-dioxane concentration was still several orders-of-magnitude higher than potentially-applicable criteria. Within the downgradient plume, diffusion contributed to higher concentrations and enhanced penetration of 1,4-dioxane into the low-permeability zones relative to 1,1,1-TCA. In the second scenario, elevated 1,4-dioxane concentrations were predicted at a site impacted by migration of a weak source from an upgradient site. Plume cutoff was beneficial because it could be implemented in time to prevent further loading of the low-permeability zone at the downgradient site. Overall, this study documented that 1,4-dioxane within transmissive portions of the source zone is quickly depleted due to characteristics that favor both diffusion-based storage and groundwater transport, leaving little mass to treat using conventional means. Furthermore, the results highlight the differences between 1,4-dioxane and chlorinated solvent source zones, suggesting that back diffusion of 1,4-dioxane mass may be serving as the dominant long-term "secondary source" at many contaminated sites that must be managed using alternative approaches.

Image analysis procedure for studying Back-Diffusion phenomena from low-permeability layers in laboratory tests

In this study, the long-term tailing derived from the storage process of contaminants in low-permeability zones is investigated. The release from these areas in the groundwater can be considered a long-term source that often undermines remediation efforts. An Image Analysis technique is used to analyze the process and evaluate the concentrations of a tracer at different points of the test section. Furthermore, the diffusive flux from the low-permeability lenses is determined. To validate the proposed technique, the results are compared with samples, and the diffusive fluxes resulting from the low-permeability zones of the reconstructed aquifer are compared with a theoretical approach.

Functionalization of micromodels with kaolinite for investigation of low salinity oil-recovery processes

Sandstone formations are ubiquitous in both aquifers and petroleum reservoirs, of which clay is a major constituent. The release of clay particles from pore surfaces as a result of reduced injection fluid salinity can greatly modify the recovery of hydrocarbons from subsurface formations by shifting the wettability properties of the rock. In this paper we demonstrate a microfluidic approach whereby kaolinite is deposited into a two-dimensional microfluidic network (micromodel) to enable direct pore-scale, real-time visualization of fluid-solid interactions with representative pore-geometry and realistic surface interactions between the reservoir fluids and the formation rock. Structural characterization of deposited kaolinite particles agrees well with natural modes of occurrence in Berea sandstones; hence, the clay deposition method developed in this work is validated. Specifically, more than 90% of the deposited clay particles formed pore-lining structures and the remainder formed pore bridging structures. Further, regions of highly concentrated clay deposition likely leading to so-called Dalmatian wetting properties were found throughout the micromodel. Two post-deposition treatments are described whereby clay is adhered to the silicon surface reversibly and irreversibly resulting in microfluidic systems that are amenable to studies on (i) the fundamental mechanisms governing the increased oil recovery during low salinity waterfloods and (ii) the effect of a mixed-wet surface on oil recovery, respectively. The reversibly functionalized platform is used to determine the conditions at which stably adhered clay particles detach. Specifically, injection brine salinity below 6000 ppm of NaCl induced kaolinite particle release from the silicon surface. Furthermore, when applied to an aged system with crude oil, the low salinity waterflood recovered an additional 14% of the original oil in place compared to waterflooding with the formation brine.