[1] This paper presents an idealized morphodynamic model to predict river dune evolution. The flow field is solved in a vertical plane assuming hydrostatic pressure conditions. The sediment transport is computed using a Meyer-Peter-Müller type of equation, including gravitational bed slope effects and a critical bed shear stress. To avoid the necessity of modeling the complex flow inside the flow separation zone, we follow an approach similar to one used earlier to simulate the dynamics of wind-blown desert dunes. In case of flow separation, the separation streamline acts as an artificial bed and sediment avalanches down the leeside distributing evenly on the leeside at the angle of repose. Model results show that bed slope effects play a crucial role in the determination of the fastest-growing wavelength from a linear analysis. Flow separation is shown to be crucial to take into account if the dune lee exceeds a certain threshold slope. If flow separation is not included, dune shapes are incorrectly predicted and the dune height saturates at an early stage of bed form evolution, yielding an underprediction of dune height and time to equilibrium. The local bed slope at the dune crest plays a critical role for obtaining an equilibrium dune height. The simulation model is able to predict the main characteristics of dune evolution, such as dune asymmetry, dune growth, and saturation at a certain dune height. Dune dimensions, migration rates, and times to equilibrium compare reasonably well to various data sets.

[1] Flow separation plays a key role in the development of dunes, and modeling the complicated flow behavior inside the flow separation zone requires much computational effort. To make a first step toward modeling dune development at reasonable temporal and spatial scales, a parameterization of the shape of the flow separation zone over twodimensional dunes is proposed herein, in order to avoid modeling the complex flow inside the flow separation zone. Flow separation behind dunes, with an angle-of-repose slip face, is characterized by a large circulating leeside eddy, where a separation streamline forms the upper boundary of the recirculating eddy. Experimental data of turbulent flow over two-dimensional subaqueous bed forms are used to parameterize this separation streamline. The bed forms have various heights and height to length ratios, and a wide range of flow conditions is analyzed. This paper shows that the shape of the flow separation zone can be approximated by a third-order polynomial as a function of the distance away from the flow separation point. The coefficients of the polynomial can be estimated, independent of flow conditions, on the basis of bed form shape at the flow separation point and a constant angle of the separation streamline at the flow reattachment point.

This paper presents an approach to incorporate time-dependent dune evolution in the determination of bed roughness coefficients applied in hydraulic models. Dune roughness is calculated by using the process-based dune evolution model of Paarlberg et al. (2009) and the empirical dune roughness predictor of Van Rijn (1984). The approach is illustrated by applying it to a river of simple geometry in the 1-D hydraulic model SOBEK for two different flood wave shapes. Calculated dune heights clearly show a dependency on rate of change in discharge with time: dunes grow to larger heights for a flood wave with a smaller rate of change. Bed roughness coefficients computed using the new approach can be up to 10% higher than roughness coefficients based on calibration, with the largest differences at low flows. As a result of this larger bed roughness, computed water depths can be up to 15% larger at low flow. The new approach helps to reduce uncertainties in bed roughness coefficients of flow models, especially for river systems with strong variations in discharge with time.