Two‐dimensional and three‐dimensional computational models in hydrodynamic and morphodynamic reconstructions of a river bend: sensitivity and functionality

This study assesses hydrodynamic and morphodynamic model sensitivity and functionality in a curved channel. The sensitivity of a depth‐averaged model to user‐defined parameters (grain size, roughness, transverse bed slope effect, transport relations and secondary flow) is tested. According to the sensitivity analysis, grain size, transverse bed slope effect and sediment transport relations are critical to simulated meander bend morphodynamics. The parametrization of grain size has the most remarkable effect: field‐based grain size parametrization is necessary in a successful morphodynamic reconstruction of a meander bend. The roughness parametrization method affects the distribution of flow velocities and therefore also morphodynamics. The combined effect of various parameters needs further research. Two‐dimensional (2D) and three‐dimensional (3D) reconstructions of a natural meander bend during a flood event are assessed against field measurements of acoustic Doppler current profiler and multi‐temporal mobile laser scanning data. The depth‐averaged velocities are simulated satisfactorily (differences from acoustic Doppler current profiler velocities 5–14%) in both 2D and 3D simulations, but the advantage of the 3D hydrodynamic model is unquestionable because of its ability to model vertical and near‐bed flows. The measured and modelled near‐bed flow, however, differed notably from each other's, the reason of which was left open for future research. It was challenging to model flow direction beyond the apex. The 3D flow features, which also affected the distribution of the bed shear stress, seem not to have much effect on the predicted morphodynamics: the 2D and 3D morphodynamic reconstructions over the point bar resembled each other closely. Although common features between the modelled and measured morphological changes were also found, some specific changes that occurred were not evident in the simulation results. Our results show that short‐term, sub‐bend scale morphodynamic processes of a natural meander bend are challenging to model, which implies that they are affected by factors that have been neglected in the simulations. The modelling of short‐term morphodynamics in natural curved channel is a challenge that requires further study. Copyright © 2014 John Wiley & Sons, Ltd.     

A hydro‐economic modelling framework for flood damage estimation and the role of riparian vegetation

A modelling framework for the quick estimate of flood inundation and the resultant damages is developed in this paper. The model, called the flood economic impact analysis system (FEIAS), can be applied to a river reach of any hydrogeological river basin. For the development of the integrated modelling framework, three models were employed: (1) a modelling scheme based on the Hydrological Simulation Program FORTRAN model that was developed for any geomorphological river basin, (2) a river flow/floodplain model, and (3) a flood loss estimation model. The first sub‐model of the flood economic impact analysis system simulates the hydrological processes for extended periods of time, and its output is used as input to a second component, the river/floodplain model. The hydraulic model MIKE 11 (quasi‐2D) is the river/floodplain model employed in this study. The simulated flood parameters from the hydraulic model MIKE 11 (quasi‐2D) are passed, at the end of each time step, to a third component, the flood loss model for the estimation of flood damage. In the present work, emphasis was given to the seasonal variation of Manning's coefficient (n), which is an important parameter for the determination of the flood inundation in hydraulic modelling. High values of Manning's coefficient for a channel indicate high flow resistance. The riparian vegetation can have a large impact on channel resistance. The modelling framework developed in this paper was used to investigate the role of riparian vegetation in reducing flood damage. Moreover, it was used to investigate the influence of cutting riparian vegetation scenarios on the flow characteristics. The proposed framework was applied to the downstream part of the Koiliaris River basin in Crete, Greece, and was tested and validated with historical data. Copyright © 2012 John Wiley & Sons, Ltd.     

Numerical simulation of overbank processes in topographically complex floodplain environments

This article presents results from an investigation of the hydraulic characteristics of overbank flows on topographically‐complex natural river floodplains. A two‐dimensional hydraulic model that solves the depth‐averaged shallow water form of the Navier–Stokes equations is used to simulate an overbank flow event within a multiple channel reach of the River Culm, Devon, UK. Parameterization of channel and floodplain roughness by the model is evaluated using monitored records of main channel water level and point measurements of floodplain flow depth and unit discharge. Modelled inundation extents and sequences are assessed using maps of actual inundation patterns obtained using a Global Positioning System, observational evidence and ground photographs. Simulation results suggest a two‐phase model of flooding at the site, which seems likely to be representative of natural floodplains in general. Comparison of these results with previous research demonstrates the complexity of overbank flows on natural river floodplains and highlights the limitations of laboratory flumes as an analogue for these environments. Despite this complexity, frequency distributions of simulated depth, velocity and unit discharge data closely follow a simple gamma distribution model, and are described by a shape parameter (α) that exhibits clear systematic trends with changing discharge and floodplain roughness. Such statistical approaches have the potential to provide the basis for computationally efficient flood routing and overbank sedimentation models. Copyright © 2002 John Wiley & Sons, Ltd.     

A coupled vertically integrated model to describe lateral exchanges between surface and subsurface in large alluvial floodplains with a fully penetrating river

This paper presents a vertically averaged model for studying water and solute exchanges between a large river and its adjacent alluvial aquifer. The hydraulic model couples horizontal 2D Saint Venant equations for river flow and a 2D Dupuit equation for aquifer flow. The dynamic coupling between river and aquifer is provided by continuity of fluxes and water level elevation between the two domains. Equations are solved simultaneously by linking the two hydrological system matrices in a single global matrix in order to ensure the continuity conditions between river and aquifer and to accurately model two‐way coupling between these two domains. The model is applied to a large reach (about 36 km²) of the Garonne River (south‐western France) and its floodplain, including an instrumented site in a meander. Simulated hydraulic heads are compared with experimental measurements on the Garonne River and aquifer in the floodplain. Model verification includes comparisons for one point sampling date (27 piezometers, 30 March 2000) and for hydraulic heads variations measured continuously over 5 months (5 piezometers, 1 January to 1 June 2000). The model accurately reproduces the strong hydraulic connections between the Garonne River and its aquifer, which are confirmed by the simultaneous variation of the water level in the river and in piezometers located near the river bank. The simulations also confirmed that the model is able to reproduce groundwater flow dynamics during flood events. Given these results, the hydraulic model was coupled with a solute‐transport component, based on advection‐dispersion equations, to investigate the theoretical dynamics of a conservative tracer over 5 years throughout the 36 km² reach studied. Meanders were shown to favour exchanges between river and aquifer, and although the tracer was diluted in the river, the contamination moved downstream from the injection plots and affected both river banks. Copyright © 2008 John Wiley & Sons, Ltd.     

A two‐dimensional morphodynamic model of gravel‐bed river with floodplain vegetation

Artificially straight river channels tend to be unstable, and ultimately develop into river meanders through bank erosion and point‐bar deposition. In this paper account is taken of the effects of riparian and floodplain vegetation on bank strength, floodplain flow resistance, shear stress partitioning, and bedload transport. This is incorporated into an existing 2D hydrodynamic‐morphological model. By applying the new model to an initially straight and single‐threaded channel, the way that its planform and cross‐sectional geometry evolve for different hydraulic and floodplain vegetation conditions is demonstrated. The results show the formation and upstream migration of gravel bars, confluence scouring and the development of meandering and braiding channel patterns. In cases where the channel becomes unstable, the instability grows out of bar formation. The resulting braiding patterns are similar to analytical results. The formation of a transition configuration requires a strong influence from vegetation. Copyright © 2010 John Wiley & Sons, Ltd.     

Investigating the influence of minor hydraulic structures on modeling flood events in lowland areas

Flood inundation models have been recognized to be a valuable tool to reproduce flow dynamics in a given area and support decision‐making processes on flood management measures. In many cases, in the simulation of flood events, only the main river channel and the associated structures are represented within the model. However, during flood events involving lowland areas, the minor drainage network – and the associated hydraulic structures – may have an important role in conveying flow and determining which areas will be flooded. The objective of this study is to investigate whether – and to what extent – small hydraulic structures in drainage networks have an influence in flooding on lowland areas. The case study for this research is the 1990 flood event which occurred in the lowland plain of the Reno River, in Northern Italy. The study area is mainly used for agricultural purposes and has a drainage system with several small bridges and culverts. The influence of the minor hydraulic structures on flood dynamics was analyzed through a combined use of one‐dimensional (1D) and two‐dimensional (2D) hydraulic models. First, a number of detailed and simplified approaches to represent hydraulic structures in the computational grids were analyzed by means of the HECRAS 1D model. Second, these approaches were implemented and tested in several 2D simulations of the flood event. The simulated inundation extents and flood levels were then compared with the observed data and with each other. The analysis of results showed that simplified schematizations were sufficient to obtain good model predictions of peak inundation extent and flood levels, at least for the present case study. Moreover, the influence of the structures on the peak flood inundation extent and flood levels was found to be limited, whereas it showed to be more significant during the drainage phase of the flood. Copyright © 2013 John Wiley & Sons, Ltd.     

Hydraulic and geomorphic processes in an overbank flood along a meandering,gravel‐bed river: implications for chute formation

Hydraulic interactions between rivers and floodplains produce off‐channel chutes, the presence of which influences the routing of water and sediment and thus the planform evolution of meandering rivers. Detailed studies of the hydrologic exchanges between channels and floodplains are usually conducted in laboratory facilities, and studies documenting chute development are generally limited to qualitative observations. In this study, we use a reconstructed, gravel‐bedded, meandering river as a field laboratory for studying these mechanisms at a realistic scale. Using an integrated field and modeling approach, we quantified the flow exchanges between the river channel and its floodplain during an overbank flood, and identified locations where flow had the capacity to erode floodplain chutes. Hydraulic measurements and modeling indicated high rates of flow exchange between the channel and floodplain, with flow rapidly decelerating as water was decanted from the channel onto the floodplain due to the frictional drag provided by substrate and vegetation. Peak shear stresses were greatest downstream of the maxima in bend curvature, along the concave bank, where terrestrial LiDAR scans indicate initial floodplain chute formation. A second chute has developed across the convex bank of a meander bend, in a location where sediment accretion, point bar development and plant colonization have created divergent flow paths between the main channel and floodplain. In both cases, the off‐channel chutes are evolving slowly during infrequent floods due to the coarse nature of the floodplain, though rapid chute formation would be more likely in finer‐grained floodplains. The controls on chute formation at these locations include the flood magnitude, river curvature, floodplain gradient, erodibility of the floodplain sediment, and the flow resistance provided by riparian vegetation. Copyright © 2015 John Wiley & Sons, Ltd.     

Influence of junction angle on three‐dimensional flow structure and bed morphology at confluent meander bends during different hydrological conditions

Recent field and modeling investigations have examined the fluvial dynamics of confluent meander bends where a straight tributary channel enters a meandering river at the apex of a bend with a 90° junction angle. Past work on confluences with asymmetrical and symmetrical planforms has shown that the angle of tributary entry has a strong influence on mutual deflection of confluent flows and the spatial extent of confluence hydrodynamic and morphodynamic features. This paper examines three‐dimensional flow structure and bed morphology for incoming flows with high and low momentum‐flux ratios at two large, natural confluent meander bends that have different tributary entry angles. At the high‐angle (90°) confluent meander bend, mutual deflection of converging flows abruptly turns fluid from the lateral tributary into the downstream channel and flow in the main river is deflected away from the outer bank of the bend by a bar that extends downstream of the junction corner along the inner bank of the tributary. Two counter‐rotating helical cells inherited from upstream flow curvature flank the mixing interface, which overlies a central pool. A large influx of sediment to the confluence from a meander cutoff immediately upstream has produced substantial morphologic change during large, tributary‐dominant discharge events, resulting in displacement of the pool inward and substantial erosion of the point bar in the main channel. In contrast, flow deflection is less pronounced at the low‐angle (36°) confluent meander bend, where the converging flows are nearly parallel to one another upon entering the confluence. A large helical cell imparted from upstream flow curvature in the main river occupies most of the downstream channel for prevailing low momentum‐flux ratio conditions and a weak counter‐rotating cell forms during infrequent tributary‐dominant flow events. Bed morphology remains relatively stable and does not exhibit extensive scour that often occurs at confluences with concordant beds. Copyright © 2014 John Wiley & Sons, Ltd.     

Evaluating the effect of scale in flood inundation modelling in urban environments

Cellular‐based approaches for flood inundation modelling have been extensively calibrated and evaluated for the prediction of flood flows on rural river reaches. However, there has only been limited application of these approaches to urban environments, where the need for flood management is greatest. Practical application of two‐dimensional (2D) flood inundation models is often limited by computation time and processing power on standard desktop PCs when attempting to resolve flows on the high‐resolution grids necessary to replicate urban features. Consequently, it is necessary to evaluate the effectiveness of coarse grids to represent flood flows through urban environments. To examine these effects, LISFLOOD‐FP, a 2D storage cell model, is applied to hypothetical flooding scenarios in Greenfields, Glasgow. Grid resampling techniques in GIS software packages are evaluated and a bilinear gridding technique appears to provide the most accurate and physically intuitive results. A gridding method maintaining sharp elevation changes at building interfaces and neighbouring land is presented and estimates of the discretization noise associated with the coarse resolution grids suggest little improvement over current gridding methods. The variation in model results from the friction sensitivity analysis suggests a non‐stationary response to Manning's n with changing model resolution. Model results suggests that a coarse resolution model for urban applications is limited by the representation of urban media in coarse model grids. Furthermore, critical length scales related to building dimensions and building separation distances exist in urban areas that determine maximum possible grid resolutions for hydraulic models of urban flooding. Copyright ©, 2008 John Wiley & Sons, Ltd.     

Morphological changes on meander point bars associated with flow structure at different discharges

This study investigates the impact of flow structure of different discharges on meander point bar morphology. We carried out mobile and terrestrial laser scanning campaigns before and after a flood on two sandy‐bed point bars. Between the scans, the flow structure was examined using an Acoustic Doppler Current Profiler at three flow stages. The results indicated that a meander point bar both affects and in turn, is itself modified by the flow at different discharges. The lower flow stages also have a significant effect on point bar morphology, especially on deposition over the bar head. Secondary circulation is responsible for scroll bar formation on the point bar margin beyond the apex. Flow separation at the inner bank, by contrast, does not require secondary circulation, but is dependent on flow depth over the point bar. A sudden increase in depth beyond the point bar top causes decreased stream power over the bar tail. The flow separation and decreased stream power cause a slow flow zone and net deposition over point bar tail. The backwater effect, if evident, may strengthen the process. Thus, filling over the bar tail seems generic for point bars and independent on secondary flow. Chutes and chute bars, scroll bars, bar head filling and bar platform filling, by contrast, require special fluvio‐morphological circumstances discussed in this paper. Whilst this paper confirms that the three‐dimensional flow structure has a major effect on point bar morphology, the flow structure seems to depend on how the point bar affects the flow trajectory which, in turn, depends upon the flow stage. Finally, the shape of the bend and the grain size distribution control the impacts of the flow structure, leading to divergent morphologies of point bars with certain generic features. Copyright © 2012 John Wiley & Sons, Ltd.     

Flow separation at the inner (convex) and outer (concave) banks of constant‐width and widening open‐channel bends

There is a paucity of data and insight in the mechanisms of, and controls on flow separation and recirculation at natural sharply‐curved river bends. Herein we report on successful laboratory experiments that elucidate flow structure in one constant‐width bend and a second bend with an outer‐bank widening. The experiments were performed with both a flat immobile gravel bed and mobile sand bed with dominant bedload sediment transport. In the constant‐width bend with immobile bed, a zone of mainly horizontal flow separation (vertical rotational axis) formed at the inner bank that did not contain detectable flow recirculation, and an outer‐bank cell of secondary flow with streamwise oriented rotational axis. Surprisingly, the bend with widening at the outer bank and immobile bed did not lead to a transverse expansion of the flow. Rather, flow in the outer‐bank widening weakly recirculated around a vertical axis and hardly interacted with the inner part of the bend, which behaved as a constant‐width bend. In the mobile bed experiment, downstream of the bend apex a pronounced depositional bar developed at the inside of the bend and pronounced scour occurred at the outside. Moreover the deformed bed promoted flow separation over the bar, including return currents. In the constant‐width bend, the topographic steering impeded the generation of an outer‐bank cell of secondary flow. In the bend with outer‐bank widening, the topographic steering induced an outward expansion of the flow, whereby the major part of the discharge was conveyed in the central part of the widening section. Flow in the outer‐bank widening was highly three dimensional and included return currents near the bottom. In conclusion, the experiments elucidated three distinct processes of flow separation common in sharp bends: flow separation at the inner bank, an outer‐bank cell of secondary flow, and flow separation and recirculation in an outer‐bank widening. Copyright © 2012 John Wiley & Sons, Ltd.     

Flood hydrograph attenuation induced by a reservoir system: analysis with a distributed rainfall‐runoff model

Spatially distributed hydrologic models can be effectively utilized for flood event simulation over basins where a complex system of reservoirs affecting the natural flow regime is present. Flood peak attenuation through mountain reservoirs can, in fact, mitigate the impact of major floods in flood‐prone areas of the lower river valley. Assessment of this effect for a complex reservoir system is performed with a spatially distributed hydrologic model where the surface runoff formation and the hydraulic routing through each reservoir and the river system are performed at a fine spatial and time resolution. The Toce River basin is presented as a case study, because of the presence of 14 active hydroelectric dams that affect the natural flow regime. A recent extreme flood event is simulated using a multi‐realization kriging method for modelling the spatial distribution of rainfall. A sensitivity analysis of the key elements of the distributed hydrologic model is also performed. The flood hydrograph attenuation is assessed. Several possible reservoir storage conditions are used to characterize the initial condition of each reservoir. The results demonstrate how a distributed hydrologic model can contribute to defining strategies for reservoir management in flood mitigation. Copyright © 2004 John Wiley & Sons, Ltd.     

A comparison of one‐ and two‐dimensional approaches to modelling flood inundation over complex upland floodplains

A much understudied aspect of flood inundation is examined, i.e. upland environments with topographically complex floodplains. Although the presence of high‐resolution topographic data (e.g. lidar) has improved the quality of river flood inundation predictions, the optimum dimensionality of hydraulic models for this purpose has yet to be fully evaluated for situations of both topographic and topological (i.e. the connectivity of floodplain features) complexity. In this paper, we present the comparison of three treatments of upland flood inundation using: (a) a one‐dimensional (1D) model (HEC‐RAS v. 3·1·2) with the domain defined as series of extended cross‐sections; (b) the same 1D model, but with the floodplain defined by a series of storage cells, hydraulically connected to the main river channel and other storage cells on the floodplain according to floodplain topological characteristics; (c) a two‐dimensional (2D) diffusion wave treatment, again with explicit representation of floodplain structural features. The necessary topographic and topological data were derived using lidar and Ordnance Survey Landline data. The three models were tested on a 6 km upland reach of the River Wharfe, UK. The models were assessed by comparison with measured inundation extent. The results showed that both the extended cross‐section and the storage cell 1D modes were conceptually problematic. They also resulted in poorer model predictions, requiring incorrect parameterization of the main river to floodplain flux in order to approach anything like the level of agreement observed when the 2D diffusion wave treatment was assessed. We conclude that a coupled 1D–2D treatment is likely to provide the best modelling approach, with currently available technology, for complex floodplain configurations. Copyright © 2007 John Wiley & Sons, Ltd.     

Evaluation of inundation areas resulting from the diversion of an extreme discharge towards the sea: case study in Tabasco,Mexico

This investigation comprises the hydraulic characterisation of a river located in the Mexican State of Tabasco, including the performance of its flood plain under the action of an extreme river discharge. This is done through the combination of a high‐quality validation dataset, remote sensing information, and a standard 2D numerical model. The dataset was collected during an intensive field campaign that took place in August 2009. In particular, in situ measurements of river discharge, bathymetry, water level, and velocities through a whole tidal cycle are employed along with multi‐spectral satellite imagery. The purpose of this study is twofold. Firstly, the integrated approach comprising the combination of a 2D hydrodynamic model, high‐quality in situ measurements and satellite imagery reduce the uncertainty in the model parameterisation and results. Secondly, it is shown that freely available sources of information, such as the Shuttle Radar Topographic Mission (SRTM) data can be processed and utilized in 2D hydraulic models. This is particularly important in countries where high‐resolution elevation data is not yet available. It is demonstrated that the selected approach is useful when the study of possible consequences in a flood plain induced by an extreme flood discharge are sought. Copyright © 2011 John Wiley & Sons, Ltd.     

Floodplain friction parameterization in two‐dimensional river flood models using vegetation heights derived from airborne scanning laser altimetry

Two‐dimensional (2‐D) hydraulic models are currently at the forefront of research into river flood inundation prediction. Airborne scanning laser altimetry is an important new data source that can provide such models with spatially distributed floodplain topography together with vegetation heights for parameterization of model friction. The paper investigates how vegetation height data can be used to realize the currently unexploited potential of 2‐D flood models to specify a friction factor at each node of the finite element model mesh. The only vegetation attribute required in the estimation of floodplain node friction factors is vegetation height. Different sets of flow resistance equations are used to model channel sediment, short vegetation, and tall and intermediate vegetation. The scheme was tested in a modelling study of a flood event that occurred on the River Severn, UK, in October 1998. A synthetic aperture radar image acquired during the flood provided an observed flood extent against which to validate the predicted extent. The modelled flood extent using variable friction was found to agree with the observed extent almost everywhere within the model domain. The variable‐friction model has the considerable advantage that it makes unnecessary the unphysical fitting of floodplain and channel friction factors required in the traditional approach to model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

HYDRAULIC ANALYSES FOR A NEW BRIDGE OVER THE PARANA RIVER, ARGENTINA

1 INTRODUCTION The Parana River in South America ranks as the sixth largest in the world in terms of average discharge (22 x 103 m3/s according to Schumm and Winkley, 1994), and its Argentinean reach stretches from Puerto Iguazu on the north to the Rio de la Plata estuary in the south (Fig. 1). Its estimated average suspended load, bed load, and total load are 23 x 106, 2.2 x 106, and 25.2 x 106 metric tons per year, respectively (Paoli and Schreider, 2001.) The economic developmen…     

One‐dimensional and two‐dimensional hydrodynamic modeling derived flow properties: impacts on aquatic habitat quality predictions

Studies of the effects of hydrodynamic model dimensionality on simulated flow properties and derived quantities such as aquatic habitat quality are limited. It is important to close this knowledge gap especially now that entire river networks can be mapped at the microhabitat scale due to the advent of point‐cloud techniques. This study compares flow properties, such as depth and velocity, and aquatic habitat quality predicted from pseudo‐2D and fully 2D hydrodynamic modeling. The models are supported by high‐resolution, point‐cloud derived bathymetries, from which close‐spaced cross‐sections were extracted for the 1D modeling, of three morphologically and hydraulically different river systems. These systems range from small low‐gradient meandering pool–riffle to large steep confined plane‐bed rivers. We test the effects of 1D and 2D models on predicted hydraulic variables at cross‐sections and over the full bathymetry to quantify the differences due to model dimensionality and those from interpolation. Results show that streambed features, whose size is smaller than cross‐sectional spacing, chiefly determine the different results of 1D and 2D modeling whereas flow discharge, stream size, morphological complexity and model grid sizes have secondary effects on flow properties and habitat quality for a given species and life stage predicted from 1D and 2D modeling. In general, the differences in hydraulic variables are larger in the bathymetric than in the cross‐sectional analysis, which suggests that some errors are introduced from interpolation of spatially disaggregated simulated variables with a 1D model, instead of model dimensionality 1D or 2D. Flow property differences are larger for velocity than for water surface elevation and depth. Differences in weighted usable area (WUA) derived from 1D and 2D modeling are relatively small for low‐gradient meandering pool–riffle systems, but the differences in the spatial distribution of microhabitats can be considerable although clusters of same habitat quality are spatially comparable. Copyright © 2014 John Wiley & Sons, Ltd.     

A field test of tubino's (1991) model of alternate bar formation

This study investigates the fluvial dynamics of straight natural stream channels. In particular, this experimental field study quantitatively assesses a physically based non-linear mathematical theory of alternate bar formation under unsteady natural flow conditions within a straight alluvial stream. The study site is an artificially straightened section of the Embarras River located approximately 16 km south of Champaign, Illinois. Data were collected on channel form, gradient, alternate bar dimensions, bed sediment size and flow stage over a 2 year study period. Both linear and non-linear steady flow hydrodynamic theories suggest that alternate bars are critical to the process of meander development. But these theories do not predict bar development for unsteady flow conditions, which typically occur in natural alluvial channels. Tubino (1991) suggests that bar evolution for a flood hydrograph can be divided into three parts: (1) a period of limited bar growth during the rising stage of the flood; (2) a stage of modest bar decay near the peak of the flood; and (3) a stage of non-linear bar growth during the prolonged falling stage of the flood. Bars developed during the falling limb of a hydrograph, and exhibited sequential development rather than the uniform growth along the reach predicted by Tubino's model. As flow stage decreased, short, low, fine-grained bars were superimposed on long, high and coarser-grained bars that developed under preceding high flow stages. These results suggest that the process of bar formation in artificially straightened natural streams with heterogeneous bed material may occur under different flow conditions and in a different manner than predicted by theoretical models. Further work should focus on attempting to isolate the physical mechanisms responsible for alternate bar formation in straight natural streams with heterogeneous bed material and flashy hydrologic flow regimes.     

Modelling the effect of time‐dependent river dune evolution on bed roughness and stage

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. Copyright © 2010 John Wiley & Sons, Ltd.