Nonideal Transport and Extended Elution Tailing of PFOS in Soil

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.

Bio-accumulation and health risk assessments of per- and polyfluoroalkyl substances in wheat grains

Per- and polyfluoroalkyl substances (PFASs) have been widely detected in various food, which has attracted worldwide concern. However, the factors influencing the transfer and bio-accumulation of PFASs from soils to wheat in normal farmland, is still ambiguous. We investigated the PFASs accumulation in agricultural soils and grains from 10 cites, China, and evaluated the health risks of PFASs via wheat consumption. Our results show that ∑PFASs in soils range from 0.34 μg/kg to 1.59 μg/kg with PFOA and PFOS dominating, whilst ∑PFASs in wheats range from 2.74 to 6.01 μg/kg with PFOA, PFBA and PFHxS dominating. The lower pH conditions and high total organic carbon (TOC) could result in the higher accumulation of PFASs in soils and subsequently in wheat grains, whilst the bioaccumulation factors of PFASs increase with increasing pH conditions but not with TOC. The estimated daily intake (EDI) values of PFBA, PFOA, and PFHxS are relatively high, but data supports that ingesting wheat grains does not result in any potential risk to the human beings. Our studies provided more information about PFASs accumulation in wheat grains, and help us understand the current potential risks of PFASs in food.

Air-water interfacial collapse and rate-limited solid desorption control Perfluoroalkyl acid leaching from the vadose zone

Some Per- and polyfluoroalkyl substances (PFAS) are strongly retained in the vadose zone due to their sorption to both soils and air-water interfaces. While significant research has been dedicated to understanding equilibrium behavior for these multi-phase retention processes, leaching and desorption from aqueous film-forming foam (AFFF) impacted soils under field relevant conditions can exhibit significant deviations from equilibrium. Herein, laboratory column studies using field collected AFFF-impacted soils were employed to examine the leaching of perfluoroalkyl acids (PFAAs) under simulated rainfall conditions. The HYDRUS 1-D model was calibrated to estimate the unsaturated hydraulic properties of the soil in a layered system using multiple boundary condtions. Forward simulations of equilibrium PFAS partitioning using the HYDRUS model and simplified mass balance calculations showed good agreement with the net PFAS mass flux out of the column. However, neither were able to predict the PFAS concentrations in the leached porewater. To better understand the mechanisms controlling the leaching behavior, the HYDRUS 1-D two-site leaching model incorporating solid phase rate limitation and equilibrium air-water interfacial partitioning was employed. Three variations of the novel model incorporating different forms of equilibrium air-water interfacial partitioning were considered using built-in numerical inversion. Results of numerical inversion show that a combination of air-water interfacial collapse and rate-limited desorption from soils can better predict the unique leaching behavior exhibited by PFAAs in AFFF-impacted soils. A sensitivity analysis of the initial conditions and rate-limited desorption terms was conducted to assess the agreement of the model with measured data. The models demonstrated herein show that, under some circumstances, laboratory equilibrium partitioning data can provide a reasonable estimation of total mass leaching, but fail to account for the significant rate-limited, non-Fickian transport which affect PFAA leaching to groundwater in unsaturated soils.

PFAS Porewater concentrations in unsaturated soil: Field and laboratory comparisons inform on PFAS accumulation at air-water interfaces

Poly- and perfluoroalkyl substance (PFAS) leaching from unsaturated soils impacted with aqueous film-forming foams (AFFFs) is an environmental challenge that remains difficult to measure and predict. Complicating measurements and predictions of this process is a lack of understanding between the PFAS concentrations measured in a collected environmental unsaturated soil sample, and the PFAS concentrations measured in the corresponding porewater using field-deployed lysimeters. The applicability of bench-scale batch testing to assess this relationship also remains uncertain. In this study, field-deployed porous cup suction lysimeters were used to measure PFAS porewater concentrations in unsaturated soils at 5 AFFF-impacted sites. Field-measured PFAS porewater concentrations were compared to those measured in porewater extracted in the laboratory from collected unsaturated soil cores, and from PFAS concentrations measured in the laboratory using batch soil slurries. Results showed that, despite several years since the last AFFF release at most of the test sites, precursors were abundant in 3 out of the 5 sites. Comparison of field lysimeter results to laboratory testing suggested that the local equilibrium assumption was valid for at least 3 of the sites and conditions of this study. Surprisingly, PFAS accumulation at the air-water interface was orders of magnitude less than expected at two of the test sites, suggesting potential gaps in the understanding of PFAS accumulation at the air-water interface at AFFF-impacted sites. Finally, results herein suggest that bench-scale testing on unsaturated soils can in some cases be used to inform on PFAS in situ porewater concentrations.

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.

An integrated analytical modeling framework for determining site-specific soil screening levels for PFAS

Soils at many contaminated sites have accumulated a significant amount of per- and polyfluoroalkyl substances (PFAS) and may require remediation to mitigate leaching to groundwater. USEPA's current approaches for determining soil screening levels (SSLs) were developed for non-PFAS contaminants. Because many PFAS are interfacially-active with unique leaching behaviors in soils, the current non-PFAS-specific soil screening models may not be applicable. Following USEPA's general methodology, we develop a new modeling framework representing PFAS-specific transport processes for determining site-specific SSLs for PFAS-contaminated sites. We couple a process-based analytical model for PFAS leaching in the vadose zone and a dilution factor model for groundwater in an integrated framework. We apply the new modeling framework to two typical types of contaminated sites. Comparisons with the standard USEPA SSL approach suggest that accounting for the PFAS-specific transport processes may significantly increase the SSL for some PFAS. For the range of soil properties and groundwater recharge rates examined, while SSLs determined with the new model are less than a factor of 2 different from the standard-model values for less interfacially-active shorter-chain PFAS, they are up to two orders of magnitudes greater for more interfacially-active longer-chain PFAS. The new analytical modeling framework provides an effective tool for deriving more accurate site-specific SSLs and improving site characterization and remedial efforts at PFAS-contaminated sites.

The influence of NAPL distribution on the transport of PFOS in Co-contaminated media

The objective of this research was to examine the influence of nonuniform NAPL distribution and hydraulic-conductivity heterogeneity on the retention and transport of PFAS. Two-dimensional flow cell experiments were conducted to investigate the impact of NAPL configuration on PFOS retention. Two simplified NAPL configurations were used, one with decane present as residual and the other with greater-than-residual (pool) NAPL saturation. The measured and predicted retardation factors were 1.7 and 1.6 for the residual-NAPL flow cell. The fractions of retention associated with solid-phase sorption, NAPL-water interfacial adsorption, and NAPL absorption were 0.56, 0.43, and 0.1, respectively. The measured retardation factor was 1.4 for the pool-NAPL flow cell. Notably, the predicted retardation factor assuming no contribution from NAPL-water interfacial adsorption was 1.5. Comparison of the two values indicates a minimal contribution of NAPL-water interfacial adsorption to transport, which is attributed to the impact of by-pass flow effects that limited contact of PFOS with the NAPL present within the pool. Mathematical-model simulations were conducted to examine the impact of heterogeneity on PFOS transport in model NAPL-contaminated heterogeneous source zones. The results demonstrated that the impact of NAPL-water interface adsorption on PFAS transport in NAPL-contaminated source zones can be influenced by the magnitude of hydraulic-conductivity heterogeneity.

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.

Investigating rate-limited sorption, sorption to air-water interfaces, and colloid-facilitated transport during PFAS leaching

Various sorption processes affect leaching of per- and polyfluoroalkyl substances (PFAS) such as PFOA and PFOS. The objectives of this study are to (1) compare rate-limited leaching in column and lysimeter experiments, (2) investigate the relevance of sorption to air-water interfaces (AWI), and (3) examine colloid-facilitated transport as a process explaining early experimental breakthrough. A continuum model (CM) with two-domain sorption is used to simulate equilibrium and rate-limited sorption. A random walk particle tracking (PT) model was developed and applied to analyze complex leaching characteristics. Results show that sorption parameters derived from column experiments underestimate long-term PFOA leaching in lysimeter experiments due to early depletion, suggesting that transformation of precursors contributes to the observed long-term leaching in the lysimeters (approximately 0.003 µg/kg/d PFOA). Both models demonstrate that sorption to AWI is the dominant retention mechanism for PFOS in lysimeter experiments, with retardation due to AWI being 3 (CM) to 3.7 (PT) times higher than retardation due to solid phase sorption. Notably, despite a simplified conception of AWI sorption, the PT results are closer to the observations. The PT simulations demonstrate possible colloid-facilitated transport at early time; however, results using substance-specific varying transport parameters align better with the observations, which should be equal if colloid-facilitated transport without additional kinetics is the sole mechanism affecting early breakthrough. Possibly, rate-limited sorption to AWI is relevant during the early stages of the lysimeter experiment. Our findings demonstrate that rate-limited sorption is less relevant for long-term leaching under field conditions compared to transformation of precursors and that sorption to AWI can be the dominant retention mechanism on contaminated sites. Moreover, they highlight the potential of random walk particle tracking as a practical alternative to continuum models for estimating the relative contributions of various retention mechanisms.

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.

Transport behavior difference and transport model of long- and short-chain per- and polyfluoroalkyl substances in underground environmental media: A review

Perfluoroalkyl carboxylic acids (PFCAs) and perfluoroalkyl sulfonates (PFSAs), which are the most commonly regulated and most widely concerned per- and polyfluoroalkyl substances (PFAS) have received increasing attention on a global scale due to their amphiphilicity, stability, and long-range transport. Thus, understanding the typical PFAS transport behavior and using models to predict the evolution of PFAS contamination plumes is important for evaluating the potential risks. In this study, the effects of organic matter (OM), minerals, water saturation, and solution chemistry on the transport and retention of PFAS were investigated, and the interaction mechanism between long-chain/short-chain PFAS and the surrounding environment was analyzed. The results revealed that high content of OM/minerals, low saturation, low pH, and divalent cation had a great retardation effect on long-chain PFAS transport. The retention caused by hydrophobic interaction was the prominent mechanism for long-chain PFAS, whereas, the retention caused by electrostatic interaction was more relevant for short-chain PFAS. Additional adsorption at the air-water and nonaqueous-phase liquids (NAPL)-water interface was another potential interaction for retarding PFAS transport in the unsaturated media, which preferred to retard long-chain PFAS. Furthermore, the developing models for describing PFAS transport were investigated and summarized in detail, including the convection-dispersion equation, two-site model (TSM), continuous-distribution multi-rate model, modified-TSM, multi-process mass-transfer (MPMT) model, MPMT-1D model, MPMT-3D model, tempered one-sided stable density transport model, and a comprehensive compartment model. The research revealed PFAS transport mechanisms and provided the model tools, which supported the theoretical basis for the practical prediction of the evolution of PFAS contamination plumes.

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.

Role of Mineral-Organic Interactions in PFAS Retention by AFFF-Impacted Soil

A comprehensive, generalized approach to predict the retention of per- and polyfluoroalkyl substances (PFAS) from aqueous film-forming foam (AFFF) by a soil matrix as a function of PFAS molecular and soil physiochemical properties was developed. An AFFF with 34 major PFAS (12 anions and 22 zwitterions) was added to uncontaminated soil in one-dimensional saturated column experiments and PFAS mass retained was measured. PFAS mass retention was described using an exhaustive statistical approach to generate a poly-parameter quantitative structure-property relationship (ppQSPR). The relevant predictive properties were PFAS molar mass, mass fluorine, number of nitrogens in the PFAS molecule, poorly crystalline Fe oxides, organic carbon, and specific (BET-N2) surface area. The retention of anionic PFAS was nearly independent of soil properties and largely a function of molecular hydrophobicity, with the size of the fluorinated side chain as the main predictor. Retention of nitrogen-containing zwitterionic PFAS was related to poorly crystalline metal oxides and organic carbon content. Knowledge of the extent to which a suite of PFAS may respond to variations in soil matrix properties, as developed here, paves the way for the development of reactive transport algorithms with the ability to capture PFAS dynamics in source zones over extended time frames.

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.

Leaching of Perfluoroalkyl Acids during Unsaturated Zone Flushing at a Field Site Impacted with Aqueous Film Forming Foam

While several studies have focused on perfluoroalkyl acid (PFAA) leaching from soils, field studies evaluating the relationship between PFAA mass removal and porewater concentrations as the PFAA source becomes depleted are lacking. Herein, in situ water flushing was performed at a site historically impacted with AFFF to accelerate the leaching of PFAAs from unsaturated soils in a highly characterized field test cell. Porous cup suction lysimeters were used to assess the changes in PFAA porewater concentrations as a function of PFAA mass removal from the unsaturated soils, where flushing was intermittently paused to determine ambient PFAA porewater concentrations. Results showed that the fractional decreases in PFAA porewater concentrations during flushing exceeded the fractional decrease in PFAA mass removal from the soil. PFOS porewater concentrations decrease by 76% (with negligible rebound) compared to only a 7.4% decrease in overall PFOS mass removed from the unsaturated zone. Overall, the results observed herein suggest that, when considering soil impacts to groundwater, less stringent soil cleanup criteria than those that consider an equivalent relationship between mass removal and mass discharge may be appropriate. In addition, remedial approaches that remove only a modest fraction of the PFAA soil mass may be protective of underlying groundwater, particularly for perfluorinated sulfonates with at least six carbons.

Adsorption of PFAAs in the Vadose Zone and Implications for Long-Term Groundwater Contamination

Perfluoroalkyl acids (PFAAs) are persistent environmental contaminants that sorb to air-water and solid interfaces throughout the vadose zone. These sorption processes lead to decadal leaching of PFAS from the source zones to groundwater systems. While these processes are increasingly well understood, critical gaps exist in describing the vertically variable adsorption in the presence of vadose zone heterogeneity and methods for efficiently upscaling the laboratory observations to predict field-scale PFAA transport and retardation. In this work, we build upon fundamental theories and scalable relationships to define a semi-analytical framework for synthesizing and upscaling PFAA adsorption in heterogeneous vadose zone systems. Solid-phase and air-water interfacial adsorption are quantified mechanistically for several PFAAs and then applied to a contaminated site in Northern Wisconsin. The results highlight the dominance of air-water and organic carbon solid-phase adsorption processes in the vadose zone. Strong sorption heterogeneity─driven by depth-dependent adsorption mechanisms─produces complex spatially variable retardation profiles. We develop vadose zone retardation potentials to quantify this field-scale heterogeneity and propose vertical integration methods to upscale spatially resolved information for transport modeling. This work highlights the importance of accounting for multiscale and multiprocess heterogeneity for accurately describing and predicting the long-term fate and transport of PFAAs in the subsurface.

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.

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.

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

Determining the risk posed by PFAS leaching from soil to groundwater requires quantification of the magnitude and temporal/spatial variability of PFAS mass discharge from the vadose zone, which is governed in part by the concentrations of PFAS in soil porewater. Porewater concentrations are impacted and mediated by the properties of the PFAS and soil, multiple transport and fate processes, and site conditions. The objective of this research was to delineate the relationship between soil porewater concentrations and soil concentrations, based on a comprehensive model of PFAS mass distribution within a soil sample volume. Measured parameters representing solid-phase sorption and air-water interfacial adsorption are used to illustrate the impact of soil and PFAS properties on the distribution of representative PFAS between soil and soil porewater. Literature data reported for soil and soil porewater concentrations of several PFAS obtained from outdoor lysimeter experiments are used to test the distribution model. Soil-to-porewater concentration ratios predicted with the model compared very well to the measured concentration ratios. The nondimensional distribution coefficient that describes the distribution of PFAS mass amongst all domains within a soil sample was observed to be a function of PFAS molecular size. Numerical simulations conducted for a model fire-training source area were used to illustrate the ranges in magnitude of soil versus porewater concentrations for representative field conditions. The results of the measured and simulated data sets demonstrated the importance of air-water interfacial adsorption for the distribution of the longer-chain PFAS within soil samples. PFAS soil porewater concentrations are anticipated to range from ng/L to mg/L depending upon soil concentrations, which in turn depend upon the nature of the site.

Air-water interfacial adsorption of C4-C10 perfluorocarboxylic acids during transport in unsaturated porous media

The impact of chain length on air-water interfacial adsorption of perfluorocarboxylic acids (PFCAs) during transport in unsaturated quartz sand was investigated. Short-chain (C4-C7: PFBA, PFPeA, PFHxA, PFHpA) and long chain (C8-C10: PFOA, PFNA, PFDA) PFCAs were selected as a representative homologous series. Surface tensions were measured to characterize surface activities of the selected PFCAs. Miscible-displacement column experiments were conducted for each of the PFCAs to characterize the magnitudes of air-water interfacial adsorption under transport conditions. The transport of the long-chain PFCAs exhibited greater retardation than the short-chain PFCAs. Air-water interfacial adsorption (AWIA) was the predominant source of retention (≥63%) for the long-chain PFCAs. Conversely, AWIA contributed less to retention than did solid-phase sorption for the short-chain PFCAs, with the former contributions ranging from 4% to 40%. Direct examination of the breakthrough-curve profiles as well as mathematical-modeling results demonstrated that transport of the two longest-chain PFCAs was influenced by nonlinear AWIA, whereas that of the shorter-chain PFCAs was not. This disparate behavior is consistent with the input concentration used for the transport experiments in comparison to the respective surface activities and critical reference concentrations of the different PFCAs. Quantitative-structure/property-relationship (QSPR) analysis was applied to characterize the influence of molecular size on air-water interfacial adsorption. The logs of the air-water interfacial adsorption coefficients (K_ia) determined from the transport experiments are a monotonic function of molar volume, consistent with prior QSPR analyses of surface-tension measured values. The K_ia values determined from the transport experiments are very similar to those measured from surface-tension data, indicating that the transport experiments produced robust measurements of AWIA.

Assessment of Mobilization Potential of Per- and Polyfluoroalkyl Substances for Soil Remediation

This study investigated the mobilization of a wide range of per- and polyfluoroalkyl substances (PFASs) present in aqueous film-forming foams (AFFFs) in water-saturated soils through one-dimensional (1-D) column experiments with a view to assessing the feasibility of their remediation by soil desorption and washing. Results indicated that sorption/desorption of most of the shorter-carbon-chain PFASs (C ≤ 6) in soil reached greater than 99% rapidly─after approximately two pore volumes (PVs) and were well predicted by an equilibrium transport model, indicating that they will be readily removed by soil washing technologies. In contrast, the equilibrium model failed to predict the mobilization of longer-chain PFASs (C ≥ 7), indicating the presence of nonequilibrium sorption/desorption (confirmed by a flow interruption experiment). The actual time taken to attain 99% sorption/desorption was up to 5 times longer than predicted by the equilibrium model (e.g., ∼62 PVs versus ∼12 PVs predicted for perfluorooctane sulfonate (PFOS) in loamy sand). The increasing contribution of hydrophobic interactions over the electrostatic interactions is suggested as the main driving factor of the nonequilibrium processes. The inverse linear relationship (R² = 0.6, p < 0.0001) between the nonequilibrium mass transfer rate coefficient and the Freundlich sorption coefficient could potentially be a useful means for preliminary evaluation of potential nonequilibrium sorption/desorption of PFASs in soils.

A field study to assess the role of air-water interfacial sorption on PFAS leaching in an AFFF source area

Field-deployed lysimeters were used to measure the concentrations of poly- and perfluoroalkyl substances (PFASs) in soil porewater at a site historically impacted with aqueous film forming foam (AFFF). Samples collected over a 49-day period showed that perfluorooctane sulfonate (PFOS) and perfluorohexane sulfonate (PFHxS) were the PFASs with the highest concentrations in porewater, with concentrations of approximately 10,000 and 25,000 ng L^-1, respectively. The corresponding average mass flux to underlying groundwater observed for PFOS and PFHxS was 28,000 ± 11,000 and 92,000 ± 32,000 ng m^-2 d^-1, respectively. Employing the use of batch desorption isotherms (soil:water slurries) to determine desorption K_d values resulted in an overestimation of PFAS porewater concentrations by a factor for 1.4 to 4. However, using the desorption K_d values from the batch desorption isotherms in combination with a PFAS mass balance that incorporated PFAS sorption at the air-water interface resulted in improved predictions of the PFAS porewater concentrations. This improvement was most notable for PFOS, where inclusion of air-water interfacial sorption resulted in a 58% reduction in the predicted PFOS porewater concentration and predicted PFOS porewater concentrations that were identical (within the 95% confidence interval) to the lysimeter measured PFOS porewater concentration. Overall these results highlight the potentially important role of air-water interfacial sorption on PFAS migration in AFFF-impacted unsaturated soils in an in situ field setting.

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.

The impact of multiple-component PFAS solutions on fluid-fluid interfacial adsorption and transport of PFOS in unsaturated porous media

The objective of this research was to investigate the impact of multiple-component PFAS solutions on the retention of PFOS during transport in unsaturated porous media. Surface tensions were measured to characterize the impact of co-PFAS on the surface activity of PFOS. Miscible-displacement experiments were conducted to examine the air-water interfacial adsorption of PFOS during transport in single and multi-PFAS systems. Literature data for transport of PFOS in NAPL-water systems were also investigated for comparison. A mathematical model incorporating surfactant-induced flow, nonlinear rate-limited sorption, nonlinear rate-limited fluid-fluid interfacial adsorption, and competitive adsorption at the fluid-fluid interface was used to simulate the transport of PFOS. The results indicate that the presence of co-PFAS had no measurable impact on solid-phase sorption of PFOS during transport under the extant conditions of the experiments. Conversely, the air-water interfacial adsorption of PFOS was decreased by the presence of co-PFAS during transport under unsaturated-flow conditions for relatively high input concentrations. The multiple-component Langmuir model could not predict the competitive adsorption behavior observed during transport. Conversely, competitive interactions were not observed for transport with a lower input concentration. The results indicate that the retention and transport of individual PFAS in mixtures may in some cases be impacted by the presence of co-PFAS due to competitive fluid-fluid interfacial adsorption effects. Reduced retention due to competitive interfacial-adsorption interactions has the potential to decrease PFOS retardation during transport, thereby increasing migration rates in sources zones and enhancing groundwater-pollution risks. SYNOPSIS: The impact of PFAS mixtures on the retention and transport of PFOS in unsaturated porous media is examined with a series of experiments and mathematical modeling.

Field-Scale Demonstration of PFAS Leachability Following In Situ Soil Stabilization

A field-scale validation is summarized comparing the efficacy of commercially available stabilization amendments with the objective of mitigating per- and polyfluoroalkyl substance (PFAS) leaching from aqueous film-forming foam (AFFF)-impacted source zones. The scope of this work included bench-scale testing to evaluate multiple amendments and application concentrations to mitigate PFAS leachability and the execution of field-scale soil mixing in an AFFF-impacted fire-training area with nearly 2.5 years of post-soil mixing monitoring to validate reductions in PFAS leachability. At the bench scale, several amendments were evaluated and the selection of two amendments for field-scale evaluation was informed: FLUORO-SORB Adsorbent (FS) and RemBind (RB). Five ∼28 m³ test pits (approximately 3 m wide by 3 m long by 3 m deep) were mixed at a site using conventional construction equipment. One control test pit (Test Pit 1) included Portland cement (PC) only (5% dry weight basis). The other four test pits (Test Pits 2 through 5) compared 5 and 10% ratios (dry weight basis) of FS and RB (also with PC). Five separate monitoring events included two to three sample cores collected from each test pit for United States Environmental Protection Agency (USEPA) Method 1315 leaching assessment. After 1 year, a mass balance for each test pit was attempted comparing the total PFAS soil mass before, during, and after leach testing. Bench-scale and field-scale data were in good agreement and demonstrated >99% decrease in total PFAS leachability (mass basis; >98% mole basis) as confirmed by the total oxidizable precursor assay, strongly supporting the chemical stabilization of PFAS.

The Co-Transport of PFAS and Cr(VI) in porous media

PFAS and Cr are present at some sites as co-contaminants. The objective of this research was to investigate the co-transport behavior of per- and polyfluoroalkyl substances (PFAS) and hexavalent chromium (Cr(VI)) in porous media. Miscible-displacement experiments were conducted using two soils and an aquifer sediment with different geochemical properties. Perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) were employed as model PFAS. The retardation of PFOS was decreased in the presence of Cr(VI). Conversely, the transport and retardation of PFOA was not affected by the presence of Cr(VI). The reduction of PFOS retardation caused by Cr(VI) is likely due to sorption competition for both organic-carbon and inorganic (metal-oxides and clay minerals) domains. The relative contributions of the three soil constituents to PFOS sorption and the potential for competition between PFOS and Cr(VI) is a function of the geochemical composition of the porous media (i.e., organic carbon, metal-oxides and clay minerals). The PFAS had minimal impact on the retention and transport of Cr(VI). To our knowledge, the results presented herein represent the first reported data for PFOS and Cr(VI) co-transport in porous media. The results of this study indicate that the presence of Cr(VI) has the potential to increase the migration potential of PFOS in soil and groundwater, which should be considered when characterizing electroplating facilities, leather tanning facilities, and other co-contaminated sites.

Environmental Sources, Chemistry, Fate, and Transport of Per- and Polyfluoroalkyl Substances: State of the Science, Key Knowledge Gaps, and Recommendations Presented at the August 2019 SETAC Focus Topic Meeting

A Society of Environmental Toxicology and Chemistry (SETAC) Focused Topic Meeting (FTM) on the environmental management of per- and polyfluoroalkyl substances (PFAS) convened during August 2019 in Durham, North Carolina (USA). Experts from around the globe were brought together to critically evaluate new and emerging information on PFAS including chemistry, fate, transport, exposure, and toxicity. After plenary presentations, breakout groups were established and tasked to identify and adjudicate via panel discussions overarching conclusions and relevant data gaps. The present review is one in a series and summarizes outcomes of presentations and breakout discussions related to (1) primary sources and pathways in the environment, (2) sorption and transport in porous media, (3) precursor transformation, (4) practical approaches to the assessment of source zones, (5) standard and novel analytical methods with implications for environmental forensics and site management, and (6) classification and grouping from multiple perspectives. Outcomes illustrate that PFAS classification will continue to be a challenge, and additional pressing needs include increased availability of analytical standards and methods for assessment of PFAS and fate and transport, including precursor transformation. Although the state of the science is sufficient to support a degree of site-specific and flexible risk management, effective source prioritization tools, predictive fate and transport models, and improved and standardized analytical methods are needed to guide broader policies and best management practices. Environ Toxicol Chem 2021;40:3234-3260. © 2021 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Release of Per- and Polyfluoroalkyl Substances from Aqueous Film-Forming Foam Impacted Soils

Per- and polyfluoroalkyl substances (PFASs) are highly mobile in the saturated subsurface, yet aqueous film-forming foam (AFFF)-impacted source zones appear to be long lasting PFAS reservoirs. This study examined the release of over one hundred anionic and zwitterionic PFASs from two AFFF-impacted surface soils under saturated conditions with packed soil columns. Perfluoroalkyl acids (PFAAs) were released more rapidly than their polyfluorinated precursors, while anionic PFASs that were present in partially uncharged states were released more slowly than PFASs that were present entirely as anions, as were zwitterionic PFASs with terminal cationic functional groups when compared with analogous zwitterions with only anionic terminal groups. Nonideal transport was observed in both per- and polyfluorinated classes, as soil column effluent concentrations of slowly released PFASs increased by up to 107-fold with sustained artificial groundwater flow. A flow-interruption experiment suggested the influence of rate-limited desorption on diverse PFAS classes, including PFAAs with as few as four perfluorinated carbons. These results suggest that during infiltration the slow, rate-limited desorption of anionic and zwitterionic PFAA precursors may result in these compounds comprising an increasingly large fraction of the remaining PFASs in AFFF-impacted surface soils.

Incidence of Pfas in soil following long-term application of class B biosolids

This field study investigated the impact of long-term land application of biosolids on PFAS presence in soils that received annual repetitive land application of Class B biosolids from 1984 to 2019. Soil samples were collected from three depths of 30.5, 91 and 183 cm below land surface. Biosolid and groundwater samples used for irrigation were also collected. Concentrations measured for 18 PFAS compounds were evaluated to assess incidence rates and potential impact on groundwater. No PFAS analytes were detected at the three sampling depths for soil samples collected from undisturbed sites with no history of agriculture, irrigation, or biosolids application (background control sites). Relatively low mean concentrations of PFAS ranging from non-detect to 1.9 μg/kg were measured in soil samples collected from sites that were used for agriculture and that received irrigation with groundwater, but never received biosolids. PFAS concentrations in soils amended with biosolids were similarly low, ranging from non-detect to a mean concentration of 4.1 μg/kg. PFOS was observed at the highest concentrations, followed by PFOA for all locations. PFAS detected in the irrigation water were also present in the soil. These results indicate that biosolids and irrigation water are both important sources of PFAS present in the soils for all of the study sites. Not all PFAS detected in the biosolids were detected in the soil. Very long chain PFAS present in the biosolids were not detected or were detected at very low levels for soil, suggesting potential preferential retention within the biosolids. The precursor NMeFOSAA was present at the second highest concentrations in the biosolids but not detected in soil, indicating possible occurrence of transformation reactions. The total PFAS soil concentrations exhibited significant attenuation with depth, with a mean attenuation of 73% at the 183 cm depth. Monotonically decreasing concentrations with depth were observed for the longer-chain PFAS.

Perfluoroalkyl and polyfluoroalkyl substances (PFAS) in groundwater at a reclaimed water recharge facility

Reclaimed water is becoming an increasingly important source of water in arid regions worldwide. In the City of Tucson, Arizona, reclaimed water comprises approximately 10% of the annual water supply. It is used to recharge the local aquifer, create surface flow in the Santa Cruz River, and irrigate parks, golf courses, and recreational fields. In December 2018, concentrations of perfluoroalkyl and polyfluoroalkyl substances (PFAS) an order of magnitude higher than the EPA lifetime health advisory of 70 ppt were discovered in the city's reclaimed water system. The PFAS were also detected in the Sweetwater Recharge Facility (SRF), adjacent to the Santa Cruz River, where reclaimed water is stored in the alluvial aquifer. PFAS have gained national attention as contaminants of emerging concern because of their widespread occurrence, toxicological impact to humans, and persistence in the environment. However, relatively little is known about their fate and transport in managed aquifer recharge systems. Results from this study show that PFAS in the SRF likely originated from the city's retired wastewater treatment facility, while lower PFAS concentrations are observed in the treated wastewater provided by the city's new treatment facility. Moreover, the combined PFOS and PFOA concentrations appear to be correlated to rising and falling groundwater levels, indicating that PFAS are likely trapped in the vadose zone and transported to the alluvial aquifer during managed aquifer recharge events.

Release of soil colloids during flow interruption increases the pore-water PFAS concentration in saturated soil

Groundwater flow through aquifer soils or packed bed systems can fluctuate for various reasons, which could affect the concentration of natural colloids and per- and polyfluoroalkyl substances (PFAS) in the pore water. In such cases, PFAS concentration could either decrease due to matrix diffusion of PFAS or increase by the detachment of colloids carrying PFAS. Yet, the effect of flow fluctuation on PFAS transport or release in porous media has not been examined. To examine the relative importance of either process, we interrupted the flow during an injection of groundwater spiked with perfluorobutanoic acid (PFBA), perfluorooctanoic acid (PFOA), and bromide as conservative tracer through clay-rich soil, so that diffusive transport would be prominent during flow interruption. After flow interruption, the PFAS concentration did not decrease indicating an insignificant contribution of matrix diffusion. The concentration increased, potentially due to enhanced release of colloid-associated PFAS. Analysis of samples before and after flow interruption by particle size analysis and SEM confirmed an increase in soil colloid concentration after the flow interruption. XRD analysis of soil and the colloids proved that PFAS were associated with specific sites of the colloids. Due to a higher affinity of PFOA to soil colloids, the total PFOA concentration in the effluent samples increased more than PFBA after the flow interruption process. The results indicate that colloids may have a disproportionally higher role in the transport of PFAS in conditions that release colloids from porous media. Thus, fluctuations in groundwater flow can increase this colloid facilitated mobility of PFAS.

The effectiveness of PFAS management options on groundwater quality in contaminated land using numerical modelling

The concentrations of PFOS and PFHxS in a contaminated site in regional New South Wales, Australia, were above the human health screening values for industrial land use. In this study, the effects of different management options on the quality of groundwater were investigated through numerical modelling. At first, a complete transfer model including the main features of advection, dispersion, adsorption and decay, was developed to simulate the long-term migration of PFOS from topsoil subjected to full climate interaction for 100 years. The sorption isotherm characteristics of the contaminated soil were determined from chemical analysis using LC/MS equipment. The model results were confirmed by PFOS values measured from a monitoring well in the proximity of the site. The model showed that PFOS values in groundwater increased gradually and exceeded the guideline values for drinking water. Three management options were suggested: a do-nothing approach, cut and replacement, and immobilisation of the topsoil up to 2 m depth. The numerical models showed that although all these strategies reduced the PFOS level in the groundwater significantly, the values were still higher than the guideline values for drinking water. This was because PFOS migrated in the ground beyond the site location. The cut and replacement and immobilisation strategies ensured that the PFOS values were lower than the guideline values for soil screening, but PFOS levels in the groundwater were not necessarily lower than the guideline values for drinking water after a long time.

Ideal versus Nonideal Transport of PFAS in Unsaturated Porous Media

Per- and poly-fluoroalkyl substances (PFAS) adsorb at air-water interfaces during transport in unsaturated porous media. This can cause surfactant-induced flow and enhanced retention that is a function of concentration, which complicates characterization and modeling of PFAS transport under unsaturated conditions. The influence of surfactant-induced flow and nonlinear air-water interfacial adsorption (AWIA) on PFAS transport was investigated with a series of miscible-displacement transport experiments conducted with a several-log range in input concentrations. Perfluorooctane sulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and ammonium perfluoro 2-methyl-3-oxahexanoate (GenX) were used as model PFAS. The results were interpreted in terms of critical reference concentrations associated with PFAS surface activities and their relationship to the relevancy of transport processes such as surfactant-induced flow and nonlinear AWIA for concentration ranges of interest. Analysis of the measured transport behavior of PFAS under unsaturated-flow conditions demonstrated that AWIA was linear when the input concentration was sufficiently below the critical reference concentration. This includes the absence of significant arrival-front self-sharpening and extended elution tailing of the breakthrough curves, as well as the similarity of retardation factors measured for a wide range of input concentrations. Independently-predicted simulations produced with a comprehensive flow and transport model that accounts for transient variably-saturated flow, surfactant-induced flow, nonlinear rate-limited solid-phase sorption, and nonlinear rate-limited AWIA provided excellent predictions of the measured transport. A series of simulations was conducted with the model to test the specific impact of various processes potentially influencing PFOS transport. The simulation results showed that surfactant-induced flow was negligible and that AWIA was effectively linear when the input concentration was sufficiently below the critical reference concentration. PFAS retention associated with AWIA can be considered to be ideal in such cases, thereby supporting the use of simplified mathematical models. Conversely, apparent nonideal transport behavior was observed for experiments conducted with input concentrations similar to or greater than the critical reference concentration.

Impact of a Hydrocarbon Surfactant on the Retention and Transport of Perfluorooctanoic Acid in Saturated and Unsaturated Porous Media

The transport and retention behavior of perfluorooctanoic acid (PFOA) in the presence of a hydrocarbon surfactant under saturated and unsaturated conditions was investigated. Miscible-displacement transport experiments were conducted at different PFOA and sodium dodecyl sulfate (SDS) input ratios to determine the impact of SDS on PFOA adsorption at solid-water and air-water interfaces. A numerical flow and transport model was employed to simulate the experiments. The PFOA breakthrough curves for unsaturated conditions exhibited greater retardation compared to those for saturated conditions in all cases, owing to air-water interfacial adsorption. The retardation factor for PFOA with a low concentration of SDS (PFOA-SDS ratio of 10:1) was similar to that for PFOA without SDS under unsaturated conditions. Conversely, retardation was greater in the presence of higher levels of SDS (1:1 and 1:10) with retardation factors increasing from 2.4 to 2.9 and 3.6 under unsaturated conditions due to enhanced adsorption at the solid-water and air-water interfaces. The low concentration of SDS had no measurable impact on PFOA air-water interfacial adsorption coefficients (K_ia) determined from the transport experiments. The presence of SDS at the higher PFOA-SDS concentration ratios increased the surface activity of PFOA, with transport-determined K_ia values increased by 27 and 139%, respectively. The model provided very good independently predicted simulations of the measured breakthrough curves and showed that PFOA and SDS experienced various degrees of differential transport during the experiments. These results have implications for the characterization and modeling of poly-fluoroalkyl substances (PFAS) migration potential at sites wherein PFAS and hydrocarbon surfactants co-occur.

Transport of PFOS in aquifer sediment: Transport behavior and a distributed-sorption model

The objectives of this research were to examine the transport of perfluorooctane sulfonic acid (PFOS) in aquifer sediment comprising different geochemical properties, and to compare the behavior to that observed for PFOS transport in soil and sand. PFOS retardation was relatively low for transport in all aquifer media. The PFOS breakthrough curves were asymmetrical and exhibited extensive concentration tailing, indicating that sorption/desorption was significantly nonideal. The results of model simulations indicated that rate-limited sorption/desorption was the primary cause of the nonideal PFOS transport. Comparison of PFOS transport in aquifer media to data reported for PFOS transport in two soils and a quartz sand showed that PFOS exhibited more extensive elution tailing for the soils, likely reflecting differences in the relative contributions of various media constituents to sorption. A three-component distributed-sorption model was developed that accounted for contributions from soil organic carbon, metal oxides, and silt + clay fraction. The model produced very good predictions of K_d for the five media with lower soil organic‑carbon contents (≤0.1%). Soil organic carbon was estimated to contribute 19-42% of the total sorption for all media except the sand, to which it contributed ~100%. The contribution of silt + clay ranged from 51 to 80% for all media except the sand. The only medium for which the contribution of metal-oxides was significant is Hanford, with an estimated contribution of 15%. Overall, the results of the study indicate that sorption of PFOS by these aquifer media comprised contributions from multiple soil constituents.

Simulating PFAS adsorption kinetics, adsorption isotherms, and nonideal transport in saturated soil with tempered one-sided stable density (TOSD) based models

Reliable quantification of per- and polyfluoroalkyl substances (PFAS) adsorption and mobility in geomedia provides critical information (i.e., evaluation and prediction) for risk characterization and mitigation strategy development. Given the limited PFAS data available and various competing theories for modeling pollutant kinetics, it is indispensable to better understand and quantify the adsorption and transport of PFAS in geomedia using generalized models built upon a consistent physical theory. This study proposed a universal physical law (called the tempered stable law) in PFAS adsorption/transport by interpreting PFAS adsorption kinetics and nonideal transport as a nonequilibrium process dominated by adsorption/desorption with multiple rates following the tempered one-sided stable density (TOSD) distribution. This universal TOSD function led to novel TOSD-based models which were then tested by successfully simulating PFAS adsorption kinetics, adsorption isotherms, and nonideal transport data reported in the literature. Model comparisons and extensions were also discussed to further check the feasibility of the TOSD models and their adaptability to capture PFAS transport in more complex geomedia at all scales.

The Case for Direct Measures of Soil-to-Groundwater Contaminant Mass Discharge at AFFF-Impacted Sites

Many entities around the world are initiating massive field campaigns to characterize the environmental distribution of per- and polyfluoroalkyl substances (PFAS), particularly at aqueous film-forming foam (AFFF) impacted sites where historic point-source discharges occurred at the ground surface. Concurrently, many regulatory agencies are publishing criteria used in practice to define the "nature and extent" of PFAS-impacted environmental media. Specific to the soil-to-groundwater transport pathway protective of the groundwater ingestion end point, these soil criteria (or screening values) are to date exclusively based on the traditional approach used for hydrophobic organics with a number of simplifying assumptions. Research has clearly contradicted these assumptions, yet alternative methodologies have yet to emerge from the literature. This Perspective provides a rationale for why alternative approaches are critically necessary and proposes the application of lysimetry as a practical solution to accurately assess PFAS transport within unsaturated vadose zone soils. Ultimately, there is an urgent need to justify soil remediation on the basis of groundwater protection and to prioritize remedial efforts in order to optimize limited fiscal resources.

Contribution of Nonaqueous-Phase Liquids to the Retention and Transport of Per and Polyfluoroalkyl Substances (PFAS) in Porous Media

Per and polyfluoroalkyl substances (PFAS) cocontamination with nonaqueous-phase organic liquids (NAPLs) has been observed or suspected at various sites, particularly at fire-training areas at which aqueous film-forming foams (AFFFs) were applied. The objectives of this study are to (1) delineate the relative significance of specific PFAS-NAPL processes on PFAS retention, including partitioning into the bulk NAPL phase and adsorption to the NAPL-water interface; (2) investigate the influence of NAPL properties, saturation, and mass-transfer constraints on PFAS retention; and (3) determine whether PFAS may impact NAPL distribution through mobilization or dissolution. Perfluorooctanesulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) are used as representative PFAS, and trichloroethene (TCE) and decane are used as representative NAPLs. NAPL-water interfacial adsorption was quantified with NAPL-water interfacial-tension measurements; partitioning into NAPL was quantified with batch experiments, and retardation factors (R) in the absence and presence of residual NAPL were determined with miscible-displacement transport experiments. R values increased in the presence of residual NAPL, with adsorption to the NAPL-water interface accounting for as much as ∼77% of retention and solid-phase adsorption also significantly contributing to retention. Additionally, this study provides the first QSPR analysis focused on NAPL-water interfacial adsorption coefficients, with results consistent with those from previous air-water studies. Lastly, this initial investigation into PFAS impacts on NAPL behavior determined that PFOS/PFOA are unlikely to enhance solubilization or mobilization of NAPL under the conditions present at many AFFF legacy sites.

A critical review of modeling Poly- and Perfluoroalkyl Substances (PFAS) in the soil-water environment

Due to their health effects and the recalcitrant nature of their CF bonds, Poly- and Perfluoroalkyl Substances (PFAS) are widely investigated for their distribution, remediation, and toxicology in ecosystems. However, very few studies have focused on modeling PFAS in the soil-water environment. In this review, we summarized the recent development in PFAS modeling for various chemical, physical, and biological processes, including sorption, volatilization, degradation, bioaccumulation, and transport. PFAS sorption is kinetic in nature with sorption equilibrium commonly quantified by either a linear, the Freundlich, or the Langmuir isotherms. Volatilization of PFAS depends on carbon chain length and ionization status and has been simulated by a two-layer diffusion process across the air water interface. First-order kinetics is commonly used for physical, chemical, and biological degradation processes. Uptake by plants and other biota can be passive and/or active. As surfactants, PFAS have a tendency to be sorbed or concentrated on air-water or non-aqueous phase liquid (NAPL)-water interfaces, where the same three isotherms for soil sorption are adopted. PFAS transport in the soil-water environment is simulated by solving the convection-dispersion equation (CDE) that is coupled to PFAS sorption, phase transfer, as well as physical, chemical, and biological transformations. As the physicochemical properties and concentration vary greatly among the potentially thousands of PFAS species in the environment, systematic efforts are needed to identify models and model parameters to simulate their fate, transport, and response to remediation techniques. Since many process formulations are empirical in nature, mechanistic approaches are needed to further the understanding of PFAS-soil-water-plant interactions so that the model parameters are less site dependent and more predictive in simulating PFAS remediation efficiency.

Examining the robustness and concentration dependency of PFAS air-water and NAPL-water interfacial adsorption coefficients

Determining robust values for the air-water or NAPL-water interfacial adsorption coefficient, KIA, is key to characterizing and modeling PFAS transport and fate in several environmental systems. Direct, high-resolution measurements of surfactant adsorption at the fluid-fluid interface were aggregated from the literature. This data set was used to examine the accuracy and applicability of Γ and KIA measurements determined for three PFAS from transport experiments and surface-tension data. The transport-measured Γ and KIA data were observed to be fully consistent with the directly-measured data. Specifically, Γ values for the two methods were entirely coincident in the region of overlapping concentrations, which spanned ~4 orders-of-magnitude. Furthermore, the two data sets adhered to an identical Γ-C profile. These results conclusively demonstrate the accuracy of the transport-measured values. Γ and KIA values determined from the application of the Gibbs adsorption equation to measured surface-tension data were fully consistent with the directly-measured and transport-measured data sets, demonstrating their applicability for representing PFAS transport in environmental systems. The directly-measured data were used to examine the concentration dependency of KIA values, absent the potential confounding effects associated with the use of surface-tension or transport-measured data. The directly-measured data clearly demonstrate that KIA attains a constant, maximum limit at lower concentrations. Two separate analyses of the transport-measured data both produced observations of constant KIA values at lower concentrations, consistent with the directly-measured data. These outcomes are discussed in terms of surface activities, relative surface coverages, and critical concentrations.

Effects of ionic strength and cation type on the transport of perfluorooctanoic acid (PFOA) in unsaturated sand porous media

Current understanding of perfluorooctanoic acid (PFOA) transport in unsaturated porous media is still limited with significant variability in solution chemistry. Column experiments were conducted to systematically evaluate the impacts of ionic strength (1.5-30 mM) and cation type (Na⁺ and Ca²⁺) on PFOA transport in unsaturated quartz sand. The results showed that an increase in ionic strength (1.5-30 mM) led to greater PFOA retardation in unsaturated columns. Meanwhile, Ca²⁺ caused more PFOA retardation than Na⁺ at the same unsaturated conditions. These findings were supported by bubble column experiments, which indicated greater PFOA adsorption at the air-water interface with increasing ionic strength or in the presence of Ca²⁺ in comparison to Na⁺. Furthermore, the air-water interfacial (AWI) adsorption coefficients calculated from surface tension isotherms also increased with increasing ionic strength or in the presence of Ca²⁺ in comparison to Na⁺. These results clearly confirm that higher ionic strength or cation valence significantly promoted PFOA adsorption at the air-water interface, and thus caused greater PFOA retardation during transport in unsaturated porous media. This work points out the importance of considering solution ionic strength and cation type in assessing the transport behavior of PFOA in unsaturated porous media.

Column versus batch methods for measuring PFOS and PFOA sorption to geomedia

The objective of this study is to compare the consistency between column and batch experiment methods for measuring solid-phase sorption coefficients and isotherms for per and polyfluoroalkyl substances (PFAS). Perfluorooctane sulfonic acid (PFOS) and perfluorooctanoic acid (PFOA) are used as representative PFAS, and experiments are conducted with three natural porous media with differing geochemical properties. Column-derived sorption isotherms are generated by conducting multiple experiments with different input concentrations (multi-C₀ method) or employing elution-front integration wherein the entire isotherm is determined from a single breakthrough curve (BTC) elution front. The isotherms generated with the multi-C₀ column method compared remarkably well to the batch isotherms over an aqueous concentration range of 3-4 orders of magnitude. Specifically, the 95% confidence intervals for the individual isotherm variables overlapped, producing statistically identical regressions. The elution-front integration isotherms generally agreed with the batch isotherms, but exhibited noise and systematic deviation at lower concentrations in some cases. All data sets were well described by the Freundlich isotherm model. Freundlich N values ranged from 0.75 to 0.81 for PFOS and was 0.87 for PFOA and are consistent with values reported in the literature for different geomedia. The results of this study indicate that column and batch experiments can measure consistent sorption isotherms and sorption coefficients for PFOS and PFOA when robust experimental setup and data analysis are implemented.

Influences of Chemical Properties, Soil Properties, and Solution pH on Soil-Water Partitioning Coefficients of Per- and Polyfluoroalkyl Substances (PFASs)

The aim of this study was to assess the soil-water partitioning behavior of a wider range of per- and polyfluoroalkyl substances (PFASs) onto soils covering diverse soil properties. The PFASs studied include perfluoroalkyl carboxylates (PFCAs), perfluoroalkane sulfonates (PFSAs), fluorotelomer sulfonates (FTSs), nonionic perfluoroalkane sulfonamides (FASAs), cyclic PFAS (PFEtCHxS), per- and polyfluoroalkyl ether acids (GenX, ADONA, 9Cl-PF3ONS), and three aqueous film-forming foam (AFFF)-related zwitterionic PFASs (AmPr-FHxSA, TAmPr-FHxSA, 6:2 FTSA-PrB). Soil-water partitioning coefficients (log K_d values) of the PFASs ranged from less than zero to approximately three, were chain-length-dependent, and were significantly linearly related to molecular weight (MW) for PFASs with MW > 350 g/mol (R² = 0.94, p < 0.0001). Across all soils, the K_d values of all short-chain PFASs (≤5 -CF2- moieties) were similar and varied less (<0.5 log units) compared to long-chain PFASs (>0.5 to 1.5  log units) and zwitterions AmPr- and TAmPr-FHxSA (∼1.5 to 2 log units). Multiple soil properties described sorption of PFASs better than any single property. The effects of soil properties on sorption were different for anionic, nonionic, and zwitterionic PFASs. Solution pH could change both PFAS speciation and soil chemistry affecting surface complexation and electrostatic processes. The K_d values of all PFASs increased when solution pH decreased from approximately eight to three. Short-chain PFASs were less sensitive to solution pH than long-chain PFASs. The results indicate the complex interactions of PFASs with soil surfaces and the need to consider both PFAS type and soil properties to describe mobility in the environment.

PFAS concentrations in soils: Background levels versus contaminated sites

Per- and polyfluoroalkyl substances (PFAS) are contaminants of critical concern due to their persistence, widespread distribution in the environment, and potential human-health impacts. In this work, published studies of PFAS concentrations in soils were compiled from the literature. These data were combined with results obtained from a large curated database of PFAS soil concentrations for contaminated sites. In aggregate, the compiled data set comprises >30,000 samples collected from >2500 sites distributed throughout the world. Data were collected for three types of sites- background sites, primary-source sites (fire-training areas, manufacturing plants), and secondary-source sites (biosolids application, irrigation water use). The aggregated soil-survey reports comprise samples collected from all continents, and from a large variety of locations in both urban and rural regions. PFAS were present in soil at almost every site tested. Low but measurable concentrations were observed even in remote regions far from potential PFOS sources. Concentrations reported for PFAS-contaminated sites were generally orders-of-magnitude greater than background levels, particularly for PFOS. Maximum reported PFOS concentrations ranged upwards of several hundred mg/kg. Analysis of depth profiles indicates significant retention of PFAS in the vadose zone over decadal timeframes and the occurrence of leaching to groundwater. It is noteworthy that soil concentrations reported for PFAS at contaminated sites are often orders-of-magnitude higher than typical groundwater concentrations. The results of this study demonstrate that PFAS are present in soils across the globe, and indicate that soil is a significant reservoir for PFAS. A critical question of concern is the long-term migration potential to surface water, groundwater, and the atmosphere. This warrants increased focus on the transport and fate behavior of PFAS in soil and the vadose zone, in regards to both research and site investigations.

Transport of GenX in Saturated and Unsaturated Porous Media

The objective of this research was to investigate the retention and transport behavior of GenX in five natural porous media with similar median grain diameters but different geochemical properties. Surface tensions were measured to characterize surface activity. Miscible-displacement experiments were conducted under saturated conditions to characterize the magnitude of solid-phase adsorption, while unsaturated-flow experiments were conducted to examine the impact of air-water interfacial adsorption on retention and transport. The results from surface-tension measurements showed that the impact of solution composition is greater for the ammonium form of GenX than for the acid form, due to the presence of the NH₄ counterion. The breakthrough curves for the experiments conducted under saturated conditions were asymmetrical, and a solute-transport model employing a two-domain representation of nonlinear, rate-limited sorption provided reasonable simulations of the measured data. The magnitudes of solid-phase adsorption were relatively small, with the highest adsorption associated with the medium containing the greatest amount of metal oxides. The breakthrough curves for the experiments conducted under unsaturated conditions exhibited greater retardation due to the impact of adsorption at the air-water interface. The contributions of air-water interfacial adsorption to GenX retention ranged from ∼24% to ∼100%. The overall magnitudes of retardation were relatively low, with retardation factors < ∼3, indicating that GenX has significant migration potential in soil and the vadose zone. To our knowledge, the results presented herein represent the first reported data for solid-water and air-water interfacial adsorption of GenX by soil. These data should prove useful for assessing the transport and fate behavior of GenX in soil and groundwater.

Fluorinated Surfactant Adsorption on Mineral Surfaces: Implications for PFAS Fate and Transport in the Environment

Fluorinated surfactants, which fall under the class of per- and polyfluoroalkyl substances (PFAS), are amphiphilic molecules that comprise hydrophobic fluorocarbon chains and hydrophilic head-groups. Fluorinated surfactants have been utilized in many applications, e.g., fire-fighting foams, paints, household/kitchenware items, product packaging, and fabrics. These compounds then made their way into the environment, and have been detected in soil, fresh water, and seawater. From there, they can enter human bodies. Fluorinated surfactants are persistent in water and soil environments, and their adsorption onto mineral surfaces contributes to this persistence. This review examines how fluorinated surfactants adsorb onto mineral surfaces, by analyzing the thermodynamics and kinetics of adsorption, and the underlying mechanisms. Adsorption of fluorinated surfactants onto mineral surfaces can be explained by electrostatic interactions, hydrophobic interactions, hydrogen bonding, and ligand and ion exchange. The aqueous pH, varying salt or humic acid concentrations, and the surfactant chemistry can influence the adsorption of fluorinated surfactants onto mineral surfaces. Further research is needed on fluorinated surfactant adsorbent materials to treat drinking water, and on strategies that can modulate the fate of these compounds in specific environmental locations.

Importance of Al/Fe oxyhydroxide coating and ionic strength in perfluorooctanoic acid (PFOA) transport in saturated porous media

Understanding subsurface transport of per- and polyfluoroalkyl substances (PFASs) is of critical importance for the benign use and risk management of PFASs. As one of the most commonly found PFASs, perfluorooctanoic acid (PFOA) is used as a representative PFAS and water-saturated column experiments were conducted to investigate the effect of Al/Fe oxyhydroxide coating and ionic strength on its transport at an environmentally relevant PFOA concentration (6.8 μg L^-1). Our results showed a clear increase in PFOA retardation in Al/Fe oxyhydroxide coated sand (retardation factor: Al: 1.87-5.58, Fe: 1.28-4.05) than those in uncoated sand (1.00-1.05), due to the stronger electrostatic attraction between anionic PFOA and Al/Fe oxyhydroxide coated sand surface. Notably, Al oxyhydroxide have a more profound effect on PFOA retention and retardation than Fe oxyhydroxide. Besides, higher ionic strength significantly inhabited PFOA retention and retardation in positively charged sand, and the considerable retention of PFOA (∼90%) in deionized water than those in 1.5 mM and 30.0 mM NaCl (<10%) clearly proves the role of competitive adsorption of Cl^- on PFOA transport in positively charged sand. In contrast, higher ionic strength (0 mM-30 mM NaCl) slightly increased PFOA retardation in negatively charged sand, illustrating the dominance of electrostatic interaction. Our findings advance current knowledge to understand PFOA transport in natural media with different surface charge property under environmental PFOA concentrations.

A Mathematical Model for the Release, Transport, and Retention of Per- and Polyfluoroalkyl Substances (PFAS) in the Vadose Zone

Per- and Polyfluoroalkyl Substances (PFAS) are emerging contaminants of critical concern. As surfactants, PFAS tend to accumulate at air-water interfaces and may stay in the vadose zone for long times before contaminating groundwater. Yet not well understood, the extent of retention in the vadose zone has critical implications for risk management and remediation strategies. We present the first mathematical model that accounts for surfactant-induced flow and solid-phase and air-water interfacial adsorption. We apply the model to simulate PFOS (a PFAS compound of primary concern) transport in the vadose zone at a model fire-training area site impacted by Aqueous Film-Forming Foam (AFFF). Air-water interfacial adsorption is shown to have a significant impact-amplified by the low water content due to gravity drainage-total retardation factors range from 233 to 1355 for the sand and 146 to 792 for the soil used in the study. The simulations illustrate it can take several decades or longer for PFOS to reach groundwater. Counterintuitively, the lower water content in the sand-due to stronger drainage and weaker capillary retention-leads to retardation factors greater than for the soil. Also, most PFOS is adsorbed at air-water interfaces with only 1-2% in the aqueous phase. The implications include 1) fine-texture materials could have lower retardation factors than sand due to higher retained water content, 2) soil PFAS concentrations are likely to be orders of magnitude higher than those in groundwater at source zones. Both implications are consistent with recent field observations at hundreds of AFFF-impacted sites.

Simulating PFAS transport influenced by rate-limited multi-process retention

The transport of per- and poly-fluoroalkyl substances (PFAS) in the vadose zone is complicated by the fact that multiple mass-transfer processes can contribute to their retention and retardation. In addition, PFAS transport at some sites can be further complicated by the presence of organic immiscible liquids (OIL). Mass-transfer processes are inherently rate limited and, therefore, have the potential to cause nonideal transport of PFAS. The objectives of this research were to: (1) develop a solute-transport model that explicitly accounts for multiple retention processes, including adsorption at air-water and OIL-water interfaces, adsorption by the solid phase, and diffusive mass-transfer between advective and nonadvective domains, and (2) apply the model to measured transport data to delineate which processes are rate limited and contribute to observed nonideal transport. Breakthrough curves for transport of two PFAS and one hydrocarbon surfactant in sand obtained from prior miscible-displacement experiments exhibited nonideal transport. The multiprocess model effectively simulated the measured transport data. The results of the analyses indicate that adsorption at the air-water and OIL-water interface can generally be treated as effectively instantaneous for transport in porous media. The rate limitations associated with solid-phase adsorption and diffusive mass transfer between advective and nonadvective domains were of greater significance.