Dendritic, i.e., tree-like, river networks are ubiquitous features on Earth’s landscapes; however, how and why river networks organize themselves into this form are incompletely understood. A branching pattern has been argued to be an optimal state. Therefore, we should expect models of river evolution to drastically reorganize (suboptimal) purely nondendritic networks into (more optimal) dendritic networks. To date, current physically based models of river basin evolution are incapable of achieving this result without substantial allogenic forcing. Here, we present a model that does indeed accomplish massive drainage reorganization. The key feature in our model is basin-wide lateral incision of bedrock channels. The addition of this submodel allows for channels to laterally migrate, which generates river capture events and drainage migration. An important factor in the model that dictates the rate and frequency of drainage network reorganization is the ratio of two parameters, the lateral and vertical rock erodibility constants. In addition, our model is unique from others because its simulations approach a dynamic steady state. At a dynamic steady state, drainage networks persistently reorganize instead of approaching a stable configuration. Our model results suggest that lateral bedrock incision processes can drive major drainage reorganization and explain apparent long-lived transience in more »
- Award ID(s):
- 1833025
- Publication Date:
- NSF-PAR ID:
- 10221459
- Journal Name:
- Proceedings of the National Academy of Sciences
- Volume:
- 118
- Issue:
- 16
- Page Range or eLocation-ID:
- Article No. e2015770118
- ISSN:
- 0027-8424
- Publisher:
- Proceedings of the National Academy of Sciences
- Sponsoring Org:
- National Science Foundation
More Like this
-
The depth to unweathered bedrock beneath landscapes influences subsurface runoff paths, erosional processes, moisture availability to biota, and water flux to the atmosphere. Here we propose a quantitative model to predict the vertical extent of weathered rock underlying soil-mantled hillslopes. We hypothesize that once fresh bedrock, saturated with nearly stagnant fluid, is advected into the near surface through uplift and erosion, channel incision produces a lateral head gradient within the fresh bedrock inducing drainage toward the channel. Drainage of the fresh bedrock causes weathering through drying and permits the introduction of atmospheric and biotically controlled acids and oxidants such that the boundary between weathered and unweathered bedrock is set by the uppermost elevation of undrained fresh bedrock, Z b . The slow drainage of fresh bedrock exerts a “bottom up” control on the advance of the weathering front. The thickness of the weathered zone is calculated as the difference between the predicted topographic surface profile (driven by erosion) and the predicted groundwater profile (driven by drainage of fresh bedrock). For the steady-state, soil-mantled case, a coupled analytical solution arises in which both profiles are driven by channel incision. The model predicts a thickening of the weathered zone upslope and, consequently,more »
-
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.
-
Valla, Pierre (Ed.)Abstract Over the past few decades, tectonic geomorphology has been widely implemented to constrain spatial and temporal patterns of fault slip, especially where existing geologic or geodetic data are poor. We apply this practice along the eastern margin of Bull Mountain, Southwest Montana, where 15 transient channels are eroding into the flat, upstream relict landscape in response to an ongoing period of increased base level fall along the Western North Boulder fault. We aim to improve constraints on the spatial and temporal slip rates across the Western North Boulder fault zone by applying channel morphometrics, cosmogenic erosion rates, bedrock characteristics, and calibrated reproductions of the modern river profiles using a 1-dimensional stream power incision model that undergoes a change in the rate of base level fall. We perform over 104 base level fall simulations to explore a wide range of fault slip dynamics and stream power parameters. Our best fit simulations suggest that the Western North Boulder fault started as individual fault segments along the middle to southern regions of Bull Mountain that nucleated around 6.2 to 2.5 Ma, respectively. This was followed by the nucleation of fault segments in the northern region around 1.5 to 0.4 Ma. We recreate themore »
-
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)more »
-
Abstract. The width of valleys and channels affects the hydrology, ecology,and geomorphic functionality of drainage networks. In many studies, thewidth of valleys and/or channels (W) is estimated as a power-law function ofthe drainage area (A), W=kcAd. However, in fluvial systemsthat experience drainage reorganization, abrupt changes in drainage areadistribution can result in valley or channel widths that are disproportionalto their drainage areas. Such disproportionality may be more distinguishedin valleys than in channels due to a longer adjustment timescale forvalleys. Therefore, the valley width–area scaling in reorganized drainagesis expected to deviate from that of drainages that did not experiencereorganization. To explore the effect of reorganization on valley width–drainage areascaling, we studied 12 valley sections in the Negev desert, Israel,categorized into undisturbed, beheaded, and reversed valleys. We found thatthe values of the drainage area exponents, d, are lower in the beheadedvalleys relative to undisturbed valleys but remain positive. Reversedvalleys, in contrast, are characterized by negative d exponents, indicatingvalley narrowing with increasing drainage area. In the reversed category, wealso explored the independent effect of channel slope (S) through theequation W=kbAbSc, which yieldednegative and overall similar values for b and c. A detailed study in one reversed valley section shows that the valleynarrows downstream, whereasmore »