skip to main content

Title: A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimensionxalmost disappear when the vertical dimensionzis expressed as a fractionz/H(x) of the laterally flowing system thicknessH(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.

more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Hydrological Processes
Page Range / eLocation ID:
p. 466-475
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.

    more » « less
  2. Abstract

    Water is known to play an essential role in initiating and maintaining subsurface weathering reactions. However, the interaction between the weathering front and the water table is unclear and under intense debate. Here, we present a high‐fidelity, 3D image of a variably saturated weathering front beneath a granite terrain in the Laramie Range, Wyoming, constructed using full‐3D ambient‐noise adjoint tomography and calibrated with data from an extensive drilling and hydraulic well testing effort. The imaged weathering front between saprolite and weathered bedrock is overall shallower than the water table under ridge but deeper than water table under valleys. We propose that downward‐advancing weathering front coevolves with water table in a positive cycle that gradually flattens the water table, enhances the rate of groundwater drainage, and exposes underlying bedrock to weathering. As a result, we expect this cycle to become “sluggish” with time as water table gradient decreases.

    more » « less
  3. Abstract

    Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first‐order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated‐soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.

    more » « less
  4. Abstract

    Deeper flows through bedrock in mountain watersheds could be important, but lack of data to characterize bedrock properties limits understanding. To address data scarcity, we combine a previously published integrated hydrologic model of a snow‐dominated, headwater basin of the Colorado River with a new method for dating baseflow age using dissolved gas tracers SF6, CFC‐113, N2, and Ar. The original flow model predicts the majority of groundwater flow through shallow alluvium (<8 m) sitting on top of less permeable bedrock. The water moves too quickly and is unable to reproduce observed SF6concentrations. To match gas data, bedrock permeability is increased to allow a larger fraction of deeper and older groundwater flow (median 112 m). The updated hydrologic model indicates interannual variability in baseflow age (3–12 years) is controlled by the volume of seasonal interflow and tightly coupled to snow accumulation and monsoon rain. Deeper groundwater flow remains stable (11.7 ± 0.7 years) as a function mean historical recharge to bedrock hydraulic conductivity (R/K). A sensitivity analysis suggests that increasing bedrock K effectively moves this alpine basin away from its original conceptualization of hyperlocalized groundwater flow (high R/K) with groundwater age insensitive to changes in water inputs. Instead, this basin is situated close to the precipitation threshold defining recharge controlled groundwater flow conditions (low R/K) in which groundwater age increases with small reductions in precipitation. Work stresses the need to explore alternative methods characterizing bedrock properties in mountain basins to better quantify deeper groundwater flow and predict their hydrologic response to change.

    more » « less
  5. Mineral weathering is an important soil-forming process driven by the interplay of water, organisms, solution chemistry, and mineralogy. The influence of hillslope-scale patterns of water flux on mineral weathering in soils is still not well understood, particularly in humid postglacial soils, which commonly harbor abundant weath- erable primary minerals. Previous work in these settings showed the importance of lateral hydrologic patterns to hillslope-scale pedogenesis. In this study, we hypothesized that there is a corresponding relationship between hydrologically driven pedogenesis and chemical weathering in podzols in the White Mountains of New Hamp- shire, USA. We tested this hypothesis by quantifying the depletion of plagioclase in the fine fraction (≤2 mm) of closely spaced, similar-age podzols along a gradient in topography and depth to bedrock that controls lateral water flow. Along this gradient, laterally developed podzols formed through frequent, episodic flushing by up- slope groundwater, and vertically developed podzols formed through characteristic vertical infiltration. We estimated the depletion of plagioclase-bound elements within the upper mineral horizons of podzols using mass transfer coefficients (τ) and quantified plagioclase losses directly through electron microscopy and microprobe analysis. Elemental depletion was significantly more pronounced in the upslope lateral eluvial (E horizon- dominant) podzols relative to lateral illuvial (B horizon-dominant) and vertical (containing both E and B hori- zons) podzols downslope, with median Na losses of ~74 %, ~56 %, and ~40 %, respectively. When comparing genetic E horizons, Na and Al were significantly more depleted in laterally developed podzols relative to vertically developed podzols. Microprobe analysis revealed that ~74 % of the plagioclase was weathered from the mineral pool of lateral eluvial podzols, compared to ~39 % and ~23 % for lateral illuvial podzols and vertically developed podzols, respectively. Despite this intense weathering, plagioclase remains the second most abundant mineral in soil thin sections. These findings confirm that the concept of soil development as occurring vertically does not accurately characterize soils in topographically complex regions. Our work improves the current understanding of pedogenesis by identifying distinct, short-scale gradients in mineral weathering shaped by local patterns of hydrology and topography. 
    more » « less