skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: REWTCrunch: A Modeling Framework for Vegetation Induced Reactive Zone Processes in the Critical Zone
Abstract Vegetation optimizes its geochemical environment for resource management via root exudation. We refer to the soil zone where biogeochemical behavior is significantly influenced, directly or indirectly, by root processes as the vegetation induced reactive zone (VIRZ). Root exudates react with VIRZ soil substrates creating temporally variable chemical environments through depth that extend below the rooting zone, impacting weathering, and releasing solutes and gases. We present a new framework, REWTCrunch, to capture VIRZ dynamics by integrating three modeling advances: the multicomponent reactive transport model CrunchFlow, the root exudation model REWT, and the multilayer canopy‐root ecohydrologic model MLCan. REWTCrunch's high‐resolution, process‐based simulation of root exudation, and the transport and transformation of carbon (C) and nutrients according to mass‐balanced and charge‐balanced reaction networks gives new insight into vertically resolved root‐soil‐microbe‐water interactions and their influence on solute fluxes at a daily timescale. We benchmark REWT and CrunchFlow, illustrate coupling mechanisms, and present REWTCrunch simulations for an agricultural site in the US Midwest. Results demonstrate root‐sourced reactive C can augment or reduce solute concentrations in the soil by several orders of magnitude. Silicate weathering products illustrate after‐harvest effects of plant C inputs in leaching patterns. Calcium simulations reveal the development of a stable weathering front. Aluminum concentrations are particularly responsive to root‐sourced reactivity, and analysis of leaching concentration versus leaching flux indicates hysteresis behavior. REWTCrunch significant improves our ability to simulate the link between root processes and soil biogeochemistry, thereby filling an important gap in the numerical simulation of root processes, weathering, and long‐term soil health.  more » « less
Award ID(s):
2012850 1835543
PAR ID:
10363350
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Biogeosciences
Volume:
127
Issue:
2
ISSN:
2169-8953
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Lithium isotope ratios (δ7Li) of rivers are increasingly serving as a diagnostic of the balance between chemical and physical weathering contributions to overall landscape denudation rates. Here, we show that intermediate weathering intensities and highly enriched stream δ7Li values typically associated with lowland floodplains can also describe small upland watersheds subject to cool, wet climates. This behavior is revealed by stream δ7Li between +22.4 and +23.5‰ within a Critical Zone observatory located in the Cévennes region of southern France, where dilute stream solute concentrations and significant atmospheric deposition otherwise mask evidence of incongruence. The water‐rock reaction pathways underlying this behavior are quantified through a multicomponent, isotope‐enabled reactive transport model. Using geochemical characterization of soil profiles, bedrock, and long‐term stream samples as constraints, we evolve the simulation from an initially unweathered granite to a steady state weathering profile which reflects the balance between (a) fluid infiltration and drainage and (b) bedrock uplift and soil erosion. Enriched stream δ7Li occurs because Li is strongly incorporated into actively precipitating secondary clay phases beyond what prior laboratory experiments have suggested. Chemical weathering incongruence is maintained despite relatively slow reaction rates and moderate clay accumulation due to a combination of two factors. First, reactive primary mineral phases persist across the weathering profile and effectively “shield” the secondary clays from resolubilization due to their greater solubility. Second, the clays accumulating in the near‐surface profile are relatively mature weathering byproducts. These factors promote characteristically low total dissolved solute export from the catchment despite significant input of exogenous dust. 
    more » « less
  2. Abstract. Watersheds are the fundamental Earth surface functioning units that connect the land to aquatic systems. Many watershed-scale models represent hydrological processes but not biogeochemical reactive transport processes. This has limited our capability to understand and predict solute export, water chemistry and quality, and Earth system response to changing climate and anthropogenic conditions. Here we present a recently developed BioRT-Flux-PIHM (BioRT hereafter) v1.0, a watershed-scale biogeochemical reactive transport model. The model augments the previously developed RT-Flux-PIHM that integrates land-surface interactions, surface hydrology, and abiotic geochemical reactions. It enables the simulation of (1) shallow and deep-water partitioning to represent surface runoff, shallow soil water, and deeper groundwater and of (2) biotic processes including plant uptake, soil respiration, and nutrient transformation. The reactive transport part of the code has been verified against the widely used reactive transport code CrunchTope. BioRT-Flux-PIHM v1.0 has recently been applied in multiple watersheds under diverse climate, vegetation, and geological conditions. This paper briefly introduces the governing equations and model structure with a focus on new aspects of the model. It also showcases one hydrology example that simulates shallow and deep-water interactions and two biogeochemical examples relevant to nitrate and dissolved organic carbon (DOC). These examples are illustrated in two simulation modes of complexity. One is the spatially lumped mode (i.e., two land cells connected by one river segment) that focuses on processes and average behavior of a watershed. Another is the spatially distributed mode (i.e., hundreds of cells) that includes details of topography, land cover, and soil properties. Whereas the spatially lumped mode represents averaged properties and processes and temporal variations, the spatially distributed mode can be used to understand the impacts of spatial structure and identify hot spots of biogeochemical reactions. The model can be used to mechanistically understand coupled hydrological and biogeochemical processes under gradients of climate, vegetation, geology, and land use conditions. 
    more » « less
  3. Abstract Chemical weathering in mountain critical zones controls river chemistry and regulates long‐term climate. Mountain landscapes contain diverse landforms created by geomorphic processes, including landslides, glacial moraines, and rock glaciers. These landforms generate unique flowpaths and water‐rock interactions that modify water chemistry as precipitation is transformed to streamflow. Variations in lithology and vegetation also strongly control water chemistry. Prior work has shown that landslides generate increased dissolved solute concentrations in rapidly uplifting mountains. However, there is still uncertainty regarding the magnitude which different geomorphic processes and land cover variations influence solute chemistry across tectonic and climatic regimes. We measured ion concentrations in surface water from areas that drain a variety of landforms and across land cover gradients in the East River watershed, a tributary of the Colorado River. Our results show that landslides produce higher solute concentrations than background values measured in streams draining soil‐mantled hillslopes and that elevated concentrations persist centuries to millennia after landslide occurrence. Channels with active bedrock incision also generate high solute concentrations, whereas solute concentrations in waters draining moraines and rock glaciers are comparable to background values. Solute fluxes from landslides and areas of bedrock incision are 1.6–1.8 times greater than nearby soil‐mantled hillslopes. Carbonic acid weathering dominates surface water samples from watersheds with greater vegetation coverage. Geomorphically enhanced weathering generates hotspots for net CO2release or sequestration, depending on lithology, that are 1.5–3.5 times greater than background values, which has implications for understanding links among surface processes, chemical weathering, and carbon cycle dynamics in alpine watersheds. 
    more » « less
  4. Climate models project changing patterns of precipitation and increases in temperature that modify soil moisture dynamics. Land use and changing frequency and intensity of precipitation can induce changes in soil structure and rooting abundances at timescales shorter than commonly considered. Soil structure is a critical ecosystem that governs water flow through soil profiles and across landscapes, and can influence weathering rates and thus solute release and soil development. We hypothesize that the altered soil structure and modification of rooting depth distributions linked to land use change can influence soil solute concentrations, and that those shifts in solute release are dependent on patterns of precipitation. We installed suction lysimeters to collect soil water for ~3 y in two grassland regions with distinct mean annual precipitation (800 mm y-1, 1100 mm y-1) in native prairie, agriculture, and post-agriculture land uses at depths of 10, 40, and 120 cm. We linked solute concentrations to soil moisture, aggregate-size distribution, pore geometry, and rooting depth distributions to assess how land use change and the altered rooting abundance it imposes can modify soil structure and hydrologic fluxes, and to infer how soil weathering can shift deep in the subsurface. We reveal how soil moisture residence time and the soil pore network can govern solute production, and the importance of precipitation and thus of soil moisture accumulation over growing seasons for mineral weathering and solute production. Specifically, we find that the solubility potential of multiple weathering products and organic carbon increases with precipitation, dominance of relatively small aggregates at the surface, and fewer coarse roots. Enhanced solute concentrations at depth may also reflect transport down-profile. Our findings reveal unintended consequences of land use change that influence important hydrologic dynamics and nutrient cycling in the vadose zone and how deeply and how persistently unexpected consequences of changes in land cover can propagate. 
    more » « less
  5. Abstract Projections of future conditions within the critical zone—earthcasts—can be used to understand the potential effects of changes in climate on processes affecting landscapes. We are developing an approach to earthcast how weathering will change in the future using scenarios of climate change. As a first step here, we use the earthcasting approach to model aspect‐related effects on soil water chemistry and weathering on hillsides in a well‐studied east‐west trending watershed (Shale Hills, Pennsylvania, USA). We completed model simulations of solute chemistry in soil water with and without the effect of aspect for comparison to catchment observations. With aspect included, aqueous weathering fluxes were higher on the sunny side of the catchment. But the effect of aspect on temperature (0.8 °C warmer soil on sunny side) and recharge (100 mm/year larger on shaded side) alone did not explain the magnitude of the observed higher weathering fluxes on the sunny side. Modeled aspect‐related differences in weathering fluxes only approach field observations when we incorporated the measured differences in clay content observed in augered soils on the two hillslopes. We also had to include a biolifting module to accurately describe cation concentrations in soil water versus depth. Biolifting lowered some mineral dissolution rates while accelerating kaolinite precipitation. These short‐duration simulations also highlighted that the inherited differences in particle size on the two sides of the catchment might in themselves be explained by weathering under different microclimates caused by aspect—over longer durations than simulated with our models. 
    more » « less