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.

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

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

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

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

Estimation of Transport Parameters of Perfluoroalkyl Acids (PFAAs) in Unsaturated Porous Media: Critical Experimental and Modeling Improvements.

Predicting the transport of perfluoroalkyl acids (PFAAs) in the vadose zone is critically important for PFAA site cleanup and risk mitigation. PFAAs exhibit several unusual and poorly understood transport behaviors, including partitioning to the air-water interface, which is currently the subject of debate. This study develops a novel use of quasi-saturated (residual air saturation) column experiments to estimate chemical partitioning parameters of both linear and branched perfluorooctane sulfonate (PFOS) in unsaturated soils. The ratio of linear-to-branched air-water interfacial partitioning constants for all six experiments was 1.62 ± 0.24, indicating significantly greater partitioning of linear PFOS isomers at the air-water interface. Standard breakthrough curve analysis and numerical inversion of HYDRUS models support the application of a Freundlich isotherm for PFOS air-water interfacial partitioning below a critical reference concentration (CRC). Data from this study and previously reported unsaturated column data on perfluorooctanoate (PFOA) were reevaluated to examine unsaturated systems for transport nonidealities. This reanalysis suggests both transport nonidealities and Freundlich isotherm behavior for PFOA below the CRC using drainage-based column methods, contrary to the assertions of the original authors. Finally, a combined Freundlich-Langmuir isotherm was proposed to describe PFAA air-water interfacial partitioning across the full range of relevant PFAA concentrations.

Subsurface transport potential of perfluoroalkyl acids (PFAAs): Column experiments and modeling.

Transport of ten perfluoroalkyl acids (PFAAs) was studied with one-dimensional (1-D) saturated column experiments using four soil types with an organic carbon fraction (foc) range of ~0-0.045. Columns were operated under conditions relevant to aqueous film-forming foam (AFFF)-impacted fire protection training areas to determine the ability of equilibrium transport parameters to describe 1-D PFAA transport, if rate-limited sorption influences PFAA transport, and if kinetic parameters can be used to evaluate factors causing rate-limited sorption. Results of initial screening of PFAA breakthrough found that over half of the breakthrough curves deviated from equilibrium transport and merited further investigation. Subsequent analysis showed that, in many cases, these deviations could be accounted for by considering the range of applicable equilibrium Kd values (i.e. based on standard deviation) applicable to the solid phase. Thus, transport of the majority of PFAAs in 3 soils with foc of 0-0.017 was not impacted by rate-limited sorption. Further, low sorption led to transport that was essentially simultaneous for the majority of PFAAs in these porous media. Exceptions were observed for long-chain PFAAs, and also in a fourth soil with foc of 0.045, which indicated the potential for rate-limited sorption to impact transport in some scenarios. Subsequent flow interruption experiments isolating kinetic behavior confirmed rate-limited sorption caused nonequilibrium transport. Linear free energy relationships (LFERs) developed in previous work to predict the inverse relationship between mass transfer coefficients (k) and sorption parameters (i.e., Kd) were used to estimate values of k for PFAAs in this study. Resulting k values were 10-3 to 10-8 h-1, consistent with previously measured kinetic parameters for other polar and anionic compounds. Models incorporating estimated k values resulted in improved predictions of breakthrough observed in nonequilibrium scenarios (R2 0.83-0.98), but k values will require further validation prior to broader application. This work illustrates rate-limited sorption considerations are needed to describe 1-D column saturated transport for some PFAAs and solid phases. At field scales, subsurface heterogeneity and PFAA precursor transformation may be equally or even more important in determining saturated PFAA transport, but kinetic parameters in this study may help to determine relative contributions of rate-limited sorption to overall transport.