skip to main content

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

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
Author(s) / Creator(s):
 ;  ;  ;  ;  ;  ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract

    We present a statistical model of soil and rock weathering in deep profiles to expand the capacity to assess weathering to heterogeneous bedrock types, which are common at the Earth's surface. We developed the Weathering Trends (WT) model by extending the fractional mass change calculation (tau) of the geochemical mass balance model in two important ways. First, WT log transforms the elemental ratio data, to discern the log‐linear patterns that naturally develop from thermodynamic and kinetic laws of chemistry. Second, WT statistically fits log‐transformed element concentration ratio data – log(cj/ci), the only depth‐varying term in tau – as a function of depth to determine characteristic depths of transitions in weathering processes, along with confidence intervals. With no prior assumptions, WT estimates average parent material composition, average composition of the upper weathered zone and mean fractional mass change of each element over the entire weathering profile. WT displays the mean shape of weathering profiles of log‐transformed geochemical data bounded by calculated confidence intervals. We share the WT model code as an open‐source R package ( The WT model was designed to interpret two 21 m cores from the Laurels Schist bedrock in the Christina River Basin Critical Zone Observatory in the Pennsylvania Piedmont, where our morphological and elemental data provided inconclusive estimates of bedrock depth. The WT model differentiated between rock variability and weathering to delineate the maximum extent of weathering at 12.3 m (CI 95% [9.2, 21.3]) in Ridge Well 1 and 7.2 m (CI 95% [4.3, 13.0]) in Interfluve Well 2. The water table was 5–8 m below fresh rock at Ridge Well 1, but at the same depth as fresh rock at the lower elevation interfluve. We assess statistical approaches to identify the best immobile element for use in WT and tau calculations. Copyright © 2017 John Wiley & Sons, Ltd.

    more » « less
  4. 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
  5. 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 collectingin situmeasurements. 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