The porous near‐surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long‐term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope‐aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north‐facing slope, regolith on the south‐facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 John Wiley & Sons, Ltd.
This content will become publicly available on July 1, 2024
Understanding how soil thickness and bedrock weathering vary across ridge and valley topography is needed to constrain the flowpaths of water and sediment production within a landscape. Here, we investigate saprolite and weathered bedrock properties across a ridge‐valley system in the Northern California Coast Ranges, USA, where topography varies with slope aspect such that north‐facing slopes have thicker soils and are more densely vegetated than south‐facing slopes. We use active source seismic refraction surveys to extend observations made in boreholes to the hillslope scale. Seismic velocity models across several ridges capture a high velocity gradient zone (from 1,000 to 2,500 m/s) located ∼4–13 m below ridgetops that coincides with transitions in material strength and chemical depletion observed in boreholes. Comparing this transition depth across multiple north‐ and south‐facing slopes, we find that the thickness of saprolite does not vary with slope aspects. Additionally, seismic survey lines perpendicular and parallel to bedding planes reveal weathering profiles that thicken upslope and taper downslope to channels. Using a rock physics model incorporating seismic velocity, we estimate the total porosity of the saprolite and find that inherited fractures contribute a substantial amount of pore space in the upper 6 m, and the lateral porosity structure varies strongly with hillslope position. The aspect‐independent weathering structure suggests that the contemporary critical zone structure at Rancho Venada is a legacy of past climate and vegetation conditions.
more » « less- Award ID(s):
- 2012616
- NSF-PAR ID:
- 10469131
- Publisher / Repository:
- AGU
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 128
- Issue:
- 7
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
Abstract -
more » « less
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
Abstract To investigate how bedrock transforms to soil, we mapped the topography of the interface demarcating onset of weathering under an east‐west trending shale watershed in the Valley and Ridge province in the USA Using wave equation travel‐time tomography from a seismic array of >4,000 geophones, we obtained a 3D P‐wave velocity (Vp) model that resolves structures ∼20 m below land surface (mbls). The depth of mobile soil and the onset of dissolution of chlorite roughly match Vp = 600 m/s and Vp = 2,700 m/s, respectively. Chlorite dissolution initiates porosity growth in the shale matrix. Depth to the 2,700 m/s contour is greater under the N‐ as compared to S‐facing hillslopes and under sub‐planar as compared to concave‐up land surfaces. Broadly, the geometries of the ‘soil’ and ‘chlorite’ Vp contours are consistent with the calculated potential for shear fracture opening under weak regional compression. However, this calculated fracture potential does not consistently explain observations related to N‐ versus S‐facing aspect nor fracture density observed by borehole televiewer. Apparently, regional compression is only a secondary influence on Vp: the primary driver of P‐wave slowing in the upper layers of this catchment is topographic control of reactive water flowpaths and their integrated effects on weathering. The Vp result is best explained as the long‐term integrated effect of groundwater flow‐induced geochemical weathering of shale in response to climate‐driven patterns of micro‐ and macro‐topography.
-
more » « less
Abstract Shallow bedrock strength controls both landslide hazard and the rate and form of erosion, yet regional patterns in near‐surface mechanical properties are rarely known quantitatively due to the challenge in collecting
in situ measurements. Here we present seismic and geomechanical characterizations of the shallow subsurface across the central Himalayan Range in Nepal. By pairing widely distributed 1D shear wave velocity surveys and engineering outcrop descriptions per the Geological Strength Index classification system, we evaluate landscape‐scale patterns in near‐surface mechanical characteristics and their relation to environmental factors known to affect rock strength. We find that shallow bedrock strength is more dependent on the degree of chemical and physical weathering, rather than the mineral and textural differences between the metamorphic lithologies found in the central Himalaya. Furthermore, weathering varies systematically with topography. Bedrock ridge top sites are highly weathered and have S‐wave seismic velocities and shear strength characteristics that are more typical of soils, whereas sites near valley bottoms tend to be less weathered and characterized by high S‐wave velocities and shear strength estimates typical of rock. Weathering on hillslopes is significantly more variable, resulting in S‐wave velocities that range between the ridge and channel endmembers. We speculate that variability in the hillslope environment may be partly explained by the episodic nature of mass wasting, which clears away weathered material where landslide scars are recent. These results underscore the mechanical heterogeneity in the shallow subsurface and highlight the need to account for variable bedrock weathering when estimating strength parameters for regional landslide hazard analysis. -
more » « less
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.
