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

Numerical modeling of PFAS movement through the vadose zone: Influence of plant water uptake and soil organic carbon distribution

In this study, we investigated the effects of soil organic carbon (SOC) distribution and water uptake by plant roots on PFAS movement in the vadose zone with a deep groundwater table under temperate, humid climate conditions. Two series of numerical simulations were performed with the HYDRUS computer code, representing the leaching of historical PFOS contamination and the infiltration of water contaminated with PFOA, respectively. We considered soil profiles with three distributions of SOC (no SOC, realistic SOC distribution decreasing with depth, and uniform SOC equal to the content measured in topsoil), three root distributions (bare soil, grassland, and forest), and three soil textures (sand, sandy loam, and loam). The SOC distribution had a profound impact on the velocity of PFOS movement. The apparent retardation factor for realistic SOC distribution was twice as large as for the scenario with no SOC and more than three times smaller than for the scenario with uniformly high SOC content. We also showed that the root distribution in soil profoundly impacts the simulations of PFAS migration through soil. Including the root zone significantly slows down the movement of PFAS, primarily due to increased evapotranspiration and reduced downward water flux. Another effect of water uptake by plant roots is an increase of PFAS concentrations in soil water (evapo-concentration). The evapo-concentration and the slowdown of PFAS movement due to root water uptake are more significant in fine-textured soils than in sand.

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.

PFAS release from the subsurface and capillary fringe during managed aquifer recharge

Managed aquifer recharge (MAR) is a sustainable way of harvesting groundwater in water-stressed urbanized areas, where reclaimed wastewater or stormwater is applied on a large basin to infiltrate water into the groundwater aquifer naturally. This process could rapidly fluctuate the water table and move the capillary fringe boundary, and the change in flow dynamic and associated geochemical changes could trigger the release of sequestered pollutants, including per- and polyfluoroalkyl substances (PFAS), also known as 'forever chemicals', from the subsurface and capillary fringe. Yet, the potential of PFAS release from the subsurface and capillary zone during recharge events when the water table rapidly fluctuates has not been evaluated. This study uses laboratory column experiments to simulate PFAS release from pre-contaminated subsurface and capillary fringe during groundwater table fluctuation. The results reveal that the groundwater level fluctuations during MAR increased the release of perfluorobutanesulfonic acid (PFBS) and perfluorooctanesulfonic acid (PFOS) from the capillary fringe, but the fraction released depended on PFAS type and their association with soil colloids. A higher proportion of PFOS in column effluent was found to be associated with particles, while a greater portion of released PFBS was in a free or dissolved state. The direction of water table fluctuation did not affect the release of PFAS in this study. A lack of change in the concentration of bromide, a conservative tracer, during flow interruption, indicates that diffusion of PFAS through reconnected pores during water table rise had an insignificant effect on PFAS release. Overall, this study provides insights into how PFAS can be released from the subsurface and capillary fringe during managed aquifer recharge when the groundwater level is expected to fluctuate quickly.

Per- and polyfluoroalkyl substances (PFAS) in final treated solids (Biosolids) from 190 Michigan wastewater treatment plants

Trends in concentration, distribution, and variability of per- and polyfluoroalkyl substances (PFAS) in biosolids are characterized using an extensive dataset of 350 samples from 190 wastewater treatment plants (WWTPs) across Michigan. All samples are comprised of final treated solids generated at the end of the wastewater treatment process. Concentrations of both individual and Σ24 PFAS are lognormally distributed, with Σ24 PFAS concentrations ranging from 1-3200 ng/g and averaging 108 ± 277 ng/g dry wt. PFAS with carboxyl and sulfonic functional groups comprise 29% and 71% of Σ24 PFAS concentrations, respectively, on average. Primary sample variability in concentration is associated with long-chain PFAS with higher tendency for partitioning to biosolids. Short-chain carboxylic compounds, most notably PFHxA, are responsible for secondary concentration variability. Usage of FTSA and PFBS replacements to long-chain sulfonic compounds also contributes to variance in biosolids concentrations. Sulfonamide precursor compounds as a collective group are detected at a similar frequency as PFOS and often have higher concentrations. Trends in PFAS enrichment for individual PFAS vary at least 3 orders-of-magnitude and generally increase with compound hydrophobicity; however, partitioning of PFAS onto solids in WWTPs is a complex process not easily described nor constrained using experimentally-derived partitioning coefficients.

Predicting the impact and duration of persistent and mobile organic compounds in groundwater systems using a contaminant mass discharge approach

This study investigated methods for predicting the duration and impact on groundwater quality from persistent and mobile organic compounds (PMOCs) at a drinking water well field affected by multiple contaminant sources. The fungicide metabolite N,N-dimethylsulfamide (DMS), which frequently occurs above the Danish groundwater quality criterion (0.1 μg/L), was used as an example. By combining contaminant mass discharge (CMD) estimations, modeling, and groundwater dating, a number of important discoveries were made. The current center of contaminant mass was located near the source area. The CMD at the well field was predicted to peak in 2040, and an effect from the investigated sources on groundwater quality could be expected until the end of the 21st century. A discrepancy in the current CMD at the well field and the estimated arrival time from the studied source area suggested an additional pesticide source, which has not yet been thoroughly investigated. The presence of the unknown source was supported by model simulations, producing an improved mass balance after inclusion of a contaminant source closer to the well field. The approach applied here was capable of predicting the duration and impact of DMS contamination at a well field at catchment scale. It furthermore shows potential for identification and quantification of the contribution from individual sources, and is also applicable for other PMOCs. Predicting the duration of the release and impact of contaminant sources on abstraction wells is highly valuable for water resources management and authorities responsible for contaminant risk assessment, remediation, and long-term planning at water utilities.

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.

The occurrence, spatial distribution, and well-depth dependence of PFASs in groundwater from a reclaimed water irrigation area

With the growing development of modern agriculture and industry, groundwater is facing more and more complex contaminants. One such contaminant is per- and polyfluoroalkyl substances (PFASs), which pose a potential risk to human health, particularly for those who rely on groundwater as their primary source of drinking water. In this study, we conducted a comprehensive investigation on the occurrence, spatial distribution, and source apportionment of PFASs in shallow (<60 m) and deep (>80 m) groundwater samples from a reclaimed water irrigation area in Beijing's suburbs. Our results showed that the average total PFAS concentration (∑10PFAS) for all samples was 10.55 ± 7.77 ng/L, ranging from 1.05 to 34.28 ng/L. The dominant congeners were PFBA, PFOA, and PFBS. No significant linear relationship was observed between PFAS concentrations and the well depth. However, the averaged ΣPFASs in groundwater were highest in the uppermost layer and declined sharply to a few ng/L in the deep aquifer below 80 m. PFASs showed elevated concentration in shallow aquifers in 9 out of 11 paired wells, indicating an overall descending trend of PFASs with increasing aquifer depth. The spatial distribution of PFASs was highly heterogeneous and showed different patterns in shallow and deep groundwater, which may be related to the complicated attenuation behavior of PFAS compounds when they transport and diffuse through overlapping aquifer layers. The influence of the landfill on groundwater PFASs was most pronounced within a 5 km radius. Source apportionment results indicated that reclaimed water irrigation is the main non-point source of PFASs in shallow groundwater. In contrast, deep groundwater is primarily subject to point sources and lateral recharge flow. This investigation of PFASs in shallow and deep wells provides a foundation for further exploration of PFASs transportation and risk prevention in regions where groundwater is a major water resource for domestic and industrial development.

Immobilization of per- and polyfluoroalkyl substances (PFAS): Comparison of leaching behavior by three different leaching tests

The evaluation of PFAS immobilization performance in laboratory experiments, especially the long-term stability, is a challenge. To contribute to the development of adequate experimental procedures, the impact of experimental conditions on the leaching behavior was studied. Three experiments on different scales were compared: batch, saturated column, and variably saturated laboratory lysimeter experiments. The Infinite Sink (IS) test - a batch test with repeated sampling - was applied for PFAS for the first time. Soil from an agricultural field amended with paper-fiber biosolids polluted with various perfluoroalkyl acids (PFAAs; 655 μg/kg ∑¹⁸PFAAs) and polyfluorinated precursors (1.4 mg/kg ∑¹⁸precursors) was used as the primary material (N-1). Two types of PFAS immobilization agents were tested: treatment with activated carbon-based additives (soil mixtures: R-1 and R-2), and solidification with cement and bentonite (R-3). In all experiments, a chain-length dependent immobilization efficacy is observed. In R-3, the leaching of short-chain PFAAs was enhanced relative to N-1. In column and lysimeter experiments with R-1 and R-2, delayed breakthrough of short-chain PFAAs (C4) occurred (> 90 days; in column experiments at liquid-to-solid ratio (LS) > 30 L/kg) with similar temporal leaching rates suggesting that leaching in these cases was a kinetically controlled process. Observed differences between column and lysimeter experiments may be attributed to varying saturation conditions. In IS experiments, PFAS desorption from N-1, R-1, and R-2 is higher than in the column experiments (N-1: +44 %; R-1: +280 %; R-2: +162 %), desorption of short-chain PFAS occurred predominantly in the initial phase (< 14 days). Our findings demonstrate that sufficient operating times are essential in percolation experiments, e.g., in column experiments >100 days and LS > 30 L/kg. IS experiments may provide a faster estimate for nonpermanent immobilization. The comparison of experimental data from various experiments is beneficial to evaluate PFAS immobilization and to interpret leaching characteristics.