An official website of the United States government
Abstract Knowing little about how porosity and permeability are distributed at depth, we commonly develop models of groundwater by treating the subsurface as a homogeneous black box even though porosity and permeability vary with depth. One reason for this depth variation is that infiltrating meteoric water reacts with minerals to affect porosity in localized zones called reaction fronts. We are beginning to learn to map and model these fronts beneath headwater catchments and show how they are distributed. The subsurface landscapes defined by these fronts lie subparallel to the soil‐air interface but with lower relief. They can be situated above, below, or at the water table. These subsurface landscapes of reaction are important because porosity developed from weathering can control subsurface water storage. In addition, porosity often changes at the weathering fronts, and when this affects permeability significantly, the front can act like a valve that re‐orients water flowing through the subsurface. We explore controls on the positions of reaction fronts under headwater landscapes by accounting for the timescales of erosion, chemical equilibration, and solute transport. One strong control on the landscape of subsurface reaction is the land surface geometry, which is in turn a function of the erosion rate. In addition, the reaction fronts, like the water table, are strongly affected by the lithology and water infiltration rate. We hypothesize that relationships among the land surface, reaction fronts, and the water table are controlled by feedbacks that can push landscapes towards an ‘ideal hill’. In this steady state, reaction‐front valves partition water volumes into shallow and deep flowpaths. These flows dissolve low‐ and high‐solubility minerals, respectively, allowing their reaction fronts to advance at the erosion rate. This conceptualization could inform better models of subsurface porosity and permeability, replacing the black box.
more »
« less
Climate driven coevolution of weathering profiles and hillslope topography generates dramatic differences in critical zone architecture
Abstract Considerable debate revolves around the relative importance of rock type, tectonics, and climate in creating the architecture of the critical zone. We demonstrate the importance of climate and in particular the rate of water recharge to the subsurface, using numerical models that incorporate hydrologic flowpaths, chemical weathering, and geomorphic rules for soil production and transport. We track alterations in both solid phase (plagioclase to clay) and water chemistry along hydrologic flowpaths that include lateral flow beneath the water table. To isolate the role of recharge, we simulate dry and wet cases and prescribe identical landscape evolution rules. The weathering patterns that develop differ dramatically beneath the resulting parabolic interfluves. In the dry case, incomplete weathering is shallow and surface parallel, whereas in the wet case, intense weathering occurs to depths approximating the base of the bounding channels, well below the water table. Exploration of intermediate cases reveals that the weathering state of the subsurface is strongly governed by the ratio of the rate of advance of the weathering front itself controlled by the water input rate, and the rate of erosion of the landscape. The system transitions between these end‐member behaviours rather abruptly at a weathering front speed ‐ erosion rate ratio of approximately 1. Although there are undoubtedly direct roles for tectonics and rock type in critical zone architecture, and yet more likely feedbacks between these and climate, we show here that differences in hillslope‐scale weathering patterns can be strongly controlled by climate.
more »
« less
- Award ID(s):
- 1331828
- PAR ID:
- 10573832
- Publisher / Repository:
- John Wiley & Sons
- Date Published:
- Journal Name:
- Hydrological Processes
- Volume:
- 33
- Issue:
- 1
- ISSN:
- 0885-6087
- Page Range / eLocation ID:
- 4 to 19
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The creation of fractures in bedrock dictates water movement through the critical zone, controlling weathering, vadose zone water storage, and groundwater recharge. However, quantifying connections between fracturing, water flow, and chemical weathering remains challenging because of limited access to the deep critical zone. Here we overcome this challenge by coupling measurements from borehole drilling, groundwater monitoring, and seismic refraction surveys in the central California Coast Range. Our results show that the subsurface is highly fractured, which may be driven by the regional geologic and tectonic setting. The pervasively fractured rock facilitates infiltration of meteoric water down to a water table that aligns with oxidation in exhumed rock cores and is coincident with the adjacent intermittent first‐order stream channel. This work highlights the need to incorporate deep water flow and weathering due to pervasive fracturing into models of catchment water balances and critical zone weathering, especially in tectonically active landscapes.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
-
Chemical weathering influences many aspects of the Earth system, including biogeochemical cycling, climate, and ecosystem function. Physical erosion influences chemical weathering rates by setting the supply of fresh minerals to the critical zone. Vegetation also influences chemical weathering rates, both by physical processes that expose mineral surfaces and via production of acids that contribute to mineral dissolution. However, the role of vegetation in setting surface process rates in different landscapes is unclear. Here we use 10Be and geochemical mass balance to quantify soil production, physical erosion, and chemical weathering rates in a landscape where a migrating drainage divide separates catchments with an order-of magnitude contrast in erosion rates and where vegetation spans temperate rainforest, tussock grassland, and unvegetated alpine ecosystems in the western Southern Alps of New Zealand. Soil production, physical erosion, and chemical weathering rates are significantly higher on the rapidly eroding versus the slowly eroding side of the drainage divide. However, chemical weathering intensity does not vary significantly across the divide or as a function of vegetation type. Soil production rates are correlated with ridgetop curvature, and ridgetops are more convex on the rapidly eroding side of the divide, where soil mineral residence times are lowest. Hence our findings suggest fluvially-driven erosion rates control soil production and soil chemical weathering rates by influencing the relationship between hillslope topography and mineral residence times. In the western Southern Alps, soil production and chemical weathering rates are more strongly mediated by physical rock breakdown driven by landscape response to tectonics, than by vegetation.more » « less
-
Abstract The hydrologic dynamics and geomorphic evolution of watersheds are intimately coupled—runoff generation and water storage are controlled by topography and properties of the surface and subsurface, while also affecting the evolution of those properties over geologic time. However, the large disparity between their timescales has made it difficult to examine interdependent controls on emergent hydrogeomorphic properties, such as hillslope length, drainage density, and extent of surface saturation. In this study, we develop a new model coupling hydrology and landscape evolution to explore how runoff generation affects long‐term catchment evolution, and analyze numerical results using a nondimensional scaling framework. We focus on hydrologic processes dominating in humid climates where storm runoff primarily arises from shallow subsurface flow and from precipitation on saturated areas. The model solves hydraulic groundwater equations to predict the water‐table elevation given prescribed, constant groundwater recharge. Water in excess of the subsurface capacity for transport becomes overland flow, which generates shear stress on the surface and may detach and transport sediment. This affects the landscape form that in turn affects runoff generation. We show that (a) four dimensionless parameters describe the possible steady state landscapes that coevolve under steady recharge; (b) hillslope length increases with increasing transmissivity relative to the recharge rate; (c) three topographic metrics—steepness index, Laplacian curvature, and topographic index—together provide a basis for interpreting landscapes that have coevolved with runoff generated via shallow subsurface flow. Finally we discuss the possibilities and limitations for quantitative comparisons between the model results and real landscapes.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1002/hyp.13307
