Vadose Research Articles

Modeling the influence of forest vegetation and climate change on the long-term performance of a cover with capillary barrier effects used to control acid mine drainage: the Lorraine case study

ABSTRACT The Lorraine mine site in Témiscamingue region (Québec, Canada) was reclaimed with a cover with capillary barrier effects (CCBE) to limit acid mine drainage generation. The main objective of this study is to predict the long-term hydrogeological behaviour of the CCBE at Lorraine site, including the influence of vegetation and climate change (CC). Numerical simulations showed that performance of the CCBE is: (1) not affected by forest vegetation and CC in the saturated zone in the north of the Lorraine site and (2) more sensitive to vegetation in the unsaturated zone in the south where the phreatic level is low.

Modulation of evapotranspiration and stream runoff by weathered bedrock in arid and semi-arid mountains

Earth's Critical Zone exhibits remarkable heterogeneity and complexity. Hence, further investigation is required to examine the composition of Earth's Critical Zone as well as the diverse eco-hydrological patterns they exhibit under varying climatic and geological circumstances. This exploration should primarily be conducted through the investigation and experiments of the hillslope unit, where the topography and weathered bedrock are representative, with particular emphasis on semi-arid regions where water resources serve as the primary limiting factor. Here, we have determined that the structure of the weathering profile displays systematic variation across the topography and heterogeneous landscape on uninterrupted slopes. Differences in the structure of the subsurface critical zone led to differences in its water storage capacity at the same time. Runoff in alpine shrubs and forests was dominated by subsurface runoff, and grassland was dominated by surface runoff. In the alpine shrub immediately adjacent to the watershed, an estimated quantity of 129 mm of water is stored within the unsaturated zone of the soil, serving as exchange water to replenish moisture in the underlying bedrock. In contrast to alpine shrubs, an estimated quantity of 62.7 mm of water originates from the unsaturated zone of soil and weathered bedrock in the forest. However, approximately 21.1 mm of moisture is unavailable to plants. The soil water storage in grasslands exhibits a decline throughout the growing season, with a subsequent augmentation occurring solely after substantial precipitation events exceeding 20 mm. In wet years, dynamic storage predominantly manifests as groundwater saturation throughout the entire ground and high subsurface runoff. In dry years, the limited runoff response indicates that the catchment's dynamic water storage primarily comprises “indirect” water storage, which predominantly resides within the soil, saprolite, and weathered rock below the “field capacity”, subsequently being released into the atmosphere through evapotranspiration.

Laboratory and numerical simulations of infiltration process and solute transport in karst vadose zone

Groundwater in karst systems is a vital resource, which is often highly vulnerable to contaminants that infiltrate through the vadose zone. Understanding the processes of water and contaminant infiltration in the vadose zone is extremely important for protecting and managing karst water resources. A laboratory-scale physical model based on a typical conceptual model of the karst vadose zone was constructed to characterize the influence of the rainfall intensity and internal structure on the infiltration process and solute transport in the karst vadose zone under controlled conditions. Five factors, including the rainfall intensity, the surface slope, the degree of karstification, the thickness of the transfer zone, and the existence of epikarst, were considered in laboratory experiments. A corresponding lumped-parameter model was subsequently developed to further analyze and investigate the water and solute infiltration processes in the laboratory scale physical model. Recession- and breakthrough curves generated from the lab experiments and the water and solute loss from each reservoir generated from the simulation results were analyzed. The downward trend of breakthrough curves' peak value and the variation of solute loss's proportion in three rainfall events demonstrate that both the increase in thickness and karstification degree will intensify the buffer effect of the karst vadose zone, and this buffer effect is mainly affected by the slow flow system which can be inferred from the increasing percentage of cumulative discharge of slow reservoir (qS). With increasing karstification degree, the recession coefficient (α) exhibits opposite trends depending on the existence of epikarst or not, indicating the karstification degree has a more significant effect on hydraulic conductivity (K) when an epikarst is present while affecting effective porosity (φ) more noticeably in the absence of epikarst. The rainfall intensity affects the results of both experimental and numerical simulation most significantly. The percentage of discharge and solutes flowing from the slow reservoir (qS and msoil_S) have negative correlations with rainfall intensity, suggesting a greater influence on fast flow with increasing rainfall intensity.

Multivariate and scale-dependent controls of deep soil carbon after afforestation in a typical loess-covered region

Afforestation is beneficial to improving soil carbon pools. However, due to the lack of deep databases, the variations in soil carbon and the combined effects of multiple factors after afforestation have yet to be adequately explored in >1 m deep soils, especially in areas with deep-rooted plants and thick vadose zones. This study examined the multivariate controls of soil organic carbon (SOC) and inorganic carbon (SIC) in 0–18 m deep under farmland, grassland, willow, and poplar in loess deposits. The novelty of this study is that the factors concurrently affecting deep soil carbon were investigated by multiwavelet coherence and structural equation models. On average, the SOC density (53.1 ± 5.0 kg m−2) was only 12% of SIC density (425.4 ± 13.8 kg m−2), with depth-dependent variations under different land use types. In the soil profiles, the variations in SOC were more obvious in the 0–6 m layer, while SIC variations were mainly observed in the 6–12 m layer. Compared with farmland (SOC: 17.0 kg m−2; SIC: 122.9 kg m−2), the plantation of deciduous poplar (SOC: 28.5 kg m−2; SIC: 144.2 kg m−2) increased the SOC and SIC density within the 0–6 m layer (p < 0.05), but grassland and evergreen willow impacted SOC and SIC density insignificantly. The wavelet coherence analysis showed that, at the large scale (>4 m), SOC and SIC intensities were affected by total nitrogen-magnetic susceptibility and magnetic susceptibility-water content, respectively. The structural equation model further identified that SOC density was directly controlled by total nitrogen (path coefficient = 0.64) and indirectly affected by magnetic susceptibility (path coefficient = 0.36). Further, SOC stimulated the SIC deposition by improving water conservation and electrical conductivity. This study provides new insights into afforestation-induced deep carbon cycles, which have crucial implications for forest management and enhancing ecosystem sustainability in arid regions.

Spatial variations in groundwater hydrochemistry, sources, and controls across catchments on Chinese Loess Plateau

Study regionThis study was carried out at the Chinese Loess Plateau (CLP), which is located in the middle reach of the Yellow River in the northern China. Study focusUnderstanding groundwater hydrochemistry and its controlling factors in regions like CLP with thick unsaturated zones is challenging due to limited accessibility. The primary objective of this study was to investigate groundwater hydrochemistry, sources, and controls, and assess groundwater suitability for drinking and irrigation purposes across 14 major catchments in CLP. New hydrological insights for the regionOur results showed that the composition of groundwater dissolved solutes is primarily influenced by weathering rock, particularly silicate rock. However, human activities, especially in agriculture and urban areas, also have a notable impact on groundwater solutes. The presence of elevated nitrate in groundwater in various regions underscores the urgency of addressing long-term nitrogen pollution. Environmental changes such as climate change, erosion rates, and land-use patterns are key factors shaping groundwater quality. As air temperature decreases and erosion intensifies, along with more intensive land use, ion concentrations in groundwater tend to rise. Note that, approximately 91% of groundwater in CLP is suitable for direct drinking, but ensuring proper drainage is crucial to prevent sodium damage during irrigation. These findings enhance our understanding of the factors influencing groundwater solutes and are essential for effectively managing groundwater quality in CLP.

Aeolian desert dust as a primary estuarine sediment source: Fine sediment transport and dynamics in the northern Arabian/Persian gulf

Fluvial discharge and freshwater runoff are recognized as the main sediment sources in many estuaries worldwide. Precipitation is, however, limited in arid and semi-arid areas, allowing for high aeolian dust activity, and inducing limited river runoff. As such, aeolian dust could contribute to the estuarine sediment sources in such areas. Uniquely, the Northern Arabian/Persian Gulf (NG) is continuously loaded with vast amounts of aeolian desert dust annually. This study focuses on the role of desert dust in the distribution of suspended sediment in the Northern Arabian/Persian Gulf (NG), employing in-situ measurements of physical parameters (collected in 2015–2016), satellite imagery, and a flexible mesh hydrodynamic model (Delft3D-FM) coupled with a sediment transport model (Delft-WAQ). The model showed reasonable performance in representing the complex system of NG. The study revealed that more than 90% of the sediment along the NG is aeolian dust, with the remaining being fluvial. Owing to salinity-induced density gradients, estuarine circulation keeps the fluvial sediment confined to the estuarine zone of the Shatt Al-Arab even during the events of increased freshwater discharges. In the NG most of the suspended sediments have an aeolian origin. They were found in shallow zones within the 5 m isobath, including the estuary of Shatt Al-Arab where they are mixed with the fluvial sediment fractions. Tides were found to play a key role in sediment mixing and horizontal dispersion in the shallow area of the NG most of the time, while the Shamal increased the transport of fine sediment to offshore areas (seaward). Based on current knowledge, this study represents a rare instance where estuarine sediment dynamics are primarily governed by aeolian deposits from a nearby desert.

Wflow_sbm v0.7.3, a spatially distributed hydrological model: from global data to local applications

Abstract. The wflow_sbm hydrological model, recently released by Deltares, as part of the Wflow.jl (v0.7.3) modelling framework, is being used to better understand and potentially address multiple operational and water resource planning challenges from a catchment scale to national scale to continental and global scale. Wflow.jl is a free and open-source distributed hydrological modelling framework written in the Julia programming language. The development of wflow_sbm, the model structure, equations and functionalities are described in detail, including example applications of wflow_sbm. The wflow_sbm model aims to strike a balance between low-resolution, low-complexity and high-resolution, high-complexity hydrological models. Most wflow_sbm parameters are based on physical characteristics or processes, and at the same time wflow_sbm has a runtime performance well suited for large-scale high-resolution model applications. Wflow_sbm models can be set a priori for any catchment with the Python tool HydroMT-Wflow based on globally available datasets and through the use of point-scale (pedo)transfer functions and suitable upscaling rules and generally result in a satisfactory (0.4 ≥ Kling–Gupta efficiency (KGE) &lt; 0.7) to good (KGE ≥ 0.7) performance for discharge a priori (without further tuning). Wflow_sbm includes relevant hydrological processes such as glacier and snow processes, evapotranspiration processes, unsaturated zone dynamics, (shallow) groundwater, and surface flow routing including lakes and reservoirs. Further planned developments include improvements on the computational efficiency and flexibility of the routing scheme, implementation of a water demand and allocation module for water resource modelling, the addition of a deep groundwater concept, and computational efficiency improvements through for example distributed computing and graphics processing unit (GPU) acceleration.

Volume change characteristic of compacted bentonite during hydration upon high suction change: Experimental study and constitutive simulation

The assessment and predicting of volume change characteristic in compacted bentonite during hydration upon high suction change are crucial considerations in the design and construction of geological repositories. However, the mechanism for volume change behavior of bentonite upon suction change path is not yet well understood. Existing elastoplastic models show some limitations for describing hydration deformation of compacted bentonite upon high suction change path. In this work, suction-controlled hydration deformation tests were conducted on compacted GMZ01 bentonite under constant vertical stresses. The various forms of void ratio-suction curves observed suggested that the swelling-collapse behavior during hydration was influenced by the over-consolidation ratio (OCR). Furthermore, it was observed that each hydration deformation curve (e-lns) exhibited a saturation suction value. When the suction decreased to the saturation suction, the bentonite reached a state where further hydration was no longer possible, resulting in the absence of additional swelling or collapse deformation. The unique relationship between saturation suction and stress applied, i.e. critical saturated state (CSS) curve was proved to existed, which dividing the stress space (s-p plane) into unsaturated and saturated zones. A model framework with the concept of CSS curve was introduced and the suction-induced volume change equation was determined based on the tests results. The simulation results of hydration deformation tests demonstrated excellent agreement with the experimental data. The test results in this work provide the important mechanical parameters for the design and construction of geological repositories. Meanwhile, the proposed model could describe volume change characteristic of compacted bentonite during hydration upon high suction change, showing great significance for the safety assessment of repositories.

Rising Damp: The Rising of Moist Earth-Air in the Mogao Grottoes, China

ABSTRACT Usually, the term ‘rising damp’ refers to capillary water rising from the ground which may damage architectural heritage. In this work, the essence of rising damp in extremely arid regions and its driving force are revealed based on experiments monitoring the relative humidity (RH) and atmospheric pressure (AP) in the Dunhuang Mogao Grottoes. The air in the vadose zone, the unsaturated region between ground level and the top of the water table, is here referred to as ‘earth-air’. When the AP rises, the earth-air is compressed, and atmospheric air enters into the soil. Then, when the AP drops, moist earth-air expands and rises into the structure. The RH in the soil is thus negatively correlated with the AP, yielding a correlation coefficient of up to –0.94. Under the action of this long-term dry–wet alternation, the salt present in the building near ground undergoes repeated cycles of crystallization and dissolution, resulting in efflorescence and a deterioration zone. Therefore, the deterioration due to rising damp in extremely arid regions is caused by the rising of moist earth-air rather than capillarity. The height to which it rises is directly proportional to the amplitude of the daily AP fluctuation and thickness of the vadose zone, exhibiting a bimodal fluctuation pattern on a daily scale. The discovery of this mechanism of rising damp provides a scientific basis for preventive conservation interventions.

Effects of microtopography on patterns and dynamics of groundwater–surface water interactions

Groundwater–surface water interaction is greatly affected by surface topography and agricultural activities in managed agricultural regions. Such an interaction increases the complexity of hierarchical groundwater flow systems. However, the combined effects of microtopography and agricultural activities on groundwater flow patterns remain unclear. In this study, an integrated groundwater–surface water model was used to quantify the impact of microtopographic variations on groundwater flow patterns in an agricultural region in northwest China. Numerical experiments were conducted under different scenarios representing various degrees of microtopography. The results indicated that microtopography plays a pivotal role in the transfer of signals from spatiotemporally varying surface processes to the groundwater flow. Without consideration of the microtopography, the stationary fractional Gaussian noise–type infiltration processes in the vadose zone–saturated zone interface may be incorrectly simulated as non-stationary fractional Brownian motion. Furthermore, fractal analysis reveals that the scaling exponents of the pressure head are highly sensitive within the vadose zone and local flow subsystems. Microtopography plays an important role in shaping the groundwater response time (GRT) pattern, which results in a smaller average GRT but a larger core range. Contrarily, neglecting microtopography may lead to overestimation of the water table ratio (WTR), resulting in erroneous identification of groundwater patterns as topography-controlled types. This study highlights the importance of microtopography in the formation and evolution of nested groundwater flow systems and provides guidance for groundwater–surface water interaction simulations in managed regions.

Reduction in Apparent Permeability Owing to Surface Precipitation of Solutes by Drying Process and Its Effect on Geological Disposal

Disposal tunnels in geological repositories are ventilated continuously for over 50 years until their closure. Under these conditions, an unsaturated zone of mixed liquid and gas phases forms around the tunnels. Moreover, drying is assumed to progress from the host rock to the tunnels. To understand these drying processes, this study investigated the migration and precipitation of solutes via capillary forces during drying in packed columns using silica sand or glass beads as packed layers and X-ray CT analysis. In addition, the apparent permeability of a column packed with silica sand containing precipitation was examined using a flow experiment. The results indicate that the precipitation and accumulation of solutes were significant near the drying surfaces of the columns. The apparent mass transfer coefficient at a relatively early stage of the drying process indicates that the migration rate of solutes depends strongly on the capillary forces during the drying process. Furthermore, the apparent permeability of the columns with precipitation decreased significantly. These indicate that the precipitation and accumulation of solutes with drying in the groundwater reduce the porosity and permeability, and the advection of groundwater around the repository may be suppressed.

Co-transport behavior of Am(III) and natural colloids in the vadose zone sediments

Americium (III) (Am(III)) in the natural environment is considered immobile due to its low solubility, strong adsorption, and high affinity to solid surfaces. However, the presence of natural colloids may carry Am(III) transport for long distance. The individual and co-transport behaviors of Am(III) and natural colloids through the unsaturated packed columns were investigated under the influence of pH, electrolyte concentration, velocity, Am(III) concentration and natural colloids concentration. Under all experimental conditions, Am(III) individual transport construct sight breakthrough curves (BTCs, CAm/C0 < 3%), but the presence of natural colloids increased the BTCs plateau of Am(III) significantly (30% < CAm/C0 < 80%), indicating that the colloids were able to promote Am(III) transport in the unsaturated porous media. DLVO theoretical calculations reveal that the increased pH and decreased electrolyte concentration lead to a rase in electrostatic repulsion, and the natural colloids tend to be dispersed and stabilized, which facilitates elution. In addition to this, the increase of velocity and colloids concentration will lead to greater breakthrough of natural colloids. The non-equilibrium two-site model and the two-site kinetic retention model well-described the BTCs of Am(III) and natural colloids, respectively. This study provide new insights into the behavior of natural colloids carrying the Am(III) into aquifers through the vadose zone sediments.

High-temperature oxidation behavior of TA15 aerospace titanium alloy at 500 °C and 800 °C

This study investigates the high-temperature oxidation behavior and performance of TA15 Ti alloy subjected to 500 and 800 °C service conditions. The oxidation process, microstructural evolution, mechanical property changes, and comprehensive performance were analyzed. At 500 °C, the alloy developed a distinctive blue-colored surface primarily composed of Ti6O, indicating a slow growth rate of the oxide layer. In contrast, the 800 °C condition led to severe oxidation with notable TiO2 and Al2O3 phases and extensive peeling of the oxidation layer. Microstructural analysis revealed that at 500 °C, there was minimal surface oxidation, and the microstructure showed a shallow heat-affected zone. At 800 °C, repeated oxidation and peeling cycles resulted in grain growth and a significant decrease in the β-phase proportion. Mechanical property assessments indicated an increase in micro Vickers hardness at both temperatures, with the 800 °C sample exhibiting higher values. However, tensile strength decreased significantly at 800 °C due to the severe oxidation layer peeling, reducing the effective thickness during tensile testing. In summary, TA15 demonstrated good high-temperature service performance at 500 °C, while at 800 °C, severe oxidation and peeling adversely affected the overall performance. This research provides insights into the differential high-temperature behavior of TA15, crucial for applications in elevated temperature environments.

Assessing Light Non-Aqueous Phase Liquids in the Subsurface Using the Soil Gas Rn Deficit Technique: A Literature Overview of Field Studies

222Radon (Rn) was proposed in the late 1990s as a naturally occurring tracer for light non-aqueous phase liquids (LNAPLs) in the subsurface, due to its preferential partitioning behavior in the non-aqueous phase, resulting in a reduction in Rn activities in areas with LNAPLs in the subsurface compared to unimpacted areas (Rn deficit). The Rn deficit technique emerged as a cost-effective, non-invasive, and sustainable method to rapidly identify and quantify LNAPLs, for the characterization and monitoring of contaminated sites. This paper presents an overview of the technique and its field applications, with a specific focus on the use of the method in the vadose zone based on soil gas Rn measurements. Although various configurations have shown favorable outcomes, limitations persist in the application of the soil gas Rn deficit technique. Deep LNAPL contamination, soil matrix heterogeneity, and temporal variations in Rn emissions pose challenges to quantitative evaluations of LNAPL contamination. Recognizing these factors is crucial for site-specific assessments. This review aims to highlight both the strengths and limitations of the method, providing insights into potential areas for future research while acknowledging the positive outcomes achieved in different configurations over the past decades.

Diffusion capability of Cs-137 in compacted soil and bentonite under groundwater-saturated condition

Several researches have been studied the diffusion coefficient (Da) of radionuclide in soil samples in varies location around the world; however, the effect of cations in groundwater as a medium of saturation for Cs-137 diffusion in soil and bentonite remains to be accomplished. Compacted soil and bentonite layers are utilized in the construction of radioactive waste disposal sites, which serve to prevent the release of Cs-137 from the repository to the unsaturated zone. The samples used were soil from the prospective disposal site in SP4—Serpong Nuclear Center (SNC) and bentonite taken from Santrijaya, Tasikmalaya. The diffusion coefficient (Da) is used as indicator of Cesium-137 diffusion in the samples. Experimental measurements were conducted on a diffusion column unit that varied in diffusion time and density conditions to determine diffusion coefficients. The samples were cut into slices and analyzed using Multichannel Analyzer (MCA) unit to measure Cs-137 radioactivity. The measurement results were then incorporated into Fick's Law equation to obtain the value of the diffusion coefficient (Da) of Cs-137 in soil and bentonite. The experimental results showed that the diffusion coefficients of samples with groundwater simulation are in the range of 10–12 m2 s−1 for the local soil and bentonite samples, respectively. Compare to the similar experiment that was conducted previously, the results indicated that the groundwater content could slightly impact the increase in the diffusion coefficient value of Cs-137.

Consideration of Vadose Zone Moisture Dynamics in Remediation of PFAS‐Impacted Soils

AbstractOilfield brines are soluble, non‐volatile, and recalcitrant contaminants that are vulnerable to capillary effects and evapoconcentration and have been known to recontaminate clean fill soils emplaced after excavations. This problem may also occur for some poly‐ and perfluoroalkyl substances (PFAS) at PFAS sites. In this paper, we consider the soil moisture characteristic curve, wetting front dynamics, and the moisture redistribution cycle to understand the potential behavior of different PFAS in the vadose zone. Drawing on experience with remediation of soils impacted by brine spills, we provide recommendations on the design and implementation of soil excavation and cover remedies such as capillary barriers to prevent recontamination.

Soil water dynamics and deep percolation in an agricultural experimental area of the North China Plain over the past 50 years: Based on field monitoring and numerical modeling

The unregulated irrigation systems used in the late 20th century have led to increasingly severe deep percolation (DP) in the agricultural irrigation areas of the North China Plain. This has become an important factor limiting the efficient utilization of water resources and sustainable environmental development in these irrigation areas. However, the thick vadose zone is hydrodynamically exceptionally complex. The soil hydrological cycle is constantly changing under the influence of major climate change and human activity, thereby causing changes in DP that are difficult to quantify accurately. Here, the Luancheng Agricultural Irrigation District in North China was selected for a continuous 20-year in situ experiment. Soil-water dynamics were monitored using neutron probes and tensiometers, to determine the complete annual soil-water cycle and the hydrodynamic properties of the thick vadose zone irrigation district. For 1971–2021, DP was simulated using the HYDRUS-1D model and was verified by fitting observed values. Soil water content (SWC) exhibited similar trends in years that differed in terms of the amounts of irrigation and precipitation. The 0–100 cm soil layer was significantly affected by precipitation and other factors, and recharge >60 mm/d caused percolation. DP occurred mostly after irrigation or during the period of intensive precipitation in July–October. The maximum percolation rate was 16.9 mm/d under the present irrigation method. The main factors leading to DP were soil water storage capacity (R2 = 0.86) and precipitation (R2 = 0.54). Under the evolution of irrigation measures in the last 50 years, the average DP has gradually decreased from 574.2 mm (1971–1990) to 435.5 mm (2005–2021). However, a substantial amount of precipitation and irrigation water infiltrated the soil and percolated into the deep soil layer without being utilized by the crop. Therefore, there is an urgent need to consider measures to reduce DP to improve water-use efficiency in agriculture.

Nature-based solution enhances resilience to flooding and catalyzes multi-benefits in coastal cities in the Global South

Coastal cities are facing a rise in groundwater levels induced by sea level rise, further triggering saturation excess flooding where groundwater levels reach the topographic surface or reduce the storage capacity of the soil, thus stressing the existing infrastructure. Lowering groundwater levels is a priority for sustaining the long-term livelihood of coastal cities. In the absence of studies assessing the possibility of using tree-planting as a measure of alleviating saturation excess flooding in the context of rising groundwater levels, the multi-benefit nature of tree-planting programs as sustainable Nature-based solutions (NBSs) in coastal cities in the Global South is discussed. In environments where groundwater is shallow, trees uptake groundwater or reduce groundwater recharge, thereby contributing to lower groundwater levels and increasing the unsaturated zone thickness, further reducing the risk of saturation excess flooding. Tree-planting programs represent long-term solutions sustained by environmental factors that are complementary to conventional engineering solutions. The multi-benefit nature of such NBSs and the expected positive environmental, economic, and social outcomes make them particularly promising. Wide social acceptance was identified as crucial for the long-term success of any tree-planting program, as the social factor plays a major role in addressing most weaknesses and threats of the solution. In the case of Nouakchott City (Mauritania), where a rise in groundwater levels has led to permanent saturation excess flooding, a tree-planting program has the potential to lower the groundwater levels, thereby reducing flooding during the rainy season.

Inclusion of bedrock vadose zone in dynamic global vegetation models is key for simulating vegetation structure and function

Abstract. Across many upland environments, soils are thin and plant roots extend into fractured and weathered bedrock where moisture and nutrients can be obtained. Root water extraction from unsaturated weathered bedrock is widespread and, in many environments, can explain gradients in vegetation community composition, transpiration, and plant sensitivity to climate. Despite increasing recognition of its importance, the “rock moisture” reservoir is rarely incorporated into vegetation and Earth system models. Here, we address this weakness in a widely used dynamic global vegetation model (DGVM; LPJ-GUESS). First, we use a water flux-tracking deficit approach to more accurately parameterize plant-accessible water storage capacity across the contiguous United States, which critically includes the water in bedrock below depths typically prescribed by soil databases. Secondly, we exploit field-based knowledge of contrasting plant-available water storage capacity in weathered bedrock across two bedrock types in the Northern California Coast Ranges as a detailed case study. For the case study in Northern California, climate and soil water storage capacity are similar at the two study areas, but the site with thick weathered bedrock and ample rock moisture supports a temperate mixed broadleaf–needleleaf evergreen forest, whereas the site with thin weathered bedrock and limited rock moisture supports an oak savanna. The distinct biomes, seasonality and magnitude of transpiration and primary productivity, and baseflow magnitudes only emerge from the DGVM when a new and simple subsurface storage structure and hydrology scheme is parameterized with storage capacities extending beyond the soil into the bedrock. Across the contiguous United States, the updated hydrology and subsurface storage improve annual evapotranspiration estimates as compared to satellite-derived products, particularly in seasonally dry regions. Specifically, the updated hydrology and subsurface storage allow for enhanced evapotranspiration through the dry season that better matches actual evapotranspiration patterns. While we made changes to both the subsurface water storage capacity and the hydrology, the most important impacts on model performance derive from changes to the subsurface water storage capacity. Our findings highlight the importance of rock moisture in explaining and predicting vegetation structure and function, particularly in seasonally dry climates. These findings motivate efforts to better incorporate the rock moisture reservoir into vegetation, climate, and landscape evolution models.

Modeling water flow and volumetric water content in a degraded peat comparing unimodal with bimodal porosity and flux with pressure head boundary condition

AbstractDegraded peatlands release large amounts of greenhouse gases. The development of effective mitigation and management measures requires an understanding of relevant site‐specific biogeochemical and hydraulic processes. However, the simulation of water fluxes and vadose zone state variables of degrading peatlands relies on proper process description, parameterization of hydraulic functions, and representation of the boundary conditions. The objective of this study was to analyze the effects of unimodal versus bimodal soil hydraulic functions and pressure head versus flux‐type lower boundary conditions (LBCs) on the calculated hydraulic characteristics of a degraded peat profile. HYDRUS‐1D was used to study the hydraulic flow dynamics parameterized with data from a weighable groundwater lysimeter for the period from May 1 to December 31, 2019. Simulations comparing uni‐ and bimodal hydraulic functions showed only minor differences. Simulations of soil water pressure at a depth of 30 cm using a flux‐type LBC (RMSE: 27 cm, where RMSE is root mean square error) performed better than simulations using a pressure head LBC (RMSE: 48 cm). The pressure head LBC performed better at simulating volumetric water contents in 30‐cm depth than the flux LBC variant (RMSE: 0.05 vs. 0.09 cm3 cm−3). For specific site conditions with a shallow, fluctuating groundwater table and temporary air entrapment, the choice of LBC was important for a more accurate simulation of soil water fluxes and volumetric water content.