PFAS concentrations in soil versus soil porewater: Mass distributions and the impact of adsorption at air-water interfaces

High-Resolution Depth-Discrete Analysis of PFAS Distribution and Leaching for a Vadose-Zone Source at an AFFF-Impacted Site

The long-term leaching of polyfluoroalkyl substances (PFAS) within the vadose zone of an AFFF application site for which the depth to groundwater is approximately 100 m was investigated by characterizing the vertical distribution of PFAS in a high spatial resolution. The great majority (99%) of PFAS mass resides in the upper 3 m of the vadose zone. The depths to which each PFAS migrated, quantified by moment analysis, is an inverse function of molar volume, demonstrating chromatographic separation. The PFAS were operationally categorized into three chain-length groups based on the three general patterns of retention observed. The longest-chain (>∼335 cm3/mol molar volume) PFAS remained within the uppermost section of the core, exhibiting minimal leaching. Conversely, the shortest-chain (<∼220 cm3/mol) PFAS accumulated at the bottom of the interval, which coincides with the onset of a calcic horizon. PFAS with intermediate-chain lengths were distributed along the length of the core, exhibiting differential magnitudes of leaching. The minimal or differential leaching observed for the longest- and intermediate-chain-length PFAS, respectively, demonstrates that retention processes significantly impacted migration. The accumulation of shorter-chain PFAS at the bottom of the core is hypothesized to result from limited deep infiltration and potential-enhanced retention associated with the calcic horizon.

PFAS occurrence and distribution in yard waste compost indicate potential volatile loss, downward migration, and transformation

We discovered high concentrations of PFAS (18.53 ± 1.5 μg kg-1) in yard waste compost, a compost type widely acceptable to the public. Seventeen out of forty targeted PFAS, belonging to six PFAS classes were detected in yard waste compost, with PFCAs (13.51 ± 0.99 μg kg-1) and PFSAs (4.13 ± 0.19 μg kg-1) being the dominant classes, comprising approximately 72.5% and 22.1% of the total measured PFAS. Both short-chain PFAS, such as PFBA, PFHxA, and PFBS, and long-chain PFAS, such as PFOA and PFOS, were prevalent in all the tested yard waste compost samples. We also discovered the co-occurrence of PFAS with low-density polyethylene (LDPE) and polyethylene terephthalate (PET) plastics. Total PFAS concentrations in LDPE and PET separated from incoming yard waste were 7.41 ± 0.41 μg kg-1 and 1.35 ± 0.1 μg kg-1, which increased to 8.66 ± 0.81 μg kg-1 in LDPE and 5.44 ± 0.56 μg kg-1 in PET separated from compost. An idle mature compost pile revealed a clear vertical distribution of PFAS, with the total PFAS concentrations at the surface level approximately 58.9-63.2% lower than the 2 ft level. This difference might be attributed to the volatile loss of short-chain PFCAs, PFAS's downward movement with moisture, and aerobic transformations of precursor PFAS at the surface.

Modeling 1-D aqueous film forming foam transport through the vadose zone under realistic site and release conditions

Aqueous film forming foams (AFFFs) have been used to extinguish fires since the 1960s, leading to widespread subsurface contamination by per- and polyfluoroalkyl substances (PFAS), an essential component of AFFF. This study presents 1-D simulations of PFAS migration in the vadose zone resulting from AFFF releases. Simulation scenarios used soil profiles from three US Air Force (USAF) installations, encompassing a range of climatic conditions and hydrogeologic environments. A three-component mixture, representative of major constituents of AFFF, facilitated the exploration of competitive and synergistic effects of co-constituents on PFAS migration. To accurately capture unsaturated transport of PFAS in porous media, the model considers (1) surfactant-induced flow, (2) non-linear sorption to the solid phase, (3) competitive accumulation at the air-water interface, and (4) the moisture-dependence of the air-water interfacial area. Defined PFAS releases were consistent with fire training exercises, emergency responses, and accidental spills of record. Simulation results illustrate the importance of hydrogeologic, climatic, geochemical, and AFFF release conditions on PFAS transport and retention. Comparison of field observations and model simulations for Ellsworth AFB indicate that much of the PFOA and PFOS mass is associated with the air-water interface and the solid phase, which limits their migration potential in the vadose zone. Results also show that rates of migration in the aqueous phase are largely controlled by hydrogeologic properties, including recharge rates and hydraulic conductivity. AFFF spill scenarios varying in volume, concentration, and frequency reveal the importance of release characteristics in determining rates of PFAS migration and concentration peaks. Variability is attributed to non-linear sorption processes, where, contrary to simple linear partitioning formulations, transport is strongly affected by the concentration of PFAS species. Simulations also demonstrate the importance of modeling the AFFF as a mixture since competitive interfacial accumulation effects are shown to enhance the mobility of less surface-active PFAS compounds.

Per- and polyfluoroalkyl substance (PFAS) mixtures induce gut microbiota dysbiosis and metabolic disruption in silkworm (Bombyx mori L.)

Mixed legacy and emerging per- and polyfluoroalkyl substances (PFASs) are commonly found in soil and dust; however, the potential toxicity of PFAS mixtures (mPFASs) in insects is unknown. Using 16S rRNA gene sequencing and transcriptome sequencing (RNA-Seq), we evaluated the adverse effects of mPFASs on silkworms, a typical lepidopteran insect. After exposure to mPFASs, the silkworm midgut was enriched with high levels of PFASs, which induced histopathological changes. The composition of the midgut microbiota was significantly affected by mPFAS exposure, and functional predictions revealed significant disruption of some metabolic pathways. RNA-seq analysis revealed that mPFASs significantly changed the transcription profiles. Functional enrichment analysis of the differentially expressed genes also revealed that biological processes related to metabolic pathways and the digestive system were significantly affected, similar to the results of the gut microbiota analysis, suggesting that mPFAS exposure had an adverse effect on the metabolic function of silkworms and may further affect their normal growth. Finally, the significant correlation between abundance changes in the gut microbiota and metabolism/digestion-related genes further highlighted the role of the gut microbiota in mPFAS-related processes affecting the metabolic functions of silkworms. To our knowledge, this study is the first to evaluate the toxic effects of mPFASs in insects and provide basic data for further PFAS toxicity investigations in insects and comprehensive ecological risk assessments of mPFASs.

PFAS release from the subsurface and capillary fringe during managed aquifer recharge

Managed aquifer recharge (MAR) is a sustainable way of harvesting groundwater in water-stressed urbanized areas, where reclaimed wastewater or stormwater is applied on a large basin to infiltrate water into the groundwater aquifer naturally. This process could rapidly fluctuate the water table and move the capillary fringe boundary, and the change in flow dynamic and associated geochemical changes could trigger the release of sequestered pollutants, including per- and polyfluoroalkyl substances (PFAS), also known as 'forever chemicals', from the subsurface and capillary fringe. Yet, the potential of PFAS release from the subsurface and capillary zone during recharge events when the water table rapidly fluctuates has not been evaluated. This study uses laboratory column experiments to simulate PFAS release from pre-contaminated subsurface and capillary fringe during groundwater table fluctuation. The results reveal that the groundwater level fluctuations during MAR increased the release of perfluorobutanesulfonic acid (PFBS) and perfluorooctanesulfonic acid (PFOS) from the capillary fringe, but the fraction released depended on PFAS type and their association with soil colloids. A higher proportion of PFOS in column effluent was found to be associated with particles, while a greater portion of released PFBS was in a free or dissolved state. The direction of water table fluctuation did not affect the release of PFAS in this study. A lack of change in the concentration of bromide, a conservative tracer, during flow interruption, indicates that diffusion of PFAS through reconnected pores during water table rise had an insignificant effect on PFAS release. Overall, this study provides insights into how PFAS can be released from the subsurface and capillary fringe during managed aquifer recharge when the groundwater level is expected to fluctuate quickly.

A systematic investigation of single solute, binary and ternary PFAS transport in water-saturated soil using batch and 1-dimensional column studies: Focus on mixture effects

This study aimed to investigate the transport and release of per- and polyfluoroalkyl substances (PFAS), as single solutes and binary and ternary mixtures, and associated competitive sorption effects in water-saturated soil. Batch sorption isotherm and desorption, and one-dimensional miscible displacement studies were conducted. For the batch study, the mixtures exhibited extensive sorption isotherm nonlinearity at aqueous concentrations exceeding 20 µg/L. At and above this threshold, competitive effects significantly decreased PFAS sorption, mostly affecting perfluorooctanoic acid (PFOA) and perfluorohexane sulfonate (PFHxS). Importantly, mixture effects exacerbated isotherm nonlinearity and may increase the leaching of PFAS in subsurface soil and groundwater. Further, up to 100% desorption occurred for single solutes and mixtures, indicating that the studied PFAS were weakly sorbed. For the column study, at influent concentrations (21 - 27 µg/L, depending on PFAS) near the threshold, PFOA and PFHxS breakthrough curves (BTC) generally exhibited equilibrium (nonlinear) transport, whereas perfluorooctane sulfonate (PFOS) exhibited nonequilibrium transport, with minimal or no mixture effects. Nonequilibrium transport of PFOS was driven by rate-limited sorption, especially as flow interruption tests confirmed the absence of physical nonequilibrium. The sorption distribution coefficients (Kd) from moment and frontal analyses, and 2-site modelling of the BTC, were consistent with the batch-derived Kd, although comparatively smaller. Such discrepancies may limit the applicability of batch-derived Kd values for predictive transport modelling purposes. Overall, understanding mixture impacts may aid effective predictive modelling of PFAS transport and leaching, especially in aqueous film forming foam (AFFF)-source zone areas associated with elevated PFAS concentrations. At low or environmental PFAS concentrations, mixture effects can be expected to be play a minor role in influencing PFAS transport.

Revising the EPA Dilution-Attenuation Soil Screening Model for PFAS

Per and polyfluoroalkyl substances (PFAS) have been shown to be ubiquitous in the environment, and one issue of critical concern is the leaching of PFAS from soil to groundwater. The risk posed by contaminants present in soil is often assessed in terms of the anticipated impact to groundwater through the determination of soil screening levels (SSLs). The U.S. Environmental Protection Agency (EPA) established a soil screening model for determining SSLs. However, the model does not consider the unique retention properties of PFAS and, consequently, the SSLs established with the model may not represent the actual levels that are protective of groundwater quality. The objective of this work is to revise the standard EPA SSL model to reflect the unique properties and associated retention behavior of PFAS. Specifically, the distribution parameter used to convert soil porewater concentrations to soil concentrations is revised to account for adsorption at the air-water interface. Example calculations conducted for PFOS and PFOA illustrate the contrasting SSLs obtained with the revised and standard models. A comparison of distribution parameters calculated for a series of PFAS of different chain length shows that the significance of air-water interfacial adsorption can vary greatly as a function of the specific PFAS. Therefore, the difference between SSLs calculated with the revised versus standard models will vary as a function of the specific PFAS, with greater differences typically observed for longer-chain PFAS. It is anticipated that this revised model will be useful for developing improved SSLs that can be used to enhance site investigations and management for PFAS-impacted sites.

Long-term removal of perfluoroalkyl substances via activated carbon process for general advanced treatment purposes

Granular or powdered activated carbon (GAC/PAC) processes are installed in full-scale drinking water treatment plants (DWTPs) to reduce disinfection byproduct precursors, odor, ammonia, and pesticides. This study investigated the ability of GAC/PAC processes in 23 DWTPs to remove per- and polyfluoroalkyl substances (PFASs). In the GAC process, filter breakthrough of perfluoroalkyl carboxylic acids (PFCAs) occurred faster as the PFCA chain length is decreased. During periods of high water temperatures (20-29 °C), the effluent concentration of two short-chain PFCAs (C4 and C5) surpassed that of the influent after the throughput reached 5,000-7,500 bed volumes (equivalent to 2-3 months) due to desorption. However, such desorption was not observed during periods of low water temperatures (5-19 °C). Meanwhile, long-chain PFCAs were consistently removed, as the GAC was replaced before breakthrough became noticeable. PFAS removal deteriorated at a remarkably fast rate after a partial breakthrough of several tens of percent. Biological activated carbon was proved ineffective in removing PFASs due to its diminished adsorption capacity after long-term use. The PAC process, however, exhibited a slight decrease in PFCA residual (10%) at higher water temperatures (15-30 °C). The PAC dose required for a certain residual ratio was lower with an increase in the hydrophobicity of PFAS; C8-PFCA only required 20 mg/L of PAC for 50% removal, while C4-PFCA required a significantly higher dose of 100-700 mg/L. Consequently, the activated carbon process, which removes organic contaminants in surface water, was inadequate in removing PFASs, particularly those with short chains. Thus, it is recommended that GAC filters be replaced more frequently (within two months) for short-chain PFAS removal. Further, the adsorption performance of PAC must be enhanced.

Determining air-water interfacial areas for the retention and transport of PFAS and other interfacially active solutes in unsaturated porous media

The objective of this work was to determine the methods that produce the most representative measurements and estimations of air-water interfacial area specifically for the retention and transport of PFAS and other interfacially active solutes in unsaturated porous media. Published data sets of air-water interfacial areas obtained with multiple measurement and prediction methods were compared for paired sets of porous media comprising similar median grain diameters, but one with solid-surface roughness (sand) and one without roughness (glass beads). All interfacial areas produced with the multiple diverse methods were coincident for the glass beads, providing validation of the aqueous interfacial tracer-test methods. The results of this and other benchmarking analyses demonstrated that the differences in interfacial areas measured for sands and soil by different methods are not due to errors or artifacts in the methods but rather the result of method-dependent differential contributions of solid-surface roughness. The contributions of roughness to interfacial areas measured by interfacial tracer-test methods were quantified and shown to be consistent with prior theoretical and experiment-based investigations of air-water interface configurations on rough solid surfaces. Three new methods for estimating air-water interfacial areas were developed, one based on the scaling of thermodynamic-determined values and the other two comprising empirical correlations incorporating grain diameter or NBET solid surface area. All three were developed based on measured aqueous interfacial tracer-test data. The three new and three existing estimation methods was tested using independent data sets of PFAS retention and transport. The results showed that the method based on treating air-water interfaces as smooth surfaces as well as the standard thermodynamic method produced inaccurate air-water interfacial areas that failed to reproduce the multiple measured PFAS retention and transport data sets. In contrast, the new estimation methods produced interfacial areas that accurately represented air-water interfacial adsorption of PFAS and associated retention and transport. The measurement and estimation of air-water interfacial areas for field-scale applications is discussed in light of these results.

Immobilization of per- and polyfluoroalkyl substances (PFAS): Comparison of leaching behavior by three different leaching tests

The evaluation of PFAS immobilization performance in laboratory experiments, especially the long-term stability, is a challenge. To contribute to the development of adequate experimental procedures, the impact of experimental conditions on the leaching behavior was studied. Three experiments on different scales were compared: batch, saturated column, and variably saturated laboratory lysimeter experiments. The Infinite Sink (IS) test - a batch test with repeated sampling - was applied for PFAS for the first time. Soil from an agricultural field amended with paper-fiber biosolids polluted with various perfluoroalkyl acids (PFAAs; 655 μg/kg ∑¹⁸PFAAs) and polyfluorinated precursors (1.4 mg/kg ∑¹⁸precursors) was used as the primary material (N-1). Two types of PFAS immobilization agents were tested: treatment with activated carbon-based additives (soil mixtures: R-1 and R-2), and solidification with cement and bentonite (R-3). In all experiments, a chain-length dependent immobilization efficacy is observed. In R-3, the leaching of short-chain PFAAs was enhanced relative to N-1. In column and lysimeter experiments with R-1 and R-2, delayed breakthrough of short-chain PFAAs (C4) occurred (> 90 days; in column experiments at liquid-to-solid ratio (LS) > 30 L/kg) with similar temporal leaching rates suggesting that leaching in these cases was a kinetically controlled process. Observed differences between column and lysimeter experiments may be attributed to varying saturation conditions. In IS experiments, PFAS desorption from N-1, R-1, and R-2 is higher than in the column experiments (N-1: +44 %; R-1: +280 %; R-2: +162 %), desorption of short-chain PFAS occurred predominantly in the initial phase (< 14 days). Our findings demonstrate that sufficient operating times are essential in percolation experiments, e.g., in column experiments >100 days and LS > 30 L/kg. IS experiments may provide a faster estimate for nonpermanent immobilization. The comparison of experimental data from various experiments is beneficial to evaluate PFAS immobilization and to interpret leaching characteristics.

Predicting Interfacial Tension and Adsorption at Fluid-Fluid Interfaces for Mixtures of PFAS and/or Hydrocarbon Surfactants

Many per- and polyfluoroalkyl substances (PFAS) are surface-active and adsorb at fluid-fluid interfaces. The interfacial adsorption controls PFAS transport in multiple environmental systems, including leaching through soils, accumulation in aerosols, and treatment methods such as foam fractionation. Most PFAS contamination sites comprise mixtures of PFAS as well as hydrocarbon surfactants, which complicates their adsorption behaviors. We present a mathematical model for predicting interfacial tension and adsorption at fluid-fluid interfaces for multicomponent PFAS and hydrocarbon surfactants. The model is derived from simplifying a prior advanced thermodynamic-based model and applies to nonionic and ionic mixtures of the same charge sign with swamping electrolytes. The only required model inputs are the single-component Szyszkowski parameters obtained for the individual components. We validate the model using literature interfacial tension data of air-water and NAPL (non-aqueous phase liquid)-water interfaces covering a wide range of multicomponent PFAS and hydrocarbon surfactants. Application of the model to representative porewater PFAS concentrations in the vadose zone suggests competitive adsorption can significantly reduce PFAS retention (up to 7 times) at some highly contaminated sites. The multicomponent model can be readily incorporated into transport models to simulate the migration of mixtures of PFAS and/or hydrocarbon surfactants in the environment.

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.

Centurial Persistence of Forever Chemicals at Military Fire Training Sites

Drinking water contamination by per- and polyfluoroalkyl substances (PFAS) is widespread near more than 300 United States (U.S.) military bases that used aqueous film-forming foams (AFFF) for fire training and firefighting activities. Much of the PFAS at these sites consist of precursors that can transform into terminal compounds of known health concern but are omitted from standard analytical methods. Here, we estimate the expected duration and contribution of precursor biotransformation to groundwater PFAS contamination at an AFFF-contaminated military base on Cape Cod, Massachusetts, United States, by optimizing a geochemical box model using measured PFAS concentrations from a multidecadal time series of groundwater and a soil survey in the source zone. A toolbox of analytical techniques used to reconstruct the mass budget of PFAS showed that precursors accounted for 46 ± 8% of the extractable organofluorine (a proxy for total PFAS) across years. Terminal PFAS still exceed regulatory limits by 2000-fold decades after AFFF use ceased. Measurements and numerical modeling show that sulfonamido precursors are retained in the vadose zone and their slow biotransformation into perfluoroalkyl sulfonates (half-life > 66 yr) sustains groundwater concentrations of perfluorobutane sulfonate (PFBS) and perfluorohexane sulfonate (PFHxS). The estimated PFAS reservoir in the vadose zone and modeled flux into groundwater suggest PFAS contamination above regulatory guidelines will persist for centuries without remediation.

Predicting Concentration- and Ionic-Strength-Dependent Air-Water Interfacial Partitioning Parameters of PFASs Using Quantitative Structure-Property Relationships (QSPRs)

Air-water interfacial retention of poly- and perfluoroalkyl substances (PFASs) is increasingly recognized as an important environmental process. Herein, column transport experiments were used to measure air-water interfacial partitioning values for several perfluoroalkyl ethers and for PFASs derived from aqueous film-forming foam, while batch experiments were used to determine equilibrium K_ia data for compounds exhibiting evidence of rate-limited partitioning. Experimental results suggest a Freundlich isotherm best describes PFAS air-water partitioning at environmentally relevant concentrations (10¹-10⁶ ng/L). A multiparameter regression analysis for K_ia prediction was performed for the 15 PFASs for which equilibrium K_ia values were determined, assessing 246 possible combinations of 8 physicochemical and system properties. Quantitative structure-property relationships (QSPRs) based on three to four parameters provided predictions of high accuracy without model overparameterization. Two QSPRs (R² values of 0.92 and 0.83) were developed using an assumed average Freundlich n value of 0.65 and validated across a range of relevant concentrations for perfluorooctane sulfonate (PFOS), perfluorooctanoate (PFOA), and hexafluoropropylene oxide-dimer acid (i.e., GenX). A mass action model was further modified to account for the changing ionic strength on PFAS air-water interfacial sorption. The final result was two distinct QSPRs for estimating PFAS air-water interfacial partitioning across a range of aqueous concentrations and ionic strengths.

Assessment of PFAS in collocated soil and porewater samples at an AFFF-impacted source zone: Field-scale validation of suction lysimeters

Per- and polyfluoroalkyl substances (PFAS) measurable in soil porewater authoritatively represent the mobile mass fraction critical to accurate assessment of leaching from source zones. This study evaluated PFAS occurrence in lysimeter-collected porewater samples for two depth intervals at a decades-old aqueous film-forming foam (AFFF)-impacted field site quarterly for a year. Notably, site-wide Log₁₀ (∑PFAS) concentrations did not significantly differ among sampling events despite highly variable sample yields due to a heterogeneous and dynamic soil moisture regime. However, Log₁₀ (∑PFAS) concentrations were significantly higher in the shallow interval concordant with higher mean soil concentrations and higher total organic carbon (TOC) reflecting net retention, which is supported by soil-to-groundwater annual mass discharge estimates less than 0.2% of the total source mass for any given PFAS. Interestingly, PFAS-specific Log₁₀ (soil-to-porewater ratios) significantly increased with soil concentration in both depth intervals contrary to concentration dependence resulting from the saturation of sorption sites potentially implicating self-assembly as an additional operative retention mechanism. Overall, these data validate the use of suction lysimeters for short-term site characterization deployments and emphasize the importance of in situ porewater samples for interrogating PFAS transport within source zones.

Model-based identification of vadose zone controls on PFAS mobility under semi-arid climate conditions

Contamination through per-and poly-fluoroalkyl substances (PFAS) have occurred globally in soil and groundwater systems at military, airport and industrial sites due to the often decades-long periodic application of firefighting foams. At PFAS contaminated sites, the unsaturated soil horizon often serves as a long-term source for sustained PFAS contamination for both groundwater and surface water runoff. An understanding of the processes controlling future mass loading rates to the saturated zone from these source zones is imperative to design efficient remediation measures. In the present study, hydrochemical data from a site where PFAS transport was observed as a result of the decades-long application of AFFF were used to develop and evaluate conceptual and numerical models that determine PFAS mobility across the vadose zone under realistic field-scale conditions. The simulation results demonstrate that the climate-driven physical flow processes within the vadose zone exert a dominating control on the retention of PFAS. Prolonged periods of evapotranspiration exceeding rainfall under the semi-arid conditions trigger periods of upward flux and evapoconcentration, leading to the observed persistence of PFAS compounds in the upper ca. 2 metres of the vadose zone, despite cessation of AFFF application to soils since more than a decade. Physico-chemical retention mechanisms, namely sorption to the air-water interface (AWI) and sediment surfaces, contribute further to PFAS retention. The simulations demonstrate how PFAS downward transport is effectively confined to short periods following discrete rain events when soils display a high degree of saturation. During these periods, AWI sorption is at a minimum. In addition, high PFAS concentrations measured and simulated below the source zone reduce the effect of the AWI further due to a decrease in surface tension associated with elevated PFAS concentrations. Consequently, time-integrated PFAS migration and retardation illuminates that the field-relevant PFAS transport rates are predominantly controlled by the physical flow processes with a lower relative importance of AWI and sediment sorption adding to PFAS retention.

Global distributions, source-type dependencies, and concentration ranges of per- and polyfluoroalkyl substances in groundwater

A meta-analysis was conducted of published literature reporting concentrations of per- and polyfluoroalkyl substances (PFAS) in groundwater for sites distributed in 20 countries across the globe. Data for >35 PFAS were aggregated from 96 reports published from 1999 to 2021. The final data set comprises approximately 21,000 data points after removal of time-series and duplicate samples as well as non-detects. The reported concentrations range over many orders of magnitude, from ng/L to mg/L levels. Distinct differences in concentration ranges are observed between sites located within or near sources versus those that are not. Perfluorooctanoic acid (PFOA), ranging from <0.03 ng/L to ~7 mg/L, and perfluorooctanesulfonic acid (PFOS), ranging from 0.01 ng/L to ~5 mg/L, were the two most reported PFAS. The highest PFAS concentration in groundwater is ~15 mg/L reported for the replacement-PFAS 6:2 fluorotelomer sulfonate (6:2 FTS). Maximum reported groundwater concentrations for PFOA and PFOS were compared to concentrations reported for soils, surface waters, marine waters, and precipitation. Soil concentrations are generally significantly higher than those reported for the other media. This accrues to soil being the primary entry point for PFAS release into the environment for many sites, as well as the generally significantly greater retention capacity of soil compared to the other media. The presence of PFAS has been reported for all media in all regions tested, including areas that are far removed from specific PFAS sources. This gives rise to the existence of a "background" concentration of PFAS that must be accounted for in both regional and site-specific risk assessments. The presence of this background is a reflection of the large-scale use of PFAS, their general recalcitrance, and the action of long-range transport processes that distribute PFAS across regional and global scales. This ubiquitous distribution has the potential to significantly impact the quality and availability of water resources in many regions. In addition, the pervasive presence of PFAS in the environment engenders concerns for impacts to ecosystem and human health.