We explore how rock properties and channel morphology vary with rock type in Last Chance Canyon, Guadalupe Mountains, New Mexico, USA. The rocks here are composed of horizontally to near-horizontally interbedded carbonate and sandstone. This study focuses on first- and second-order channel sections, where the streams have a lower channel steepness index (ksn) upstream and transition to higher ksn values downstream. We hypothesize that differences in bed thickness and rock strength influence ksn values, both locally by influencing bulk bedrock strength and also nonlocally through the production of coarse sediment. We collected discontinuity intensity data (the length of bedding planes and fractures per unit area), Schmidt hammer rebound measurements, and measured the largest boulder at every 12.2 m elevation contour to test this hypothesis. Bedrock and boulder mineralogy were determined using a lab-based carbonate dissolution method. High-resolution orthomosaics and digital surface models (DSMs) were generated from drone and ground-based photogrammetry. The orthomosaics were used to map channel sections with exposed bedrock. The United States Geological Survey (USGS) 10 m digital elevation models (DEMs) were used to measure channel slope and hillslope relief. We find that discontinuity intensity is negatively correlated with Schmidt hammer rebound values in sandstone bedrock. Channel steepness tends to be higher where reaches are primarily incising through more thickly bedded carbonate bedrock and lower where more thinly bedded sandstone is exposed. Bedrock properties also influence channel morphology indirectly, through coarse sediment input from adjacent hillslopes. Thickly bedded rock layers on hillslopes erode to contribute larger colluvial sediment to adjacent channels, and these reaches have higher ksn values. Larger and more competent carbonate sediment armors both the carbonate and the more erodible sandstone and reduces steepness contrasts across rock types. We interpret that in the relatively steep, high-level ksn downstream channel sections, the slope is primarily controlled by the coarse alluvial cover. We further posit that the upstream low-level ksn reaches have a base level that is fixed by the steep downstream reaches, resulting in a stable configuration, where channel slopes have adjusted to lithologic differences and/or sediment armor.
more »
« less
Sediment Cover Modulates Landscape Erosion Patterns and Channel Steepness in Layered Rocks: Insights From the SPACE Model
Abstract Erosional perturbations from changes in climate or tectonics are recorded in the profiles of bedrock rivers, but these signals can be challenging to unravel in settings with non‐uniform lithology. In layered rocks, the surface lithology at a given location varies through time as erosion exposes different layers of rock. Recent modeling studies have used the Stream Power Model (SPM) to highlight complex variations in erosion rates that arise in bedrock rivers incising through layered rocks. However, these studies do not capture the effects of coarse sediment cover on channel evolution. We use the “Stream Power with Alluvium Conservation and Entrainment” (SPACE) model to explore how sediment cover influences landscape evolution and modulates the topographic expression of erodibility contrasts in horizontally layered rocks. We simulate river evolution through alternating layers of hard and soft rock over million‐year timescales with a constant and uniform uplift rate. Compared to the SPM, model runs with sediment cover have systematically higher channel steepness values in soft rock layers and lower channel steepness values in hard rock layers. As more sediment accumulates, the contrast in steepness between the two rock types decreases. Effective bedrock erodibilities back‐calculated assuming the SPM are strongly influenced by sediment cover. We also find that sediment cover can significantly increase total relief and timescales of adjustment toward landscape‐averaged steady‐state topography and erosion rates.
more »
« less
- Award ID(s):
- 1918351
- PAR ID:
- 10531278
- Publisher / Repository:
- DOI PREFIX: 10.1029
- Date Published:
- Journal Name:
- Journal of Geophysical Research: Earth Surface
- Volume:
- 129
- Issue:
- 7
- ISSN:
- 2169-9003
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract We present 17 new 10Be erosion rates from southern Peru sampled across an extreme orographic rainfall gradient. Using a rainfall-weighted variant of the normalized channel steepness index, ksnQ, we show that channel steepness values, and thus topography, are adjusted to spatially varying rainfall. Rocks with similar physical characteristics define distinct relationships between ksnQ and erosion rate (E), suggesting ksnQ is also resolving lithologic variations in erodibility. However, substantial uncertainty exists in parameters describing these relationships. By combining our new data with 38 published erosion rates from Peru and Bolivia, we collapse the range of compatible parameter values and resolve robust, nonlinear ksnQ–E relationships suggestive of important influences of erosional thresholds, rock properties, sediment characteristics, and temporal runoff variability. In contrast, neither climatic nor lithologic effects are clear using the traditional channel steepness metric, ksn. Our results highlight that accounting for spatial rainfall variations is essential for disentangling the multiple influences of climate, lithology, and tectonics common in mountain landscapes, which is a necessary first step toward greater understanding of how these landscapes evolve.more » « less
-
Abstract. Landscape morphology reflects drivers such as tectonicsand climate but is also modulated by underlying rock properties. Whilegeomorphologists may attempt to quantify the influence of rock strengththrough direct comparisons of landscape morphology and rock strengthmetrics, recent work has shown that the contact migration resulting from the presence of mixed lithologies may hinder such an approach. Indeed, this work counterintuitively suggests that channel slopes within weaker units can sometimes be higher than channel slopes within stronger units. Here, we expand upon previous work with 1-D stream power numerical models in which we have created a system for quantifying contact migration over time. Although previous studies have developed theories for bedrock rivers incising through layered stratigraphy, we can now scrutinize these theories with contact migration rates measured in our models. Our results show that previously developed theory is generally robust and that contact migration rates reflect the pattern of kinematic wave speed across the profile. Furthermore, we have developed and tested a new approach for estimating kinematic wave speeds. This approach utilizes channel steepness, a known base-level fall rate, and contact dips. Importantly, we demonstrate how this new approach can be combined with previous work to estimate erodibility values. We demonstrate this approach by accurately estimating the erodibility values used in our numerical models. After this demonstration, we use our approach to estimate erodibility values for a stream near Hanksville, UT. Because we show in our numerical models that one can estimate the erodibility of the unit with lower steepness, the erodibilities we estimate for this stream in Utah are likely representative of mudstone and/or siltstone. The methods we have developed can be applied to streams with temporally constant base-level fall, opening new avenues of research within the field of geomorphology.more » « less
-
Abstract Prior numerical modeling work has suggested that incision into sub‐horizontal layered stratigraphy with variable erodibility induces non‐uniform erosion rates even if base‐level fall is steady and sustained. Erosion rates of cliff bands formed in the stronger rocks in a stratigraphic sequence can greatly exceed the rate of base‐level fall. Where quartz in downstream sediment is sourced primarily from the stronger, cliff‐forming units, erosion rates estimated from concentrations of cosmogenic beryllium‐10 (10Be) in detrital sediment will reflect the locally high erosion rates in retreating cliff bands. We derive theoretical relationships for threshold hillslopes and channels described by the stream‐power incision model as a quantitative guide to the potential magnitude of this amplification of10Be‐derived erosion rates above the rate of base‐level fall. Our analyses predict that the degree of erosion rate amplification is a function of bedding dip and either the ratio of rock erodibility in alternating strong and weak layers in the channel network, or the ratio of cliff to intervening‐slope gradient on threshold hillslopes. We test our predictions in the cliff‐and‐bench landscape of the Grand Staircase in southern Utah, USA. We show that detrital cosmogenic erosion rates in this landscape are significantly higher (median 300 m/Ma) than the base‐level fall rate (~75 m/Ma) determined from the incision rate of a trunk stream into a ~0.6 Ma basalt flow emplaced along a 16 km reach of the channel. We infer a 3–6‐fold range in rock strength from near‐surface P‐wave velocity measurements. The approximately four‐fold difference between the median10Be‐derived erosion rate and the long‐term rate of base‐level fall is consistent with our model and the observation that the stronger, cliff‐forming lithologies in this landscape are the primary source of quartz in detrital sediments. © 2020 John Wiley & Sons, Ltd.more » « less
-
Abstract In most landscape evolution models, extreme rainfall enhances river incision. In steep landscapes, however, these events trigger landslides that can buffer incision via increased sediment delivery and aggradation. We quantify landslide sediment aggradation and erosional buffering with a natural experiment in southern Taiwan where a northward gradient in tectonic activity drives increasing landscape steepness. We find that landscape response to extreme rainfall during the 2009 typhoon Morakot varied along this gradient, where steep areas experienced widespread channel sediment aggradation of >10 m and less steep areas did not noticeably aggrade. We model sediment export to estimate a sediment removal timeline and find that steep, tectonically active areas with the most aggradation may take centuries to resume bedrock incision. Expected sediment cover duration reflects tectonic uplift. We find that despite high stream power, sediment cover may keep steep channels from eroding bedrock for up to half of a given time period. This work highlights the importance of dynamic sediment cover in landscape evolution and suggests a mechanism by which erosional efficiency in tectonically active landscapes may decrease as landscape steepness increases.more » « less
An official website of the United States government
