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

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.

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.

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.

New simplified design methods for engineering barriers around contaminated sites with Cauchy boundaries

Since the 1980s, low-permeability slurry trench cutoff walls have been widely constructed as barriers to retard the migration of contaminants. The thickness of the cutoff walls is a key determinant of the wall service life. Through a series of theoretical derivations, simplified methods for determining the flux limit and concentration limit were proposed to determine the thickness of cutoff walls for contaminated sites with constant pollutant flux. The relative errors of both the flux-based and concentration-based methods increase as the breakthrough criterion of the ratio between the specified limit concentration of the contaminant to the source concentration (C*) and the ratio of the limited value of contaminant flux to the constant source flux (F*) increases, with a given Peclet number PL. The maximum relative error reaches 4% and 6% when C* and F* are both 0.1, which covers most practical situations in cutoff wall design. Good agreements of wall thickness were obtained between the proposed simplified methods and analytical solutions via a clear example. The proposed method can efficiently simplify the design process of cutoff walls with high accuracy, providing a basis for containing contaminated sites.

Repeated Aqueous Film-Forming Foams Applications: Impacts on Polyfluoroalkyl Substances Retention in Saturated Soil

Historical practices at firefighter-training areas involved repeated aqueous film-forming foams (AFFFs) applications, resulting in source zones characterized by high concentrations of perfluoroalkyl and polyfluoroalkyl substances (PFAS). Repeated applications of AFFF composed of 14 anionic and 23 zwitterionic perfluoroalkyl substances (PFAS) were conducted on a single one-dimensional saturated soil column to quantify PFAS retention. An electrofluorination-based (3M) Milspec AFFF, which was above the mixture's critical micelle concentration (CMC), was at application strength (3%, v/v). Retention and retardation of PFAS mass increased with each successive AFFF addition, although the PFAS concentration profiles for subsequent applications differed from the initial. Greater degree of mass retention and retardation correlated with longer PFAS carbon-fluorine chain length and charged-headgroup type and as a function of AFFF application number. Anionic PFAS were increasingly retained with each subsequent AFFF application, while zwitterionic PFAS exhibited an alternating pattern of sorption and desorption. Surfactant-surfactant adsorption and competition during repeat AFFF applications that are at concentrations above the CMC resulted in adsorbed PFAS from the first application, changing the nature of the soil surface with preferential sorption of anionic PFAS and release of zwitterionic PFAS due to competitive elution. Applying a polyparameter quantitative structure-property relationship developed to describe sorption of AFFF-derived PFAS to uncontaminated, saturated soil was attempted for our experimental conditions. The model had been derived for data where AFFF is below the apparent CMC and our experimental conditions that included the presence of mixed micelles (aggregates consisting of different kinds of surfactants that exhibit characteristics properties different from micelles composed of a single surfactant) resulted in overall PFAS mass retained by an average of 27.3% ± 2.7% (standard error) above the predicted values. The correlation was significantly improved by adding a "micelle parameter" to account for cases where the applied AFFF was above the apparent CMC. Our results highlight the importance of interactions between the AFFF components that can only be investigated by employing complex PFAS mixtures at concentrations present in actual AFFF at application strength, which are above their apparent CMC. In firefighter-training areas (AFFF source zones), competitive desorption of PFAS may result in downgradient PFAS retention when desorbed PFAS become resorbed to uncontaminated soil.

Potential impact of bacteria on the transport of PFAS in porous media

The transport and fate of per- and poly-fluoroalkyl substances (PFAS) in soil and groundwater is a topic of critical concern. A number of factors and processes may influence the transport and fate of PFAS in porous media. One factor that has received minimal attention to date is the impact of bacteria on the retention and transport of PFAS, which is the focus of this current study. The first part of this work comprised a critical review of prior studies to delineate observed PFAS-bacteria interactions and to summarize the mechanisms of PFAS sorption and retention by bacteria. Retention of PFAS by bacteria can occur through sorption onto cell surfaces and/or by incorporation into the cell interior. Factors such as the molecular structure of PFAS, solution chemistry, and bacterial species can affect the magnitude of PFAS sorption. The influence of bacteria on the retention and transport of PFAS was investigated in the second part of the study with a series of batch and miscible-displacement experiments. Batch experiments were conducted using Gram-negative Pseudomonas aeruginosa and Gram-positive Bacillus subtilis to quantify the sorption of perfluorooctane sulfonic acid (PFOS). The results indicated that both bacteria showed strong adsorption of PFOS, with no significant difference in adsorption capacity. Miscible-displacement experiments were then conducted to examine the retention and transport of PFOS in both untreated sand and sand inoculated with Pseudomonas aeruginosa or Bacillus subtilis for 1 and 3 days. The transport of PFOS exhibited greater retardation for the experiments with inoculated sand. Furthermore, the enhanced sorption was greater for the 3-day inoculation compared to the 1-day, indicating that biomass is an important factor affecting PFOS transport. A mathematical model representing transport with nonlinear and rate-limited sorption successfully simulated the observed PFOS transport. This study highlights the need for future studies to evaluate the effect of bacteria on the transport of PFAS in soil and groundwater.

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.

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.

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.

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.

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.

Modelling the sorption behaviour of perfluoroalkyl carboxylates and perfluoroalkane sulfonates in soils

A simple parametric model was developed to predict the sorption of perfluoroalkyl substances (PFASs) in soils. Initially, sorption and desorption solid-liquid distribution coefficients (K_d and K_d,des respectively) of eight PFASs (five perfluoroalkyl carboxylates, PFCAs, and three perfluoroalkane sulfonates, PFSAs) in seven soils with organic carbon (OC) content ranging from 1.6 to 41% were quantified using batch experiments. The information obtained helped to fill the gaps in a literature-based database of K_d values of PFASs, which was lacking data on soils with high OC content. The overall dataset finally comprised 435 entries. Normalized sorption coefficients for the soil OC and mineral fraction contents (K_OC and K_MIN respectively) were deduced for each PFAS by correlating the corresponding K_d values obtained under a wide range of experimental conditions with the fraction of organic carbon (f_OC) of the soils. Furthermore, the sorption mechanisms in each phase were shown to depend mainly on PFAS chain length. The dependence of K_OC and K_MIN values on PFAS chain length defined the basic equations to construct the model for predicting PFAS sorption, applicable to both PFCAs and PFSAs with chain lengths ranging from 3 to 11 fluorinated carbons. The validation of the proposed model confirmed its ability to predict the K_d of PFASs based only on the soil OC and silt+clay contents and PFAS chain length. Therefore, it can be used in the first stages of a risk assessment process aiming at estimating the potential mobility of PFASs in soils after a contamination event. SYNOPSIS: This study develops a new parametric model to predict the sorption of perfluoroalkyl substances (PFASs) in soils.

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.

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.