Strategies for Managing Risk due to Back Diffusion

Managing the remediation strategy of contaminated megasites using field-scale calibration of geo-electrical imaging with chemical monitoring

Groundwater contamination is a threat to drinking water resources and ecosystems. Remediation by injection of chemical reagents into the aquifer may be preferred to excavation to reduce cost and environmental footprint. Yet, successful remediation requires complete contact between contamination and reagents. Subsurface heterogeneities are often responsible for diffusion into low-permeable zones, which may inhibit this contact. Monitoring the spatial distribution of injected reagents over time is crucial to achieve complete interaction. Source zone contamination at megasites is particularly challenging to remediate and monitor due to the massive scale and mixture of contaminants. Source zone remediation at Kærgård Plantation megasite (Denmark) is monitored here, with a new methodology, using high-resolution cross-borehole electrical resistivity tomography (XB-ERT) imaging calibrated by chemical analyses on groundwater samples. At this site, high levels of toxic non-aqueous phase liquids (NAPL) are targeted by in-situ chemical oxidation using activated persulfate. It may take numerous injection points with extensive injection campaigns to distribute reagents, which requires an understanding of how reagent may transport within the aquifer. A geophysical (XB-ERT) monitoring network of unprecedented size was installed to identify untreated zones and help manage the remediation strategy. The combination of spatially continuous geophysical information with discrete but precise chemical information, allowed detailed monitoring of sulfate distribution, produced during persulfate activation. Untreated zones identified in the first remediation campaign were resolved in the second campaign. The monitoring allowed adjusting the number of injection screens and the injection strategy from one campaign to the next, which resulted in better persulfate distribution and contaminant degradation in the second campaign. Furthermore, geophysical transects repeated over the timespan of a remediation campaign allowed high-resolution time-lapse imaging of reagent transport, which could in the future improve the predictability of transport models, compared to only using on a-priori assumptions of the hydraulic conductivity field.

Processes governing treatment of contaminants in low permeability zones

Herein we embrace the premise that aquifers are commonly composed of transmissive and low k (permeability) zones. Contaminants stored and subsequently released from low k zones sustain aqueous phase plumes for problematic periods. Processes governing the occurrence and treatment of contaminants in low k zones are advanced via conceptual models and a laboratory tank study. A two-dimensional sand tank with interbedded low k clay layers is flushed for 92 days with water spiked with 100 mg/L fluorescein, a proxy for chlorinated solvent contamination, and 67 mg/L bromide, a conservative tracer. Given active sources, fluorescein and bromide diffuses into the clay layers. Subsequently, the tank is flushed with water for 38 days. Water only flushing illustrates how the release of contaminants stored in low k zones sustains downgradient plumes. Next, an alkaline persulfate solution (40,000 mg/L persulfate at pH 11) is delivered to the tank for eight days. A fiber optic cable, placed on the glass wall of the sand tank, and a spectrometer with an ultraviolet light source are used to track depletion of fluorescein in transmissive sand and low k clay zones through time. Lastly, the tank is flushed with water only for 69 days to evaluate the efficacy of treatment with respect to mitigating releases from low k zones. Results indicate that flushing the tank with an alkaline persulfate solution, at a laboratory-scale, was effective in depleting fluorescein in both the transmissive and low k zones. Novelly, results capture concurrent transport of reactants and contaminants in domains governed by advection in transmissive zones and diffusion in low k zones.

Contaminant Back Diffusion from Low-Conductivity Matrices: Case Studies of Remedial Strategies

Recalcitrant groundwater contamination is a common problem at hazardous waste sites worldwide. Groundwater contamination persists despite decades of remediation efforts at many sites because contaminants sorbed or dissolved within low-conductivity zones can back diffuse into high-conductivity zones, and therefore act as a continuing source of contamination to flowing groundwater. A review of the available literature on remediation of plume persistence due to back diffusion was conducted, and four sites were selected as case studies. Remediation at the sites included pump and treat, enhanced bioremediation, and thermal treatment. Our review highlights that a relatively small number of sites have been studied in sufficient detail to fully evaluate remediation of back diffusion; however, three general conclusions can be made based on the review. First, it is difficult to assess the significance of back diffusion without sufficient data to distinguish between multiple factors contributing to contaminant rebound and plume persistence. Second, high-resolution vertical samples are decidedly valuable for back diffusion assessment but are generally lacking in post-treatment assessments. Third, complete contaminant mass removal from back diffusion sources may not always be possible. Partial contaminant mass removal may nonetheless have potential benefits, similar to partial mass removal from primary DNAPL source zones.

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.

Modeling a well-characterized perfluorooctane sulfonate (PFOS) source and plume using the REMChlor-MD model to account for matrix diffusion

Two of the most important retention processes for per- and polyfluoroalkyl substances (PFAS) in groundwater likely are sorption and matrix diffusion. The objective of this study was to model concentration and mass discharge of one PFAS, perfluorooctane sulfonate (PFOS), with matrix diffusion processes incorporated using data from a highly chemically- and geologically-characterized site. When matrix diffusion is incorporated into the REMChlor-MD model for PFOS at this research site, it easily reproduces the field data for three key metrics (concentration, mass discharge, and total mass). However, the no-matrix diffusion model produced a much poorer match. Additionally, after about 40 years of groundwater transport, field data and the REMChlor-MD model both showed the majority (80%) of the measured PFOS mass that exited the source zones was located in downgradient low permeability zones due to matrix diffusion. As such, most of the PFOS mass is not available to immediately migrate downgradient via advection in the more permeable sands at this site, which has important implications for monitored natural attenuation (MNA). Plume expansion over the next 50 years is forecasted to be limited, from a 350-m plume length in 2017 to 550 m in 2070, as matrix diffusion will attenuate groundwater plumes by slowing their expansion. This phenomenon is important for constituents that do not degrade, such as PFOS, compared to those susceptible to degradation. Overall, this work shows that matrix diffusion is a relevant process in environmental PFAS persistence and slows the rate of plume expansion over time.

Impact of matrix diffusion on the migration of groundwater plumes for Perfluoroalkyl acids (PFAAs) and other non-degradable compounds

Groundwater fate and transport modeling results demonstrate that matrix diffusion plays a role in attenuating the expansion of groundwater plumes of "non-degrading" or highly recalcitrant compounds. This is especially significant for systems where preferred destructive attenuation processes, such as biological and abiotic degradation, are weak or ineffective for plume control. Under these conditions, models of nondestructive physical attenuation processes, traditionally dispersion or sorption, do not demonstrate sufficient plume control unless matrix diffusion is considered. Matrix diffusion has been shown to be a notable emergent impact of geological heterogeneity, typically associated with back diffusion and extending remediation timeframes through concentration tailing of the trailing edge of a plume. However, less attention has been placed on evaluating how matrix diffusion can serve as an attenuation mechanism for the leading edge of a plume of non-degrading compounds like perfluoroalkyl acids (PFAAs), including perfluorooctane sulfonate (PFOS). In this study, the REMChlor-MD model was parametrically applied to a generic unconsolidated and heterogeneous geologic site with a constant PFOS source and no degradation of PFOS in the downgradient edge of the plume. Low levels of mechanical dispersion and retardation were used in the model for three different geologic heterogeneity cases ranging from no matrix diffusion (e.g., sand only) to considerable matrix diffusion using low permeability ("low-k") layers/lenses and/or aquitards. Our analysis shows that, in theory, many non-degrading plumes may expand for significant time periods before dispersion alone would eventually stabilize the plume; however, matrix diffusion can significantly slow the rate and degree of this migration. For one 100-year travel time scenario, consideration of matrix diffusion results in a simulated PFOS plume length that is over 80% shorter than the plume length simulated without matrix diffusion. Although many non-degrading plumes may continue to slowly expand over time, matrix diffusion resulted in lower concentrations and smaller plume footprints. Modeling multiple hydrogeologic settings showed that the effect of matrix diffusion is more significant in transmissive zones containing multiple low-k lenses/layers than transmissive zones underlain and overlain by low-k aquitards. This study found that at sites with significant matrix diffusion, groundwater plumes will be shorter, will expand more slowly, and may be amenable to a physical, retention-based, Monitored Natural Attenuation (MNA) paradigm. In this case, a small "Plume Assimilative Capacity Zone" in front of the existing plume could be reserved for slow, de minimus, future expansion of a non-degrading plume. If potential receptors are protected in this scenario, then this approach is similar to allowances for expanding plumes under some existing environmental regulatory programs. Accounting for matrix diffusion may support new strategic approaches and alternative paradigms for remediation even for sites and conditions with "non-degrading" constituents such as PFAAs, metals/metalloids, and radionuclides.

On the interplay between electromigration and electroosmosis during electrokinetic transport in heterogeneous porous media

Electrokinetic techniques represent a valuable approach to enhance solute transport, reactant delivery and contaminant degradation in complex environmental matrices, such as contaminated soil and groundwater, and have a great potential for the remediation of many organic and inorganic pollutants. This study investigates the complex interplay between the key electrokinetic transport mechanisms, electromigration and electroosmosis, in physically heterogeneous porous media and its impact on tracer distribution, reactant mixing and degradation efficiency. We perform experiments in a multidimensional setup, considering different types of heterogeneities, injected tracers and reactants, as well as background electrolyte pore water with different chemical composition and pH. We show that EK transport is significantly affected by the physical heterogeneities, due to the interaction between electrokinetic and hydraulic processes, and by the pore water chemistry that plays a key role on the magnitude and spatial distribution of electroosmotic fluxes. The latter affect the overall transport of charged and non-charged species, including the migration velocity of injected plumes, their spatial patterns, spreading and mixing with the background groundwater, and the extent of degradation and the spatio-temporal evolution of reactive zones in the heterogeneous porous media. Process-based numerical modeling allowed us to interpret the experimental observations and to disentangle the coupled effects of physical, chemical and electrostatic processes in the multidimensional, heterogeneous setups. Besides elucidating the mechanisms controlling electrokinetic transport, the results of this study have also important implications for practical field implementation of EK approaches in intrinsically heterogeneous subsurface systems.

Evaluation of the Bioelectrochemical Approach and Different Electron Donors for Biological Trichloroethylene Reductive Dechlorination

Trichloroethylene (TCE) and more in general chlorinated aliphatic hydrocarbons (CAHs) can be removed from a contaminated matrix thanks to microorganisms able to perform the reductive dechlorination reaction (RD). Due to the lack of electron donors in the contaminated matrix, CAHs' reductive dechlorination can be stimulated by fermentable organic substrates, which slowly release molecular hydrogen through their fermentation. In this paper, three different electron donors constituted by lactate, hydrogen, and a biocathode of a bioelectrochemical cell have been studied in TCE dechlorination batch experiments. The batch reactors evaluated in terms of reductive dechlorination rate and utilization efficiency of the electron donor reported that the bio-electrochemical system (BES) showed a lower RD rate with respect of lactate reactor (51 ± 9 µeq/d compared to 98 ± 4 µeq/d), while the direct utilization of molecular hydrogen gave a significantly lower RD rate (19 ± 8 µeq/d), due to hydrogen low solubility in liquid media. The study also gives a comparative evaluation of the different electron donors showing the capability of the bioelectrochemical system to reach comparable efficiencies with a fermentable substrate without the use of other chemicals, 10.7 ± 3.3% for BES with respect of 3.5 ± 0.2% for the lactate-fed batch reactor. This study shows the BES capability of being an alternative at classic remediation approaches.

3D dynamic model empowering the knowledge of the decontamination mechanisms and controlling the complex remediation strategy of a contaminated industrial site

Knowledge of the geology and hydrogeology of the polluted site emblematize a key requirement for environmental remediation, through assembling and synthesizing findings from various sources of physical evidence. In an increasingly virtual era, digital and geo-referenced metadata may serve as tools for collecting, merging, matching, and understanding multi-source information. The main goal of this paper is to emphasize the significance of a 3D hydrogeochemical model to the portrayal and the understanding of contamination dynamics and decontamination mechanisms at a highly contaminated industrial site. Some remediation measures are active on-site, due to the evidence-based presence of chlorinated solvents in groundwater. These are attributable to a slow-release source of pollutants in the saturated zone associated with very low permeability sediments. Therefore, in this research, a new technique for the remediation of secondary sources of dense non-aqueous phase liquid (DNAPL) contamination was investigated for the first time on a full-scale application. The combination of groundwater circulation wells (IEG-GCW®) and a continuous electron donor production device was set up to boost in situ bioremediation (ISB). A multi-phase approach was followed handling and releasing data during various remediation stages, from site characterization via pilot testing to full-scale remediation, thus allowing users to monitor, analyze, and manipulate information in 3D space-time. Multi-source and multi-temporal scenarios reveal the impact of ongoing hydraulic dynamics and depict the decontamination mechanisms in response to the interventions implemented over time, by quantifying the overall performance of the adopted strategies in terms of removal of secondary sources of pollution still active at the site.