skip to main content


Title: Quantifying Depth‐Dependent Seismic Anisotropy in the Critical Zone Enhanced by Weathering of a Piedmont Schist
Abstract

Weathering processes weaken and break apart rock, freeing nutrients and enhancing permeability through the subsurface. To better understand these processes, it is useful to constrain physical properties of materials derived from weathering within the critical zone. Foliated rocks exhibit permeability, strength and seismic anisotropy–the former two bear hydrological and geomorphological consequences while the latter is geophysically quantifiable. Each of these types of anisotropy are related to rock fabric (fractures and foliation); thus, characterizing weathering‐dependent changes in rock fabric with depth may have a range of implications (e.g., landslide susceptibility, groundwater modeling, and landscape evolution). To better understand how weathering effects rock fabric, we quantify seismic anisotropy in saprolite and weathered bedrock within two catchments underlain by the Precambrian Loch Raven schist, located in Oregon Ridge Park, MD. Using circular geophone arrays and perpendicular seismic refraction profiles, anisotropy versus depth functions are created for material 0–25 m below ground surface (bgs). We find that anisotropy is relatively low (0%–15%) in the deepest material sampled (12–25 m bgs) but becomes more pronounced (29%–33%) at depths corresponding with saprolite and highly weathered bedrock (5–12 m bgs). At shallow soil depths (0–5 m bgs), material is seismically isotropic, indicating that mixing processes have destroyed parent fabric. Therefore, in situ weathering and anisotropy appear to be correlated, suggesting that in‐place weathering amplifies the intrinsic anisotropy of bedrock.

 
more » « less
Award ID(s):
1945431 2012073 1654194 2012353 2012316
NSF-PAR ID:
10360063
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Volume:
126
Issue:
10
ISSN:
2169-9003
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    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
  2. Abstract

    As bedrock weathers to regolith – defined here as weathered rock, saprolite, and soil – porosity grows, guides fluid flow, and liberates nutrients from minerals. Though vital to terrestrial life, the processes that transform bedrock into soil are poorly understood, especially in deep regolith, where direct observations are difficult. A 65-m-deep borehole in the Calhoun Critical Zone Observatory, South Carolina, provides unusual access to a complete weathering profile in an Appalachian granitoid. Co-located geophysical and geochemical datasets in the borehole show a remarkably consistent picture of linked chemical and physical weathering processes, acting over a 38-m-thick regolith divided into three layers: soil; porous, highly weathered saprolite; and weathered, fractured bedrock. The data document that major minerals (plagioclase and biotite) commence to weather at 38 m depth, 20 m below the base of saprolite, in deep, weathered rock where physical, chemical and optical properties abruptly change. The transition from saprolite to weathered bedrock is more gradational, over a depth range of 11–18 m. Chemical weathering increases steadily upward in the weathered bedrock, with intervals of more intense weathering along fractures, documenting the combined influence of time, reactive fluid transport, and the opening of fractures as rock is exhumed and transformed near Earth’s surface.

     
    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. Fractures in Earth's critical zone influence groundwater flow and storage and promote chemical weathering. Fractured materials are difficult to characterize on large spatial scales because they contain fractures that span a range of sizes, have complex spatial distributions, and are often inaccessible. Therefore, geophysical characterizations of the critical zone depend on the scale of measurements and on the response of the medium to impulses at that scale. Using P-wave velocities collected at two scales, we show that seismic velocities in the fractured bedrock layer of the critical zone are scale-dependent. The smaller-scale velocities, derived from sonic logs with a dominant wavelength of ~0.3 m, show substantial vertical and lateral heterogeneity in the fractured rock, with sonic velocities varying by 2,000 m/s over short lateral distances (~20 m), indicating strong spatial variations in fracture density. In contrast, the larger-scale velocities, derived from seismic refraction surveys with a dominant wavelength of ~50 m, are notably slower than the sonic velocities (a difference of ~3,000 m/s) and lack lateral heterogeneity. We show that this discrepancy is a consequence of contrasting measurement scales between the two methods; in other words, the contrast is not an artifact but rather information—the signature of a fractured medium (weathered/fractured bedrock) when probed at vastly different scales. We explore the sample volumes of each measurement and show that surface refraction velocities provide reliable estimates of critical zone thickness but are relatively insensitive to lateral changes in fracture density at scales of a few tens of meters. At depth, converging refraction and sonic velocities likely indicate the top of unweathered bedrock, indicative of material with similar fracture density across scales. 
    more » « less
  5. Subduction of the very young (<15 Myr old) oceanic lithosphere of the Nazca plate in central to southern Colombia is observationally related to an unusually high and unusually variable amount of intermediate (>50 km) depth seismicity. From 2010 through 2019 89% of central and southern Colombia’s 11,466 intermediate depth events occurred between 3.5°N and 5.5°N, highlighting these unusual characteristics of the young slab. In addition, morphologic complexity and possible tears characterize the Nazca slab in Colombia and complicate mantle flow in the region. Prior SKS-phase shear-wave splitting results indicate sub-slab anisotropy is dominated by plate motion parallel-to-subparallel orientations in the region, suggesting the young slab has entrained a relatively thick portion of the sub-slab mantle. These observations suggest the subduction of young lithosphere has significant effects on both the overlying and underlying asthenosphere in the Colombia subduction zone. Here we use more than 10 years of data to calculate receiver functions for the Red Sismológica Nacional de Colombia’s network of broadband seismometers. These receiver functions allow us to tie these prior observations of the Colombia subduction zone to distinct, structural features of the slab. We find that the region of high seismicity corresponds to a low seismic velocity feature along the top of the subducting plate between 3.5°N and 5.5°N that is not present to the south. Moderately elevated P-wave velocity to S-wave velocity ratios are also observed within the slab in the north. This feature likely represents hydrated slab mantle and/or uneclogitized oceanic crust extending to a deeper depth in the north of the region which may provide fluids to drive slab seismicity. We further find evidence for a thick layer of material along the slab’s lithosphere-asthenosphere boundary characterized by spatially variable anisotropy. This feature likely represents entrained asthenosphere at the base of the plate sheared by both the overlying plate and complex flow related to proposed slab tears just north and south of the study region. These observations highlight how structural observations provide key contextual constraints on short-term (seismogenic) and long-term (anisotropic fabric) dynamic processes in the Colombia subduction zone. Plain-language Summary The Nazca oceanic plate is very young (<15 million years old) where it is pulled or subducted beneath the South America plate in central and southern Colombia. Earthquakes occurring in the subducted Nazca plate at depths greater than 50 km are nearly 9x more common in central Colombia than in southern Colombia. The subducted Nazca plate also has a complex shape in this region and may have been torn both in northern Colombia and to the south near the Colombia-Ecuador border. The slow flow of mantle rock beneath the subducted plate is believed to be affected by this and earlier studies have inferred this flow is mostly in the same direction as the subducting plate's motion. We have used 10+ years of data to calculate receiver functions, which can detect changes in the velocity of seismic waves at the top and bottom of the subducted plate to investigate these features. We found that the Nazca plate is either hydrated or has rocks with lower seismic velocities at its top in the central part of Colombia where earthquakes are common. We also find that a thick layer of mantle rock at the base of the subducted plate has been sheared. 
    more » « less