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.

Experimental upscaling analyses for a surfactant-enhanced in-situ chemical oxidation (S-ISCO) remediation design.

Surfactant-enhanced in-situ chemical oxidation (S-ISCO) is an emerging innovative remediation technology for the treatment of dense non-aqueous phase liquids (DNAPLs). S-ISCO combines the solubilization of contaminants by means of surfactants with the chemical oxidation by an oxidizing agent, thus, potentially increasing the efficiency of the state-of-the-art ISCO technique. Scientific investigations are needed to enable the technology transfer for potential field applications based on the development of a remediation design under well-defined boundary conditions. For this purpose, experimental upscaling analyses were performed using the special infrastructure of the research facility for subsurface remediation (VEGAS). Batch tests showed that oxidation of the selected surfactant E-Mulse 3® (EM3) by activated persulfate (Na-PS) reduced the solubilization of the model contaminants 1,4-DCB, naphthalene, and PCE. As a consequence, the processes of contaminant solubilization and degradation were temporally and spatially separated in the developed remediation design. A proof of concept was provided by performing an S-ISCO medium-scale experiment (100 cm length, 70 cm height, 12.5 cm width), with 1,2-DCB as model DNAPL contaminant to be treated. A groundwater circulation well (GCW) was used to inject a 60 g/L Na-PS solution and to effectively mix the reagents. Sampling of the experiment's outflow and the soil material after treatment showed that neither rebound effects nor residual mass loadings on the soil material could be detected after termination of the S-ISCO treatment. To further evaluate the S-ISCO remediation design under field-like conditions, a large-scale S-ISCO experiment was conducted (6 m length, 3 m height, 1 m width), allowing for an extensive sampling campaign to monitor relevant processes. An efficient contaminant removal from the former source zone could be reached by surfactant solubilization, decreasing contaminant levels from initially over 2000 mg/L 1,2-DCB to final concentrations below 5 mg/L 1,2-DCB. The heterogeneously distributed contaminant degradation, implemented by a three-filter GCW, was attributed to density-induced migration processes that impeded an optimal reaction zone. A density-dependent numerical transport could qualitatively match the observations. By comparing different simulation scenarios, an adapted operation of the GCW was established that provides for a more efficient distribution of the density-influenced oxidant injection.

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 ∑18PFAAs) and polyfluorinated precursors (1.4 mg/kg ∑18precursors) 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.

An Open, Object-Based Framework for Generating Anisotropy in Sedimentary Subsurface Models.

The spatial distribution of hydraulic properties in the subsurface controls groundwater flow and solute transport. However, many approaches to modeling these distributions do not produce geologically realistic results and/or do not model the anisotropy of hydraulic conductivity caused by bedding structures in sedimentary deposits. We have developed a flexible object-based package for simulating hydraulic properties in the subsurface-the Hydrogeological Virtual Realities (HyVR) simulation package. This implements a hierarchical modeling framework that takes into account geological rules about stratigraphic bounding surfaces and the geometry of specific sedimentary structures to generate realistic aquifer models, including full hydraulic-conductivity tensors. The HyVR simulation package can create outputs suitable for standard groundwater modeling tools (e.g., MODFLOW), is written in Python, an open-source programming language, and is openly available at an online repository. This paper presents an overview of the underlying modeling principles and computational methods, as well as an example simulation based on the Macrodispersion Experiment site in Columbus, Mississippi. Our simulation package can currently simulate porous media that mimic geological conceptual models in fluvial depositional environments, and that include fine-scale heterogeneity in distributed hydraulic parameter fields. The simulation results allow qualitative geological conceptual models to be converted into digital subsurface models that can be used in quantitative numerical flow-and-transport simulations, with the aim of improving our understanding of the influence of geological realism on groundwater flow and solute transport.

Estimating a Representative Value and Proportion of True Zeros for Censored Analytical Data with Applications to Contaminated Site Assessment.

This paper demonstrates a maximum likelihood (ML)-based approach to derive representative ("best guess") contaminant concentrations from data with censored values (e.g., less than the detection limit). The method represents an advancement over existing techniques because it is capable of estimating the proportion of measurements that are true zeros and incorporating varying levels of censorship (e.g., sample specific detection limits, changes through time in method detection). The ability of the method to estimate the proportion of true zeros is validated using precipitation data. The stability and flexibility of the method are demonstrated with stochastic simulation, a sensitivity analysis, and unbiasedness analysis including varying numbers of significant digits. A key aspect of this paper is the application of the statistical analysis to real site rock core contaminant concentration data collected within a plume at two locations using high resolution depth-discrete sampling. Comparison of the representative values for concentrations at each location along the plume center-line shows a larger number of true zeros and generally lower concentrations at the downgradient location according to the conceptual site model, leading to improved estimates of attenuation with distance and/or time and associated confidence; this is not achievable using deterministic methods. The practical relevance of the proposed method is that it provides an improved basis for evaluating change (spatial, temporal, or both) in environmental systems.