A major control on bedrock incision is the interaction between alluvial cover and erosive mobile grains. The extent of alluvial cover is typically predicted as a function of relative sediment flux (sediment supply rate over bed load transport capacity, q_bs/q_bc), yet little is known about how the bed roughness affects the alluvial cover. We performed field experiments with various flow discharges, sediment supply rates, grain sizes, and bed surface topographies. We then developed physically based models for estimating the threshold of sediment movement and the extent of alluvial cover, so as to include the effect of roughness change. The results for the threshold of sediment movement and the extent of alluvial cover obtained from our models show reasonable agreement with the results of the field experiments. We explored the sensitivity of the models to variations in sediment supply and bedrock relative roughness (bedrock hydraulic roughness height over grain size, k_sb/d). The results suggest the following: (1) a larger relative roughness yields a greater dimensionless critical shear stress required for initial sediment motion; (2) at a given sediment supply rate, the extent of alluvial cover is larger when the relative roughness is larger; (3) when the sediment supply rate and the relative roughness are small, throughput bed load moves over (and can abrade) a purely bedrock channel with no alluvial cover; and (4) the critical value of sediment supply rate below which throughput bed load transport occurs increases with decreasing relative roughness. The experimental results and analysis provide a framework for treating the (a) incisional morphodynamics of purely bedrock rivers by throughput bed load with no alluvial cover, (b) incisional/alluvial morphodynamics of mixed bedrock‐alluvial rivers, and (c) purely alluvial morphodynamics, as well as the transition between these states.

Results from computational morphodynamics modeling of coupled flow–bed–sediment systems are described for 10 applications as a review of recent advances in the field. Each of these applications is drawn from solvers included in the public‐domain International River Interface Cooperative (iRIC) software package. For mesoscale river features such as bars, predictions of alternate and higher mode river bars are shown for flows with equilibrium sediment supply and for a single case of oversupplied sediment. For microscale bed features such as bedforms, computational results are shown for the development and evolution of two‐dimensional bedforms using a simple closure‐based two‐dimensional model, for two‐ and three‐dimensional ripples and dunes using a three‐dimensional large‐eddy simulation flow model coupled to a physics‐based particle transport model, and for the development of bed streaks using a three‐dimensional unsteady Reynolds‐averaged Navier–Stokes solver with a simple sediment‐transport treatment. Finally, macroscale or channel evolution treatments are used to examine the temporal development of meandering channels, a failure model for cantilevered banks, the effect of bank vegetation on channel width, the development of channel networks in tidal systems, and the evolution of bedrock channels. In all examples, computational morphodynamics results from iRIC solvers compare well to observations of natural bed morphology. For each of the three scales investigated here, brief suggestions for future work and potential research directions are offered. © 2019 The Authors Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd

The 1-D saltation–abrasion model of channel bedrock incision of Sklar and Dietrich (2004), in which the erosion rate is buffered by the surface area fraction of bedrock covered by alluvium, was a major advance over models that treat river erosion as a function of bed slope and drainage area. Their model is, however, limited because it calculates bed cover in terms of bedload sediment supply rather than local bedload transport. It implicitly assumes that as sediment supply from upstream changes, the transport rate adjusts instantaneously everywhere downstream to match. This assumption is not valid in general, and thus can give rise to unphysical consequences. Here we present a unified morphodynamic formulation of both channel incision and alluviation that specifically tracks the spatiotemporal variation in both bedload transport and alluvial thickness. It does so by relating the bedrock cover fraction to the ratio of alluvium thickness to bedrock macro-roughness, rather than to the ratio of bedload supply rate to capacity bedload transport. The new formulation (MRSAA) predicts waves of alluviation and rarification, in addition to bedrock erosion. Embedded in it are three physical processes: alluvial diffusion, fast downstream advection of alluvial disturbances, and slow upstream migration of incisional disturbances. Solutions of this formulation over a fixed bed are used to demonstrate the stripping of an initial alluvial cover, the emplacement of alluvial cover over an initially bare bed and the advection–diffusion of a sediment pulse over an alluvial bed. A solution for alluvial–incisional interaction in a channel with a basement undergoing net rock uplift shows how an impulsive increase in sediment supply can quickly and completely bury the bedrock under thick alluvium, thus blocking bedrock erosion. As the river responds to rock uplift or base level fall, the transition point separating an alluvial reach upstream from an alluvial–bedrock reach downstream migrates upstream in the form of a "hidden knickpoint". A tectonically more complex case of rock uplift subject to a localized zone of subsidence (graben) yields a steady-state solution that is not attainable with the original saltation–abrasion model. A solution for the case of bedrock–alluvial coevolution upstream of an alluviated river mouth illustrates how the bedrock surface can be progressively buried not far below the alluvium. Because the model tracks the spatiotemporal variation in both bedload transport and alluvial thickness, it is applicable to the study of the incisional response of a river subject to temporally varying sediment supply. It thus has the potential to capture the response of an alluvial–bedrock river to massive impulsive sediment inputs associated with landslides or debris flows.

In this study, we conducted multiple physical experiments to estimate the efficacy and spatial pattern of erosion by abrading sediment moving through a simple U‐shaped channel bend with erodible bed and banks. The experiments showed that in the bend, lateral abrasion followed a monotonically increasing linear relationship with sediment feed rate. However, vertical incision had a more complex relation with the sediment feed rate, with an initial increase in abrasion as the feed rate increased followed by a decrease in abrasion of the bed as cover effects became dominant at higher feed rates. Bank erosion was large in places where the width and the lateral slope of the point bar were relatively large. On the other hand, in places where the width of the point bar was smaller, the bedrock bed was eroded primarily along the boundary of the point bar, resulting in a bedrock bench near the outer bank.