skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 10:00 PM ET on Friday, December 8 until 2:00 AM ET on Saturday, December 9 due to maintenance. We apologize for the inconvenience.

Title: Resistance Is Not Futile: Grain Resistance Controls on Observed Critical Shields Stress Variations

Estimates of the onset of sediment motion are integral for flood protection and river management but are often highly inaccurate. The critical shear stress (τ*c) for grain entrainment is often assumed constant, but measured values can vary by almost an order of magnitude between rivers. Such variations are typically explained by differences in measurement methodology, grain size distributions, or flow hydraulics, whereas grain resistance to motion is largely assumed to be constant. We demonstrate that grain resistance varies strongly with the bed structure, which is encapsulated by the particle height above surrounding sediment (protrusion,p) and intergranular friction (ϕf). We incorporate these parameters into a novel theory that correctly predicts resisting forces estimated in the laboratory, field, and a numerical model. Our theory challenges existing models, which significantly overestimate bed mobility. In our theory, small changes inpandϕfcan induce large changes inτ*cwithout needing to invoke variations in measurement methods or grain size. A data compilation also reveals that scatter in empirical values ofτ*ccan be partly explained by differences inpbetween rivers. Therefore, spatial and temporal variations in bed structure can partly explain the deviation ofτ*cfrom an assumed constant value. Given that bed structure is known to vary with applied shear stresses and upstream sediment supply, we conclude that a constantτ*cis unlikely. Values ofτ*care not interchangeable between streams, or even through time in a given stream, because they are encoded with the channel history.

more » « less
Award ID(s):
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
DOI PREFIX: 10.1029
Date Published:
Journal Name:
Journal of Geophysical Research: Earth Surface
Page Range / eLocation ID:
p. 3308-3322
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Steep, boulder bed streams often contain sediment patches, which are areas of the bed with relatively well‐defined boundaries that are occupied by distinct grain size distributions (GSD). In sediment mixtures, the underlying GSD affects the critical Shields stress for a given grain size, which is commonly predicted using hiding functions. Hiding functions may vary with reach‐wide bed GSD, but the effect of local GSD on relative sediment mobility between sediment patches is poorly understood. We explore the effects of patch‐scale GSD on sediment mobility using tracer particles combined with local shear stresses to develop hiding functions for different patch classes within a steep stream. Hiding functions for all tested patch classes were similar, which indicates that the same hiding function can be used for different patches. However, the critical Shields stress for a given grain size generally decreased with lower patch median grain size (D50) suggesting that patches control the relative mobility of each size through both the underlying GSD and local shear stresses. The effects of the underlying GSD partly depend on grain protrusion, which we measured for all grain sizes present on each patch class. Protrusion was generally greater for larger grains regardless of patch class, but for a given grain size, protrusion was increased with smaller patchD50. For a given grain size, higher protrusion results in greater applied fluid forces and reduced resisting forces to partly explain our lower critical Shields stresses in finer patches. Patches therefore can importantly modulate relative sediment mobility through bed structure and may need to be included in reach‐scale sediment transport and channel stability estimates.

    more » « less
  2. Abstract

    To explore the causes of history‐dependent sediment transport in rivers, we use a 19‐year record of coarse sediment transport from a steep channel in Switzerland. We observe a strong dependence of the threshold for sediment motion (τc) on the magnitude of previous flows for prior shear stresses ranging from 104 to 340 Pa, resulting in seasonally increasingτcfor 10 of 19 years. This stabilization occurs with and without measureable bedload transport, suggesting that small‐scale riverbed rearrangement increasesτc. Following large transport events (>340 Pa), this history dependence is disrupted. Bedload tracers suggest that significant reorganization of the bed erases memory of previous flows. We suggest that the magnitude of past flows controls the organization of the bed, which then modifiesτc, paralleling the evolution of granular media under shear. Our results support the use of a state function to better predict variability in bedload sediment transport rates.

    more » « less
  3. Abstract

    Flow resistance in mountain streams is important for assessing flooding hazard and quantifying sediment transport and bedrock incision in upland landscapes. In such settings, flow resistance is sensitive to grain‐scale roughness, which has traditionally been characterized by particle size distributions derived from laborious point counts of streambed sediment. Developing a general framework for rapid quantification of resistance in mountain streams is still a challenge. Here we present a semi‐automated workflow that combines millimeter‐ to centimeter‐scale structure‐from‐motion (SfM) photogrammetry surveys of bed topography and computational fluid dynamics (CFD) simulations to better evaluate surface roughness and rapidly quantify flow resistance in mountain streams. The workflow was applied to three field sites of gravel, cobble, and boulder‐bedded channels with a wide range of grain size, sorting, and shape. Large‐eddy simulations with body‐fitted meshes generated from SfM photogrammetry‐derived surfaces were performed to quantify flow resistance. The analysis of bed microtopography using a second‐order structure function identified three scaling regimes that corresponded to important roughness length scales and surface complexity contributing to flow resistance. The standard deviationσzof detrended streambed elevation normalized by water depth, as a proxy for the vertical roughness length scale, emerges as the primary control on flow resistance and is furthermore tied to the characteristic length scale of rough surface‐generated vortices. Horizontal length scales and surface complexity are secondary controls on flow resistance. A new resistance predictor linking water depth and vertical roughness scale, i.e. H/σz, is proposed based on the comparison betweenσzand the characteristic length scale of vortex shedding. In addition, representing streambeds using digital elevation models (DEM) is appropriate for well‐sorted streambeds, but not for poorly sorted ones under shallow and medium flow depth conditions due to the missing local overhanging features captured by fully 3D meshes which modulate local pressure gradient and thus bulk flow separation and pressure distribution. An appraisal of the mesh resolution effect on flow resistance shows that the SfM photogrammetry data resolution and the optimal CFD mesh size should be about 1/7 to 1/14 of the standard deviation of bed elevation. © 2019 John Wiley & Sons, Ltd.

    more » « less
  4. Abstract

    The impacts of aquatic vegetation on bed load transport rate and bedform characteristics were quantified using flume measurements with model emergent vegetation. First, a model for predicting the turbulent kinetic energy,kt, in vegetated channels from channel average velocityUand vegetation volume fractionϕwas validated for mobile sediment beds. Second, using data from several studies, the predictedktwas shown to be a good predictor of bed load transport rate,Qs, allowingQsto be predicted fromUandϕfor vegetated channels. The control ofQsbyktwas explained by statistics of individual grain motion recorded by a camera, which showed that the number of sediment grains in motion per bed area was correlated withkt. Third, ripples were observed and characterized in channels with and without model vegetation. For low vegetation solid volume fraction (ϕ ≤ 0.012), the ripple wavelength was constrained by stem spacing. However, at higher vegetation solid volume fraction (ϕ=0.025), distinct ripples were not observed, suggesting a transition to sheet flow, which is sediment transport over a plane bed without the formation of bedforms. The fraction of the bed load flux carried by migrating ripples decreased with increasingϕ, again suggesting that vegetation facilitated the formation of sheet flow.

    more » « less
  5. Abstract

    We present estimates of line-of-sight distortion fields derived from the 95 and 150 GHz data taken by BICEP2, BICEP3, and the Keck Array up to the 2018 observing season, leading to cosmological constraints and a study of instrumental and astrophysical systematics. Cosmological constraints are derived from three of the distortion fields concerning gravitational lensing from large-scale structure, polarization rotation from magnetic fields or an axion-like field, and the screening effect of patchy reionization. We measure an amplitude of the lensing power spectrumALϕϕ=0.95±0.20. We constrain polarization rotation, expressed as the coupling constant of a Chern–Simons electromagnetic termgaγ≤ 2.6 × 10−2/HI, whereHIis the inflationary Hubble parameter, and an amplitude of primordial magnetic fields smoothed over 1 MpcB1Mpc≤ 6.6 nG at 95 GHz. We constrain the rms of optical depth fluctuations in a simple “crinkly surface” model of patchy reionization, findingAτ< 0.19 (2σ) for the coherence scale ofLc= 100. We show that all of the distortion fields of the 95 and 150 GHz polarization maps are consistent with simulations including lensed ΛCDM, dust, and noise, with no evidence for instrumental systematics. In some cases, theEBandTBquadratic estimators presented here are more sensitive than our previous map-based null tests at identifying and rejecting spuriousB-modes that might arise from instrumental effects. Finally, we verify that the standard deprojection filtering in the BICEP/Keck data processing is effective at removing temperature to polarization leakage.

    more » « less