Sediment transport through self-adjusting, bedrock-walled waterfall plunge pools

Many waterfalls have deep plunge pools that are often partially or fully filled with sediment. Sediment fill may control plunge-pool bedrock erosion rates, partially determine habitat availability for aquatic organisms, and affect sediment routing and debris flow initiation. Currently, there exists no mechanistic model to describe sediment transport through waterfall plunge pools. Here we develop an analytical model to predict steady-state plunge-pool depth and sediment-transport capacity by combining existing jet theory with sediment transport mechanics. Our model predicts plunge-pool sediment-transport capacity increases with increasing river discharge, flow velocity, and waterfall drop height and decreases with increasing plunge-pool depth, radius, and grain size. We tested the model using flume experiments under varying waterfall and plunge-pool geometries, flow hydraulics, and sediment size. The model and experiments show that through morphodynamic feedbacks, plunge pools aggrade to reach shallower equilibrium pool depths in response to increases in imposed sediment supply. Our theory for steady-state pool depth matches the experiments with an R^2 value of 0.8, with discrepancies likely due to model simplifications of the hydraulics and sediment transport. Analysis of 75 waterfalls suggests that the water depths in natural plunge pools are strongly influenced by upstream sediment supply, and our model provides a mass-conserving framework to predict sediment and water storage in waterfall plunge pools for sediment routing, habitat assessment, and bedrock erosion modeling

A Mechanistic Model of Waterfall Plunge Pool Erosion into Bedrock

Landscapes often respond to changes in climate and tectonics through the formation and upstream propagation of knickzones composed of waterfalls. Little work has been done on the mechanics of waterfall erosion, and instead most landscape-scale models neglect waterfalls or use rules for river erosion, such as stream power, that may not be applicable to waterfalls. Here we develop a physically based model to predict waterfall plunge pool erosion into rock by abrasion from particle impacts and test the model against flume experiments. Both the model and experiments show that evolving plunge pools have initially high vertical erosion rates due to energetic particle impacts, and erosion slows and eventually ceases as pools deepen and deposition protects the pool floor from further erosion. Lateral erosion can continue after deposition on the pool floor, but it occurs at slow rates that become negligible as pools widen. Our work points to the importance of vertical drilling of successive plunge pools to drive upstream knickzone propagation in homogenous rock, rather than the classic mechanism of headwall undercutting. For a series of vertically drilling waterfalls, we find that upstream knickzone propagation is faster under higher combined water and sediment fluxes and for knickzones composed of many waterfalls that are closely spaced. Our model differs significantly from stream-power-based erosion rules in that steeper knickzones can retreat faster or more slowly depending on the number and spacing of waterfalls within a knickzone, which has implications for interpreting climatic and tectonic history through analysis of river longitudinal profiles

A Mechanistic Model of Waterfall Plunge Pool Erosion into Bedrock

Landscapes often respond to changes in climate and tectonics through the formation and upstream propagation of knickzones composed of waterfalls. Little work has been done on the mechanics of waterfall erosion, and instead most landscape-scale models neglect waterfalls or use rules for river erosion, such as stream power, that may not be applicable to waterfalls. Here we develop a physically based model to predict waterfall plunge pool erosion into rock by abrasion from particle impacts and test the model against flume experiments. Both the model and experiments show that evolving plunge pools have initially high vertical erosion rates due to energetic particle impacts, and erosion slows and eventually ceases as pools deepen and deposition protects the pool floor from further erosion. Lateral erosion can continue after deposition on the pool floor, but it occurs at slow rates that become negligible as pools widen. Our work points to the importance of vertical drilling of successive plunge pools to drive upstream knickzone propagation in homogenous rock, rather than the classic mechanism of headwall undercutting. For a series of vertically drilling waterfalls, we find that upstream knickzone propagation is faster under higher combined water and sediment fluxes and for knickzones composed of many waterfalls that are closely spaced. Our model differs significantly from stream-power-based erosion rules in that steeper knickzones can retreat faster or more slowly depending on the number and spacing of waterfalls within a knickzone, which has implications for interpreting climatic and tectonic history through analysis of river longitudinal profiles

Self-formed waterfall plunge pools in homogeneous rock

Waterfalls are ubiquitous, and their upstream propagation can set the pace of landscape evolution, yet no experimental studies have examined waterfall plunge pool erosion in homogeneous rock. We performed laboratory experiments, using synthetic foam as a bedrock simulant, to produce self-formed waterfall plunge pools via particle impact abrasion. Plunge pool vertical incision exceeded lateral erosion by approximately tenfold until pools deepened to the point that the supplied sediment could not be evacuated and deposition armored the pool bedrock floor. Lateral erosion of plunge pool sidewalls continued after sediment deposition, but primarily at the downstream pool wall, which might lead to undermining of the plunge pool lip, sediment evacuation, and continued vertical pool floor incision in natural streams. Undercutting of the upstream pool wall was absent, and our results suggest that vertical drilling of successive plunge pools is a more efficient waterfall retreat mechanism than the classic model of headwall undercutting and collapse in homogeneous rock

Deltaic deposits at Aeolis Dorsa: Sedimentary evidence for a standing body of water on the northern plains of Mars

A fundamental long-standing question regarding Mars history is whether the flat and low-lying northern plains ever hosted an ocean. The best opportunity to solve this problem is provided by stratigraphic observations of sedimentary deposits onlapping the crustal dichotomy. Here, we use high-resolution imagery and topography to analyze a branching network of inverted channel and channel lobe deposits in the Aeolis Dorsa region, just north of the dichotomy boundary. Observations of stacked, cross-cutting channel bodies and stratal geometries indicate that these landforms represent exhumed distributary channel deposits. Observations of depositional trunk feeder channel bodies, a lack of evidence for past topographic confinement, channel avulsions at similar elevations, and the presence of a strong break in dip slope between topset and foreset beds suggest that this distributary system was most likely a delta, rather than an alluvial fan or submarine fan. Sediment transport calculations using both measured and derived channel geometries indicate a minimum delta deposition time on the order of 400  years. The location of this delta within a thick and widespread clastic wedge abutting the crustal dichotomy boundary, unconfined by any observable craters, suggests a standing body of water potentially 105 km2 in extent or greater and is spatially consistent with hypotheses for a northern ocean

Influence of bed patchiness, slope, grain hiding, and form drag on gravel mobilization in very steep streams

Steep streams are a major portion of channel networks and provide a link to transport sediment from hillslopes to lower gradient rivers. Despite their importance, key unknowns remain, perhaps foremost of which is evaluating in steep streams empirical laws for fluvial sediment transport developed for low-gradient rivers. To address this knowledge gap, we painted sediment in situ over 3 years to monitor incipient sediment motion and sediment-patch development in five small (drainage areas of 0.04–2 km^2) and steep (slopes of 5–37%) tributaries of Elder Creek, California, United States. We found that channel beds organized into size-sorted sediment patches which displayed active fluvial transport of gravel annually, consistent year-to-year patch median grain sizes, partial transport of bed material, and significantly higher values of critical Shields stress for incipient sediment motion compared to that observed for lower gradient rivers. The high critical Shields stresses (up to ≈0.5 for the median grain size) agree within a factor of ~3 to theoretical predictions which account for slope-dependent hydraulics, grain hiding, and sediment patches. For grains of approximately the same size as the roughness length scale, slope-dependent hydraulics and bed patchiness are the dominant controls on critical Shields stress values, while grain hiding is important for grains larger or smaller than the roughness length scale. Form drag exists in our monitored tributaries but has a smaller influence than the above effects. Our field observations show fluvial processes contribute to sediment mobilization in steep channels which are often considered to be dominated by debris flows

A model for fire‐induced sediment yield by dry ravel in steep landscapes

Sediment flux from hillslopes to channels commonly increases following wildfires, with implications for the carbon cycle, river habitats, and debris-flow hazards. Although much of this material is transported via dry ravel, existing ravel models are not applicable to hillslopes with gradients greater than the angle of repose, which can constitute the majority of mountainous terrain. To fill this knowledge gap, we develop a continuity model for sediment storage by vegetation dams on steep hillslopes to predict sediment yields following wildfire. The maximum volume of sediment stored prior to wildfire is set to be a function of vegetation density, the capacity of plants to impound sediment, and the contributing hillslope area. Time is required after fire to establish vegetation and replenish hillslope sediment storage, which introduces vegetation regrowth rate, soil production rate, and fire recurrence interval as important variables that affect ravel yield. Model results for the San Gabriel Mountains, California, predict that sediment yield can increase by several orders of magnitude following fire. These results are consistent with field data of ravel yield (~30 mm per contributing area of hillslope in 5 months) we collected following the 2009 Station Fire, as well as postfire sediment flux recorded by 93 debris basins. In contrast to previous work, our model shows that heightened postfire sediment yields can be explained by a change in hillslope sediment storage independent of major changes in the soil production rate and landscape form over geomorphic timescales

Fluvial organic carbon cycling regulated by sediment transit time and mineral protection

Rivers transfer terrestrial organic carbon (OC) from mountains to ocean basins, playing a key role in the global carbon cycle. During fluvial transit, OC may be oxidized and emitted to the atmosphere as CO2 or preserved and transported to downstream depositional sinks. The balance between oxidation and preservation determines the amount of particulate OC (POC) that can be buried long term, but the factors regulating this balance are poorly constrained. Here, we quantify the effects of fluvial transit on POC fluxes along an ~1,300 km lowland channel with no tributaries. We show that sediment transit time and mineral protection regulate the magnitude and rate of POC oxidation, respectively. Using a simple turnover model, we estimate that annual POC oxidation is a small percentage of the POC delivered to the river. Modelling shows that lateral erosion into POC-rich floodplains can increase POC fluxes to downstream basins, thereby offsetting POC oxidation. Consequently, rivers with high channel mobility can enhance CO2 drawdown while management practices that stabilize river channels may reduce the potential for CO2 drawdown

Fluvial organic carbon composition regulated by seasonal variability in lowland river migration and water discharge

Identifying drivers of seasonal variations in fluvial particulate organic carbon (POC) composition can aid sediment provenance and biogeochemical cycling studies. We evaluate seasonal changes in POC composition in the Río Bermejo, Argentina, a lowland river running ∼1,270 km without tributaries. Weekly POC concentration and isotopic composition from 2016 to 2018 show that during the wet season, increased lateral channel migration generates an influx of 13C-enriched and 14C-enriched floodplain-sourced material, overprinting the 13C-depleted and 14C-depleted headwater signature that is observed during the dry season. These findings demonstrate how channel morphodynamics can drive variability of POC composition in lowland rivers, and may modulate the composition of POC preserved in sedimentary archives