Spatial variations in surface sediment structure in riffle–pool sequences: a preliminary test of the Differential Sediment Entrainment Hypothesis (DSEH)

Riffle–pool sequences are maintained through the preferential entrainment of sediment grains from pools rather than riffles. This preferential entrainment has been attributed to a reversal in the magnitude of velocity and shear stress under high flows; however the Differential Sediment Entrainment Hypothesis (DSEH) postulates that differential entrainment can instead result from spatial sedimentological contrasts. Here we use a novel suite of in situ grain‐scale field measurements from a riffle–pool sequence to parameterize a physically‐based model of grain entrainment. Field measurements include pivoting angles, lift forces and high resolution digital elevation models (DEMs) acquired using terrestrial laser scanning, from which particle exposure, protrusion and surface roughness were derived. The entrainment model results show that grains in pools have a lower critical entrainment shear stress than grains in either pool exits or riffles. This is because pool grains have looser packing, hence greater exposure and lower pivoting angles. Conversely, riffle and pool exit grains have denser packing, lower exposure and higher pivoting angles. A cohesive matrix further stabilizes pool exit grains. The resulting predictions of critical entrainment shear stress for grains in different subunits are compared with spatial patterns of bed shear stress derived from a two‐dimensional computational fluid dynamics (CFD) model of the reach. The CFD model predicts that, under bankfull conditions, pools experience lower shear stresses than riffles and pool exits. However, the difference in sediment entrainment shear stress is sufficiently large that sediment in pools is still more likely to be entrained than sediment in pool exits or riffles, resulting in differential entrainment under bankfull flows. Significantly, this differential entrainment does not require a reversal in flow velocities or shear stress, suggesting that sedimentological contrasts alone may be sufficient for the maintenance of riffle–pool sequences. This finding has implications for the prediction of sediment transport and the morphological evolution of gravel‐bed rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Turbulence characteristics of flows passing through a tetrahedron frame in a smooth open-channel

The turbulence characteristics of flows passing through a tetrahedron frame were investigated by using a 2-dimensional fiber-optic laser Doppler velocimeter (2-D FLDV). Experiments for uniform flows with different bed slopes under both submerged and un-submerged conditions were carried out in a re-circulating flume with glass side walls. The experimental bed was a smooth fixed bed. It was observed that with the tetrahedron frame the mean longitudinal velocity decrease in the retardation zone. However, both the longitudinal and the vertical turbulence intensities are larger than those for the undisturbed approach flow. The tetrahedron frame may reduce the probability of sediment entrainment by retarding the flow and reducing the boundary shear stress. In addition, it may induce sediment deposition in a sediment laden flow by changing the flow direction and increasing the energy dissipation.     

An analysis of the characteristics of rough bed turbulent shear stresses in an open channel

Entrainment of sediment particles from channel beds into the channel flow is influenced by the characteristics of the flow turbulence which produces stochastic shear stress fluctuations at the bed. Recent studies of the structure of turbulent flow has recognized the importance of bursting processes as important mechanisms for the transfer of momentum into the laminar boundary layer. Of these processes, the sweep event has been recognized as the most important bursting event for entrainment of sediment particles as it imposes forces in the direction of the flow resulting in movement of particles by rolling, sliding and occasionally saltating. Similarly, the ejection event has been recognized as important for sediment transport since these events maintain the sediment particles in suspension. In this study, the characteristics of bursting processes and, in particular, the sweep event were investigated in a flume with a rough bed. The instantaneous velocity fluctuations of the flow were measured in two-dimensions using a small electromagnetic velocity meter and the turbulent shear stresses were determined from these velocity fluctuations. It was found that the shear stress applied to the sediment particles on the bed resulting from sweep events depends on the magnitude of the turbulent shear stress and its probability distribution. A statistical analysis of the experimental data was undertaken and it was found necessary to apply a Box-Cox transformation to transform the data into a normally distributed sample. This enabled determination of the mean shear stress, angle of action and standard error of estimate for sweep and ejection events. These instantaneous shear stresses were found to be greater than the mean flow shear stress and for the sweep event to be approximately 40 percent greater near the channel bed. Results from this analysis suggest that the critical shear stress determined from Shield's diagram is not sufficient to predict the initiation of motion due to its use of the temporal mean shear stress. It is suggested that initiation of particle motion, but not continuous motion, can occur earlier than suggested by Shield's diagram due to the higher shear stresses imposed on the particles by the stochastic shear stresses resulting from turbulence within the flow.     

Aerodynamic entrainment thresholds and dislodgement rates on impervious and permeable beds

In the literature it has been suggested that on permeable, granular beds, both the threshold and rate of aerodynamic entrainment may be affected significantly by seepage flows into and out of the bed induced by fluctuating pressures in the overlying turbulent boundary layer. Using a range of grain sizes and flow conditions, the series of laboratory experiments reported here compares directly the aerodynamic entrainment of loose grains overlying fixed permeable sediment beds with that occuring over fixed impervious beds. For a given granular material, no significant differences in entrainment dynamics on the two types of bed were observed and in the range of flow conditions examined both the threshold shear velocity (U _*T) and the aerodynamic entrainment coefficient (k) were found to be independent of bed permeability.     

Transport of moderately sorted gravel at low bed shear stresses: The role of fine sediment infiltration

A reliable estimation of sediment transport in gravel‐bed streams is important for various practical engineering and biological studies (e.g., channel stability design, bed degradation/aggradation, restoration of spawning habitat). In the present work, we report original laboratory experiments investigating the transport of gravel particles at low bed shear stresses. The laboratory tests were conducted under unsteady flow conditions inducing low bed shear stresses, with detailed monitoring of the bed topography using a laser scanner. Effects of bed surface arrangements were documented by testing loose and packed bed configurations. Effects of fine sediments were examined by testing beds with sand, artificial fine sand or cohesive silt infiltrated in the gravel matrix. Analysis of the experimental data revealed that the transport of gravel particles depends upon the bed arrangement, the bed material properties (e.g., size and shape, consolidation index, permeability) and the concentration of fine sediments within the surface layer of moving grains. This concentration is directly related to the distribution of fine particles within the gravel matrix (i.e., bottom‐up infiltration or bridging) and their transport mode (i.e., bedload or suspended load). Compared to loose beds, the mobility of gravel is reduced for packed beds and for beds clogged from the bottom up with cohesive fine sediments; in both cases, the bed shear stress for gravel entrainment increases by about 12%. On the other hand, the mobility of gravel increases significantly (bed shear stress for particle motion decreasing up to 40%) for beds clogged at the surface by non‐cohesive sand particles. Copyright © 2017 John Wiley & Sons, Ltd.     

Bed scour by debris flows: experimental investigation of effects of debris‐flow composition

Debris flows can grow greatly in size by entrainment of bed material, enhancing their runout and hazardous impact. Here, we experimentally investigate the effects of debris‐flow composition on the amount and spatial patterns of bed scour and erosion downstream of a fixed to erodible bed transition. The experimental debris flows were observed to entrain bed particles both grain by grain and en masse, and the majority of entrainment was observed to occur during passage of the flow front. The spatial bed scour patterns are highly variable, but large‐scale patterns are largely similar over 22.5–35° channel slopes for debris flows of similar composition. Scour depth is generally largest slightly downstream of the fixed to erodible bed transition, except for clay‐rich debris flows, which cause a relatively uniform scour pattern. The spatial variability in the scour depth decreases with increasing water, gravel (= grain size) and clay fraction. Basal scour depth increases with channel slope, flow velocity, flow depth, discharge and shear stress in our experiments, whereas there is no correlation with grain collisional stress. The strongest correlation is between basal scour and shear stress and discharge. There are substantial differences in the scour caused by different types of debris flows. In general, mean and maximum scour depths become larger with increasing water fraction and grain size, and decrease with increasing clay content. However, the erodibility of coarse‐grained experimental debris flows (gravel fraction = 0.64) is similar on a wide range of channel slopes, flow depths, flow velocities, discharges and shear stresses. This probably relates to the relatively large influence of grain‐collisional stress to the total bed stress in these flows (30–50%). The relative effect of grain‐collisional stress is low in the other experimental debris flows (<5%), causing erosion to be largely controlled by basal shear stress. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of bedload transport in laminar flows

Bedload transport generally depends on the bed shear stress and Reynolds number. Many studies conducted for the condition of turbulent flows have revealed the dependence of the transport rate on the bed shear stress, while knowledge of the Reynolds number effect on the transport rate is very limited. As an extreme case to reflect the viscous effect on sediment transport, sediment transport in laminar flows is considered in this paper. A stochastic approach is adopted to explore how the transport rate can be associated with characteristics of laminar flows. First, the probability of erosion in the absence of turbulence is assumed to depend only on the randomness of bed particles. The probability is then applied to formulate the sediment transport rate, of which the derivation is made largely based on Einstein’s bedload theory. The theoretical result indicates that the dimensionless transport rate for laminar flows is dependent on the dimensionless shear stress and dimensionless particle diameter or the shear Reynolds number. Comparisons are finally made between the derived formula and an empirical correlation available in the literature.     

EXPERIMENTAL STUDY ON INCIPIENT MOTION OF SEDIMENT PARTICLES ON GENERALIZED SLOPING FLUVIAL BEDS

1 INTRODUCTIONWhen water flows over a fluvial bed, hydro-dynandc force induced by the flow is acting on thesediment particles lying on the bed. A further increase in flow velocity results in an increase in themagnitude of this fOrce; and sediment particles begin to move if a situation is eventu8lly reached whenthe hydro-dynandc force exceeds a certain critical value. This initial movement of sediment pallicles istermed inciPient motion. The erosion and sedimentation of nuvial beds can be…     

Three‐dimensional flow structure and patterns of bed shear stress in an evolving compound meander bend

Compound meander bends with multiple lobes of maximum curvature are common in actively evolving lowland rivers. Interaction among spatial patterns of mean flow, turbulence, bed morphology, bank failures and channel migration in compound bends is poorly understood. In this paper, acoustic Doppler current profiler (ADCP) measurements of the three‐dimensional (3D) flow velocities in a compound bend are examined to evaluate the influence of channel curvature and hydrologic variability on the structure of flow within the bend. Flow structure at various flow stages is related to changes in bed morphology over the study timeframe. Increases in local curvature within the upstream lobe of the bend reduce outer bank velocities at morphologically significant flows, creating a region that protects the bank from high momentum flow and high bed shear stresses. The dimensionless radius of curvature in the upstream lobe is one‐third less than that of the downstream lobe, with average bank erosion rates less than half of the erosion rates for the downstream lobe. Higher bank erosion rates within the downstream lobe correspond to the shift in a core of high velocity and bed shear stresses toward the outer bank as flow moves through the two lobes. These erosion patterns provide a mechanism for continued migration of the downstream lobe in the near future. Bed material size distributions within the bend correspond to spatial patterns of bed shear stress magnitudes, indicating that bed material sorting within the bend is governed by bed shear stress. Results suggest that patterns of flow, sediment entrainment, and planform evolution in compound meander bends are more complex than in simple meander bends. Moreover, interactions among local influences on the flow, such as woody debris, local topographic steering, and locally high curvature, tend to cause compound bends to evolve toward increasing planform complexity over time rather than stable configurations. Copyright © 2016 John Wiley & Sons, Ltd.     

Development of bedforms due to waves blocked by a counter current

Experiments are conducted in a laboratory flume on the propagation of a surface wave against unidirectional flow with a sediment bed. This article presents the spatial variation of bedforms induced by the wave-blocking phenomenon by a suitably tuned uniform fluid flow and a counter-propagating wave. The occurrence of wave-blocking is confirmed by finding a critical wave frequency in a particular flow discharge in which the waves are effectively blocked and is established using the linear dispersion relation. The purpose of this work is to identify wave-blocking and its influence on the development of bedforms over the sediment bed. Interestingly bedform signatures are observed at a transition of bedforms in three zones, with asymmetric ripples having a steeper slope downstream face induced by the incoming current, followed by flat sand bars beneath the wave-blocking zone and more symmetric ripples below the wave-dominated region at the downstream. This phenomenon suggests that the sediment bed is segmented into three different regions of bed geometry along the flow. The deviations of mean flows, Reynolds stresses, turbulent kinetic energy, and power spectral density due to the wave-blocking phenomenon are presented along the non-uniform flow over sediment bed. The bottom shear stress, bed roughness and stochastic nature of the bedform features are also discussed. The results are of relevance to engineers and geoscientists concerned with contemporary process as well as those interested in the interpretation of palaeoenvironmental conditions from fossil bedforms. © 2019 John Wiley & Sons, Ltd.     

A probabilistic model for sediment entrainment:The role of bed irregularity

A generalized probabilistic model is developed in this study to predict sediment entrainment under the incipient motion, rolling, and pickup modes. A novelty of the proposed model is that it incorporates in its formulation the probability density function of the bed shear stress, instead of the near-bed velocity fluctuations, to account for the effects of both flow turbulence and bed surface irregularity on sediment entrainment. The proposed model incorporates in its formulation the collective effects of three para-meters describing bed surface irregularity, namely the relative roughness, the volumetric fraction and relative position of sediment particles within the active layer. Another key feature of the model is that it provides a criterion for estimating the lift and drag coefficients jointly based on the recognition that lift and drag forces acting on sediment particles are interdependent and vary with particle protrusion and packing density. The model was validated using laboratory data of both fine and coarse sediment and was compared with previously published models. The study results show that all the examined models perform adequately for the fine sediment data, where the sediment particles have more uniform gra-dation and relative roughness is not a factor. The proposed model was particularly suited for the coarse sediment data, where the increased bed irregularity was captured by the new parameters introduced in the model formulation. As a result, the proposed model yielded smaller prediction errors and physically acceptable values for the lift coefficient compared to the other models in case of the coarse sediment data.     

Formation and erosion of sediment cover in an experimental bedrock‐alluvial channel

Sediment grains in a bedrock‐alluvial river will be deposited within or adjacent to a sediment patch, or as isolated grains on the bedrock surface. Previous analysis of grain geometry has demonstrated that these arrangements produce significant differences in grain entrainment shear stress. However, this analysis neglected potential interactions between the sediment patches, local hydraulics and grain entrainment. We present a series of flume experiments that measure the influence of sediment patches on grain entrainment. The flume had a planar bed with roughness that was much smaller than the diameters of the mobile grains. In each experiment sediment was added either as individual grains or as a single sediment pulse. Flow was then increased until the sediment was entrained. Analysis of the experiments demonstrates that: (1) for individual grains, coarse grains are entrained at a higher discharge than fine grains; (2) once sediment patches are present, the different in entrainment discharge between coarse and fine grains is greatly reduced; (3) the sheltering effect of patches also increases the entrainment discharge of isolated grains; (4) entire sediment patches break‐up and are eroded quickly, rather than through progressive grain‐by‐grain erosion; (5) as discharge increases there is some tendency for patches to become more elongate and flow‐aligned, and more randomly distributed across the bed. One implication of this research is that the critical shear stress in bedrock‐alluvial channels will be a function of the extent of the sediment cover. Another is that the influence of sediment patches equalizes critical shear stresses between different grain sizes and grain locations, meaning that these factors may not need to be accounted for. Further research is needed to quantify interactions between sediment patches, grain entrainment and local hydraulics on rougher bedrock surfaces, and under different types of sediment supply. Copyright © 2016 John Wiley & Sons, Ltd.     

EXPERIMENTAL STUDY ON BED SHEAR STRESS AND SEDIMENT TRANSPORT IN A MEANDERING LABORATORY CHANNEL

1 INTRODUCTION River erosion is a complex phenomenon. The rate of bank retreat is determined by flow, bed topography, sediment transport, bank properties, and water quality. Prediction of future river planform changes and the knowledge of river erosion and river meandering are required for land use planning in alluvial river valleys and determining locations for bridges and hydraulic structures. The control of riverbank erosion requires prediction of flow and bed features in a meanderin…     

Experimental investigation on local shear stress and turbulence intensities over a rough non-uniform bed with and without sediment using 2D particle image velocimetry

The current work focuses on locally resolving velocities,turbulence,and shear stresses over a rough bed with locally non-uniform character.A nonporous subsurface layer and fixed interfacial sublayer of gravel and sand were water-worked to a nature-like bed form and additionally sealed in a hydraulic flume.Two-dimensional Particle Image Velocimetry(2 D-PIV) was applied in the vertical plane of the experimental flume axis.Runs with clear water and weak sediment transport were done under slightly supercritical flow to ensure sediment transport conditions without formation of considerable sediment deposits or dunes.The study design included analyzing the double-averaged flow parameters of the entire measurement domain and investigating the flow development at 14 consecutive vertical subsections.Local geometrical variabilities as well the presence of sediment were mainly reflected in the vertical velocity component.Whereas the vertical velocity decreased over the entire depth in presence of sediment transport,the streamwise velocity profile was reduced only within the interfacial sublayer.In the region with decelerating flow conditions,however,the streamwise velocity profile systematically increased along the entire depth extent.The increase in the main velocity(reduction of flow resistance)correlated with a decrease of the turbulent shear and main normal stresses.Therefore,effects of rough bed smoothening and drag force reduction were experimentally documented within the interfacial sublayer due to mobile sediment.Moreover,the current study leads to the conclusion that in nonuniform flows the maximum Reynolds stress values are a better predictor for the bed shear stress than the linearly extrapolated Reynolds stress profile.This is an important finding because,in natural flows,uniform conditions are rare.     

Examining the physical components of boundary shear stress for water‐worked gravel deposits

It is argued in this commentary that, in order to understand better the physical mechanisms that generate boundary shear stress over water‐worked gravel beds, flow velocity data should be re‐evaluated by spatial averaging the Reynolds equations to produce time‐ and space‐averaged (double‐averaged) momentum equations. A series of laboratory experiments were conducted in which the flow velocities were measured using a PIV system over two water‐worked gravel deposits. Combined with detailed data on the bed surface topography and vertical porosity, the physical components of shear stress were obtained. This enabled the various momentum transfer mechanisms present above, within and at the interface of a porous, fluvial deposit, to be quantified. This included the examination of the relevant contributions of temporal and spatial fluctuations in velocity and surface drag to the overall momentum transfer. It is demonstrated that double‐averaging represents a logical framework for assessing the fluid forces responsible for sediment entrainment and for investigating intragravel flow and sediment–water interface exchange mechanisms within the roughness layer in water‐worked gravel deposits. By considering the physical components of shear stress and their relative sizes it was possible to provide a physically based explanation for existing observations of enhanced mobility of gravel–sand mixtures and the transfer of solutes into porous, gravel deposits. This analysis reveals the importance of obtaining co‐located, high quality spatial data on the flow field and bed surface topography in order to gain a physical understanding of the mechanisms which generate boundary shear stress. Copyright © 2010 John Wiley & Sons, Ltd.     

Assessing entrainment of bed material in a debris‐flow event: a theoretical approach incorporating Monte Carlo method

Entrainment of underlying bed sediment by a debris flow can significantly increase the debris‐flow magnitude. To study this phenomenon, a theoretical approach to assessing bed‐sediment entrainment is presented. The approach is based on a static approximation that bed‐sediment entrainment occurs when the shearing stress of the flow is sufficiently high to overcome the basal resistance of the bed sediment. In order to delineate erodible zones in a channel, we analyze the critical condition of this static equilibrium model, and subsequently propose a new concept of a critical line to detect the entrainment reaches in a channel. Considering the spatial and temporal uncertainties of the input parameter, the approach is further incorporated within a Monte Carlo method, and the distribution of entrainment zones and post‐entrainment volumes can be analyzed. This approach is illustrated by back‐analysis of the 2010 Yohutagawa debris‐flow event, Japan. Results from 10 000 trials of Monte Carlo simulation are compared with the in situ surveys. It is shown that the present approach can be satisfactorily used to delineate erodible zones and estimate possible entrainment volume of the event. Discussion regarding the sensitivities and limitations of the approach concludes the paper. Copyright © 2015 John Wiley & Sons, Ltd.     

Deriving formulas for an unsteady virtual velocity of bedload tracers

In a flume experiment with steady flow conditions, H. A. Einstein recognised the transport of bedload particles as consisting of steps of rolling, sliding, or saltation with intermittent rest periods, and introduced the concept of an average, ‘virtual’ transport velocity. This virtual velocity then has also been derived from tracer studies in the field by dividing the travelled distance of a tracer by the duration of competent flow. As a consequence, the virtual velocity in the field is represented by one single value only, despite the unsteady flow variables. Tracer measurements in a river have not been yet used to express transport velocity as a direct function of these actual variables, and insights from tracer measurements into the processes of sediment transport remain limited. In particular, the unsteady conditions for bedload in the field have impeded the derivation of sediment transport characteristics as determined from laboratory experiments, as well as the transfer of laboratory insights to a field setting. We introduce a method of data regression for the derivation of an ‘unsteady’ virtual velocity from repeated surveys of tracer positions. The regression program called graVel (provided as supplementary material) relates the integral of an excess flow variable term to measured travel distances, yielding the most probable threshold value for entrainment and the coefficient of linear and non‐linear formulas. An extended regression allows additional fitting of the exponent in non‐linear formulas. Application to published tracer data from the Mameyes River, Puerto Rico, shows that the unsteady virtual velocity is more likely governed by non‐linear relations to excess Shields stress, similar to bedload transport, than by relations linking the particle velocity linearly to excess shear velocity. Partial agreements with non‐dimensional results derived from the larger, non‐wadeable Mur River encourage the establishment of a generalised formula for the unsteady virtual velocity. Copyright © 2017 John Wiley & Sons, Ltd.     

Surficial sediment erodibility from time-series measurements of suspended sediment concentrations: development and validation

Numerical models of fine sediment transport depend on different approaches to parameterize the erosion properties of surficial sediment strata. These properties, namely the critical shear stress for erosion and the erosion rate coefficient, are crucial for reproducing the short-term and long-term sediment dynamics of the system. Methods to parameterize these properties involve either specialized laboratory measurements on sediment samples or optimization by model calibration. Based on observations of regular patterns in the variation of suspended sediment concentrations (SSC) over the tidal cycle in a small, narrow estuary, an alternate approach, referred to as the entrainment flux method, for quantifying the erosion properties of surficial bed strata is formulated and applied. The results of this method are shown to be analogous to the erosion data used to formulate the standard linear erosion formulation developed by various authors. The erosion properties inferred from the entrainment flux method are also compared to direct measurements of erodibility on sediment samples from the same site using the Gust microcosm apparatus. The favorable comparison of the two approaches suggests that the entrainment flux method can be used to infer and quantify the erodibility of surficial sediment strata in similar small and narrow estuaries. This method has certain advantages, chiefly its ease of implementation and the fact that it uses SSC time series which would typically be expected to be available for the study of or for model application at a given site. Guidelines for selecting the appropriate dataset for the application of the method are also presented.     

The impact of inter-flood duration on non-cohesive sediment bed stability

Limited field and flume data suggests that both uniform and graded beds appear to progressively stabilize when subjected to inter-flood flows as characterized by the absence of active bedload transport. Previous work has shown that the degree of bed stabilization scales with duration of inter-flood flow, however, the sensitivity of this response to bed surface grain size distribution has not been explored. This article presents the first detailed comparison of the dependence of graded bed stability on inter-flood flow duration. Sixty discrete experiments, including repetitions, were undertaken using three grain size distributions of identical D₅₀ (4.8 mm); near-uniform (σ_g = 1.13), unimodal (σ_g = 1.63) and bimodal (σ_g = 2.08). Each bed was conditioned for between 0 (benchmark) and 960 minutes by an antecedent shear stress below the entrainment threshold of the bed (τ*_c50). The degree of bed stabilization was determined by measuring changes to critical entrainment thresholds and bedload flux characteristics. Results show that (i) increasing inter-flood duration from 0 to 960 minutes increases the average threshold shear stress of the D₅₀ by up to 18%; (ii) bedload transport rates were reduced by up to 90% as inter-flood duration increased from 0 to 960 minutes; (iii) the rate of response to changes in inter-flood duration in both critical shear stress and bedload transport rate is non-linear and is inversely proportional to antecedent duration; (iv) there is a grade dependent response to changes in critical shear stress where the magnitude of response in uniform beds is up to twice that of the graded beds; and (v) there is a grade dependent response to changes in bedload transport rate where the bimodal bed is most responsive in terms of the magnitude of change. These advances underpin the development of more accurate predictions of both entrainment thresholds and bedload flux timing and magnitude, as well as having implications for the management of environmental flow design. © 2019 John Wiley & Sons, Ltd. © 2019 John Wiley & Sons, Ltd.     

Modelling dam-break flows over mobile beds using a 2D coupled approach

Dam-break flows usually propagate along rivers and floodplains, where the processes of fluid flow, sediment transport and bed evolution are closely linked. However, the majority of existing two-dimensional (2D) models used to simulate dam-break flows are only applicable to fixed beds. Details are given in this paper of the development of a 2D morphodynamic model for predicting dam-break flows over mobile beds. In this model, the common 2D shallow water equations are modified, so that the effects of sediment concentrations and bed evolution on the flood wave propagation can be considered. These equations are used together with the non-equilibrium transport equations for graded sediments and the equation of bed evolution. The governing equations are solved using a matrix method, thus the hydrodynamic, sediment transport and morphological processes can be jointly solved. The model employs an unstructured finite volume algorithm, with an approximate Riemann solver, based on the Roe-MUSCL scheme. A predictor–corrector scheme is used in time stepping, leading to a second-order accurate solution in both time and space. In addition, the model considers the adjustment process of bed material composition during the morphological evolution process. The model was first verified against results from existing numerical models and laboratory experiments. It was then used to simulate dam-break flows over a fixed bed and a mobile bed to examine the differences in the predicted flood wave speed and depth. The effects of bed material size distributions on the flood flow and bed evolution were also investigated. The results indicate that there is a great difference between the dam-break flow predictions made over a fixed bed and a mobile bed. At the initial stage of a dam-break flow, the rate of bed evolution could be comparable to that of water depth change. Therefore, it is often necessary to employ the turbid water governing equations using a coupled approach for simulating dam-break flows.