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

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

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 (Knw < 0.1) and adsorption on 60-80 mesh Ottawa sand (Kd < 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.

A similarity solution for reaction front propagation in a fracture-matrix system.

We propose a new composite similarity variable, based on which a similarity solution is derived for reaction front propagation in fracture-matrix systems. The similarity solution neglects diffusion/dispersion within the fracture and assumes the existence of a sharp reaction front in the rock matrix. The reaction front location in the rock matrix is shown to follow a linear decrease with distance along the fracture. The reaction front propagation along the fracture is shown to scale like diffusion (i.e. as the square root of time). The similarity solution using the composite similarity variable appears to be applicable to a broad class of reactive transport problems involving mineral reactions in fracture-matrix systems. It also reproduces the solutions for non-reactive solute and heat transport when diffusion/dispersion/conduction are neglected in the fracture. We compared our similarity solution against numerical simulations for nonlinear reactive transport of an aqueous species with a mineral in the rock matrix. The similarity solutions agree very well with the numerical solutions, especially at later times when diffusion limitations are more pronounced.This article is part of the themed issue 'Energy and the subsurface'.