Seulgi Moon

Geophysical imaging reveals topographic stress control of bedrock weathering.

Bedrock weathering runs to the hills Fractures in bedrock drive the breakdown of rock into soil. Soil makes observations of bedrock processes challenging. St. Clair et al. combined a three-dimensional stress model with geophysical measurements to show that bedrock erosion rates mirror changes in topography (see the Perspective by Anderson). Seismic reflection and electromagnetic profiles allowed mapping of the bedrock fracture density. The profiles mirror changes in surface elevation and thus provide a way to study the critical zone between rock and soil. Science, this issue p. 534; see also p. 506 Geophysical survey data and stress modeling connect surface topography to Earth’s critical zone. [Also see Perspective by Anderson] Bedrock fracture systems facilitate weathering, allowing fresh mineral surfaces to interact with corrosive waters and biota from Earth’s surface, while simultaneously promoting drainage of chemically equilibrated fluids. We show that topographic perturbations to regional stress fields explain bedrock fracture distributions, as revealed by seismic velocity and electrical resistivity surveys from three landscapes. The base of the fracture-rich zone mirrors surface topography where the ratio of horizontal compressive tectonic stresses to near-surface gravitational stresses is relatively large, and it parallels the surface topography where the ratio is relatively small. Three-dimensional stress calculations predict these results, suggesting that tectonic stresses interact with topography to influence bedrock disaggregation, groundwater flow, chemical weathering, and the depth of the “critical zone” in which many biogeochemical processes occur.

The chemical, mechanical, and hydrological evolution of weathering granitoid

Surprisingly few studies connect the chemical, mechanical, and hydrological evolution of rock as it weathers to saprolite and soil. We assess this coevolution in granodiorite from Monterey Peninsula, California, by measuring changes in bulk chemistry, mineralogy, volumetric strain, the oxidation state of Fe in biotite crystals, tensile strength, abrasion rate, connected porosity, and hydraulic conductivity in samples covering a range of weathering grades. We identify the oxidative dissolution of biotite as the key chemical reaction because of the volumetric expansion that accompanies formation of altered biotite and precipitation of ferrihydrite. We show how the associated accumulation of elastic strain produces an energy density that is sufficient to support rock fracturing over length scales equivalent to constituent crystals. The resulting intragranular and intergranular cracking profoundly reduces tensile strength and increases the abrasion rate, connected porosity, and hydraulic conductivity of the rock matrix. These changes increase the rate of plagioclase weathering, and ultimately the rock disintegrates into grus and clay. Major changes in rock properties can occur with only minor element leaching, and the threshold behavior of weathering that arises from the coevolution of chemical, hydrological, and mechanical properties may be difficult to capture using simplified weathering models that fail to incorporate these properties. Our results, which combine the mechanical and hydrological evolution of weathering rock with more common measurements of chemical changes, should help to more accurately model the effects of, and mechanical and hydrological feedbacks upon, chemical weathering of rock.

Critical Zone Structure Under a Granite Ridge Inferred From Drilling and Three‐Dimensional Seismic Refraction Data

1Department of Geology and Geophysics, University of Wyoming, Laramie, Wyoming, USA 2Wyoming Center for Environmental Hydrology and Geophysics, University of Wyoming, Laramie, Wyoming, USA 3Department of Earth Sciences, Dickinson College, Carlisle, Pennsylvania, USA 4Jackson School of Geosciences, University of Texas at Austin, Austin, Texas, USA 5Earth, Planetary, and Space Sciences, University of California Los Angeles, Los Angeles, California, USA 6Idaho National Laboratory, Idaho Falls, Idaho, USA 7CSIRO Land and Water, PMB 2, Glen Osmond, Adelaide, SA 5064, Australia

Concentration–discharge patterns of weathering products from global rivers

Quantifying the functional relationships relating river discharge and weathering products places key constraints on the negative feedback between the silicate weathering and climate. In this study we analyze the concentration–discharge relationships of weathering products from global rivers using previously compiled time-series datasets for concentrations and discharge from global rivers. To analyze the nature of the covariation between specific discharge and concentrations, we use both a power law equation and a recently developed solute production equation. The solute production equation allows us to quantify weathering efficiency, or the resistance to dilution at high runoff, via the Damköhler coefficient. These results are also compared to those derived using average concentration–discharge pairs. Both the power law exponent and the Damköhler coefficient increase and asymptote as catchments exhibit increasingly chemostatic behavior, resulting in an inverse relationship between the two parameters. We also show that using the distribution of average concentration–discharge pairs from global rivers, rather than fitting concentration–discharge relationships for each individual river, underestimates global median weathering efficiency by up to a factor of ~10. This study demonstrates the utility of long time-series sampling of global rivers to elucidate controlling processes needed to quantify patterns in global silicate weathering rates.

Recent topographic evolution and erosion of the deglaciated Washington Cascades inferred from a stochastic landscape evolution model

In this study, we model postglacial surface processes and examine the evolution of the topography and denudation rates within the deglaciated Washington Cascades to understand the controls on and time scales of landscape response to changes in the surface process regime after deglaciation. The postglacial adjustment of this landscape is modeled using a geomorphic-transport-law-based numerical model that includes processes of river incision, hillslope diffusion, and stochastic landslides. The surface lowering due to landslides is parameterized using a physically based slope stability model coupled to a stochastic model of the generation of landslides. The model parameters of river incision and stochastic landslides are calibrated based on the rates and distribution of thousand-year-time scale denudation rates measured from cosmogenic 10Be isotopes. The probability distributions of those model parameters calculated based on a Bayesian inversion scheme show comparable ranges from previous studies in similar rock types and climatic conditions. The magnitude of landslide denudation rates is determined by failure density (similar to landslide frequency), whereas precipitation and slopes affect the spatial variation in landslide denudation rates. Simulation results show that postglacial denudation rates decay over time and take longer than 100 kyr to reach time-invariant rates. Over time, the landslides in the model consume the steep slopes characteristic of deglaciated landscapes. This response time scale is on the order of or longer than glacial/interglacial cycles, suggesting that frequent climatic perturbations during the Quaternary may produce a significant and prolonged impact on denudation and topography.