skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Rock properties and sediment caliber govern bedrock river morphology across the Taiwan Central Range
Feedbacks between surface and deep Earth processes in collisional mountain belts depend on how erosion and topographic relief vary in space and time. One outstanding unknown lies in how rock strength influences bedrock river morphology and thus mountain relief. Here, we quantify boulder cover and channel morphology using uncrewed aerial vehicle surveys along 30 kilometers of bedrock-bound river corridors throughout the Taiwan Central Range where regional gradients in rock properties relate to tectonic history. We find that boulder size systematically increases with increasing metamorphic grade and depth of exhumation. Boulder size correlates with reach-scale channel steepness but does not explain observations of highly variable channel width. Transport thresholds indicate that rivers are adjusted to mobilize boulders and are well in excess of the threshold to transport gravel and cobbles, as previously assumed. The linkage between metamorphic history, boulder size, and channel steepness reveals how rock properties can influence feedbacks between tectonics and topography throughout the life span of a mountain range.  more » « less
Award ID(s):
1933203 2231170
PAR ID:
10502479
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
American Association for the Advancement of Science
Date Published:
Journal Name:
Science Advances
Volume:
9
Issue:
46
ISSN:
2375-2548
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. 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
  3. Abstract River profiles are shaped by climatic and tectonic history, lithology, and internal feedbacks between flow hydraulics, sediment transport and erosion. In steep channels, waterfalls may self‐form without changes in external forcing (i.e., autogenic formation) and erode at rates faster or slower than an equivalent channel without waterfalls. We use a 1‐D numerical model to investigate how self‐formed waterfalls alter the morphology of bedrock river longitudinal profiles. We modify the standard stream power model to include a slope threshold above which waterfalls spontaneously form and a rate constant allowing waterfalls to erode faster or slower than other fluvial processes. Using this model, we explore how waterfall formation alters both steady state and transient longitudinal profile forms. Our model predicts that fast waterfalls create km‐scale reaches in a dynamic equilibrium with channel slope held approximately constant at the threshold slope for waterfall formation, while slow waterfalls can create local channel slope maxima at the location of slow waterfall development. Furthermore, slow waterfall profiles integrate past base level histories, leading to multiple possible profile forms, even at steady‐state. Consistency between our model predictions and field observations of waterfall‐rich rivers in the Kings and Kaweah drainages in the southern Sierra Nevada, California, supports the hypothesis that waterfall formation can modulate river profiles in nature. Our findings may help identify how bedrock channels are influenced by waterfall erosion and aid in distinguishing between signatures of external and internal perturbations, thereby strengthening our ability to interpret past climate and tectonic changes from river longitudinal profiles. 
    more » « less
  4. 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
  5. Although rates of fluvial incision across the Colorado Plateau are known reasonably well, rate variability through time and its controlling processes are still poorly understood. We used boulder armored benches from the Teasdale-Torrey lowlands reach of the Fremont River in the northwestern Colorado Plateau (Utah, USA) as temporal markers to determine regional incision rates and explore controls on rate variability. Bench gravels are sourced from Tertiary volcanic rocks capping nearby Boulder and Thousand Lakes Mountains. The sedimentology of bench deposits suggests that most form from mass movement with later fluvial reworking. Volcanic boulders are tougher than the local sedimentary bedrock, which promotes boulder armoring and topographic inversion. Thirty-seven boulder cosmogenic 3He exposure ages from 11 different benches range from >600 ka to ca. 100 ka. Soil carbonate stages from two benches are in good agreement with surface exposure ages. Averaged Fremont River and tributary incision rates determined from bench exposure ages are 32% faster for tributaries off of Thousand Lakes Mountain (0.41 m/k.y.) than tributaries off of Boulder Mountain (0.28 m/k.y.). This difference in incision rate may be due to Laramideage structures limiting incision for the tributaries that drain Boulder Mountain and extensive Pleistocene ice caps on Boulder Mountain creating a wider and thicker boulder armor slowing incision. 
    more » « less