Complex behavior of petroleum hydrocarbons in vadose zone: A holistic analysis using unsaturated soil columns

Explainable machine learning for arsenic remobilization potential in the vadose zone: Leveraging readily available soil properties

The vadose zone acts as a natural buffer that prevents contaminants such as arsenic (As) from contaminating groundwater resources. Despite its capability to retain As, our previous studies revealed that a substantial amount of As could be remobilized from soil under repeated wet-dry conditions. Overlooking this might underestimate the potential risk of groundwater contamination. This study quantified the remobilization of As in the vadose zone and developed a prediction model based on soil properties. 22 unsaturated soil columns were used to simulate vadose zones with varying soil properties. Repeated wet-dry cycles were conducted upon the As-retaining soil columns. Consequently, 13.9-150.6 mg/kg of As was remobilized from the columns, which corresponds to 37.0-74.6 % of initially retained As. From the experimental results, a machine learning model using a random forest algorithm was established to predict the potential for As remobilization based on readily accessible soil properties, including organic matter (OM) content, iron (Fe) content, uniformity coefficient, D30, and bulk density. Shapley additive explanation analyses revealed the interrelated effects of multiple soil properties. D30, which is inter-related with Fe content, exhibited the highest contribution to As remobilization, followed by OM content, which was partially mediated by bulk density.

Spatial distribution characteristics and degradation mechanism of microorganisms in n-hexadecane contaminated vadose zone

The damage caused by petroleum hydrocarbon pollution to soil and groundwater environment is becoming increasingly significant. The vadose zone is the only way for petroleum hydrocarbon pollutants to leak from surface into groundwater. The spatial distribution characteristics of indigenous microorganisms in vadose zone, considering presence of capillary zones, have rarely been reported. To explore the spatial distribution characteristics of indigenous microorganisms in vadose zone contaminated by petroleum hydrocarbons, a one-dimensional column migration experiment was conducted using n-hexadecane as characteristic pollutant. Soil samples were collected periodically from different heights during experiment. Corresponding environmental factors were monitored online. The microbial community structure and spatial distribution characteristics of the cumulative relative abundance were systematically analyzed using 16S rRNA sequencing. In addition, the microbial degradation mechanism of n-hexadecane was analyzed using metabolomics. The results showed that presence of capillary zone had a strong retarding effect on n-hexadecane infiltration. Leaked pollutants were mainly concentrated in areas with strong capillary action. Infiltration and displacement of NAPL-phase pollutants were major driving force for change in moisture content (θ) and electric conductivity (EC) in vadose zone. The degradation by microorganisms results in a downward trend in potential of hydrogen (pH) and oxidation-reduction potential (ORP). Five petroleum hydrocarbon-degrading bacterial phyla and 11 degradable straight-chain alkane bacterial genera were detected. Microbial degradation was strong in the area near edge of capillary zone and locations of pollutant accumulation. Mainly Sphingomonas and Nocardioides bacteria were involved in microbial degradation of n-hexadecane. Single-end oxidation involved microbial degradation of n-hexadecane (C16H34). The oxygen consumed, hexadecanoic acid (C16H32O2) produced during this process, and release of hydrogen ions (H+) were the driving factors for reduction of ORP and pH. The vadose zone in this study considered presence of capillary zone, which was more in line with actual contaminated site conditions compared with previous studies. This study systematically elucidated vertical distribution characteristics of petroleum hydrocarbon pollutants and spatiotemporal variation characteristics of indigenous microorganisms in vadose zone considered presence of capillary zone. In addition, the n-hexadecane degradation mechanism was elucidated using metabolomics. This study provides theoretical support for development of natural attenuation remediation measures for petroleum-hydrocarbon-contaminated soil and groundwater.

Groundwater-level fluctuation effects on petroleum hydrocarbons in vadose zones and their potential risks: Laboratory studies

Despite the role of the vadose zone protecting groundwater from contamination, the non-stationarity in this zone makes it difficult to predict the behavior of petroleum hydrocarbons (PH) therein. In laboratory soil columns with sandy and sandy loam soils, we simulated a vadose zone subjected to repeated groundwater-level fluctuation (GLF) to evaluate the behavior of PH under hydrodynamic conditions. The GLF vertically redistributed the PH, the extent of which was pronounced in the sandy soil with a high initial concentration due to the enhanced transport of the immiscible PH through the larger pores. The frequency of GLF did not show a substantial effect on the extent of PH redistribution but largely affected their attenuation. The greater GLF hindered PH volatilization by maintaining a high degree of water saturation, while the subsequent development of a local anaerobic regime inhibited biodegradation, which was more apparent in the sandy loam. Finally, a specific potential risk index was introduced to quantitatively compare the potential risk of PH contamination in different vadose zones exposed to GLF. Overall, the sandy soil contaminated with the higher total PH (TPH) concentration showed markedly higher potential risk indices (i.e., 18.4-29.0%), while the ones comprised of the sandy loam showed 0.6-4.9%, which increased under the greater number of GLF cycles.

Evaluation of the impact of transition from porous to fractured rock media on 3D field-scale DNAPLs contamination

A 3D high-resolution subsurface characteristic (HSC) numerical model to assess migration and distribution of subsurface DNAPLs was developed. Diverse field data, including lithologic, hydrogeologic, petrophysical, and fracture information from both in situ observations and laboratory experiments were utilized for realistic model representation. For the first time, the model integrates hydrogeologic characteristics of both porous (unconsolidated soil (US) and weathered rock (WR)) and fractured rock (FR) media distinctly affecting DNAPLs migration. This allowed for capturing DNAPLs behavior within US, WR, and FR as well as at the boundary between the media, simultaneously. In the 3D HSC model, hypothetical 100-year DNAPLs contamination was simulated, quantitatively analyzing its spatiotemporal distributions by momentum analyses. Twelve sensitivity scenarios examined the impact of WR and FR characteristics on DNAPLs migration, delineating significant roles of WR. DNAPLs primarily resided in WR due to low permeability and limited penetration into FR through sparse inlet fractures. The permeability anisotropy in WR was most influential to determine the DNAPLs fate, surpassing the impacts of FR characteristics, including rock matrix permeability, fracture aperture size, and fracture + rock mean porosity. This study first attempted to apply the field-data-based multiple geological media concept in the DNAPLs prediction model. Consequently, the field-scale effects of WR and media transitions, which have been often overlooked in evaluating DNAPLs contamination, were underscored.