An official website of the United States government
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
A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport
Abstract 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):
- 1654194
- PAR ID:
- 10453984
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 33
- Issue:
- 4
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- p. 466-475
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
Volcanic provinces are among the most active but least well understood landscapes on Earth. Here, we show that the central Cascade arc, USA, exhibits systematic spatial covariation of topography and hydrology that are linked to aging volcanic bedrock, suggesting systematic controls on landscape evolution. At the Cascade crest, a locus of Quaternary volcanism, water circulates deeply through the upper 1 km of crust but transitions to shallow and dominantly horizontal flow as rocks age away from the arc front. We argue that this spatial pattern reflects a temporal state shift in the deep Critical Zone. Chemical weathering at depth, surface particulate deposition, and tectonic forcing drive landscapes away from an initial state with minimal topographic dissection, large vertical hydraulic conductivity, abundant lakes, and muted hydrographs toward a state of deep fluvial dissection, small vertical hydraulic conductivity, few lakes, and flashy hydrographs. This state shift has major implications for regional water resources. Drill hole temperature profiles imply at least km of active groundwater currently stored at the Cascade Range crest, with discharge variability a strong function of bedrock age. Deeply circulating groundwater also impacts volcanism, and Holocene High Cascades eruptions reflect explosive magma–water interactions that increase regional volcanic hazard potential. We propose that a Critical Zone state shift drives volcanic landscape evolution in wet climates and represents a framework for understanding interconnected solid earth dynamics and climate in these terrains.more » « less
-
The depth to unweathered bedrock beneath landscapes influences subsurface runoff paths, erosional processes, moisture availability to biota, and water flux to the atmosphere. Here we propose a quantitative model to predict the vertical extent of weathered rock underlying soil-mantled hillslopes. We hypothesize that once fresh bedrock, saturated with nearly stagnant fluid, is advected into the near surface through uplift and erosion, channel incision produces a lateral head gradient within the fresh bedrock inducing drainage toward the channel. Drainage of the fresh bedrock causes weathering through drying and permits the introduction of atmospheric and biotically controlled acids and oxidants such that the boundary between weathered and unweathered bedrock is set by the uppermost elevation of undrained fresh bedrock, Z b . The slow drainage of fresh bedrock exerts a “bottom up” control on the advance of the weathering front. The thickness of the weathered zone is calculated as the difference between the predicted topographic surface profile (driven by erosion) and the predicted groundwater profile (driven by drainage of fresh bedrock). For the steady-state, soil-mantled case, a coupled analytical solution arises in which both profiles are driven by channel incision. The model predicts a thickening of the weathered zone upslope and, consequently, a progressive upslope increase in the residence time of bedrock in the weathered zone. Two nondimensional numbers corresponding to the mean hillslope gradient and mean groundwater-table gradient emerge and their ratio defines the proportion of the hillslope relief that is unweathered. Field data from three field sites are consistent with model predictions.more » « less
-
We investigated how runoff-to-groundwater partitioning changes as a function of substrate age and degree of regolith development in the Island of Hawai’i, by modeling watershed-scale hydrodynamic properties for a series of volcanic catchments of different substrate age developed under different climates. In the younger catchments, rainfall infiltrates directly into the groundwater system and surface runoff is minimal, consisting of ephemeral streams flowing on the scale of hours to days. The older catchments show increasing surface runoff, with deeper incision and perennial discharge. We hypothesize that watershed-scale hydrodynamic properties change as a function of their weathering history—the convolution of time and climate: as surfaces age and become increasingly weathered, hydraulic conductivity is reduced, leading to increased runoff-to-recharge ratios. To test this relationship, we calculated both saturated hydraulic conductivity (k) and aquifer thickness (D) using recession flow analysis. We show that the average k in the younger catchments can be between 3 to 6 orders of magnitude larger than in older catchments, whereas modeled D increases with age. Ephemeral streams with zero baseflow at daily timescales cannot be evaluated using the same method. Instead, we calculated the recession constant for two contiguous catchments developed on young ash or lava deposits of different ages. Increasing bedrock age results in slower recession response in these ephemeral streams, which is consistent with decreasing hydraulic conductivity. Our results highlight the role of the weathering history in determining the evolution of watershed-scale hydrologic properties in volcanic catchments.more » « less
