{A numerical scheme is developed to simulate the transport of natural gravel. Starting with computerized tomographic (CT) scans of natural grains, our method approximates the shapes of these grains by “gluing” spheres of different sizes together with overlaps. The conglomerated spheres move using a Discrete Element Method (DEM) which is coupled with a Lattice Boltzmann Method (LBM) fluid solver, forming the first complete workflow from particle shape measurement to high resolution simulations with hundreds of distinct shapes. The simulations are quantitatively benchmarked by flume experiments. The numerical tool is used to further validate a recently proposed modified sediment transport relation, which takes particle shape effects into account, including the competition between hydrodynamic drag and material friction. Unlike a physical experiment, our simulations allow us to vary the hydrodynamic drag coefficient of the natural gravel independently of the material friction. Our studies support the modified sediment transport relation. The simulations also provide insights on the particle-level kinematics, such as particle orientations, in the bedload transport process. Particles below the bed surface prefer to orient with their shortest axes perpendicular to the bed surface, but the tendency goes down as the packing fraction decreases far from the bed surface. The particles rotate freely in the dilute particle flow regime. }

In deltas and estuaries throughout the world, a fluvial-to-tidal transition zone (FTTZ) exists where both the river discharge and the tidal motion drive the flow. It is unclear how bedform characteristics are impacted by changes in tidal flow strength, and how this is reflected in the hydraulic roughness. To understand bedform geometry and variability in the FTTZ and possible impacts on hydraulic roughness, we assess dune variability from multibeam bathymetric surveys, and we use a calibrated 2D hydrodynamic model (Delft3D-FM) of a sand-bedded lowland river (Fraser River, Canada). We focus on a period of low river discharge during which tidal impact is strong. We find that the fluvial-tidal to tidal regime change is not directly reflected in dune height, but local patterns of increasing and decreasing dune height are present. The calibrated model is able to predict local patterns of dune heights using tidally-averaged values of bed shear stress. However, the spatially variable dune morphology hampers local dune height predictions. The fluvial-to-tidal regime change is reflected in dune shape, where dunes have lower leeside angles and are more symmetrical in the tidal regime. Those tidal effects do not significantly impact the reach-scale roughness, and predicted dune roughness using dune height and length is similar to the dune roughness inferred from model calibration. Hydraulic model performance with a calibrated, constant roughness is not improved by implementing dune-derived bed roughness. Instead, large-scale river morphology may explain differences in model roughness and corresponding estimates from dune predictors.

Near the threshold of grain motion, sediment transport is “on-off” intermittent, characterized by large but rare bursts separated by long periods of low transport. Without models that can predict the presence of intermittency, measurements of average sediment flux can be in error by up to an order of magnitude. Despite its known presence and impact, it is not clear whether on-off intermittency arises from the grain activity (the number of moving grains) or grain velocities, which together determine the sediment flux. We use laboratory flume experiments to show that the on-off intermittency has its origins in the velocity distributions of grains that move by rolling along the bed, whereas grain activity is not on-off intermittent. Improved predictions of sediment flux require that the types of intermittency we identify be incorporated into stochastic models of sediment flux. Their recognition opens the door to physically based uncertainty estimates of time-averaged sediment flux.

Dunes dominate the bed of sandy rivers and they respond to flow by changing shape and size, modifying flow and sediment transport dynamics of rivers. Our understanding and ability to predict dune adaptation, particularly dune growth and decay, remains incomplete. Here we investigate dune growth from an initial flatbed in a laboratory setting by continuously mapping the 3D bed topography using a line laser scanner combined with a 3D camera. High-resolution profiles of flow velocity and sediment concentration providing both bedload and suspended sediment fluxes were obtained by deploying Acoustic Concentration and Velocity Profiler technology. Our analysis reveals that the magnitude of the dune slipface angle, which determines flow separation and controls turbulence production, adjusts to the imposed flow at time scales similar to the evolution of dune height and length. The initiation of a flow separation zone intensifies through scour, and results in acceleration of the dune growth. Gradients in sediment transport and the rate of dune growth are inherently linked to spatial variations in slipface angles. During dune growth, the slipface angle evolves differently than the ratio of dune height to length, which immediately reaches its equilibrium value after dune initiation.

Sediment transport in rivers near the threshold of grain motion is characterized by rare but large transport events. This intermittency makes it difficult to relate average sediment flux to average flow conditions, or to define an unambiguous threshold for grain entrainment. Although intermittent sediment transport can be observed and characterized, its origins are unclear. In this study we investigate bedload sediment transport near the threshold of grain motion in an experimental flume to examine the origins of intermittency. We apply image-processing techniques to high-speed video of grains in a narrow flume, which allows us to track individual particles and measure statistics of particle motion. Bedload sediment transport near the threshold of grain motion is very low, allowing us to approximate the time evolution of the sediment flux via a polynomial expansion, including a linear growth rate and a nonlinear term which saturates the growth. We introduce a noisy coefficient to the linear growth rate term (“multiplicative noise”), rather than adding the noise to the equation, to model the inherent fluctuations in the system. We demonstrate that multiplicative noise near the threshold of grain motion can account for the observed intermittency. We use analytical results from bifurcation theory in the presence of multiplicative noise to analyze our experimental results, quantifying the noise responsible for the intermittency and calculating the critical shear stress for grain entrainment in a novel way that is consistent with the physics of grain motion at low transport stages.

Bedrock rivers get wider by lateral erosion. Lateral erosion is widely thought to occur when the bed is covered by alluvium, which deflects the downstream transport of bedload particles into channel walls. Here we develop a model for lateral bedrock erosion by bedload particle impacts. The lateral erosion rate is the product of the volume eroded per particle impact and the impact rate on the wall. The volume eroded per particle impact is modelled by tracking the motion of bedload particles from collision with roughness elements to impacts on the wall. The impact rate on the wall is zero if the bedload particle deflected by roughness elements cannot reach the wall. Otherwise, the impact rate on the wall is the same with that on roughness elements. The model further incorporates the co-evolution of wall morphology, shear stress and erosion rate. The model predicts the undercut wall shape observed in physical experiments. The non-dimensional lateral erosion rate is used to explore how lateral erosion varies under different relative sediment supply (ratio of supply to transport capacity) and transport stage conditions. Maximum lateral erosion rates occur at high relative sediment supply rates (~ 0.7) and moderate transport stages (~10). The competition between lateral and vertical erosion is investigated by coupling the saltation-abrasion vertical erosion model with our lateral erosion model. The results suggest that vertical erosion dominates under near 75% of supply and transport stage conditions, but is outpaced by lateral erosion near the threshold for full bed coverage.