An official website of the United States government
Abstract We used seismic refraction to image the P‐wave velocity structure of a shale watershed experiencing regional compression in the Valley and Ridge Province (USA). From estimates showing strong compressional stress, we expected the depth to unweathered bedrock to mirror the hill‐valley‐hill topography (“bowtie pattern”) by analogy to seismic velocity patterns in crystalline bedrock in the North American Piedmont that also experience compression. Previous researchers used failure potentials calculated for strong compression in the Piedmont to suggest fractures are open deeper under hills than valleys to explain the “bowtie” pattern. Seismic images of the shale watershed, however, show little evidence of such a “bowtie.” Instead, they are consistent with weak (not strong) compression. This contradiction could be explained by the greater importance of infiltration‐driven weathering than fracturing in determining seismic velocities in shale compared to crystalline bedrock, or to local perturbations of the regional stress field due to lithology or structures.
more »
« less
3D Seismic Anatomy of a Watershed Reveals Climate‐Topography Coupling That Drives Water Flowpaths and Bedrock Weathering
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
- PAR ID:
- 10372470
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 126
- Issue:
- 12
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The rate of chemical weathering has been observed to increase with the rate of physical erosion in published comparisons of many catchments, but the mechanisms that couple these processes are not well understood. We investigated this question by exam- ining the chemical weathering and porosity profiles from catchments developed on marine shale located in Pennsylvania, USA (Susquehanna Shale Hills Critical Zone Observatory, SSHCZO); California, USA (Eel River Critical Zone Observatory, ERC- ZO); and Taiwan (Fushan Experimental Forest). The protolith compositions, protolith porosities, and the depths of regolith at these sites are roughly similar while the catchments are characterized by large differences in erosion rate (1–3 mm yr1 in Fushan 0.2–0.4 mm yr1 in ERCZO 0.01–0.025 mm yr1 in SSHCZO). The natural experiment did not totally isolate erosion as a variable: mean annual precipitation varied along the erosion gradient (4.2 m yr1 in Fushan > 1.9 m yr1 in ERCZO > 1.1 m yr1 in SSHCZO), so the fastest eroding site experiences nearly twice the mean annual temperature of the other two. Even though erosion rates varied by about 100, the depth of pyrite and carbonate depletion (defined here as regolith thickness) is roughly the same, consistent with chemical weathering of those minerals keeping up with erosion at the three sites. These minerals were always observed to be the deepest to react, and they reacted until 100% depletion. In two of three of the catchments where borehole observations were available for ridges, these minerals weathered across narrow reaction fronts. On the other hand, for the rock-forming clay mineral chlorite, the depth interval of weathering was wide and the extent of depletion observed at the land surface decreased with increasing erosion/precipitation. Thus, chemical weathering of the clay did not keep pace with erosion rate. But perhaps the biggest difference among the shales is that in the fast-eroding sites, microfractures account for 30–60% of the total porosity while in the slow-eroding shale, dissolution could be directly related to secondary porosity. We argue that the microfractures increase the influx of oxygen at depth and decrease the size of diffusion-limited internal domains of matrix, accelerating weathering of pyrite and carbonate under high erosion-rate condi- tions. Thus, microfracturing is a process that can couple physical erosion and chemical weathering in shales.more » « less
-
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
-
Abstract Projections of future conditions within the critical zone—earthcasts—can be used to understand the potential effects of changes in climate on processes affecting landscapes. We are developing an approach to earthcast how weathering will change in the future using scenarios of climate change. As a first step here, we use the earthcasting approach to model aspect‐related effects on soil water chemistry and weathering on hillsides in a well‐studied east‐west trending watershed (Shale Hills, Pennsylvania, USA). We completed model simulations of solute chemistry in soil water with and without the effect of aspect for comparison to catchment observations. With aspect included, aqueous weathering fluxes were higher on the sunny side of the catchment. But the effect of aspect on temperature (0.8 °C warmer soil on sunny side) and recharge (100 mm/year larger on shaded side) alone did not explain the magnitude of the observed higher weathering fluxes on the sunny side. Modeled aspect‐related differences in weathering fluxes only approach field observations when we incorporated the measured differences in clay content observed in augered soils on the two hillslopes. We also had to include a biolifting module to accurately describe cation concentrations in soil water versus depth. Biolifting lowered some mineral dissolution rates while accelerating kaolinite precipitation. These short‐duration simulations also highlighted that the inherited differences in particle size on the two sides of the catchment might in themselves be explained by weathering under different microclimates caused by aspect—over longer durations than simulated with our models.more » « less
-
Abstract. Mineral specific surface area (SSA) increases as primaryminerals weather and restructure into secondary phyllosilicate, oxide, andoxyhydroxide minerals. SSA is a measurable property that captures cumulativeeffects of many physical and chemical weathering processes in a singlemeasurement and has meaningful implications for many soil processes,including water-holding capacity and nutrient availability. Here we reportour measurements of SSA and mineralogy of two 21 m deep SSA profiles attwo landscape positions, in which the emergence of a very small mass percent(<0.1 %) of secondary oxide generated 36 %–81 % of the total SSAin both drill cores. The SSA transition occurred near 3 m at bothlocations and did not coincide with the boundary of soil to weathered rock. The3 m boundary in each weathering profile coincides with the depth extentof secondary iron oxide minerals and secondary phyllosilicates. Althoughelemental depletions in both profiles extend to 7 and 10 m depth, themineralogical changes did not result in SSA increase until 3 m depth. Theemergence of secondary oxide minerals at 3 m suggests that this boundary may bethe depth extent of oxidation weathering reactions. Our results suggest thatoxidation weathering reactions may be the primary limitation in thecoevolution of both secondary silicate and secondary oxide minerals. Wevalue element depletion profiles to understand weathering, but our findingof nested weathering fronts driven by different chemical processes (e.g.,oxidation to 3 m and acid dissolution to 10 m) warrants the recognition thatelement depletion profiles are not able to identify the full set ofprocesses that occur in weathering profiles.more » « less
