skip to main content


Title: Modeling cave cross‐section evolution including sediment transport and paragenesis
Summary

Deposits within caves are often used to interpret past landscape evolution and climate conditions. However, cave passage shapes also preserve information about past conditions. Despite the usefulness of passage shape, no previous models simulate cave cross‐section evolution in a realistic manner. Here we develop a model for evolving cave passage cross‐sections using a shear stress estimation algorithm and a shear stress erosion rule. Our model qualitatively duplicates observed cave passage shapes so long as erosion rates vary with shear stress, as in the case of transport limited dissolution or mechanical erosion. This result provides further evidence that erosion rates within caves are not typically limited by surface reaction rates, even though current speleogenesis models predict surface‐rate limitation under most turbulent flow conditions. By adding sediment transport and alluviation to the model we successfully simulate paragenetic channels. Simulations duplicate the hypothesized dynamics of paragenesis, whereby: 1) the cross‐section of a phreatic passage grows until shear stress is sufficiently reduced that alluviation occurs, 2) the floor of the passage becomes armored and erosion continues on the ceiling and walls, 3) negative feedback produces an equilibrium cross‐sectional area such that shear stress is sufficient to transport incoming sediment. We derive an approximate scaling relationship that indicates that equilibrium paragenetic channel width scales with the square root of discharge, and weakly with the inverse of sediment supply. Simulations confirm this relationship and show that erosion mechanism, sediment size, and roughness are secondary controls. The inverse scaling of width with sediment supply in paragenetic channels contrasts with surface bedrock channels, which respond to larger sediment supplies by widening. Our model provides a first step in simulating cave cross‐section evolution and points to the need for a better understanding of the dominant erosion mechanisms in soluble bedrock channels. © 2020 John Wiley & Sons, Ltd.

 
more » « less
NSF-PAR ID:
10455471
Author(s) / Creator(s):
 ;  
Publisher / Repository:
Wiley Blackwell (John Wiley & Sons)
Date Published:
Journal Name:
Earth Surface Processes and Landforms
Volume:
45
Issue:
11
ISSN:
0197-9337
Page Range / eLocation ID:
p. 2588-2602
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. We propose an exploratory model to describe the morphodynamics of distributary channel network growth on river deltas. The interface between deep channels and the shallow, unchannelized delta front deposits is modeled as a moving boundary. Steady flow over the unchannelized delta front is friction dominated and modeled by Laplace's equation. Shear stress along the network boundary produces nonlinear erosion rates at the interface, causing the boundary to move and network elements (channels and branches) to form. The model was run for boundary conditions resembling the Wax Lake Delta in coastal Louisiana, 20 parameterizations of sediment transport, and 3 parameterizations of discharge. In each case, the model produced a complex channel network with channel number, width, bifurcation angle, and channel shape depending on the sediment transport formula. For reasonable sediment transport parameters and gradually increasing water discharge, the model produced network characteristics and progradation rates similar to the Wax Lake Delta. This suggests that the model contains the processes responsible for network growth, despite its abstract formulation.

     
    more » « less
  3. Abstract

    Sediment transport and channel morphology in mountainous hillslope‐coupled streams reflect a mixture of hillslope and channel processes. However, the influence of lithology on channel form and adjustment and sediment transport remains poorly understood. Patterns of channel form, grain size, and transport capacity were investigated in two gravel‐bed streams with contrasting lithology (basalt and sandstone) in the Oregon Coast Range, USA, in a region in which widespread landslides and debris flows occurred in 1996. This information was used to evaluate threshold channel conditions and channel bed adjustment since 1996. Channel geometry, slope, and valley width were measured or extracted from LiDAR and sediment textures were measured in the surface and subsurface. Similar coarsening patterns in the first few kilometres of both streams indicated strong hillslope influences, but subsequent downstream fining was lithology‐dependent. Despite these differences, surface grain size was strongly related to shear stress, such that the ratio of available to critical shear stress for motion of the median surface grain size at bankfull stage was around one over most of the surveyed lengths. This indicated hydraulic sorting of supplied sediment, independent of lithology. We infer a cycle of adjustment to sediment delivered during the 1996 flooding, from threshold conditions, to non‐alluvial characteristics, to threshold conditions in both basins. The sandstone basin can also experience complete depletion of the gravel‐size alluvium to sand size, leading to bedrock exposure because of high diminution rates. Although debris flows being more frequent in a basalt basin, this system will likely display threshold‐like characteristics over a longer period, indicating that the lithologic control on channel adjustment is driven by differences in rock competence that control grain size and available gravel for bed load transport. © 2020 John Wiley & Sons, Ltd.

     
    more » « less
  4. Abstract

    Meandering gravel‐bed rivers tend to exhibit bed surface sorting patterns with coarse particles located in pools and fine particles on bar tops. The mechanism by which these patterns emerge has been explored in sand‐bed reaches; however, for gravel‐bed meandering channels it remains poorly understood. Here we present results from a flume experiment in which bed morphology, velocity, sediment sorting patterns, and bed load transport were intensively documented. The experimental channel is 1.35 m wide, 15.2 m long, and its centerline follows a sine‐generated curve with a crossing angle of 20°. Water and sediment input were held constant throughout the experiment and measurements were collected under quasi‐equilibrium conditions. Boundary shear stress calculated from near‐bed velocity measurements indicates that in a channel with mild sinuosity, deposition of fine particles on bars is a result of divergent shear stress at the inside bend of the channel, downstream of the apex. Boundary shear stress in the upstream half of the pool was below critical for coarse particles (>8 mm), leading to an armored pool. Inward directed selective transport was responsible for winnowing of fine particles in the pool. Fine and coarse sediment followed similar trajectories through the meander bend, which contrasts earlier studies of sand‐bedded meanders where the loci of fine and coarse particles cross paths. This suggests a different sorting mechanism for gravel bends. This experiment shows that a complex interaction of quasi‐equilibrium bed topography, selective sediment transport, and secondary currents are responsible for the sorting patterns seen in gravel‐bed, meandering channels.

     
    more » « less
  5. 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