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.


This content will become publicly available on July 29, 2026

Title: A mathematical model for bedrock incision in near‐threshold gravel‐bed rivers
Gravel‐bed rivers that incise into bedrock are common worldwide. These systems have many similarities with other alluvial channels: they transport large amounts of sediment and adjust their forms in response to discharge and sediment supply. At the same time, the occurrence of bedrock incision implies behaviour that falls on a spectrum between fully detachment‐limited ‘bedrock channels’ and fully transport‐limited ‘alluvial channels’. Here, we present a mathematical model of river profile evolution that integrates bedrock erosion, gravel transport and the formation of channels whose hydraulic geometry is consistent with that of near‐threshold alluvial channels. We combine theory for five interrelated processes: bedload sediment transport in equilibrium gravel‐bed channels, channel width adjustment to flow and sediment characteristics, abrasion of bedrock by mobile sediment, plucking of bedrock and progressive loss of gravel‐sized sediment due to grain attrition. This model contributes to a growing class of models that seek to capture the dynamics of both bedrock incision and alluvial sediment transport. We demonstrate the model's ability to reproduce expected fluvial features such as inverse power law scaling between slope and area, and width and depth consistent with near‐threshold channel theory, and we discuss the role of sediment characteristics in influencing the mode of channel behaviour, erosional mechanism, channel steepness and profile concavity.  more » « less
Award ID(s):
1822062 2148762
PAR ID:
10556609
Author(s) / Creator(s):
; ;
Publisher / Repository:
Wiley
Date Published:
Journal Name:
Earth Surface Processes and Landforms
Volume:
49
Issue:
13
ISSN:
0197-9337
Page Range / eLocation ID:
4168 to 4186
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract. Landslides are the main source of sediment in most mountain ranges. Rivers then act as conveyor belts, evacuating landslide-derived sediment. Sediment dynamics are known to influence landscape evolution through interactions among landslide sediment delivery, fluvial transport and river incision into bedrock. Sediment delivery and its interaction with river incision therefore control the pace of landscape evolution and mediate relationships among tectonics, climate and erosion. Numerical landscape evolution models (LEMs) are well suited to study the interactions among these surface processes. They enable evaluation of a range of hypotheses at varying temporal and spatial scales. While many models have been used to study the dynamic interplay between tectonics, erosion and climate, the role of interactions between landslide-derived sediment and river incision has received much less attention. Here, we present HyLands, a hybrid landscape evolution model integrated within the TopoToolbox Landscape Evolution Model (TTLEM) framework. The hybrid nature of the model lies in its capacity to simulate both erosion and deposition at any place in the landscape due to fluvial bedrock incision, sediment transport, and rapid, stochastic mass wasting through landsliding. Fluvial sediment transport and bedrock incision are calculated using the recently developed Stream Power with Alluvium Conservation and Entrainment (SPACE) model. Therefore, rivers can dynamically transition from detachment-limited to transport-limited and from bedrock to bedrock–alluvial to fully alluviated states. Erosion and sediment production by landsliding are calculated using a Mohr–Coulomb stability analysis, while landslide-derived sediment is routed and deposited using a multiple-flow-direction, nonlinear deposition method. We describe and evaluate the HyLands 1.0 model using analytical solutions and observations. We first illustrate the functionality of HyLands to capture river dynamics ranging from detachment-limited to transport-limited conditions. Second, we apply the model to a portion of the Namche Barwa massif in eastern Tibet and compare simulated and observed landslide magnitude–frequency and area–volume scaling relationships. Finally, we illustrate the relevance of explicitly simulating landsliding and sediment dynamics over longer timescales for landscape evolution in general and river dynamics in particular. With HyLands we provide a new tool to understand both the long- and short-term coupling between stochastic hillslope processes, river incision and source-to-sink sediment dynamics. 
    more » « less
  2. Abstract Relatively little is known about the geomorphological characteristics of floodplain secondary channels and the potential for floodplain flows to mobilize bed material within these channels. This study examines the geomorphological characteristics (channel form, material properties, wood jams) and bed‐material mobilization potential of secondary channels on the floodplain of a meandering river in Illinois, USA. It also compares these attributes to those of the main channel. Results show that secondary channels are at most about one‐third the size of the main channel but also vary in size over distance. Channel dimensions tend to be greatest near the proximal connection of secondary channels to the main channel, suggesting that flow from the main channel is effective in producing scour where it enters secondary channels. The beds of secondary channels consist mainly of mud in contrast to sand and gravel on the bed of the main channel, implying that secondary channels do not convey bed material from the main channel onto the floodplain. Secondary channels connected to the main channel at both ends have more abundant active wood jams than those connected only at the proximal end. Flow from the main channel enters secondary channels at sub‐bankfull stages, but maximum mobilization of cohesive bed material in secondary channels only occurs during flows that exceed the average bankfull stage in the main channel. Overall, secondary channels are active conduits of flow, sediment, and large wood on floodplains and can contribute to floodplain sediment fluxes through entrainment of bed material. 
    more » « less
  3. Abstract Bedrock rivers are the pacesetters of landscape evolution in uplifting fluvial landscapes. Water discharge variability and sediment transport are important factors influencing bedrock river processes. However, little work has focused on the sensitivity of hillslope sediment supply to precipitation events and its implications on river evolution in tectonically active landscapes. We model the temporal variability of water discharge and the sensitivity of sediment supply to precipitation events as rivers evolve to equilibrium over 106model years. We explore how coupling sediment supply sensitivity with discharge variability influences rates and timing of river incision across climate regimes. We find that sediment supply sensitivity strongly impacts which water discharge events are the most important in driving river incision and modulates channel morphology. High sediment supply sensitivity focuses sediment delivery into the largest river discharge events, decreasing rates of bedrock incision during floods by orders of magnitude as rivers are inundated with new sediment that buries bedrock. The results show that the use of river incision models in which incision rates increase monotonically with increasing river discharge may not accurately capture bedrock river dynamics in all landscapes, particularly in steep landslide prone landscapes. From our modeling results, we hypothesize the presence of an upper discharge threshold for river incision at which storms transition from being incisional to depositional. Our work illustrates that sediment supply sensitivity must be accounted for to predict river evolution in dynamic landscapes. Our results have important implications for interpreting and predicting climatic and tectonic controls on landscape morphology and evolution. 
    more » « less
  4. Abstract Bedrock river width is an essential geometric parameter relevant to understanding flood hazards and gauging station rating curves, and is critical to stream power incision models and many other landscape evolution models. Obtaining bedrock river width measurements, however, typically requires extensive field campaigns that take place in rugged and steep topography where river access is often physically challenging. Although prior work has turned to measuring channel width from satellite imagery, these data present a snapshot in time, are typically limited to rivers ≥ 10–30 m wide due to the image resolution, and are physically restricted to areas devoid of vegetation. For these reasons, we are generally data limited, and the factors impacting bedrock channel width remain poorly understood. Due to these limitations, researchers often turn to assumptions of width‐scaling relationships with drainage area or discharge to estimate bedrock channel width. Here we present a new method of obtaining bedrock channel width at a desired river discharge through the incorporation of a high‐resolution bare‐earth digital elevation model (DEM) using MATLAB Topotoolbox and the HEC‐RAS river analysis system. We validate this method by comparing modeled results to US Geological Survey (USGS) field measurements at existing gauging stations, as well as field channel measurements. We show that this method can capture general characteristics of discharge rating curves and predict field‐measured channel widths within uncertainty. As high‐resolution DEMs become more available across the United States through the USGS three‐dimensional elevation program (3DEP), the future utility of this method is notable. Through developing and validating a streamlined, open‐source, and freely available workflow of channel width extraction, we hope this method can be applied to future research to improve the quantity of channel width measurements and improve our understanding of bedrock channels. 
    more » « less
  5. Abstract Meandering river floodplains often contain intermittently flooded complex channel networks. Many questions remain as to the pervasiveness, function, and evolution of these floodplain channels. In this present work, we analyzed size‐specific sediment transport potential and assessed whether the channelized floodplain of the meandering East Fork White River near Seymour, Indiana is on a net erosional or depositional trajectory. We applied a two‐dimensional hydrodynamic model and used simulated model results to estimate the largest sediment size that can be moved in suspension and as bedload at various flows for grain size classes between 4 µm and 64 mm. We developed a probabilistic method that integrates the largest sediment size that can be moved at various flows to compute an effective grain size, which we compared to measured field data. Results show that the river is capable of supplying sand to the floodplain and these floodplain channels can transport sand in suspension and gravel as bedload. This suggests that sediment supplied from the river could be transported as bedload in floodplain channels. These floodplain channels are supply limited under the current hydrologic regime and the grain size distribution of the bed surface is set by the flow conditions; thus, these floodplain channels are net erosional. Finally, our proposed method of probabilistically integrating the largest sediment size that can be moved at various flows can be used to predict the upper end of the grain size distribution in suspension and in bed material, which is applicable to floodplains as well as coastal areas. 
    more » « less