skip to main content

Title: The role of lateral erosion in the evolution of nondendritic drainage networks to dendricity and the persistence of dynamic networks

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 » landscapes on Earth.

« less
; ;
Award ID(s):
Publication Date:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range or eLocation-ID:
Article No. e2015770118
Proceedings of the National Academy of Sciences
Sponsoring Org:
National Science Foundation
More Like this
  1. 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 »a progressive upslope increase in the residence time of bedrock in the weathered zone. Two nondimensional numbers corresponding to the mean hillslope gradient and mean groundwater-table gradient emerge and their ratio defines the proportion of the hillslope relief that is unweathered. Field data from three field sites are consistent with model predictions.« less
  2. 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.
  3. 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 »evolution of the Western North Boulder fault to show that through time, these individual segments propagate at the fault tips and link together to span over 40 km, with a maximum slip of 462 m in the central portion of the fault. Fault slip rates range from 0.02 to 0.45 mm/yr along strike and are consistent with estimates for other active faults in the region. We find that the timing of fault initiation coincides well with the migration of the Yellowstone hotspot across the nearby Idaho-Montana border and thus attribute the initiation of extension to the crustal bulge from the migrating hotspot. Overall, we provide the first quantitative constraints on fault initiation and evolution of the Western North Boulder fault, perhaps the farthest north basin in the Northern Basin and Range province that such constraints exist. We show that river profiles are powerful tools for documenting the spatial and temporal patterns of normal fault evolution, especially where other geologic/geodetic methods are limited, proving to be a vital tool for accurate tectonic hazard assessments.« less
  4. 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 »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.« less
  5. 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 »the channel widens, suggesting that, ashypothesized, the channel width adjusts faster to post-reorganizationdrainage area distribution. The adjusted narrow channel dictates the widthof formative flows in the reversed valley, which contrasts with the meaningfullywider formative flows of the beheaded valley across the divide. Thisdifference results in a step change in the unit stream power between thereversed and beheaded channels, potentially leading to a “width feedback”that promotes ongoing divide migration and reorganization. Our findings demonstrate that valley width–area scaling is a potential toolfor identifying landscapes influenced by drainage reorganization. Accountingfor reorganization-specific scaling can improve estimations of erosion ratedistributions in reorganized landscapes.« less