Vadose Research Articles

Responses of deep soil moisture to direct rainfall and groundwater in the red soil critical zone: A four-stage pattern

Deep soil moisture (DSM) plays a critical role in vadose zone hydrological processes. However, the recharge processes and their underlying mechanisms remain unclear due to the scarcity of comprehensive deep hydrological data, especially in regions characterized by active interactions between surface and subsurface hydrological processes. To address this knowledge gap, we performed an event-based study using a four-year hydrological monitoring dataset, encompassing rainfall, profile soil moisture, and groundwater along a hillslope at the Red Soil Critical Zone Observatory (CZO) in China. A total of 226 rainfall events were classified into four distinct types (I-IV) using a k-means clustering algorithm. The changes in DSM and groundwater table (GWT) were dependent on the rainfall type. The GWT (0.57–0.97) played a greater role in contributing to the DSM than the direct rainfall (0.48–0.64) based on a structural equation model, which was more pronounced from upper slope to foot slope. DSM recharge pattern can be described by a four-stage conceptual model, as indicated by the relationship between DSM and GWT: a slow and linear increase (stage 1), a small plateau (stage 2), a rapid linear rise (stage 3), and finally a great plateau (stage 4). DSM at the small and great plateau were closely related with the field capacity and saturation, respectively. Rainfall type, vadose zone features, and soil physical properties jointly shaped DSM recharge patterns. These findings can offer valuable insights into understanding the deep hydrological processes of the red soil region as well as regions with similar hydroclimatic conditions.

Discovery of modern living intertidal stromatolites on Sheybarah Island, Red Sea, Saudi Arabia

Abstract Microbial carbonates, and stromatolites in particular, represent the earliest geological record of life on Earth, which dominated the planet as the sole biotic carbonate factory for almost 3 b.y., from the Archean to the late Proterozoic. Rare and sparsely scattered across the globe in the present day, modern “living” stromatolites are typically relegated to extreme environmental niches, remaining as vestiges of a prodigious microbial past. Here, we report the first discovery of living shallow-marine stromatolites in the Middle East, on Sheybarah Island, Al Wajh carbonate platform, on the NE Red Sea shelf (Saudi Arabia). We detail their regional distribution and describe their environmental conditions, internal structures, and microbial diversity. We also report the first discovery of reticulated filaments in a photic setting, associated with these stromatolites. The Sheybarah stromatolites occur in the intertidal to shallow subtidal zones along the seaward-facing beach in three depth-dependent growth forms. Their inner layers were formed by microbially mediated accretion and differential lithification of sediment grains. Compositional microbial analysis revealed the presence of a wide range of microbial life forms.

Identification of Deep and Shallow Groundwater Potential Zones using GIS Techniques; Case Study: Monaragala District, Sri Lanka

One of the most crucial steps in the management of groundwater sources sustainably is the delineation of the zones with potential for groundwater. However, in addition to taking a lot of time and money, groundwater investigations can be difficult and complex in many locations. To get sufficient results, it needed to demarcate shallow and deep groundwater potentiality separately. Sri Lankan districts have various problems due to a lack of water. With population growth and urbanization, surface water is insufficient to meet the water requirement. This study focused on groundwater potential mapping in Monaragala district covered with the total are 5708 km2 in Sri Lanka. The fuzzy overlay and weighted overlay method integrated with GIS were mainly used to determine potential zones for deep and shallow groundwater. The final results from the two overlaying methods were compared to identify most accurate locations. Six parameters, namely aquifer type, soil depth, soil group, rainfall, proximity to surface water bodies, and slope, were evaluated for the shallow groundwater studies. Proximity to lineament, geology and shallow groundwater parameters were assessed for profound groundwater potentiality. The final output resulting maps were categorized into five main categories such as extremely low, low, moderate, high and extremely high by considering obtained potential values. The deep groundwater potential map demonstrates that the high groundwater potential areas are located around the lineaments, prominently areas having mountain ranges. Therefore the effect of the lineament is very high for deep groundwater availability. High shallow groundwater potential areas are located in smaller regions from Katharagama, Thanamalwila and Siyambalanduwa. Obtained results were validated using existing 220 tube well data and 100 dug well data. 
 Keywords: Fuzzy overlay, Geology, GIS, Groundwater, Weighted overlay

Spatial distribution of the diatoms <i>Dactyliosolen fragilissimus</i> and <i>Cerataulina pelagica</i> in the autumn phytoplankton of the South-Eastern Baltic in 2020–2021

The mass development of diatom Dactyliosolen fragilissimus and Cerataulina pelagica, which are atypical of the South-East Baltic Sea water area, were registered along the Sambia peninsula coastline in the mid of October 2020. Both species were found in the coastal zone of the western and northern coast of the Sambia peninsula to a depth of 32 m, from the Baltic Strait to port Pionerskiy. In early November 2021 a mass vegetation of C. pelagica was recorded –both in the shallow coastal zone of the northern coast of the Sambia peninsula and in the deep water zone of the South-Eastern Baltic Sea. The biomass of this species reached 2.8 g/m3. Both diatom species were abundant within the entire Russian Exclusive Economic Zone. The prevailing hydrological situation during the study period probably facilitated penetration of the species into the study area by surface transport from the Arkona Basin.

Paleozoic reef formation in the Pechora syneclise and the Caspian basin – a comparative analysis

In the two marginal depressions of the East European Platform, the Pechora syneclise and the Caspian depression, reef formations are widely represented, but the stratigraphic intervals of their development are far from equal. In the Pechora syneclise, reef formation began in the second half of the Llandoverian, reached its maximum in the Frasnian and was replaced by the development of reef hills in the Famennian. The structures are represented both as asymmetric reefs framing the shallow zones at their boundary with relatively deep paleobasins and in the form of single structures within the latter. The reefs of the next global maximums, the Late Visean-Serpukhovian and Lower Permian, are developed only marginally at the boundary with the Ural paleoocean and its relic, the Pre-Urals foreland basin. In the Caspian basin, reefs of all three global maxima of development are present, and there are both asymmetric reef systems framing the shelf edges and symmetric intra-basin isolated structures. Such difference is due to the different paleogeomorphological type of the basins. The Caspian basin was sharply differentiated in depth throughout the Middle and Late Paleozoic, which led to the formation of reefs rising above the bottom of these basins. In the Pechora syneclise, depth differentiation of the basins occurred only in the late Devonian. Generally weakly dissected shallow seas of the Visiean-Serpukhovian and early Permian did not provide conditions for the formation of powerful reefs protruding above the seafloor. The latter were formed only on the edge of the Ural paleoocean, and in the Permian ‒ on the edge of the Pre-Ural foreland basin.

Exploring the artistic dimensions of soils in the vadose zone

AbstractThe rise of industrialized agriculture, coupled with a global demographic shift toward urbanization, has marginalized every‐day connections between individuals and soils. This societal shift has led to a decline in the appreciation of and cultural identity with soils. Amid a broader movement aimed at fostering soil awareness and environmental action, many artists and designers have been instrumental in bringing soils to the cultural mainstream, figuring their esthetic, social, political, and ecological dimensions. Contemporary artists place great importance on the underlying idea of animating soils through artistic practice, often transcending the physical artifact itself to assign agency to soils and the myriad beings that occupy them. Artists harness a multitude of media to articulate human‐soil relationships. From the early environmental art movements of the 1960s and 1970s to contemporary works centered on urban and industrial brownfields, soils’ multifaceted roles in CO2 transformation, water and nutrient cycling, agriculture, and as living bodies buffering against pollution have become ground for public discourse. What could this mean for the vadose zone in terms of reflecting on material flows in porous media beyond understandings of soil physics? In this paper, we draw on over 20 years of experience in studying the portrayal of soil in various arts genres to deliberate the potential of creative thinking about and thinking with the vadose zone.

KarstNSim: A graph-based method for 3D geologically-driven simulation of karst networks

In karst aquifers, groundwater flow is highly influenced by the interconnected underground cavities and conduits that form the karst network. Modeling karst flows requires the use of spatially distributed approaches accounting for these networks. Their exploration is, however, often complex, and mapping them using indirect methods such as geophysical ones has proven challenging. To overcome these limitations, stochastically simulating discrete karst networks should account for the uncertainties on conduit position and geometry. Only a few existing methods can reproduce realistic and diverse karst morphologies. In this work, we propose a new algorithm, KarstNSim, for simulating discrete karst networks, that incorporates field data to generate a range of possible karst network geometries. It relies on the computation of the shortest path between the inlets and outlets of the network with the use of an anisotropic cost function defined on an n-nearest neighbor graph conformal to geological and structural heterogeneities. This cost function represents the physico-chemical processes that govern speleogenesis – such as erosion and chemical weathering – providing simplified control over the geometry of the generated networks. Our approach reproduces the vadose-phreatic partition visible in the karst networks, by generating sub-vertical conduits in the unsaturated zone and sub-horizontal ones in the saturated part. It encompasses geological parameters such as inception surfaces, fractures, permeability, and solubility of layers, along with considering the hydrological context of recharge by assigning relative weights to the inlets. We simulate various synthetic models to demonstrate the influence of different input parameters on the spatial organization of the past and present karst flows. We also apply KarstNSim to a real case study – the Ribeaucourt system (Meuse, France) – and use available data to generate realistic karst networks. The lack of knowledge on the Ribeaucourt network geometry is alleviated by the use of multiple scenarios corresponding to different sets of weights in the cost function, resulting in a variety of network morphologies which represent the uncertainty on the real geometry. The geometrical and topological metrics computed on the simulated networks generally overlap with those of real ones, suggesting the reliability of the method.

Transport of graphene oxide in the capillary fringe: Insights from sandbox experiments and numerical simulation

Studying contaminant transport in the capillary fringe (CF), a crucial part of the vadose zone, offers insights into the mechanisms controlling pollution in soils and groundwater aquifers. This paper investigated contaminant transport in the CF by continuously injecting a conservative tracer (NaCl) and graphene oxide nanoparticle (GONP), an adsorptive contaminant, into a sandbox. After entering the CF from the unsaturated zone, both NaCl and GONP underwent lateral transport. The breakthrough curves (BTCs) for NaCl and GONP were derived from water samples collected at predetermined sampling holes. Subsequently, contaminant transport in the CF was modeled using a one-dimensional–two-dimensional (1D-2D) coupled hydrodynamic model. This model incorporated lateral dispersivity (αL = 1.198 cm) and longitudinal dispersivity (αT = 0.286 cm), calculated using a point-by-point method. The hydrodynamic dispersion coefficients obtained were then applied to the Brooks and Corey (BC) and the van Genuchten (VG) parametric models. The BC model more accurately simulated the NaCl migration compared to the VG model, leading to its application in simulating GONP transport in the CF. However, the simulated BTCs for GONP showed a lag behind the measured data, especially at high ionic strengths. This discrepancy was attributed to the variable adsorption partition coefficient of GONP under different ionic conditions. During the experiment, GONP adsorption onto the porous media's surface altered the capillary dynamics, notably increasing capillary rise height, decreasing seepage velocity, and reducing GONP dispersion. Therefore, it is necessary to consider the adsorption capacity of the contaminants in order to accurately assess their transport within the vadose zone.

Extreme precipitation accelerates nitrate leaching in the intensive agricultural region with thick unsaturated zones

Nitrate accumulation in the soil profile in the intensive agricultural region has been widely concerned in the world. However, the changes in nitrate accumulation characteristics caused by climate change, such as extremely high precipitation, are not well quantified, particularly for the regions with thick unsaturated zones. Here, we resampled the soil profiles taken in normal year (2020) after extreme precipitation year (2021) (>800 cm) in three regions in the southern Loess Plateau (LP) with three different water managements including rainfed orchards (n = 10), well-irrigated orchards (n = 4) and canal-irrigated orchards (n = 8). The accumulation amounts, peak depths, and accumulation depths of nitrate soil profiles of the different regions of two years were compared. The results showed that average nitrate accumulation in normal year at the rainfed region (800-cm depth), well-irrigated region (800-cm depth) and canal-irrigated region (1400-cm depth) were 5995 kg N ha−1, 9765 kg N ha−1, and 19,608 kg N ha−1, respectively. Compared with 2020, extreme precipitation in 2021 led to 56–91% reductions (2060–3702 kg N ha−1) in nitrate accumulation in 0–200 cm soil layer, and average nitrate leaching into the aquifer was >1390 kg N ha−1 in the canal-irrigated region. Average migration depths of nitrate peak in rainfed, well-irrigated and canal-irrigated regions were 92 cm, 115 cm, and 188 cm, respectively; as for nitrate accumulation depths, they were 10 cm, 80 cm and 108 cm, respectively. Vertically, the dried soil layer and paleosol layer (high clay content) in the canal-irrigated region significantly hindered nitrate deep migration caused by the extreme precipitation. The result highlights that extreme precipitation significantly accelerated nitrate leaching in the deep soil profiles, and future vulnerability and risk assessment studies must account for the impacts of extreme precipitation on nitrate leaching.

Modeling 1-D aqueous film forming foam transport through the vadose zone under realistic site and release conditions

Aqueous film forming foams (AFFFs) have been used to extinguish fires since the 1960s, leading to widespread subsurface contamination by per- and polyfluoroalkyl substances (PFAS), an essential component of AFFF. This study presents 1-D simulations of PFAS migration in the vadose zone resulting from AFFF releases. Simulation scenarios used soil profiles from three US Air Force (USAF) installations, encompassing a range of climatic conditions and hydrogeologic environments. A three-component mixture, representative of major constituents of AFFF, facilitated the exploration of competitive and synergistic effects of co-constituents on PFAS migration. To accurately capture unsaturated transport of PFAS in porous media, the model considers (1) surfactant-induced flow, (2) non-linear sorption to the solid phase, (3) competitive accumulation at the air-water interface, and (4) the moisture-dependence of the air-water interfacial area. Defined PFAS releases were consistent with fire training exercises, emergency responses, and accidental spills of record. Simulation results illustrate the importance of hydrogeologic, climatic, geochemical, and AFFF release conditions on PFAS transport and retention. Comparison of field observations and model simulations for Ellsworth AFB indicate that much of the PFOA and PFOS mass is associated with the air-water interface and the solid phase, which limits their migration potential in the vadose zone. Results also show that rates of migration in the aqueous phase are largely controlled by hydrogeologic properties, including recharge rates and hydraulic conductivity. AFFF spill scenarios varying in volume, concentration, and frequency reveal the importance of release characteristics in determining rates of PFAS migration and concentration peaks. Variability is attributed to non-linear sorption processes, where, contrary to simple linear partitioning formulations, transport is strongly affected by the concentration of PFAS species. Simulations also demonstrate the importance of modeling the AFFF as a mixture since competitive interfacial accumulation effects are shown to enhance the mobility of less surface-active PFAS compounds.

Expanding Karst Groundwater Tracing Techniques: Incorporating Population Genetic and Isotopic Data to Enhance Flow-Path Characterization

Karst aquifers are unique among groundwater systems because of variable permeability and flow-path organization changes resulting from dissolution processes. Over time, changes in flow-path connectivity complicate interpretations of conduit network evolution in karst hydrogeology. Natural and artificial tracer techniques have long provided critical information for protecting karst aquifers and understanding the potential impacts on ecosystems and human populations. Conventional tracer methods are useful in karst hydrogeologic studies for delineating flow paths and defining recharge, storage, and discharge properties. However, these methods only provide snapshots of the current conditions and do not provide sufficient information to understand the changes in interconnection or larger-scale evolution of flow paths in the aquifer over time. With advances in population genetics, it is possible to assess population connectivity, which may provide greater insights into complex groundwater flow paths. To assess this potential, we combined the more traditional approaches collected in this and associated studies, including artificial (dye) and natural (geochemistry, isotopes, and discharge) tracers, with the population genetic data of a groundwater crustacean to determine whether these data can provide insights into seasonal or longer changes in connections between conduits. The data collected included dye trace, hydrographs, geochemistry, and asellid isopod (Caecidotea bicrenenta) population genetics in Fern Cave, AL, USA, a 25 km-long cave system. Combined, these data show the connections between two separate flow paths during flood events as the downstream populations of isopods belonging to the same subpopulation were measured in both systems. Additionally, the sub-populations found in higher elevations of the cave suggest a highly interconnected unsaturated zone that allows for genetic movement in the vadose zone. Although upstream populations show some similarities in genetics, hydrologic barriers, in the form of large waterfalls, likely separate populations within the same stream.

Modeling evaluation of the impact of residual source material on remedial time frame at a former uranium mill site

Groundwater contamination at legacy uranium processing sites is an ongoing global challenge. Plumes at many uranium-contaminated sites are more persistent than originally predicted by groundwater modeling. Previous investigations of uranium plume persistence identified residual and secondary sources that contribute to plume longevity, but there is a remaining need to revise forecasted cleanup times using information about these ongoing sources. The purpose of this study is to investigate the quantitative impact of residual vadose zone sources of uranium on groundwater remediation time frame. This objective was approached by applying numerical uranium transport simulations and uncertainty analysis to a former uranium mill site in the southwestern United States. Information from recent site investigations provided details about the distribution and release characteristics of uranium accumulations in the vadose zone. The residual uranium characteristics were incorporated as decaying source terms in the transport model. A stochastic approach using an iterative ensemble smoother was applied for history matching, and the transport model was used to assess the impact of multiple remedial alternatives on forecasted time frame. The forecasted time frame to achieve the groundwater remediation goal for uranium by monitored natural attenuation is on the order of thousands of years, and treatment of the dissolved plume does not reduce the projected time frame. The large proportion of residual uranium mass remaining in the vadose zone and the gradual leaching rate due to the site's semiarid climate create a long-lived source that can sustain a dissolved plume for thousands of years despite an estimated 99% mass removal achieved during mill tailings disposal. Residual uranium in vadose zone sediments beneath former tailings impoundments could present comparable uranium plume persistence and remediation challenges at other legacy uranium mill sites in semiarid climates. Other remaining uranium-impacted sites are similarly complex, and forecasted remedial time frames are needed to effectively achieve compliance, manage risk, assess the benefits of additional treatment, manage and project costs, and support beneficial site reuse.

Thermal conductivity in a vadose zone: A novel spectral solution and application by the implication of significant coherence of thermal condition in the shallow soil water layer

AbstractIn the context of the heterogeneity in the unsaturated or vadose zone, accurately representing the analytical mechanisms and in‐situ water content within the soil layer poses a significant challenge. Particularly in shallow layers, thermal conditions exhibit rapid changes in response to evolving surface temperatures. This study proposes a hypothesis suggesting that the in situ heat mechanism may notably impact the soil water layer. The research introduces an innovative approach to theoretically uncover thermal conditions, including soil temperature, soil temperature gradients, and heat flux, within the shallow Quaternary gravel layer at various depths through spectral analysis of temporal observations. The study presents a stochastic inverse solution to estimate thermal conductivity by leveraging spectral analysis of soil heat flux and temperature gradients. The findings reveal that thermal conditions exhibit the most prominent periodic fluctuations during the diurnal process over a 24‐hour cycle. The soil temperature gradients and heat flux measurements at depths of 0.1, 0.3, 0.6, and 1.2 m demonstrate their ability to capture changes in soil temperature and air temperature to a certain extent within the frequency domain. Furthermore, the analysis highlights the intrinsic uncertainty and sensitivity of estimating thermal conductivity in heterogeneous soil environments. The wide variability observed in thermal conductivity values, coupled with their dependence on soil type and environmental conditions, underscores the need for careful consideration of these factors in future studies and modeling efforts. Applying the derived inverse spectral solution allows for determining thermal conductivity throughout the soil‐water system across depths ranging from 0.1 to 1.2 m. As a result, this research demonstrates the feasibility and practicality of assessing the thermal conductivity of the soil layer in conjunction with heat flux and temperature gradients through spectral analysis.

A novel hybrid magnetic material based on three-phase saturated core fault current limiter

Installing a three-phase saturated core fault current limiter (TSFCL) is an effective method to limit the increasing short-circuit current in high-voltage alternating current (HVAC) systems. However, the traditional TSFCL uses silicon steels as its working material for the iron core, resulting in poor current-limiting effects, requiring a large DC excitation current, and having a slow response speed. Nanocrystalline alloys have advantages such as a lower saturation point, higher permeability in the unsaturated zone, and a steep B–H curve, making them highly suitable for TSFCL’s applications. However, their poor capability which withstands fault hindered their large-scale commercialization. Therefore, this paper proposes a novel hybrid magnetic material based on TSFCL (HMTFCL), where the materials of the working limbs are replaced by nanocrystalline alloys. Compared to TSFCL, the DC excitation current of HMTFCL decreases by 60.5%. And, the HMTFCL reduces the peak short-circuit current by 23%. Last, the magnetic flux intensity inside the working limb of HMTFCL is reduced to zero more quickly.

P3D-BRNS v1.0.0: a three-dimensional, multiphase, multicomponent, pore-scale reactive transport modelling package for simulating biogeochemical processes in subsurface environments

Abstract. The porous microenvironment of soil offers various environmental functions which are governed by physical and reactive processes. Understanding reactive transport processes in porous media is essential for many natural systems (soils, aquifers, aquatic sediments or subsurface reservoirs) or technological processes (water treatment or ceramic and fuel cell technologies). In particular, in the vadose zone of the terrestrial subsurface the spatially and temporally varying saturation of the aqueous and the gas phase leads to systems that involve complex flow and transport processes as well as reactive transformations of chemical compounds in the porous material. To describe these interacting processes and their dynamics at the pore scale requires a well-suited modelling framework accounting for the proper description of all relevant processes at a high spatial resolution. Here we present P3D-BRNS as a new open-source modelling toolbox harnessing the core libraries of OpenFOAM and coupled externally to the Biogeochemical Reaction Network Simulator (BRNS). The native OpenFOAM volume-of-fluid solver is extended to have an improved representation of the fluid–fluid interface. The solvers are further developed to couple the reaction module which can be tailored for a specific reactive transport simulation. P3D-RBNS is benchmarked against three different flow and reactive transport processes: (1) fluid–fluid configuration in a capillary corner, (2) mass transfer across the fluid–fluid interface and (3) microbial growth with a high degree of accuracy. Our model allows for simulation of the spatio-temporal distribution of all biochemical species in the porous structure (obtained from μ-CT images), for conditions that are commonly found in the laboratory and environmental systems. With our coupled computational model, we provide a reliable and efficient tool for simulating multiphase, reactive transport in porous media.

A Conceptual Framework for Modeling Spatiotemporal Dynamics of Diesel Attenuation Capacity: A Case Study across Namyangju, South Korea

The vadose zone acts as a natural buffer against groundwater contamination, and thus, its attenuation capacity (AC) directly affects groundwater vulnerability to pollutants. A regression model from the previous study predicting the overall AC of soils against diesel was further expanded to the GIS-based overlay-index model. Among the six physicochemical parameters used in the regression model, saturation degree (SD) is notably susceptible to climatological and meteorological events. To accommodate the lack of soil SD historical data, a series of infiltration simulations were separately conducted using Phydrus code with moving boundary conditions (i.e., rainfall records). The temporal variation of SD and the resulting AC under transient conditions are captured by building a space–time cube using a temporal raster across the study area within the designated time frame (1997–2022). The emerging hot spot analysis (EHSA) tool, based on the Getis–Ord Gi* and Mann–Kendall statistics, is applied to further identify any existing pattern associated with both SD and AC in both space and time simultaneously. Under stationary conditions, AC decreases along depth and is relatively lower near water bodies. Similarly, AC cold spot trends also show up near water bodies under transient conditions. The result captures not only the trends across time but also shows the exact location where the changes happen. The proposed framework provides an efficient tool to look for locations that have a persistently low or a gradually decreasing ability to attenuate diesel over time, indicating the need for stricter management regulations from a long-term perspective.

Probing the exchange of CO2 and O2 in the shallow critical zone during weathering of marl and black shale

Abstract. Chemical weathering of sedimentary rocks can release carbon dioxide (CO2) and consume oxygen (O2) via the oxidation of petrogenic organic carbon and sulfide minerals. These pathways govern Earth's surface system and climate over geological timescales, but the present-day weathering fluxes and their environmental controls are only partly constrained due to a lack of in situ measurements. Here, we investigate the gaseous exchange of CO2 and O2 during the oxidative weathering of black shales and marls exposed in the French southern Alps. On six field trips over 1 year, we use drilled headspace chambers to measure the CO2 concentrations in the shallow critical zone and quantify CO2 fluxes in real time. Importantly, we develop a new approach to estimate the volume of rock that contributes CO2 to a chamber, and assess effective diffusive gas exchange, by first quantifying the mass of CO2 that is stored in a chamber and connected rock pores. Both rock types are characterized by similar contributing rock volumes and diffusive movement of CO2. However, CO2 emissions differed between the rock types, with yields over rock outcrop surfaces (inferred from the contributing rock volume and the local weathering depths) ranging on average between 73 and 1108 tCkm-2yr-1 for black shales and between 43 and 873 tCkm-2yr-1 for marls over the study period. Having quantified diffusive processes, chamber-based O2 concentration measurements are used to calculate O2 fluxes. The rate of O2 consumption increased with production of CO2, and with increased temperature, with an average O2:CO2 molar ratio of 10:1. If O2 consumption occurs by both rock organic carbon oxidation and carbonate dissolution coupled to sulfide oxidation, either an additional O2 sink needs to be identified or significant export of dissolved inorganic carbon occurs from the weathering zone. Together, our findings refine the tools we have to probe CO2 and O2 exchange in rocks at Earth's surface and shed new light on CO2 and O2 fluxes, their drivers, and the fate of rock-derived carbon.

Novel approach to roof rainwater harvesting and aquifer recharge in an urban environment: Dry and wet infiltration wells comparison

In urban environments there is a severe reduction of infiltration and groundwater recharge due to the existence of large impervious areas. During rain events, large volumes of water that could have recharged groundwater and surface water bodies are diverted into the municipal drainage system and lost from the freshwater storage. Moreover, extreme rain events impose high peak flows and large runoff volumes, which increase the risk of urban floods. Recent studies have suggested the use of rainwater harvesting for groundwater recharge, as a plausible solution for these challenges in dense urban environments. While the benefits of this approach are well understood, research on its practical, engineering, and hydrological aspects is relatively limited. The objective of the present study was to examine the use of infiltration wells for groundwater recharge with harvested rainwater collected from building rooftops under Mediterranean climate conditions. Two types of wells with similar hydraulic and technical properties were examined: a well that reaches the groundwater (wet well); and a well that discharges the harvested water into the unsaturated zone (dry well). Infiltration capacities of the wells were compared in controlled experiments conducted during summer months, and in operational recharge of harvested rainwater, during winter. Both dry and wet wells were found to be suitable for purposes of groundwater recharge with rooftop-harvested rainwater. Infiltration capacity of the wet well was about seven times greater than the infiltration capacity of the dry well. While the infiltration capacity of the wet well was constant throughout the entire length of the study (∼10 m3/h/m), the dry well infiltration capacity improved during winter (from 0.5 m3/h/m to 1.5 m3/h/m), a result of development of the dry well with time. Considering Tel-Aviv, Israel, as a case study for a dense modern city in a Mediterranean climate, it is demonstrated herein that the use of infiltration wells may reduce urban drainage by ∼40 %.

Correction: Presence of Hummock and Hollow Microtopography Reflects Shifting Balances of Shallow Subsidence and Root Zone Expansion Along Forested Wetland River Gradients

Correction: Presence of Hummock and Hollow Microtopography Reflects Shifting Balances of Shallow Subsidence and Root Zone Expansion Along Forested Wetland River Gradients

Genome-Resolved Metagenomics and Denitrifying Strain Isolation Reveal New Insights into Microbial Denitrification in the Deep Vadose Zone.

The heavy use of nitrogen fertilizer in intensive agricultural areas often leads to nitrate accumulation in subsurface soil and nitrate contamination in groundwater, which poses a serious risk to public health. Denitrifying microorganisms in the subsoil convert nitrate to gaseous forms of nitrogen, thereby mitigating the leaching of nitrate into groundwater. Here, we investigated denitrifying microorganisms in the deep vadose zone of a typical intensive agricultural area in China through microcosm enrichment, genome-resolved metagenomic analysis, and denitrifying bacteria isolation. A total of 1000 metagenome-assembled genomes (MAGs) were reconstructed, resulting in 98 high-quality, dereplicated MAGs that contained denitrification genes. Among them, 32 MAGs could not be taxonomically classified at the genus or species level, indicating that a broader spectrum of taxonomic groups is involved in subsoil denitrification than previously recognized. A denitrifier isolate library was constructed by using a strategy combining high-throughput and conventional cultivation techniques. Assessment of the denitrification characteristics of both the MAGs and isolates demonstrated the dominance of truncated denitrification. Functional screening revealed the highest denitrification activity in two complete denitrifiers belonging to the genus Pseudomonas. These findings greatly expand the current knowledge of the composition and function of denitrifying microorganisms in subsoils. The constructed isolate library provided the first pool of subsoil-denitrifying microorganisms that could facilitate the development of microbe-based technologies for nitrate attenuation in groundwater.