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

Laboratory validation of a simplified model for estimating equilibrium PFAS mass leaching from unsaturated soils

Modelling per- and polyfluoroalkyl substance (PFAS) fate and transport in the vadose zone is inherently more complex than in the saturated zone due to the highly transient nature and the wetting phase saturation dependent hydraulic flux associated with the vadose zone. The chemical complexity of PFAS impart multiple partitioning processes which complicate the evaluation of PFAS transport in the vadose zone. To date, simplified screening models describing PFAS leaching have been developed to determine PFAS soil cleanup criteria in the vadose zone. Recent work has presented evidence that while PFAS transport in the vadose zone is governed by several non-equilibrium mechanisms, it is possible to predict PFAS mass flux using equilibrium modelling over month to year timescales. We hypothesized that by quantifying important equilibrium partitioning and hydraulic processes, we could simplify vadose zone leaching models for assessing mass flux from the vadose zone to the underlying groundwater. A mass flux, cell-based model which accounts for important partitioning processes (solid and air-water interfacial partitioning) and transience in hydraulic processes (water flux and water content) was developed and validated herein. Column studies were conducted under simulated rainfall conditions to provide transient hydraulic and PFAS leaching data. A HYDRUS 1-D with PFAS module model was calibrated to the hydraulic conditions of the simulated rainfall columns. Forward simulations were carried out using HYDRUS and the mass balance approximation models. The HYDRUS and mass balance approximations performed nearly identically for all PFAS, and both models predicted PFAS mass leaching within a half order of magnitude of most measured data. These results suggest that readily applicable empirical models and simplified numerical models can reasonably estimate month to year scale mass flux from the vadose zone for sites without major heterogeneity or transport non-ideality considerations.

Fate, distribution, and transport dynamics of Per- and Polyfluoroalkyl Substances (PFASs) in the environment

Per- and Polyfluoroalkyl Substances (PFASs) are persistent organic pollutants with significant environmental and health impacts due to their widespread occurrence, bioaccumulation potential, and resistance to degradation. This paper comprehensively reviews current knowledge of PFAS fate and transport mechanisms by correlating PFAS leaching, retention, and movement to their physicochemical properties and environmental factors based on observing PFAS fate and transport in unsaturated zones, surface water, sediments, plants, and atmosphere. The complex and unique physiochemical properties of PFASs, such as their carbon-fluorine bonds and amphiphilic nature, determine their environmental behavior and persistence. Recent studies emphasize that concentration-dependent affinity coefficients predict the transport of diverse PFAS mixtures by considering the impact of the Air-Water Interface (AWI). These studies highlight the complex interactions that influence PFAS behavior in environmental systems and the need for refined modeling techniques to account for transport dynamics. Competitive adsorption at the AWI, influenced by PFAS physicochemical properties and environmental factors, is crucial. PFAS chain length profoundly affects PFAS volatility and mobility, i.e., longer chains show higher solid matrix adsorption, while shorter chains exhibit greater atmospheric deposition potential. Solution chemistry, encompassing pH and ionic strength, variably alters PFAS sorption behaviors. Mathematical models, such as the Leverett Thermodynamic Model (LTM) and Surface Roughness Multipliers (SRM), effectively predict PFAS retention, offering enhanced accuracy for surface-active solutes through empirical adjustments. Co-contaminants' presence influences the transport behavior of PFASs in the environment. Microbial activity alters PFAS retention, while microplastics, especially polyamide, contribute to their adsorption. These complex interactions govern PFAS fate and transport in the environment. The paper identifies critical gaps in current understanding, including the fate of PFASs, analytical challenges, ecological risk assessment methods, and the influence of episodic events on PFAS transport dynamics. This paper also investigates the research gap in refining current models and experimental approaches to predict PFAS transport accurately and enhance risk mitigation efforts. Addressing these gaps is crucial for advancing remediation strategies and regulatory frameworks to mitigate PFAS contamination effectively.

Forecasting and Hindcasting PFAS Concentrations in Groundwater Discharging to Streams near a PFAS Production Facility

Per- and polyfluoroalkyl substances (PFAS) are known to be highly persistent in groundwater, making it vital to develop new approaches to important practical questions such as the time scale for future persistence of PFAS in contaminated groundwater. In the approach presented here, groundwater from beneath streambeds was analyzed for PFAS and age-dated using SF6 and 3H/3He. The results were coupled with groundwater flux measurements in a convolution approach to estimate past and future PFAS concentrations in groundwater discharge to the streams. At our test site near the Cape Fear River (CFR) of North Carolina, PFAS were detected in groundwater up to 43 years old, suggesting that some PFAS entered groundwater immediately or shortly after fluorochemical production began at the nearby Fayetteville Works. Results are consistent with little to no retardation in groundwater for perfluoroethers such as hexafluoropropylene oxide-dimer acid (HFPO-DA) and perfluoro-2-methoxypropanoic acid (PMPA), the two most abundant PFAS, with mean concentrations of 229 and 498 ng/L, respectively. Future PFAS concentrations in groundwater discharge to streams were estimated to remain above current MCL or health advisory levels through at least 2050 or 2060 (using 3H/3He and SF6, respectively). Recent atmospheric deposition data suggest lower but non-negligible amounts of PFAS may continue to enter groundwater, likely further extending PFAS persistence in groundwater and the adjacent CFR. This approach shows promise for giving an overall perspective on persistence of PFAS in groundwater discharge from a broad contaminated area.

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.

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.

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.