Adsorption of PFAAs in the Vadose Zone and Implications for Long-Term Groundwater Contamination

The effect of heavy precipitation on the leaching of heavy metals from tropical coastal legacy tailings

The continued growth in demand for mineral resources has led to a large amount of mining wastes, which is a major challenge in the context of carbon neutrality and climate change. In this study, runoff migration, batch leaching, and column experiments were used to investigate the short-, medium-, and long-term leaching of heavy metals from legacy tailings, respectively; the cumulative metal release kinetic equations were established, and the long-term effects of tailings leaching were verified by HYDRUS-1D. In runoff migration experiments, surface dissolution of tailings and the co-migration of adsorbed soil particles by erosion were the main carriers in the early stages of leachate formation (Mn ∼ 65 mg/L and SO42- up to 2697.2 mg/L). Batch leaching tests showed that the concentration of heavy metals in soil leached by acid rain were 0.1 ∼ 22.0 μg/L for Cr, 0.7 ∼ 26.0 μg/L for Cu, 4.8 ∼ 5646.0 μg/L for Mn, 0.3 ∼ 232.4 μg/L for Ni, and 1.3 ∼ 448.0 μg/L for Zn. The results of column experiments indicated that some soluble components and metals with high mobility showed a significant decreasing trend at cumulative L/S ≤ 2. Additionally, the metals have higher leaching rates under TCLP conditions, as shown by Mn > Co > Zn > Cd > Ni > Cu > Pb > Cr. The fitting results of Langmuir equation were closer to the cumulative release of metals in the real case, and the release amounts of Mn, Zn, Co, and Ni were higher with 55, 5.84, 2.66, and 2.51 mg/kg, respectively. The water flow within tailings affects the spatial distribution of metals, which mainly exist in relatively stable chemical fractions (F3 + F4 + F5 > 90 %) after leaching. Numerical simulation verified that Mn in leachate has reached 8 mg/L at a scale of up to 100 years. The research results are expected to provide technical basis for realizing the resource utilization of tailings in the future.

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.

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.

Determining air-water interfacial areas for the retention and transport of PFAS and other interfacially active solutes in unsaturated porous media

The objective of this work was to determine the methods that produce the most representative measurements and estimations of air-water interfacial area specifically for the retention and transport of PFAS and other interfacially active solutes in unsaturated porous media. Published data sets of air-water interfacial areas obtained with multiple measurement and prediction methods were compared for paired sets of porous media comprising similar median grain diameters, but one with solid-surface roughness (sand) and one without roughness (glass beads). All interfacial areas produced with the multiple diverse methods were coincident for the glass beads, providing validation of the aqueous interfacial tracer-test methods. The results of this and other benchmarking analyses demonstrated that the differences in interfacial areas measured for sands and soil by different methods are not due to errors or artifacts in the methods but rather the result of method-dependent differential contributions of solid-surface roughness. The contributions of roughness to interfacial areas measured by interfacial tracer-test methods were quantified and shown to be consistent with prior theoretical and experiment-based investigations of air-water interface configurations on rough solid surfaces. Three new methods for estimating air-water interfacial areas were developed, one based on the scaling of thermodynamic-determined values and the other two comprising empirical correlations incorporating grain diameter or NBET solid surface area. All three were developed based on measured aqueous interfacial tracer-test data. The three new and three existing estimation methods was tested using independent data sets of PFAS retention and transport. The results showed that the method based on treating air-water interfaces as smooth surfaces as well as the standard thermodynamic method produced inaccurate air-water interfacial areas that failed to reproduce the multiple measured PFAS retention and transport data sets. In contrast, the new estimation methods produced interfacial areas that accurately represented air-water interfacial adsorption of PFAS and associated retention and transport. The measurement and estimation of air-water interfacial areas for field-scale applications is discussed in light of these results.

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.