skip to main content

Title: Catchment‐Scale Architecture of the Deep Critical Zone Revealed by Seismic Imaging

Weathering and erosion processes are crucial to Critical Zone (CZ) evolution, landscape formation and availability of natural resources. Although many of these processes take place in the deep CZ (∼10–100 m), direct information about its architecture remain scarce. Near‐surface geophysics offers a cost‐effective and minimally intrusive alternative to drilling that provides information about the physical properties of the CZ. We propose a novel workflow combining seismic measurements, petrophysical modelling and geostatistical analysis to characterize the architecture of the deep CZ at the catchment scale, on the volcanic tropical island of Basse‐Terre (Guadeloupe, France). With this original workflow, we are for the first time able to jointly produce maps of the water table and the weathering front across an entire catchment, by means of a single geophysical method. This integrated view of the CZ highlights complex weathering patterns that call for going beyond “simple” hillslope CZ models.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Geophysical Research Letters
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The belowground architecture of the critical zone (CZ) consists of soil and rock in various stages of weathering and wetness that acts as a medium for biological growth, mediates chemical reactions, and controls partitioning of hydrologic fluxes. Hydrogeophysical imaging provides unique insights into the geometries and properties of earth materials that are present in the CZ and beyond the reach of direct observation beside sparse wellbores. An improved understanding of CZ architecture can be achieved by leveraging the geophysical measurements of the subsurface. Creating categorical models of the CZ is valuable for driving hydrologic models and comparing belowground architectures between different sites to interpret weathering processes. The CZ architecture is revealed through a novel comparison of hillslopes by applying facies classification in the elastic-electric domain driven by surface-based hydrogeophysical measurements. Three pairs of hillslopes grouped according to common geologic substrates — granite, volcanic extrusive, and glacially altered — are classified by five different hydrofacies classes to reveal the relative wetness and weathering states. The hydrofacies classifications are robust to the choice of initial mean values used in the classification and noncontemporaneous timing of geophysical data acquisition. These results will lead to improved interdisciplinary models of CZ processes at various scales and to an increased ability to predict the hydrologic timing and partitioning. Beyond the hillslope scale, this enhanced capability to compare CZ architecture can also be exploited at the catchment scale with implications for improved understanding of the link between rock weathering, hydrochemical fluxes, and landscape morphology. 
    more » « less
  2. Abstract. Trees, the most successful biological power plants on earth, build and plumb the critical zone (CZ) in ways that we do not yet understand. To encourage exploration of the character and implications of interactions between trees and soil in the CZ, we propose nine hypotheses that can be tested at diverse settings. The hypotheses are roughly divided into those about the architecture (building) and those about the water (plumbing) in the CZ, but the two functions are intertwined. Depending upon one's disciplinary background, many of the nine hypotheses listed below may appear obviously true or obviously false. (1) Tree roots can only physically penetrate and biogeochemically comminute the immobile substrate underlying mobile soil where that underlying substrate is fractured or pre-weathered. (2) In settings where the thickness of weathered material, H, is large, trees primarily shape the CZ through biogeochemical reactions within the rooting zone. (3) In forested uplands, the thickness of mobile soil, h, can evolve toward a steady state because of feedbacks related to root disruption and tree throw. (4) In settings where hH and the rates of uplift and erosion are low, the uptake of phosphorus into trees is buffered by the fine-grained fraction of the soil, and the ultimate source of this phosphorus is dust. (5) In settings of limited water availability, trees maintain the highest length density of functional roots at depths where water can be extracted over most of the growing season with the least amount of energy expenditure. (6) Trees grow the majority of their roots in the zone where the most growth-limiting resource is abundant, but they also grow roots at other depths to forage for other resources and to hydraulically redistribute those resources to depths where they can be taken up more efficiently. (7) Trees rely on matrix water in the unsaturated zone that at times may have an isotopic composition distinct from the gravity-drained water that transits from the hillslope to groundwater and streamflow. (8) Mycorrhizal fungi can use matrix water directly, but trees can only use this water by accessing it indirectly through the fungi. (9) Even trees growing well above the valley floor of a catchment can directly affect stream chemistry where changes in permeability near the rooting zone promote intermittent zones of water saturation and downslope flow of water to the stream. By testing these nine hypotheses, we will generate important new cross-disciplinary insights that advance CZ science.

    more » « less
  3. Abstract

    Lithium isotopes are used to trace weathering intensity, but little is known about the processes that fractionate them in highly weathered settings, where secondary minerals play a dominant role in weathering reactions. To help fill this gap in our knowledge of Li isotope systematics, we investigated Li isotope fractionation at an andesitic catchment in Puerto Rico, where the highest rates of silicate weathering on Earth have been documented. We found the lowest δ7Li values published to date for porewater (−27‰) and bulk regolith (−38‰), representing apparent fractionations relative to parent rock of −31‰ and −42‰, respectively. We also found δ7Li values that are lower in the exchangeable fraction than in the bulk regolith or porewater, the opposite than expected from secondary mineral precipitation. We interpret these large isotopic offsets and the unusual relationships between Li pools as resulting from two distinct weathering processes at different depths in the regolith. At the bedrock‐regolith transition (9.3–8.5 m depth), secondary mineral precipitation preferentially retains the lighter6Li isotope. These minerals then dissolve further up the profile, leaching6Li from the bulk solid, with a total variation of about +50‰withinthe profile, attributable primarily to clay dissolution. Importantly, streamwater δ7Li (about +35‰) is divorced entirely from these regolith weathering processes, instead reflecting deeper weathering reactions (>9.3 m). Our work thus shows that the δ7Li of waters draining highly weathered catchments may reflect bedrock mineralogy and hydrology, rather than weathering intensity in the regolith covering the catchment.

    more » « less
  4. Abstract

    The critical zone (CZ) is the region of the Earth’s surface that extends from the bottom of the weathered bedrock to the tree canopy and is important because of its ability to store water and support ecosystems. A growing number of studies use active source shallow seismic refraction to explore and define the size and structure of the CZ across landscapes. However, measurement uncertainty and model resolution at depth are generally not evaluated, which makes the identification and interpretation of CZ features inconclusive. To reliably resolve seismic velocity with depth, we implement a Transdimensional Hierarchical Bayesian (THB) framework with reversible‐jump Markov Chain Monte Carlo to generate samples from the posterior distribution of velocity structures. We also perform 2D synthetic tests to explore how well THB traveltime inversion can resolve different subsurface velocity structures. We find that THB recovers both sharp changes in velocity as well as gradual velocity increases with depth. Furthermore, we explore the velocity structure in a series of ridge‐valley systems in northern California. The posterior velocity model shows an increasing thickness of low velocity material from channels to ridgetops along a transect parallel to bedding strike, implying a deeper weathering zone below ridgetops and hillslopes than below channels. The THB method enhances the ability to reliably image CZ structure, and the model uncertainty estimates it yields provides an objective way to interpret deep CZ structure. The method can be applied across other near‐surface studies, especially in the presence of significant surface topography.

    more » « less
  5. Abstract

    The structure of the critical zone (CZ) is a product of feedbacks among hydrologic, climatic, biotic, and chemical processes. Past research within snow‐dominated systems has shown that aspect‐dependent solar radiation inputs can produce striking differences in vegetation composition, topography, and soil depth between opposing hillslopes. However, far fewer studies have evaluated the role of microclimates on CZ development within rain‐dominated systems, especially below the soil and into weathered bedrock. To address this need, we characterized the CZ of a north‐facing and south‐facing slope within a first‐order headwater catchment located in central coast California. We combined terrain analysis of vegetation distribution and topography with soil pit characterization, geophysical surveys and hydrologic measurements between slope‐aspects. We documented denser vegetation and higher shallow soil moisture on north facing slopes, which matched previously documented observations in snow‐dominated sites. However, average topographic gradients were 24° and saprolite thickness was approximately 6 m across both hillslopes, which did not match common observations from the literature that showed widespread asymmetry in snow‐dominated systems. These results suggest that dominant processes for CZ evolution are not necessarily transferable across regions. Thus, there is a continued need to expand CZ research, especially in rain‐dominated and water‐limited systems. Here, we present two non‐exclusive mechanistic hypotheses that may explain these unexpected similarities in slope and saprolite thickness between hillslopes with opposing aspects.

    more » « less