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

Per- and polyfluoroalkyl substances (PFAS) as contaminants in groundwater resources - A comprehensive review of subsurface transport processes

Per- and polyfluorinated alkyl substances (PFAS) are persistent contaminants in the environment. An increased awareness of adverse health effects related to PFAS has further led to stricter regulations for several of these substances in e.g. drinking water in many countries. Groundwater constitutes an important source of raw water for drinking water production. A thorough understanding of PFAS subsurface fate and transport mechanisms leading to contamination of groundwater resources is therefore essential for management of raw water resources. A review of scientific literature on the subject of processes affecting subsurface PFAS fate and transport was carried out. This article compiles the current knowledge of such processes, mainly focusing on perfluoroalkyl acids (PFAA), in soil- and groundwater systems. Further, a compilation of data on transport parameters such as solubility and distribution coefficients, as well as, insight gained and conclusions drawn from the reviewed material are presented. As the use of certain fire-fighting foams has been identified as the major source of groundwater contamination in many countries, research related to this type of pollution source has been given extra focus. Uptake of PFAS in biota is outside the scope of this review. The review showed a large spread in the magnitude of distribution coefficients and solubility for individual PFAS. Also, it is clear that the influence of multiple factors makes site-specific evaluation of distribution coefficients valuable. This article aims at giving the reader a comprehensive overview of the subject, and providing a base for further work.

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.

The occurrence, spatial distribution, and well-depth dependence of PFASs in groundwater from a reclaimed water irrigation area

With the growing development of modern agriculture and industry, groundwater is facing more and more complex contaminants. One such contaminant is per- and polyfluoroalkyl substances (PFASs), which pose a potential risk to human health, particularly for those who rely on groundwater as their primary source of drinking water. In this study, we conducted a comprehensive investigation on the occurrence, spatial distribution, and source apportionment of PFASs in shallow (<60 m) and deep (>80 m) groundwater samples from a reclaimed water irrigation area in Beijing's suburbs. Our results showed that the average total PFAS concentration (∑10PFAS) for all samples was 10.55 ± 7.77 ng/L, ranging from 1.05 to 34.28 ng/L. The dominant congeners were PFBA, PFOA, and PFBS. No significant linear relationship was observed between PFAS concentrations and the well depth. However, the averaged ΣPFASs in groundwater were highest in the uppermost layer and declined sharply to a few ng/L in the deep aquifer below 80 m. PFASs showed elevated concentration in shallow aquifers in 9 out of 11 paired wells, indicating an overall descending trend of PFASs with increasing aquifer depth. The spatial distribution of PFASs was highly heterogeneous and showed different patterns in shallow and deep groundwater, which may be related to the complicated attenuation behavior of PFAS compounds when they transport and diffuse through overlapping aquifer layers. The influence of the landfill on groundwater PFASs was most pronounced within a 5 km radius. Source apportionment results indicated that reclaimed water irrigation is the main non-point source of PFASs in shallow groundwater. In contrast, deep groundwater is primarily subject to point sources and lateral recharge flow. This investigation of PFASs in shallow and deep wells provides a foundation for further exploration of PFASs transportation and risk prevention in regions where groundwater is a major water resource for domestic and industrial development.

Long-term leaching through clayey till of N,N-dimethylsulfamide, a Persistent and Mobile Organic Compound (PMOC)

Environmental pollution with Persistent and Mobile Organic Compounds (PMOC) from anthropogenic activities is an increasing cause for concern. These compounds are readily leached to groundwater aquifers and are likely to resist degradation, putting pressure on groundwater resources. Pesticides can form PMOCs upon degradation in the environment. The PMOC N,N-dimethylsulfamide (DMS) was the most frequently detected pesticide metabolite in Danish drinking water wells in 2020, although the pesticidal use of the last parent compound (tolylfluanid) ended in 2007. This study aimed to improve the understanding of the leaching of the PMOC DMS from clayey tills by combining a review of compound properties, sources and use, comprehensive field observations and numerical flow and solute transport modeling. The modeling explored the mechanisms of DMS retention during vertical transport in clayey till and the fingerprint in the underlying aquifer. The results were supported by detailed field observations at an agricultural site with strawberry production. Porewater samples were collected from clayey till to a depth of 12 m bgs by a custom designed installation method of suction cups. Groundwater sampling (249 samples) was designed to provide vertical concentration profiles at various distances from the presumed sources. The review of properties showed that the parent compounds and intermediates degrade quickly in topsoil, releasing the highly persistent and mobile DMS. We tested the effect of fractures on transport with different hydraulic apertures and a scenario without fractures by numerical modeling. The results showed that the presence of fractures can smooth the breakthrough curve below the clayey till, leading to faster breakthrough, lower maximum concentration, and several decades of prolonged leaching in simulations with the largest aperture (20 μm). The fracture-matrix interaction is a possible explanation for the observed delay of leaching from clayey till. The vertical concentration profiles in groundwater were used for identifying the sources at the field site and testing source strengths. Assigning one point source (200 μg/L) and two diffuse sources (40-50 μg/L) to the model produced vertical concentration profiles that compared well with observed field data in clayey till and the aquifer. All results were integrated into a conceptual model for the environmental fate of PMOCs in soil and groundwater. The findings of this study imply that the presence of fractures in clayey till should be considered in conceptual site models, since they can substantially prolong the leaching of PMOCs to groundwater. The integration of comprehensive field investigations and numerical modeling is key to understand the fate of PMOCs in complex field systems with different source types. Together with widespread occurrences of PMOCs in groundwater systems, the results highlight the need for improved approval procedures for pesticides and biocides which considers their persistent and mobile metabolites.

Watershed scale PFAS fate and transport model for source identification and management implications

Developing strategic plans for the remediation and mitigation of pre- and polyfluoroalkyl substances (PFAS) in soil, groundwater, and surface water requires an understanding of the fate and transport of these chemicals on a regional scale. To fill this knowledge gap, we developed a distributed hydrogeochemical model and applied it to a large-scale watershed with various point and non-point sources of a long-chain, highly persistent PFAS compound known as perfluorooctane sulfonic acid (PFOS). The results showed that the developed model could reproduce the spatiotemporal concentration of PFOS across a large and diverse watershed. Herein, our first objective was to quantify the PFOS transport from the unsaturated zone to the groundwater and surface water via leaching, surface runoff, lateral flow, and sediment transport. The second objective was to identify factors influencing PFOS release from confirmed and suspected PFAS sites and urban and agricultural areas. The modeling results show that surface runoff played a significant role in PFOS transport, with urban areas and industrial sites being major contributors. In addition, sediment transport was found to be a notable pathway for PFOS release, particularly from sites with biosolids application. Further analysis revealed the relative importance of topography, soil water retention, and water-solid adsorption factors in determining PFOS transport dynamics at the watershed scale for better source identification and targeted management.

In Situ Insight into the Availability and Desorption Kinetics of Per- and Polyfluoroalkyl Substances in Soils with Diffusive Gradients in Thin Films

The physicochemical exchange dynamics between the solid and solution phases of per- and polyfluoroalkyl substances (PFAS) in soils needs to be better understood. This study employed an in situ tool, diffusive gradients in thin films (DGT), to understand the distribution and exchange kinetics of five typical PFAS in four soils. Results show a nonlinear relationship between the PFAS masses in DGT and time, implying that PFAS were partially supplied by the solid phase in all of the soils. A dynamic model DGT-induced fluxes in soils/sediments (DIFS) was used to interpret the results and derive the distribution coefficients for the labile fraction (K_dl), response time (t_c), and adsorption/desorption rates (k₁ and k_-1). The larger labile pool size (indicated by K_dl) for the longer chain PFAS implies their higher potential availability. The shorter chain PFAS tend to have a larger t_c and relatively smaller k_-1, implying that the release of these PFAS in soils might be kinetically limited but not for more hydrophobic compounds, such as perfluorooctanesulfonic acid (PFOS), although soil properties might play an important role. K_dl ultimately controls the PFAS availability in soils, while the PFAS release from soils might be kinetically constrained (which may also hold for biota uptake), particularly for more hydrophilic PFAS.

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.

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.