Evaluating air-water and NAPL-water interfacial adsorption and retention of Perfluorocarboxylic acids within the Vadose zone

Structure, orientation, and dynamics of per- and polyfluoroalkyl substance (PFAS) surfactants at the air-water interface: Molecular-level insights

HYPOTHESIS

Understanding the intricate molecular-level details of toxic per- and polyfluoroalkyl substances (PFAS) partitioning to the air-water interface holds paramount importance in evaluating their fate and transport, as well as for finding safer alternatives for various applications, including aqueous film forming foams. The behavior of these substances at interfaces strongly depends on molecular architecture, chemistry, and concentration, which define molecular packing, self-assembly, interfacial diffusion, and the surface tension.

SIMULATIONS

Modeling of three PFAS surfactants, namely, longer-tail (perfluorooctanoate (PFOA)) and shorter-tail (perfluorobutanoate (PFBA) and 2,3,3,3-tetrafluoro-2-(heptafluoropropoxy) propanoate (GenX)) has been conducted using atomistic molecular dynamics simulations. A systematic comparison between these representative PFAS of different sizes and structure reveals factors influencing their association behavior, mechanism of surface tension reduction, and interfacial mobility as a function of surface coverage.

FINDINGS

Shorter-chain PFAS surfactants (GenX or PFBA) require lower surface coverage compared to longer chain (PFOA) PFAS to achieve the same decrease in surface tension. However, a higher concentration of GenX and PFBA is necessary in the bulk aqueous solution to achieve the same surface coverage as PFOA, due to their higher solubility in water. The PFAS molecular orientation and mobility at the interface are found to be vastly influenced by the length and architecture of the hydrophobic fluorocarbon tail. A significant ordering of the water dipole moment near the anionic headgroup is apparent at high surface concentration. A direct correlation is established between the PFAS interfacial properties and PFAS-PFAS, PFAS-counterion, and PFAS-water interactions.

PFAS fate using lysimeters during degraded soil reclamation using biosolids

Carbon- and nutrient-rich biosolids are used in agriculture and land reclamation. However, per- and polyfluoroalkyl substances (PFAS) typically present in biosolids raise concerns of PFAS leaching to groundwater and plant uptake. Here, we investigated PFAS persistence and leaching from biosolids applied to a site constructed artificially to mimic degraded soils. Treatments included biosolids and biosolids blended with mulch applied at different rates to attain either one and five times the agronomic N rate for vegetable crops and a control treatment with synthetic urea and triple superphosphate fertilizer. Leachates were collected for a 2-year period from 15-cm depth zero-tension drainage lysimeters. Soils were analyzed post biosolids application. PFAS were quantified using isotope-dilution, solid-phase extraction and liquid chromatography tandem mass spectrometry. Leachate profiles exemplified an initial high total PFAS concentration, followed by a sharp decline and subsequent small fluctuations attributed to pre-existing soil conditions and rainfall patterns. Quantifiable PFAS in leachate were proportional to biosolids application rates. Short-chain perfluoroalkyl acids (CF2 < 6) were dominant in leachate, while the percentage of longer chains homologues was higher in soils. A 43% biosolids blend with mulch resulted in 21% lower PFAS leachate concentrations even with the blend application rate being 1.5 times higher than biosolids due to the blend's lower N-content. The blending effect was more pronounced for long-chain perfluoroalkyl sulfonic acids that have a greater retention by soils and the air-water interface. Biosolids blending as a pragmatic strategy for reducing PFAS leachate concentrations may aid in the sustainable beneficial reuse of biosolids.

Influence of kinetic air-water interfacial partitioning on unsaturated transport of PFAS in sandy soils

This study investigates the impact of kinetic air-water partitioning on the transport of perfluoroalkyl substances (PFAS) within homogeneous and heterogeneous sandy vadose zones under transient unsaturated flow conditions. These experimental conditions are realistic for field behavior, where transient flow foments the continual growth and collapse of air-water interfaces (AWIs), and where layered heterogenous conditions enhance the perturbations of AWIs. Short-chain PFAS behave like conservative tracers with negligible air-water interface partitioning, whereas longer-chain PFAS demonstrate non-equilibrium retention behavior, especially in heterogeneous media. AWI partitioning kinetics were found to be important in controlling PFAS transport and mass flux, particularly during PFAS sorption to the air-water interface, which results because of the different nature and more rapid changes in AWI during drainage, wherein PFAS are moving toward the interface to achieve equilibrium, than during imbibition, where PFAS are leaving the interface to achieve equilibrium. Neglecting these kinetic AWI sorption processes can result in an underestimate of the PFAS transport velocities and mass flux reaching the water table. The presence of trapped air may also inhibit PFAS partitioning in a similar manner by causing longer diffusion paths from bulk water to a portion of the AWIs. The modified HYDRUS effectively captured the transport processes and provided an excellent match to the measured breakthrough curves. To assess relevance using realistic transient infiltration rates, simulations were conducted using precipitation data from an actual site. The results showed that accounting for kinetic AWI partitioning increases the cumulative PFOS mass flux to groundwater by a factor of 2.3 compared to equilibrium conditions, significantly impacting PFAS porewater concentrations. This difference was threefold under experimental conditions, suggesting that the importance of kinetic effects may vary significantly over the long term and under different climatic conditions or soil types, due to their strong dependence on water flux.

Environmental fate and transport of PFAS in wastewater treatment plant effluent discharged to rapid infiltration basins

Fate and transport of per- and polyfluoroalkyl substances (PFAS) in wastewater treatment plant (WWTP) effluent discharged to rapid infiltration basins (RIBs) is investigated using data from 26 WWTPs in Michigan, USA. PFAS were found to accumulate in groundwater downgradient from RIBs with median groundwater-effluent enrichment factors for ten commonly detected, terminal-form perfluoroalkyl acids (PFAAs) ranging from 1.3 to 5.2. Maximum contaminant levels for drinking water were exceeded in groundwater at all WWTPs with available PFAS data. Numerical models of unsaturated fluid flow and PFAS transport honoring RIB site properties, such as median vertical separation distance to the water table and a realistic range of area-normalized effluent fluxes, show long-chain PFAS undergo significant delays from air-water interface (AWI) adsorption, requiring up to 15 times longer to reach maximum mass flux to the saturated zone under low-flux conditions, where AWI area is 2.5 times greater. Short-chain PFAS commonly detected in effluent are only minimally affected by AWI adsorption and show little to no attenuation under high-flux conditions. The nonlinear inverse relationship between water content and AWI area highlights the important role of AWI adsorption in modulating unsaturated transport of long-chain PFAS to underlying groundwater due to the broad range of flux rates applied to RIB systems.

Transport and competitive interfacial adsorption of PFOA and PFOS in unsaturated porous media: Experiments and modeling

Among emerging contaminants, per- and polyfluoroalkyl substances (PFAS) have captured public attention based upon their environmental ubiquity and potential risks to human health. Due to their typical surface release conditions and amphiphilic properties, PFAS tend to sorb to soil and accumulate at the air-water interface within the vadose zone. These processes can result in substantial plume attenuation. Although there is a growing body of literature on vadose zone transport, few studies have explored PFAS mixture transport, particularly under conditions where nonlinear sorption processes are important. The present study aims to advance our understanding of PFAS transport in variably saturated porous media through integration of experiments and mathematical modeling. Experiments include batch studies to quantify sorption to the solid phase, interfacial tension (IFT) measurements to estimate adsorption at the air-water interface (AWI), and column studies with F-70 Ottawa sand at 100 % and ca. 50 % water saturation to explore transport mechanisms. Employed PFAS solutions encompass individual solutes and binary mixtures of perfluorooctanoic acid (PFOA) and perfluorooctane sulfonate (PFOS) at concentration levels spanning four orders of magnitude to assess competitive and nonlinear sorption at the AWI. Observations demonstrate that concentration levels and competitive effects substantially influence PFAS transport in unsaturated systems. In the presence of PFOS, PFOA experienced less retention than would be anticipated based on single-solute behavior, and effluent breakthrough curves exhibited chromatographic peaking. The presented mathematical model for simultaneous flow and transport of PFAS was able to capture experimental observations with a consistent set of parameters and minimal curve fitting. These results demonstrate the robustness of the model formulation that included rate-limited interfacial mass transfer, an extended Langmuir-Szyszkowski model for adsorption at the AWI, and a scaled Leverett thermodynamic model to predict the AWI specific area. Overall, the results of this work underscore the importance of the AWI in PFAS transport and highlight the relevance of competition effects in adsorption formulations.

Occurrence, Concentration, and Distribution of 35 PFASs and Their Precursors Retained in 20 Stormwater Biofilters

Current knowledge about the fate and transport behaviors of per- and polyfluoroalkyl substances (PFASs) in urban stormwater biofilter facilities is very limited. C5-14,16 perfluoroalkyl carboxylic acids [perfluorinated carboxylic acids (PFCAs)], C4,8,10 perfluoroalkanesulfonic acids (PFSAs), methyl-perfluorooctane sulfonamide acetic acid (MeFOSAA, a PFSA precursor), and unknown C6-8 PFCA and perfluorooctanesulfonic acid precursors were frequently found in bioretention media and forebay sediments at Σ35PFAS concentrations of <0.03-19 and 0.064-16 μg/kg-DW, respectively. Unknown C6-8 PFCA precursor concentrations were up to ten times higher than the corresponding PFCAs, especially at forebays and biofilters' top layer. No significant trend could be attributed to PFAS and precursor concentrations versus depth of filter media, though PFAS concentrations were 2-3 times higher in the upper layers on average (significant difference between the upper (0-5 cm) and deepest (35-50 cm) layer). PFASs had a similar spatial concentration distribution in each filter media (no clear difference between short- and long-chain PFASs). Commercial land use and organic matter were important factors explaining the concentration variations among the biofilters and between the sampling depths, respectively. Given the comparable PFAS accumulations in deeper and superficial layers and possible increased mobility after precursor biotransformation, designing shallow-depth, nonamended sand biofilters or maintaining only the top layer may be insufficient for stormwater PFAS management.

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.

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.

Per-and polyfluoroalkyl (PFAS) eternal pollutants: Sources, environmental impacts and treatment processes

Per- and polyfluoroalkyl substances (PFAS) have become a growing environmental concern due to their tangible impacts on human health. However, due to the large number of PFAS compounds and the analytical difficulty to identify all of them, there are still some knowledge gaps not only on their impact on human health, but also on how to manage them and achieve their effective degradation. PFAS compounds originate from man-made chemicals that are resistant to degradation because of the presence of the strong carbon-fluorine bonds in their chemical structure. This review consists of two parts. In the first part, the environmental effects of fluorinated compound contamination in water are covered with the objective to highlight how their presence in the environment adversely impacts the human health. In the second part, the focus is put on the different techniques available for the degradation and/or separation of PFAS compounds in different types of waters. Examples of removal/treatment of PFAS present in either surface or ground water are presented.

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.

Removal of per- and polyfluoroalkyl substances from wastewater via aerosol capture

The widespread use of per- and polyfluoroalkyl substances (PFAS)-containing products in numerous commercial and industrial applications has resulted in their occurrence in wastewater treatment plants (WWTPs). Herein, proof-of-concept bench-scale experiments were performed to measure the extent to which PFAS could be removed from a WWTP if aerosols generated during aeration were captured. Experiments were designed to mimic the aeration rate:water volume ratio, the water volume:surface area ratio, and aeration bubble size applicable to the full-scale aeration vessel. Results showed that substantial (75%) removal of perfluorooctane sulfonate (PFOS) was observed under these operating conditions in the bench-scale system; up to 97% PFOS removal was observed if the aeration rate was increased 3-fold. PFAS removal generally increased with increasing aerosol capture and with increasing PFAS surface activity. Analysis of semi-quantified PFAS showed that the semi-quantified PFAS accounted for approximately 93% of the identified PFAS in the raw wastewater, dominated largely by the presence of 2:2 fluorotelomer carboxylic acid (2:2 FTCA). This preliminary study suggests that aerosol capture in aeration basins has potential for mitigating PFAS in WWTPs. Further testing is needed to assess the feasibility of this approach at the field scale.

The Dynamics of Per- and Polyfluoroalkyl Substances (PFAS) at Interfaces in Porous Media: A Computational Roadmap from Nanoscale Molecular Dynamics Simulation to Macroscale Modeling

Managing and remediating perfluoroalkyl and polyfluoroalkyl substance (PFAS) contaminated sites remains challenging. The major reasons are the complexity of geological media, partly unknown dynamics of the PFAS in different phases and at fluid-fluid and fluid-solid interfaces, and the presence of cocontaminants such as nonaqueous phase liquids (NAPLs). Critical knowledge gaps exist in understanding the behavior and fate of PFAS in vadose and saturated zones and in other porous media such as concrete and asphalt. The complexity of PFAS-surface interactions warrants the use of advanced characterization and computational tools to understand and quantify nanoscale behavior of the molecules. This can then be upscaled to the microscale to develop a constitutive relationship, in particular to distinguish between surface and bulk diffusion. The dominance of surface diffusion compared to bulk diffusion results in the solutocapillary Marangoni effect, which has not been considered while investigating the fate of PFAS. Without a deep understanding of these phenomena, derivation of constitutive relationships is challenging. The current Darcy scale mass-transfer models use constitutive relationships derived from either experiments or field measurements, which makes their applicability potentially limited. Here we review current efforts and propose a roadmap for developing Darcy scale transport equations for PFAS. We find that this needs to be based on systematic upscaling of both experimental and computational studies from nano- to microscales. We highlight recent efforts to undertake molecular dynamics simulations on problems with similar levels of complexity and explore the feasibility of conducting nanoscale simulations on PFAS dynamics at the interface of fluid pairs.

Delineation of a PFOA Plume and Assessment of Data Gaps in its Conceptual Model Using PlumeSeeker™

An accurate conceptual site model (CSM) and plume-delineation at contamination sites are pre-requisites for successful remediation and for satisfying regulators and stakeholders. PlumeSeeker™ is well-suited for assessing data gaps in CSMs by using available site data and for identifying the optimal number and locations of sampling locations to delineate contaminant plumes. It is an enhancement of a university research code for plume delineation using geostatistical and stochastic modeling integrated with the groundwater modeling software MODFLOW-SURFACT™. PlumeSeeker™ increases the overall confidence in the location of the plume boundary through a variance-reduction approach that selects existing- or new monitoring wells for sampling based on minimizing the uncertainty in plume boundary and on new field information. Applicable at sites with or without existing monitoring wells, PlumeSeeker™ is particularly powerful for optimally allocating project resources (labor, well installation, and laboratory costs) between existing wells and sampling at new locations. An application of PlumeSeeker™ at Lakehurst, the naval component of Joint Base McGuire-Dix-Lakehurst in New Jersey, demonstrates how the cost of delineating the migration pathway of a perfluorooctanoic acid (PFOA) plume can be minimized by requiring only 9 new sampling locations in addition to samples from 2 existing wells for achieving a 70% reduction in plume uncertainty. In addition, the use of available site data in three different scenarios identified CSM data-gaps in the source area and in the interaction between Manapaqua Branch and groundwater, where the observed high concentration in this area could have resulted from a combination of groundwater migration and induced infiltration.

Enhanced aggregation and interfacial adsorption of an aqueous film forming foam (AFFF) in high salinity matrices

Per- and polyfluoroalkyl substances (PFAS) exist in contaminated groundwater, surface water, soil, and sediments from use of aqueous film forming foams (AFFFs). Under these conditions PFAS exhibit unusual behavior due to their surfactant properties, namely, aggregation and surface activity. Environmental factors such as salinity can affect these properties, and complicate efforts to monitor PFAS. The effect of high salinity matrices on the critical micelle concentration (CMC) of a AFFF formulation manufactured by 3M and the surface accumulation of PFAS was assessed with surface tension isotherm measurements and bench-scale experiments quantifying PFAS at the air-water interface. Conditions typical of brackish and saline waters substantially depressed the CMC of the AFFF by over 50% and increased the interfacial mass accumulation of PFAS in the AFFF mixture by up to a factor of 3, relative to values measured in ultrapure water. These results indicate that high salinity matrices increase the aggregation and surface activity of PFAS in mixtures, which are key properties affecting their 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.

Investigation of the adsorption/desorption mechanism of perfluoroalkyl substances on HLB-WAX extraction phases for microextraction

The C-F alkyl structural backbone of per- and polyfluoroalkyl substances makes this class of molecules resistant to heat and degradation, leading to their high persistence and mobility in the environment and bioaccumulation in the tissues of living organisms. In this study, 15 PFAS with an alkyl chain length from C₄ to C₁₄, currently monitored by the U.S. Environmental Protection Agency (EPA), were preconcentrated by solid-phase microextraction (SPME) and analyzed by liquid chromatography-tandem mass spectrometry. The adsorption and desorption mechanisms of PFAS onto ion-exchange extraction phases was evaluated to understand the extraction process of PFAS from various environmental matrices under different conditions. This was achieved using two SPME geometries, namely fibers and thin films. The use of thin films resulted in a twofold improvement in extraction efficiency compared to fibers, especially for the short-chain PFAS. Methanol:water (80:20, v/v) was chosen as the optimized desorption solution, with ammonium formate added to minimize carryover. Extraction time profiles for both SPME geometries showed faster equilibration with thin films (30 min) compared to fibers (90-120 min). The linear dynamic range obtained with this method using fibers and thin films ranged from 10 to 5000 ng L^-1 and 2.5-5000 ng L^-1, respectively, with acceptable accuracy (70-130%) and precision (<15%). LOD ranged within 2.5-10 ng L^-1 for fibers and 0.01-0.25 ng L^-1 for thin films. Investigating the factors affecting PFAS recovery in complex samples enabled the quantitative assessment of PFAS contamination in various environmental water samples such as seawater, melted snow and biospecimens like human plasma. A 96-SPME holder was used for validation, which is compatible with sampling in 96-well plates and ensures high throughput in the analysis of real samples. The total concentration of PFAS detected in seawater and snow was 51.3 ng L^-1 and 16.4 ng L^-1, respectively.

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.

Preferential Retention and Transport of Perfluorooctanesulfonic Acid in a Dolomite Aquifer

Per- and polyfluoroalkyl substances (PFAS) can represent a significant human health risk if present in aquifers used as a drinking water source. Accurate assessment of PFAS exposure risks requires an improved understanding of field-scale PFAS transport in groundwater. Activities at a former firefighter training site in University Park, Pennsylvania introduced perfluorooctanesulfonic acid (PFOS) to the underlying dolomite aquifer. Groundwater sampling from 2015 to 2018 delineated a PFOS plume with two concentration maxima located approximately 20 and approximately 220 m downgradient of the training site, separated by a zone of lower concentrations. We use a combination of analytical and numerical models, informed by independent measurements of aquifer porosity, hydraulic conductivity, and organic carbon content, to interpret the field observations. Our analysis demonstrates that preferential retention and transport resulting from simple heterogeneity in bedrock sorption, as caused by organic carbon (OC) content variability, provides a plausible explanation for plume separation. Dissolved PFOS partitions strongly to organic solids (high K_oc ), so even a small OC (<1 wt%) significantly retards PFOS transport, whereas zones with little to no OC allow for transport rates that approximate those of a conservative solute. Our work highlights an important consideration for modeling the groundwater transport of PFOS, and other compounds with high K_oc . In aquifers with discrete layers of varying OC, models using a uniform site-average OC will underestimate transport distances, thereby misrepresenting exposure risks for downgradient communities.

Effect of geochemical conditions on PFAS release from AFFF-impacted saturated soil columns

Per- and polyfluoroalkyl substances (PFASs) are frequently found at high concentrations in the subsurface of aqueous film forming foam (AFFF)-impacted sites. Geochemical parameters affect the release of PFASs from source area soils into groundwater but have not been extensively studied for soils that have been historically impacted with AFFF. This study investigated the effects of pH and salt concentrations on release of anionic and zwitterionic PFASs from AFFF-impacted soils in flow-through saturated columns. High pH (10) columns with elevated sodium concentrations had higher cumulative masses eluted of several PFASs compared to pH 3 and pH 7 columns with lower sodium concentrations, likely caused by changes to soil organic matter surface charge. Four PFASs (e.g. 4:2 fluorotelomer sulfonate, perfluorobutane sulfonamido acetic acid) eluted significantly earlier in both pH 3 and pH 10/high NaCl columns compared to pH 7 columns. The results of this study suggest that shifts in pH for soils located at AFFF-impacted sites - particularly raising the pH - may mobilize sorbed PFASs, specifically longer-chain and zwitterionic compounds that are typically strongly sorbed to soil.

Foam fractionation of per- and polyfluoroalkyl substances (PFASs) in landfill leachate using different cosurfactants

Foam fractionation is one solution to recover per- and polyfluoroalkyl substances (PFASs) from aqueous sources. The separation process is based on the sorption of PFASs to the air-water interface of bubbles. In many practical cases, the PFAS concentration in the polluted liquid is too low to sustain foam formation and requires the support of a cosurfactant not only to act as a collector of PFAS but also to produce and sustain foam for effective separation. However, there is a lack of information regarding the appropriate choice of cosurfactant and its quantitative effect on the interfacial partitioning of PFASs on the air bubbles. This study is directed to (i) evaluate the effectiveness of four cosurfactants with different-charged headgroups (i.e., anionic, cationic, zwitterionic and nonionic) for foam fractionation of PFASs, and (ii) estimate the air-water interfacial partitioning (K_i) of PFASs in the presence of four different types of cosurfactants. The K_i values span over 4 orders of magnitude with good correlation with PFASs molar volume. All of the cosurfactants were effective for the removal of the long chain PFASs (1.2-4 logs). The cationic and zwitterionic surfactants have oppositely charged head groups with respect to the anionic PFASs and therefore facilitate increased separation due to charge interactions. Some short chain PFASs (e.g., Perfluorobutanesulfonic acid (PFBS), Perfluoropentanesulfonic acid (PFPeS)) can be effectively removed using cationic and zwitterionic cosurfactants.

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.

Simulated leaching of PFAS from land-applied municipal biosolids at agricultural sites

Biosolids are an important resource for agricultural practice but have recently received increased focus as a potential source of per- and polyfluoroalkyl substances (PFAS) in the environment. Few studies have investigated the transport of PFAS through the unsaturated zone under conditions relevant to biosolids application sites. Herein, the unsaturated flow and transport model HYDRUS is used to evaluate the leaching of per- and polyfluoroalkyl substances (PFAS) from land-applied biosolids used in agricultural practice to determine the impacts of PFAS leaching on underlying groundwater resources. This numerical case study was based on conditions and operations at two test sites in central Illinois where biosolids were applied at agronomic rates and where PFAS contents and desorption characteristics were previously characterized. Each site possessed different vadose zone soil textural heterogeneity. Simulations were performed under actual present-day meteorological conditions and extended 150 years beyond the initial biosolids application. These long-term simulations demonstrate how soil equilibrium sorption/desorption processes within the biosolids-amended surface soils effectively control the transport rate of individual PFAS to groundwater. Air-water interfacial (AWI) adsorption, which is sometimes considered to be a significant source of PFAS retention in vadose zone soils, was observed to have minimal impacts on PFAS leaching rates within the biosolids-amended surface soils at these sites. Additionally, the impact of AWI adsorption was found to be most significant for PFAS transport within the underlying vadose zone soils when these soils were more texturally homogeneous and considerably less significant within the texturally heterogeneous soils represented herein. The results of multiple long-term simulations were used to develop an empirical equation that relates predicted maximum PFAS pore-water concentrations reaching the saturated zone with changes in PFAS concentrations in the biosolids-amended soil for various biosolids re-application events. This approach is shown to be very useful in developing site-specific PFAS soil screening levels and/or maximum leachate levels for PFAS in support of establishing best management practices (BMPs) for land application of biosolids.

Calculating PFAS interfacial adsorption as a function of salt concentration using model parameters determined from chemical structure

Per- and polyfluoroalkyl substances (PFAS) are widely-detected environmental contaminants known to concentrate at surfaces and interfaces. Many of the most commonly-detected PFAS function as ionic surfactants under environmental conditions. The interfacial behaviors of ionic surfactants, including PFAS, are strongly dependent on salt concentration and composition, with interfacial affinity potentially varying by orders of magnitude for the same compound under different conditions. The work described here presents a tool for predicting the salt-dependent adsorption of PFAS compounds based entirely on chemical structure, something of great value for predicting the real-world environmental behavior of many of the large numbers of PFAS compounds for which experimental data are not available. The approach combines two different previously-developed models, one a mass-action model designed to predict the effects of salts on interfacial adsorption of ionic PFAS (the UNSW-OU salt model), and the second a group-contribution model designed to predict interfacial adsorption of PFAS in the absence of salt based on chemical structure. The challenge of combining the two models comes from the fact that both are based on different isotherms. The salt model can produce sigmoidal isotherms under salt-limited conditions (an isotherm shape that is supported by experimental evidence), while the group-contribution model can generate Langmuir parameters from calculations based on chemical structure. Equations were derived to determine salt model isotherm parameters from Langmuir parameters (either from the group-contribution model or experimental sources) by matching surface tension curves in the vicinity of the concentration of highest second derivative. Refined group-contribution model parameters were determined based on data from an additional 40 surface tension curves to allow improved structure-based predictions for important classes of PFAS that were not sufficiently well-represented in the original model. The resulting equations provide a tool allowing quantitative predictions of PFAS behavior under realistic environmental conditions for compounds for which little or no experimental data are available.

Retention of PFOS and PFOA Mixtures by Trapped Gas Bubbles in Porous Media

The transport of per- and polyfluoroalkyl substances (PFAS) in soil and groundwater is important for site investigation, risk characterization, and remediation planning. The adsorption of PFAS at air-water interfaces has been shown to significantly contribute to PFAS retention, with subsequent effects on concentrations and the time scales of transport. In this study, column experiments were conducted to investigate the transport of perfluorooctanesulfonic acid (PFOS), perfluorooctanoic acid (PFOA), and 6:2 fluorotelomer sulfonate (6:2 FTS) individually and in binary mixtures in the presence of a trapped gas phase, using clean sands to isolate adsorption to air-water interfaces. Consistent with previous studies, the transport of PFOS, PFOA, and 6:2 FTS was retarded by adsorption at the air-water interface, with greater retention of PFOS due to its higher affinity for the air-water interface. Chromatographic separation occurred in the experiments using binary mixtures of PFOS and PFOA, with greater retention at lower influent concentrations. The mixture experiments also showed enhanced breakthrough of PFOA in the presence of PFOS, where effluent concentrations of PFOA were temporarily greater than the influent concentration prior to the breakthrough of PFOS. This enhanced breakthrough was attributed to competition between PFOS and PFOA for adsorption to the air-water interface.

Enrichment of poly- and perfluoroalkyl substances (PFAS) in the surface microlayer and foam in synthetic and natural waters

Bench-scale experiments were performed to interrogate poly- and perfluoroalkyl substance (PFAS) enrichment in the water surface microlayer (SML). In initial experiments using electrolyte-only solutions, the perfluorooctane sulfonate (PFOS) and perfluorooctane carboxylate (PFOA) enrichment in the SML were reasonably (with a factor of 2) described by the Gibbs adsorption equation coupled with a Freundlich-based interfacial adsorption model. Enrichment in the SML among perfluorinated sulfonates and perfluorinated carboxylates of varying chain lengths was proportional to their surface activity. The PFOS enrichment factor (EF), defined as the PFAS concentration in the SML divided by the concentration in the bulk water, was 18 in a 200 mg/l NaCl solution. The presence of elevated organic carbon levels in synthetic surface waters inhibited PFAS accumulation in the SML, with resulting EF values of approximately 1 for all PFAS. However, in the presence of elevated organic levels coupled with foam, PFAS enrichment in the foam was observed, with a foam EF of 25 measured for PFOS in synthetic surface waters. PFAS EF values measured in several natural surface waters without foam showed little variation among the waters tested, with PFOS EF values ranging between 6 and 10. Together, these results suggest that PFAS accumulation in the SML is largely controlled by PFAS sorption at the air-water interface for the conditions examined in this study, and the presence of foam with natural organics enhances PFAS uptake at the water surface.

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.

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.

Influence of Residual Nonaqueous-Phase Liquids (NAPLs) on the Transport and Retention of Perfluoroalkyl Substances

Per- and polyfluoralkyl substances (PFAS) are known to accumulate at interfaces, and the presence of nonaqueous-phase liquids (NAPLs) could influence the PFAS fate in the subsurface. Experimental and mathematical modeling studies were conducted to investigate the effect of a representative NAPL, tetrachloroethene (PCE), on the transport behavior of PFAS in a quartz sand. Perfluorooctanesulfonate (PFOS), perfluorononanoic acid (PFNA), a 1:1 mixture of PFOS and PFNA, and a mixture of six PFAS (PFOS, PFNA, perfluorooctanoic acid (PFOA), perfluoroheptanoic acid (PFHpA), perfluorohexanesulfonate (PFHxS), and perfluorobutanesulfonate (PFBS)) were used to assess PFAS interactions with PCE-NAPL. Batch studies indicated that PFAS partitioning into PCE-NAPL (K_nw < 0.1) and adsorption on 60-80 mesh Ottawa sand (K_d < 6 × 10^-5 L/g) were minimal. Column studies demonstrated that the presence of residual PCE-NAPL (∼16% saturation) delayed the breakthrough of PFOS and PFNA, with minimal effects on the mobility of PFBS, PFHpA, PFHxS, and PFOA. Breakthrough curves (BTCs) obtained for PFNA and PFOS alone and in mixtures were nearly identical, indicating the absence of competitive adsorption effects. A mathematical model that accounts for NAPL-water interfacial sorption accurately reproduced PFAS BTCs, providing a tool to predict PFAS fate and transport in co-contaminated subsurface environments.

Predicting the impact of salt mixtures on the air-water interfacial behavior of PFAS

Salts are known to have strong impacts on environmental behavior of per- and polyfluoroalkyl substances (PFAS) including air-water interfacial adsorption. Multivalent salts impact interfacial adsorption to a greater extent than monovalent salts. Models to make a priori predictions of PFAS interfacial adsorption in the presence of multiple salts with different ionic charges are needed given the need to predict PFAS environmental fate. This study further develops a mass-action model to predict the interfacial behavior of PFAS as a function of both salt valency and concentration. The model is validated using surface tension data for a series of monovalent and divalent salt mixtures over a wide range of ionic strengths (i.e., from no added salt to 0.5 M) as well as comparison to data from literature. This model highlights the disproportionate impact of multivalent salts on interfacial adsorption and the practical utility of the model for predicting interfacial adsorption in the presence of multiple monovalent and multivalent inorganic salts. Results suggest that failure to account for divalent salt, even when concentrations are much smaller than monovalent salt, under most environmentally relevant aqueous phase conditions will result in significant underpredictions of PFAS interfacial adsorption. Simple examples of PFAS distribution in a range of salt conditions in the vadose zone and in aerated-water treatment reactors highlight the predictive utility of the model.

Comparison of methods to estimate air-water interfacial areas for evaluating PFAS transport in the vadose zone

When performing calculations or numerical simulations for the fate and transport of PFAS and other surface-active solutes in the vadose zone, accurately representing the relationship between the area of the air-water interfaces (A_aw) as a function of water saturation (S_w), and changes in that relationship resulting from changes in soil texture, are equally important as accurately characterizing interfacial adsorption coefficients and the concentration dependence for PFAS solutes. This is true because the magnitude of the A_aw directly governs the degree of air-water interfacial adsorption, which contributes to the transport retardation of these solutes within unsaturated porous media. Herein, a well-known thermodynamic-based model for predicting the A_aw-S_w relationship is evaluated through comparisons to literature data collected using various measurement techniques for model sands and a limited number of soils using data collected from the current published literature. This predictive model, herein termed the Leverett thermodynamic model (LTM), relies on the characterization of the soil-water retention curve (SWRC) for a given soil, using the van Genuchten (VG) equation for the pressure head-vs-S_w relationship. Therefore, methods to estimate the VG equation parameters are also compared as to the A_aw-S_w relationships predicted. Comparisons suggest that the LTM provides the best estimate of the actual A_aw-S_w relationships for water containing non-surface-active solutes. Because PFAS solutes are also surface-active, A_aw measurement methods utilizing surface-active tracers are considered to provide the most accurate representation of the A_aw-S_w relationship for these solutes. Differences between A_aw-S_w relationships derived from tracer methods and the LTM are described in relation to media surface roughness effects. Based on the available literature data, a practical empirical model is proposed to adjust the LTM prediction to account for the effects of surface roughness on the magnitude of the A_aw for surface-active solutes. Finally, example retention calculations are performed to demonstrate the sensitivity of the predicted A_aw-S_w relationship on the vadose zone transport of of a representative PFAS, perfluorooctane sulfonate.

A new framework for modeling the effect of salt on interfacial adsorption of PFAS in environmental systems

Per- and polyfluoroalkyl substances (PFAS) are surface active contaminants of great environmental concern, due to their widespread historical use and their environmental persistence. Salts are known to have a profound influence on the interfacial behaviors of all ionic surfactants, including some of the most commonly detected PFAS. This work describes a new mass-action model for predicting the interfacial behavior of surfactants as a function of salt concentration. The three-parameter model is fit to interfacial tension data over a range of salt concentrations, and is then able to predict interfacial adsorption isotherms for the entire range from no added salt, up to 0.5 M added salt. The phenomenological nature of the model means that it is likely to provide more robust predictions for new systems and conditions than some of the existing empirical approaches, and the minimal number of adjustable parameters ensures that unique calibrations are possible with limited data. The model is found to be consistent with experimental data, and is bracketed by experimental values at low PFAS concentrations. Of particular interest, the model predicts the existence of sigmoidal adsorption isotherms at low salt concentrations, a deviation from isotherms calculated the commonly-used Szyszkowski equation; the observation is supported by a maximum in measured interfacial adsorption coefficient calculated from low-concentration surface tension measurements. Because adsorption affinities can vary by orders of magnitude with changing salt concentration, the ability to predict the effects of salt on adsorption is of critical importance for quantitative prediction of PFAS behavior in the environment.

Unification of surface tension isotherms of PFOA or GenX salts in electrolyte solutions by mean ionic activity

The surface tension isotherms of soluble salts of per- and polyfluoroalkyl substances (PFAS) in electrolyte solutions are typically reported as functions of the PFAS concentration. However, for univalent salts and electrolytes, the Langmuir-Szyszkowski equation is a function of the mean ionic activity a*. Using previously reported data, we show that for salts of perfluorooctanoic acid (PFOA) or hexafluoropropylene oxide dimer acid (GenX™), use of a* rather than concentration provides a unified surface tension isotherm, independent of the electrolyte concentration. This suggests that the electrolyte dependence of the isotherm arises purely from its effect on PFAS activity, rather than an intrinsic surface property. This finding has important implications for the understanding of PFAS retention in saline unsaturated soils, and for PFAS extraction from saline waters by foam fractionation.

The influence of molecular structure on PFAS adsorption at air-water interfaces in electrolyte solutions

Fluid-fluid interfacial adsorption has been demonstrated to be an important retention process for per and polyfluoroalkyl substances (PFAS) in porous media with air or non-aqueous phase liquids (NAPLs) present. The objective of this study was to characterize the influence of PFAS molecular structure on air-water interfacial adsorption in electrolyte solutions. Measured and literature-reported surface-tension data sets were aggregated to generate the largest compilation of interfacial adsorption coefficients measured in aqueous solutions comprising environmentally representative ionic strengths. The surface activities and interfacial adsorption coefficients (Ki) exhibited chain length trends, with greater surface activities and larger Ki values corresponding to longer chain length. The impact of multiple-component PFAS solutions on the surface activity of a select PFAS was a function of the respective surface activities and concentrations. Quantitative structure-property relationship analysis (QSPR) employing a single molecular descriptor (molar volume) was used successfully to characterize the impact of PFAS molecular structure on air-water interfacial adsorption. A previously reported QSPR model based on PFAS data generated for deionized-water solutions was updated to include more than 60 different PFAS, comprising all head-group types and a wide variety of tail structures. The QSPR model developed for PFAS in electrolyte solution compared favorably to the model developed for deionized water. Additionally, the magnitude of ionic strength for non-zero ionic strength systems was determined to have relatively minimal impact on interfacial adsorption coefficients. The new QSPR model is therefore anticipated to be representative for a wide variety of PFAS and for a wide range of ionic compositions.

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.

Fate and transport of per- and polyfluoroalkyl substances (PFASs) in the vadose zone

Per- and polyfluoroalkyl substances (PFASs) are a heterogeneous group of persistent organic pollutants that have been detected in various environmental compartments around the globe. Emerging research has revealed the preferential accumulation of PFASs in shallow soil horizons, particularly at sites impacted by firefighting activities, agricultural applications, and atmospheric deposition. Once in the vadose zone, PFASs can sorb to soil, accumulate at interfaces, become volatilized, be taken up in biota, or leach to the underlying aquifer. At the same time, polyfluorinated precursor species may transform into highly recalcitrant perfluoroalkyl acids, changing their chemical identity and thus transport behavior along the way. In this review, we critically discuss the current state of the knowledge and aim to interconnect the complex processes that control the fate and transport of PFASs in the vadose zone. Furthermore, we identify key challenges and future research needs. Consequently, this review may serve as an interdisciplinary guide for the risk assessment and management of PFAS-contaminated sites.

Influence of aqueous film forming foams on the solubility and mobilization of non-aqueous phase liquid contaminants in quartz sands

At sites where aqueous film forming foams (AFFFs) are used for fire suppression or training activities, interactions between dissolved foam constituents and organic liquids could alter contaminant migration in the subsurface. In this study, batch reactor and column experiments were conducted to investigate the potential for AFFF solutions to enhance the solubility and mobility of three representative non-aqueous phase liquid (NAPLs), JP-4 jet fuel, trichloroethene (TCE), and tetrachloroethene (PCE). For AFFF concentrations up to 5% wt. (50,000 mg/L), aqueous solubilities of TCE and PCE increased by less than 50%, indicating the absence of micellar solubilization. However, NAPL-water interfacial tensions were reduced to less than 1.5 mN/m and resulted in accumulation of up to 2.25 mg/m² of AFFF at the NAPL-water interface. To assess the potential for AFFF to mobilize residual (entrapped) NAPL at a field application rate of 3% wt. (30,000 mg/L), columns were packed with two size fractions of Ottawa sands (20-30 mesh and 60-80 mesh) that yielded residual NAPL saturations ranging from 11.7 to 17.6%. Following injection of 3 pore volumes of the 3% wt. AFFF solution, partial mobilization of residual NAPL was observed for PCE, TCE, and JP-4, with saturation reductions of 0.7 to 2% in 20-30 mesh and 0.3% to 1.3% in 60-80 mesh Ottawa sand. The columns were then flushed with an ultralow-IFT surfactant solution consisting of 4% wt. 1:1 Aerosol AY/OT, which resulted in nearly complete mobilization of the remaining residual NAPL. When NAPL desaturation curves were expressed in terms of the total trapping number (N_T), the threshold value of N_T required for NAPL mobilization by 3% wt. AFFF was approximately 2 × 10^-5, consistent with previous studies. These findings demonstrate that AFFF solutions can lead to partial mobilization of residual NAPL, and that the total trapping number concept can be used to predict such behavior a priori. In addition, the observed IFT reductions are sufficient to alter NAPL accumulation and redistribution behavior in the subsurface through enhanced spreading above the water table (e.g., JP-4) or penetration into lower-permeability media (e.g., TCE, PCE).

A group-contribution model for predicting the physicochemical behavior of PFAS components for understanding environmental fate

The factors controlling per- and polyfluoroalkyl substances (PFAS) environmental fate remains the subject of considerable debate and study. As surfactants, PFAS readily partition to interfaces, a property that controls their transport and fate. A group contribution model is developed to predict the extent to which PFAS partitions to the air-water interface. Langmuir adsorption and Szyszkowski equation parameters were fitted to literature air-water surface tension data for a range of PFAS and conventional hydrocarbon surfactants. This approach enabled the prediction of the impact of the hydrophilic head group, and other molecular components, on PFAS interfacial partitioning in instances when PFAS data are unavailable but analogous hydrocarbon surfactant data are available. The model was extended to predict a range of parameters (i.e., solubility, critical micelle concentration (CMC), K_D, K_oc and K_ow) that are used to predict PFAS environmental fate, including long-range PFAS transport and in multimedia models. Model predictions were consistent with laboratory and field derived parameters reported in the literature. Additionally, the proposed model can predict the impact of pH and speciation on the extent of PFAS interfacial partitioning, a potentially important feature for understanding the behaviors of some ionizable PFAS, such as fluorinated carboxylic acids. The proposed model provides a conceptually straightforward method to predict a wide range of environmental fate parameters for a wide range of PFAS. As such, the model is a powerful tool that can be used to determine parameters needed to predict PFAS environmental fate.

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.

Air-water interfacial adsorption coefficients for PFAS when present as a multi-component mixture

Surface tension isotherms and calculated air-water interfacial (AWI) adsorption data are presented for solution mixtures of per- and polyfluoroalkyl substances (PFAS), specifically a series of binary and one ternary mixtures of homologous linear perfluorocarboxylic acids (PFCAs) in a simulated groundwater, and two 8-component mixtures containing both PFCAs and linear perfluoroalkane sulfonates (PFSAs). In all cases, non-ideal competitive adsorption was observed that favored the most surface-active component(s) of the solution mixture. The multi-component extended Langmuir (EL) isotherm model was observed to accurately predict the competitive adsorption observed in the binary and ternary PFCA solution mixtures. However, the predictive utility of the EL model was observed to diminish when mixtures contained both PFCAs and PFSAs, which differ in their hydrophile structure, resulting in overpredictions and underpredictions of the AWI adsorption isotherms derived from measured data depending on the specific components present in the solution mixtures. Observations indicate that the individual component adsorptive affinities for the AWI can change in response to competitive preferential adsorption as their solution concentrations increase that is not being captured by the EL model. Our results demonstrate that alternative mathematical models are needed that support concentration dependent affinity coefficients for non-similar mixtures of PFAS, such that the transport of individual target PFAS components within a larger mixture of components can be accurately predicted across a wider range of solution concentration.

Aqueous Film-Forming Foams Exhibit Greater Interfacial Activity than PFOA, PFOS, or FOSA

Perfluoroalkyl acids spontaneously concentrate at air-water and non-aqueous phase liquid (NAPL)-water interfaces, which can influence their retention during subsurface transport. This work presents measurements of air- and NAPL-water interfacial tension for synthetic groundwater containing perfluorooctanoic acid (PFOA), perfluorooctanesulfonic acid (PFOS), perfluorooctanesulfonamide (FOSA), or aqueous film-forming foam (AFFF) formulations at concentrations ranging from 0.1 to greater than 1000 mg/L. The NAPLs tested included dodecane, tetrachloroethylene, and jet fuel. AFFF formulations were less efficient at lowering interfacial tension than PFOA, FPOS, or FOSA substances below 100 mg/L, while above 100 mg/L, these formulations were more effective, achieving tensions of less than 3 mN/m. Infiltration of solutions with such low tension could lead to mobilization of residual NAPL. Equations based on interfacial tension measurements show that concentrations of PFOA, PFOS, and FOSA at the air-water interface were from 2 to 16 times greater than at the NAPL-water interface below 100 mg/L and were 10-50 times greater for AFFF below 20 mg/L. Calculations for unsaturated soil estimate that up to 87% of PFOS mass was at the air-water interface and less than 4% at the dodecane-water interface for bulk-water concentrations below 1 mg/L.

A Modified HYDRUS Model for Simulating PFAS Transport in the Vadose Zone

The HYDRUS unsaturated flow and transport model was modified to simulate the effects of non-linear air-water interfacial (AWI) adsorption, solution surface tension-induced flow, and variable solution viscosity on the unsaturated transport of per- and polyfluoroalkyl substances (PFAS) within the vadose zone. These modifications were made and completed between March 2019 and May 2019, and were implemented into both the one-dimensional (1D) and two-dimensional (2D) versions of HYDRUS. Herein, the model modifications are described and validated against the available literature-derived PFAS transport data (i.e., 1D experimental column transport data). The results of both 1D and 2D example simulations are presented to highlight the function and utility of the model to capture the dynamic and transient nature of the temporally and spatially variable interfacial area of the AWI (Aaw) as it changes with soil moisture content (Θw) and how it affects PFAS unsaturated transport. Specifically, the simulated examples show that while AWI adsorption of PFAS can be a significant source of retention within the vadose zone, it is not always the dominant source of retention. The contribution of solid-phase sorption can be considerable in many PFAS-contaminated vadose zones. How the selection of an appropriate Aaw(Θw) function can impact PFAS transport and how both mechanisms contribute to PFAS mass flux to an underlying groundwater source is also demonstrated. Finally, the effects of soil textural heterogeneities on PFAS unsaturated transport are demonstrated in the results of both 1D and 2D example simulations.

The influence of solution chemistry on air-water interfacial adsorption and transport of PFOA in unsaturated porous media

There is great interest in the transport behavior of PFAS in the vadose zone, and the impact of leaching on groundwater contamination. Air-water interfacial adsorption is an important process for PFAS retention in unsaturated porous media, and it is influenced by many factors including solution conditions such as ionic strength. The present study employed miscible-displacement column experiments to investigate the impact of ionic strength and pH on perfluorooctanoic acid (PFOA) retardation and transport under dynamic water-flow conditions. The results showed that retardation under unsaturated conditions was affected significantly by changes in ionic strength, whereas there was minimal impact for saturated conditions. This indicates that air-water interfacial adsorption, which was a major source of retardation, was influenced significantly by changes in ionic strength while they had a minor impact on solid-phase adsorption. The impact of changes in ionic strength on the magnitude of air-water interfacial adsorption observed for the column experiments was consistent with measured surface-tension data. The impact of changes in pH was less significant compared to that of ionic strength for transport under unsaturated conditions. These results illustrate the influence of solution chemistry on PFAS adsorption and transport under unsaturated conditions. This solution-dependent behavior should be considered when characterizing PFAS transport in soils and the vadose zone.

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.

Remediation of trapped DNAPL enhanced by SDS surfactant and silica nanoparticles in heterogeneous porous media: experimental data and empirical models

The remediation of nonaqueous phase liquids (NAPLs) enhanced by surfactant and nanoparticles (NP) has been investigated in numerous studies. However, the role of NP-assisted surfactants in the dissolution process is still not well discussed. Besides, there is a lack of empirical dissolution models considering the effects of initial residual saturation S_trap, NAPL distribution, and surfactant concentration in NAPL-aqueous phase systems. In this work, micromodel experiments are conducted to quantify mass transfer coefficients for different injected aqueous phases including deionized water, SDS surfactant solutions, and NP-assisted solutions with different levels of concentrations and flow rates. Observations reveal that silica nanoparticles (SNP) can significantly enhance interphase mass transfer, while SDS surfactant reduces the mass transfer coefficient. In addition, S_trap and intrinsic interfacial area a_i, as an indicator of dense nonaqueous phase liquids (DNAPL) distribution, influence the interphase mass transfer. The a_i is also independent of DNAPL saturation S_NAPL except for S_NAPL < 7% when ganglia breakup occurs. Based on these observations, new empirical dissolution models are proposed in the presence and the absence of SDS surfactant and SNP in which a_i, S_trap, and surfactant concentrations are introduced as new parameters. The evaluated mass transfer rate coefficients using the proposed models show a significant improvement compared to available empirical models. The finding of this study might be attractive for application in field-scale simulations of surfactant-enhanced aquifer remediation (SEAR) and NP-assisted methods.

Uptake of Poly- and Perfluoroalkyl Substances at the Air-Water Interface

Bench-scale experiments were performed to assess uptake of poly- and perfluoroalkyl substances (PFAS), both single compounds and mixtures, at the air-water interface. The focus was on evaluating uptake at field-relevant PFAS concentrations (<2 × 10^-4 mol m^-3 or 0.1 mg L^-1), assessing the impacts of various PFAS mixtures, and quantifying the impacts of background NaCl concentrations. Both interfacial tension measurements and direct quantification of PFAS mass sorbed at the air-water interface in water films were used to evaluate PFAS interfacial partitioning. Results showed that a Freundlich-based model, rather than a Langmuir-based model, described perfluorooctanoic acid (PFOA) and perfluorooctanesulfonic acid (PFOS) interfacial uptake. At lower and field-relevant PFOS and PFOA concentrations, the Langmuir-based model underpredicted interfacial uptake by up to several orders of magnitude. The interfacial partition coefficient, k_aw, increased as PFAS concentrations decreased. Results also showed that the interfacial tension and interfacial uptake of PFAS mixtures were (within a factor of 2) predicted based on the single solute systems assuming ideal dilute behavior. Furthermore, the addition of NaCl at concentrations of up to 0.01 M increased PFOS uptake by less than a factor of 2 at field-relevant PFOS concentrations. The results presented herein have important implications for PFAS migration in unsaturated soils as well as for remedial technologies that rely on PFAS interfacial sorption, particularly at field-relevant PFAS concentrations.