Denitrification dynamics in unsaturated soils with different porous structures and water saturation degrees: A focus on the shift in microbial community structures

Prediction of arsenic retention in vadose zone based on empirical relationship between soil properties and segmented retardation factors

Arsenic (As) is a widespread environmental contaminant that poses a significant threat to ecosystems and human health. Although previous studies have qualitatively revealed the effects of individual soil properties on the transport and fate of As in the vadose zone, their integrated impacts remain obscure. Moreover, studies investigating the retardation factor therein, which is a key parameter for comprehending As transport in the vadose zone, are extremely limited. In this study, we investigated the interplay of soil properties with As transport and retention within the vadose zone, while focusing on the retardation factor of As. We employed steady-state unsaturated water-flow soil column experiments coupled with a mobile-immobile model and multiple linear regression analysis to elucidate the dependence of As retardation factors on the soil properties. In the mobile water zone, iron and organic matter contents emerged as the two most influential properties that impedes As mobility. Whereas, in the immobile water zone, the coefficient of uniformity and bulk density were the most influential factors that enhanced As retention. Finally, we derived an empirical equation for calculating the As retardation factors in each zone, offering a valuable tool for describing and predicting As behavior to protect the groundwater resources underneath.

Cover crop residue decomposition triggered soil oxygen depletion and promoted nitrous oxide emissions

Cover cropping is a promising strategy to improve soil health, but it may also trigger greenhouse gas emissions, especially nitrous oxide (N2O). Beyond nitrogen (N) availability, cover crop residue decomposition may accelerate heterotrophic respiration to limit soil O2 availability, hence promote N2O emissions from denitrification under sub-optimal water-filled pore space (WFPS) conditions that are typically not conducive to large N2O production. We conducted a 21-day incubation experiment to examine the effects of contrasting cover crop residue (grass vs legume) decomposition on soil O2 and biogeochemical changes to influence N2O and CO2 emissions from 15N labeled fertilized soils under 50% and 80% WFPS levels. Irrespective of cover crop type, mixing cover crop residue with N fertilizer resulted in high cumulative N2O emissions under both WFPS conditions. In the absence of cover crop residues, the N fertilizer effect of N2O was only realized under 80% WFPS, whereas it was comparable to the control under 50% WFPS. The N2O peaks under 50% WFPS coincided with soil O2 depletion and concomitant high CO2 emissions when cover crop residues were mixed with N fertilizer. While N fertilizer largely contributed to the total N2O emissions from the cover crop treatments, soil organic matter and/or cover crop residue derived N2O had a greater contribution under 50% than 80% WFPS. Our results underscore the importance of N2O emissions from cover crop-based fertilized systems under relatively lower WFPS via a mechanism of respiration-induced anoxia and highlight potential risks of underestimating N2O emissions under sole reliance on WFPS.

Influence of iron (hydr)oxide mineralogy and contents in aquifer sediments on dissolved organic carbon attenuations during aquifer storage and recovery

Aquifer storage and recovery (ASR) is a promising approach for managing water resources that enhances water quality through biogeochemical reactions occurring within aquifers. Iron (hydr)oxides, which are the predominant metallic oxides in soil, play a crucial role in degrading dissolved organic carbon (DOC), primarily through a process known as dissimilatory iron reduction (DIR). However, the efficiency of this reaction varies depending on the mineralogy and composition of the aquifer, and this understanding is essential for adequate water quality in ASR. The objective of this study is to investigate the impact of iron (hydr)oxide on acetate, as an organic carbon source, attenuation during the ASR. To achieve this, three sets of laboratory sediment columns were prepared, each containing a different type of iron (hydr)oxide minerals: ferrihydrite, goethite, and hematite. Following an acclimation period of 28 days to simulate the microcosm within an aquifer, the columns were continuously supplied with the simulated river water spiked with acetate (DOC 40-60 mg L-1), and the acetate concentration in the effluent was monitored. The result revealed that the column containing ferrihydrite achieved 97% acetate attenuation through DIR with anoxic conditions (DO < 0.1 mg L-1), while the goethite and hematite columns exhibited limited attenuation rates of 40 and 50%, respectively. Furthermore, the efficiency of acetate attenuation in the ferrihydrite columns increased with the content of ferrihydrite but experienced a rapidly declined at higher contents (3-4%), possibly due to the partial conversion of ferrihydrite to goethite as a result of the interaction between ferrihydrite and the Fe(II) produced during DIR. Additionally, an analysis of the microbial community demonstrated that microorganisms known to possess the ability to reduce iron (hydr)oxides under anaerobic conditions were abundant in the ferrihydrite columns.

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.

Stable isotope probing reveals compositional and functional shifts in active denitrifying communities along the soil profile in an intensive agricultural area

Denitrifying microbial communities in the vadose zone play an essential role in eliminating the nitrate leached from agricultural practices. This nitrate could otherwise contaminate groundwater and threaten public health. Here, we utilized stable isotope probing combined with amplicon sequencing and functional gene quantification to inspect the composition and function of heterotrophic denitrifying microorganisms along a 9-m soil profile in an intensive agricultural area. Dramatic differences in the composition of the active denitrifiers were uncovered between the surface soil and deep layers of the vadose zone. The main denitrifying bacterial taxa identified from 13C-DNA fractions were Pseudomonadaceae (Pseudomonas), Rhodocyclaceae (Azoarcus), and Burkholderiaceae in the surface soil (0-0.2 m), and were Pseudomonadaceae (Pseudomonas), Burkholderiaceae, Bacillaceae (Bacillus), and Paenibacillaceae (Ammoniphilus) in the deep layers (0.5-9.0 m). Analysis of the functional genes (nirS, nirK, and nosZ) in isotope-labeled DNA revealed an upward nos/nir ratio with increasing soil depth, which may account for the higher nitrous oxide emission potential in the surface soil, as compared to the deeper sand-rich, low organic carbon layers. This study improves our understanding of active denitrifying microbes in the vadose zone and helps in developing techniques to reduce nitrate pollution in groundwater.

Different effects of taproot and fibrous root crops on pore structure and microbial network in reclaimed soil

Understanding the effects of plant roots on the pore structure and microbial community of soil is crucial to recovery and improve soil productivity in mining areas. This study aims to assess the impact of taproot (TR) and fibrous root (FR) crops on the physicochemical properties, pore structure, and microbial communities and networks in reclaimed mine soil. Results showed that reclamation positively influenced pore structure and microbial diversity. Tillage with TR and FR crops significantly increased porosity, total pore volume, and area of mining soil (p < 0.05). Compared with TR, FR produced more macropores, mesopores, and micropores. In addition, the module group, average degree, density, and connectivity of microbial network in FR cultivated soil were higher than those in TR cultivated soil. The microbial network map showed that FR had more keystone taxa than TR, and mainly consisted of Acidobacteria and Proteobacteria. In the FR microbial network, Rhizobiales, Betaproteobacteria, and Acidobacteria_Gp11 play critical roles as module hubs and Noviherbaspirillum and Zavarzinella as connectors. Furthermore, most of the key microbes were significantly correlated (p < 0.05) with the total pore area and probably tended to live in pores >75 μm and 0.1-5 μm in size. Therefore, FR crops were more effective than TR crops in improving pore structure and enhancing the development of microbial network in reclaimed soil.

Soil pore structure mediates the effects of soil oxygen on the dynamics of greenhouse gases during wetting-drying phases

The timing and magnitude of greenhouse gas (GHG) production depend strongly on soil oxygen (O2) availability, and the soil pore geometry characteristics largely regulate O2 and moisture conditions relating to GHG biochemical processes. However, the interactions between O2 dynamics and the concentration and flux of GHGs during the soil moisture transitions under various soil pore conditions have not yet been clarified. In this study, a soil-column experiment was conducted under wetting-drying phases using three pore-structure treatments, FINE, MEDIUM, and COARSE, with 0 %, 30 %, and 50 % coarse quartz sand applied to soil, respectively. The concentrations of soil gases (O2, nitrous oxide (N2O), carbon dioxide (CO2), and methane (CH4)) were monitored at a depth of 15 cm hourly, and their surface fluxes were measured daily. Soil porosity, pore size distribution, and pore connectivity were quantified using X-ray computed microtomography. The soil O2 concentrations were found to decline sharply as soil moisture increased to the water holding capacities of 0.46, 0.41, and 0.32 cm cm-3 in the FINE, MEDIUM, and COARSE, respectively. The dynamic patterns of the O2 concentrations varied across the soil pore structures, decreasing to anaerobic in FINE (<0.01 %) and MEDIUM (0.02 %), and to hypoxic (4.42 %) in COARSE. Correspondingly, the soil N2O concentration was the highest in FINE (101 μL L-1) and the lowest in COARSE (10 μL L-1), whereas the highest surface N2O flux was observed in MEDIUM (131 μg N m-2 h-1). As soil CO2 concentrations declined, CO2 fluxes increased from FINE to MEDIUM to COARSE. Most pores of FINE, MEDIUM, and COARSE were 15-80 μm, 85-100 μm, and 105-125 μm, respectively, in terms of diameter. The X-ray CT visible (>15 μm) porosity in FINE, MEDIUM and COARSE were 0.09, 0.17, and 0.28 mm3 mm-3, respectively. The corresponding Euler-Poincaré numbers were 180,280, 76,705, and -10,604, respectively, indicating higher connectivity in COARSE than in MEDIUM or FINE. In soil dominated by small air-filled porosity which limits gas diffusion and result in low soil O2 concentration, N2O concentration was increased and CO2 flux was inhibited as the moisture content increased. The turning point in the sharp decrease in O2 concentration was found to correspond with a moisture content, and a pore diameter of 95-110 μm was associated with the critical turning point between holding water and O2 depletion in soil. These findings suggest that O2-regulated biochemical processes are key to the production and flux of GHGs, which in turn are dependent on the soil pore structure and a coupling relationship between N2O and CO2. Improved understanding of the intense effect of soil physical properties provided an empirical foundation for the future development of mechanistic prediction models for how pore-space scale processes with high temporal (hourly) resolution up to GHGs fluxes at larger spatial and temporal scales.

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

The migration of petroleum hydrocarbons in vadose zone involves complex coupled processes such as downward displacement and natural attenuation. Despite its significance in determining groundwater vulnerability to petroleum contamination and optimizing the remedial strategy, it has not been comprehensively studied in terms of overall processes under field-relevant conditions. In this study, a series of unsaturated soil column experiments were conducted by simulating subsurface diesel contamination within a vadose zone using different soil textures at different soil bulk densities and initial diesel concentrations, while partly exposing them to simulated precipitation. The results showed that the soil column with less fine fraction was favorable for the downward migration of diesel but unfavorable for its natural degradation. However, precipitation complicated the relative conductivities of multiple fluids (water, air, and diesel) through the pore network, therby decreasing diesel migration and degradation. For example, the downward migration of diesel in the SL column decreased by 8.4% under precipitation, while the overall attenuation rate dropped to almost 0.24% of its original state. Lowering bulk density (from 1.5 to 1.23 g/cm3), however, could enhance the attenuation rate presumably due to the secured void space for the incoming fluids. A high initial concentration of diesel (2%; w/w) inhibited its natural attenuation, while its influence on its vertical propagation after the precipitation was not significant. The present findings provide a mechanistic basis for approximating the behavior of petroleum hydrocarbons in a random vadose zone.