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.

An overlooked mechanism underlying the attenuated temperature response of soil heterotrophic respiration.

Biogeochemical reactions occurring in soil pore space underpin gaseous emissions measured at macroscopic scales but are difficult to quantify due to their complexity and heterogeneity. We develop a volumetric-average method to calculate aerobic respiration rates analytically from soil with microscopic soil structure represented explicitly. Soil water content in the model is the result of the volumetric-average of the microscopic processes, and it is nonlinearly coupled with temperature and other factors. Since many biogeochemical reactions are driven by oxygen (O2) which must overcome various resistances before reaching reactive microsites from the atmosphere, the volumetric-average results in negative feedback between temperature and soil respiration, with the magnitude of the feedback increasing with soil water content and substrate quality. Comparisons with various experiments show the model reproduces the variation of carbon dioxide emission from soils under different water content and temperature gradients, indicating that it captures the key microscopic processes underpinning soil respiration. We show that alongside thermal microbial adaptation, substrate heterogeneity and microbial turnover and carbon use efficiency, O2 dissolution and diffusion in water associated with soil pore space is another key explanation for the attenuated temperature response of soil respiration and should be considered in developing soil organic carbon models.

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.

Role of co-occurring competition and facilitation in plant spacing hydrodynamics in water-limited environments.

Plant performance (i.e., fecundity, growth, survival) depends on an individual's access to space and resources. At the community level, plant performance is reflected in observable vegetation patterning (i.e., spacing distance, density) often controlled by limiting resources. Resource availability is, in turn, strongly dependent on plant patterning mediated by competitive and facilitative plant-plant interactions. Co-occurring competition and facilitation has never been specifically investigated from a hydrodynamic perspective. To address this knowledge gap, and to overcome limitations of field studies, three intermediate-scale laboratory experiments were conducted using a climate-controlled wind tunnel-porous media test facility to simulate the soil-plant-atmosphere continuum. The spacing between two synthetic plants, a design consideration introduced by the authors in a recent publication, was varied between experiments; edaphic and mean atmospheric conditions were held constant. The strength of the above- and belowground plant-plant interactions changed with spacing distance, allowing the creation of a hydrodynamic conceptual model based on established ecological theories. Greatest soil water loss was observed for the experiment with the smallest spacing where competition dominated. Facilitation dominated at the intermediate spacing; little to no interactions were observed for the largest plant spacing. Results suggest that there exists an optimal spacing distance range that lowers plant environmental stress, thus improving plant performance through reduced atmospheric demand and conservation of available soil water. These findings may provide a foundation for improving our understanding of many climatological, ecohydrological, and hydrological problems pertaining to the hydrodynamics of water-limited environments where plant-plant interactions and community self-organization are important.

Effect of NAPL Source Morphology on Mass Transfer in the Vadose Zone.

The generation of vapor-phase contaminant plumes within the vadose zone is of interest for contaminated site management. Therefore, it is important to understand vapor sources such as non-aqueous-phase liquids (NAPLs) and processes that govern their volatilization. The distribution of NAPL, gas, and water phases within a source zone is expected to influence the rate of volatilization. However, the effect of this distribution morphology on volatilization has not been thoroughly quantified. Because field quantification of NAPL volatilization is often infeasible, a controlled laboratory experiment was conducted in a two-dimensional tank (28 cm × 15.5 cm × 2.5 cm) with water-wet sandy media and an emplaced trichloroethylene (TCE) source. The source was emplaced in two configurations to represent morphologies encountered in field settings: (1) NAPL pools directly exposed to the air phase and (2) NAPLs trapped in water-saturated zones that were occluded from the air phase. Airflow was passed through the tank and effluent concentrations of TCE were quantified. Models were used to analyze results, which indicated that mass transfer from directly exposed NAPL was fast and controlled by advective-dispersive-diffusive transport in the gas phase. However, sources occluded by pore water showed strong rate limitations and slower effective mass transfer. This difference is explained by diffusional resistance within the aqueous phase. Results demonstrate that vapor generation rates from a NAPL source will be influenced by the soil water content distribution within the source. The implications of the NAPL morphology on volatilization in the context of a dynamic water table or climate are discussed.

Long-term mass transfer and mixing-controlled reactions of a DNAPL plume from persistent residuals.

Understanding and being able to predict the long-term behavior of DNAPL (i.e., PCE and TCE) residuals after active remediation has ceased have become increasingly important as attention at many sites turns from aggressive remediation to monitored natural attenuation and long-term stewardship. However, plume behavior due to mass loading and reactions during these later phases is less studied as they involve large spatial and temporal scales. We apply both theoretical analysis and pore-scale simulations to investigate mass transfer from DNAPL residuals and subsequent reactions within the generated plume, and, in particular, to show the differences between early- and late-time behaviors of the plume. In the zone of entry of the DNAPL entrapment zone where the concentration boundary layer in the flowing groundwater has not fully developed, the pore-scale simulations confirm the past findings based on laboratory studies that the mass transfer increases as a power-law function of the Peclét number, and is enhanced due to reactions in the plume. Away from the entry zone and further down gradient, the long-term reactions are limited by the available additive and mixing in the porous medium, thereby behave considerably differently from the entry zone. For the reaction between the contaminant and an additive with intrinsic second-order bimolecular kinetics, the late-time reaction demonstrates a first-order decay macroscopically with respect to the mass of the limiting additive, not with respect to that of the contaminant. The late-time decay rate only depends on the intrinsic reaction rate and the solubility of the entrapped DNAPL. At the intermediate time, the additive decays exponentially with the square of time (t(2)), instead of time (t). Moreover, the intermediate decay rate also depends on the initial conditions, the spatial distribution of DNAPL residuals, and the effective dispersion coefficient.

Experimental study of the effects of DNAPL distribution on mass rebound.

The release of stored dissolved contaminants from low permeability zones contributes to plume persistence beyond the time when dense nonaqueous phase liquid (DNAPL) has completely dissolved. This is fundamental to successfully meeting acceptable low concentrations in groundwater that are driven by site-specific cleanup goals. The study goals were to assess the role of DNAPL entrapment morphology on mass storage and plume longevity. As controlled field studies are not feasible, two-dimensional (2D) test tanks were used to quantify the significance of mass loading processes from source dissolution and stored mass rebound. A simple two-layer soil domain representing a high permeable formation sand overlying a zone of lower permeability sand was used in the tests. DNAPL mass depletion through dissolution was monitored via X-ray photon attenuation, and effluent samples were used to monitor the plume. These data enabled analysis of the DNAPL distribution, the dissolved plume, and the dissolved phase distribution within the low permeability layer. Tests in an intermediate tank showed that mass storage contributes substantially to plume longevity. Detectable effluent concentrations persisted long after DNAPL depletion. The small tank results indicated that the DNAPL morphology influenced the flow field and caused distinctive transport mechanisms contributing to mass storage. Zones of high DNAPL saturation at the interface between the low and high permeability layers exhibited flow bypassing and diffusion dominated transport into the low permeability layer. In the absence of a highly saturated DNAPL zone near the soil interface the contaminant penetrated deeper into the low permeability layer caused by a combination of advection and diffusion.

Biologically enhanced mass transfer of tetrachloroethene from DNAPL in source zones: experimental evaluation and influence of pool morphology.

High-saturation pools of dense nonaqueous phase liquid (DNAPL) are long-term sources of groundwater contamination at many hazardous-waste sites. DNAPL pools consist of a high saturation zone with slow dissolution overlaid by a transition zone with lower saturations and more rapid dissolution. Effects of biological activity on pool dissolution must be understood to evaluate and implement bioremediation strategies. Bioenhanced dissolution of tetrachloroethene (PCE) in transition zones of high-saturation pools was investigated in a custom-designed 5-cm flow cell. Experiments were conducted to characterize mass transfer following DNAPL emplacement, with and without an active microbial culture capable of reductive dehalogenation. For average pool saturations < or = 0.55, mass transfer during biodegradation was enhanced by factors of 4-13, due primarily to high mass flux of PCE degradation products. However, at an average pool saturation of 0.74, mass transfer was enhanced by factors less than 1.5. Mass transfer was significantly greater from pools with an observable transition zone than without. Advective flow through multiphase transition zones enhanced dissolution and biological activity. These laboratory-scale experimental results suggest that biotechnologies may be effective remediation strategies for depletion of source zones within pool transition zones.

Effects of heterogeneity and experimental scale on the biodegradation of diesel.

Biodegradation of petroleum hydrocarbon contamination is a common method for remediating soils and groundwater. Due to complexities with field-scale studies, biodegradation rates are typically evaluated at the bench-scale in laboratory studies. However, important field conditions can be difficult to mimic in the laboratory. This study investigates three scaling factors that can impact laboratory biodegradation rates and that are frequently unaccounted for in typical laboratory experimental procedures. These factors are soil heterogeneity, morphology of petroleum hydrocarbon non-aqueous phase liquids (NAPLs) and soil moisture distribution. The effects of these factors on the biodegradation rate of diesel NAPL is tested under a variety of experimental procedures from well-mixed batch studies to four-foot static soil columns. The results indicate that a high degree of variability results from even small-scale heterogeneities. In addition, it appears that as the experimental scale increases, the measured biodegradation rates slow. The results indicate that diesel biodegradation rates derived from small-scale experiments are not necessarily representative of field-scale biodegradation rates.

The effect of entrapped nonaqueous phase liquids on tracer transport in heterogeneous porous media: laboratory experiments at the intermediate scale.

This work considers the applicability of conservative tracers for detecting high-saturation nonaqueous-phase liquid (NAPL) entrapment in heterogeneous systems. For this purpose, a series of experiments and simulations was performed using a two-dimensional heterogeneous system (10x1.2 m), which represents an intermediate scale between laboratory and field scales. Tracer tests performed prior to injecting the NAPL provide the baseline response of the heterogeneous porous medium. Two NAPL spill experiments were performed and the entrapped-NAPL saturation distribution measured in detail using a gamma-ray attenuation system. Tracer tests following each of the NAPL spills produced breakthrough curves (BTCs) reflecting the impact of entrapped NAPL on conservative transport. To evaluate significance, the impact of NAPL entrapment on the conservative-tracer breakthrough curves was compared to simulated breakthrough curve variability for different realizations of the heterogeneous distribution. Analysis of the results reveals that the NAPL entrapment has a significant impact on the temporal moments of conservative-tracer breakthrough curves.

Convergent-flow tracer tests in heterogeneous media: combined experimental-numerical analysis for determination of equivalent transport parameters.

In modeling transport within naturally heterogeneous aquifers, it is usually assumed that the transport equations valid at local scales can also be applied at larger scales. At larger scales, the heterogeneous domain is represented by an equivalent homogeneous medium. Convergent-flow tracer tests constitute one of the most frequently used field tests to estimate effective input parameters of equivalent homogeneous aquifers. Traditionally, statistical approaches applied to groundwater flow and solute transport have provided tools to estimate these equivalent parameters. These approaches are based on a number of simplifications including the assumption that the point transmissivity values follow a multilog-normal random function. Several investigators have found that this assumption may not be valid in many field cases. In order to study the applicability of the equivalent homogeneous formulation in a nontraditional stochastic field, a number of experimental and numerical studies were conducted. The results are used to determine the apparent values of porosity and dispersivity that would be obtained if convergent-flow tracer tests were conducted in a deterministically generated heterogeneous transmissivity field displaying anisotropy in the correlation structure. It is shown that in this particular heterogeneous media, apparent porosity strongly depends on connectivity rather than on transmissivity. This dependence on connectivity questions the theoretical results obtained in continuum equivalent fields to estimate effective porosity.