Morphodynamics of alternate bars in the presence of riparian vegetation

Alpine gravel-bed rivers are dynamic systems that have been subjected to many anthropic alterations in the past centuries. Riparian vegetation development on previously bare sediment bedforms has been a common adjustment, raising important management issues in terms of flood risks and biodiversity. Many of these rivers are also channelized, and as a result present a pattern of alternate bars. Considering recent advances in numerical biomorphodynamic modeling, this study aims at exploring numerically the morphodynamics of alternate bars in the presence of riparian vegetation. To this end, a dynamic vegetation module has been implemented on top of an existing morphodynamic model, accounting for ecological processes of seed dispersal, seedling recruitment, growth, and mortality. Numerical simulations have been performed on a simplified reach of a gravel-bed river with free migrating alternate bars at initial state. In this work 96 scenarios have been simulated, each representing 50 years of channel evolution, with different flood regimes characterized by various peak discharges and flood durations. Yearly peak discharge variability is explicitly modeled in 48 scenarios. Model outcomes present two possible equilibrium biomorphodynamic behaviors: stationary vegetated bars, or free migrating bars in the case of frequent vegetation removal during floods. This binary behavior holds true when the stochasticity of annual peak discharges is represented, and for a wide range of parameter values included in vegetation dynamic modeling. Transient mobility of vegetated bars is observed in specific configurations where large sediment deposits deflect the flow field, eroding bar heads. Modeled bar wavelengths are in the range of values predicted for free bars by linear bar theory, and remain far from the theoretical values of hybrid, steady bars. The shift from unvegetated migrating bars to steady vegetated bars seems to show that in these simulations vegetation constitutes a hydraulic forcing, leading to a shift from free bars to forced bars, with a final configuration largely inherited from the initial state. © 2019 John Wiley & Sons, Ltd.     

Riffle‐pool morphometry and stage‐dependant morphodynamics of a large floodplain river (Vereinigte Mulde,Sachsen‐Anhalt,Germany)

Riffle‐pool sequences are a common feature of gravel‐bed rivers. However, mechanisms of their generation and maintenance are still not fully understood. In this study a monitoring approach is employed that focuses on analysing cross‐sectional and longitudinal channel geometry of a large floodplain river (Vereinigte Mulde, Sachsen‐Anhalt, Germany) with a high temporal and spatial resolution, in order to conclude from stage‐dependant morphometric changes to riffle and pool maintaining processes. In accordance with previous authors, pool cross‐sections of the Mulde River are narrow and riffle cross‐sections are wide suggesting that they should rather be addressed as two general types of channel cross‐sections than solely as bedforms. At high flows, riffles and pools in the study reaches changed in length and height but not in position. Pools were scoured and riffles aggraded, a development which was reversed during receding flows below the threshold of 0·4Q_bf (40% bankfull discharge). An index for the longitudinal amplitude of riffle‐pool sequences, the bed undulation intensity or bedform amplitude, is introduced and proved to be highly significant as a form parameter, its first derivative as a process parameter. The process of pool scour and riffle fill is addressed as bedform maintenance or bedform accentuation. It is indicated by increasing longitudinal bed amplitudes. According to the observed dynamics of bed amplitudes, maintenance of riffle‐pool sequences lags behind discharge peaks. Maximum bed amplitudes may be reached with a delay of several days after peak discharges. Increasing bed undulation intensity is interpreted to indicate bed mobility. Post‐flood decrease of the bed undulation intensity indicates a retrograde phase when transport from pools to riffles has ceased and bed mobility is restricted to riffle tails and heads of pools. This type of transport behaviour is referred to as disconnected mobility. The comparison of two river reaches, one with undisturbed sediment supply, the other with sediment deficit, suggests that high bed undulation intensity values at low flows indicate sediment deficit and potentially channel degrading conditions. It is more generally hypothesized that channel bed undulations constitute a major component of form roughness and that increased bed amplitudes are an important feature of channel bed adjustment to sediment deficit be it temporally during late floods or permanently due to a supply limitation of bedload. Copyright © 2011 John Wiley & Sons, Ltd.     

Effects of river discharge on hyporheic exchange flows in salmon spawning areas of a large gravel‐bed river

The flow magnitude and timing from hydroelectric dams in the Snake River Basin of the Pacific north‐western US is managed in part for the benefit of salmon. The objective of this research was to evaluate the effects of Hells Canyon Dam discharge operations on hydrologic exchange flows between the river and riverbed in Snake River fall Chinook salmon spawning areas. Interactions between river water and pore water within the upper 1 m of the riverbed were quantified through the use of self‐contained temperature and water level data loggers suspended inside of piezometers. The data were recorded at 20 min intervals over a period of 200 days when the mean daily discharge was 218–605 m³ s^?1, with hourly stage changes as large as 1·9 m. Differences in head pressure between the river and riverbed were small, often within ± 2 cm. Measured temperature gradients in the riverbed indicated significant interactions between the surface and subsurface water. At the majority of sites, neither hydraulic nor temperature gradients were significantly affected by either short‐ or long‐term changes in discharge operations from Hells Canyon Dam. Only 2 of 14 study sites exhibited acute flux reversals between the river and riverbed resulting from short‐term, large magnitude changes in discharge. The findings suggest that local scale measurements may not be wholly explanatory of the hydrological exchange between the river and riverbed. The processes controlling surface water exchange at the study sites are likely to be bedform‐induced advective pumping, turbulence at the riverbed surface, and large‐scale hydraulic gradients along the longitudinal profile of the riverbed. By incorporating the knowledge of hydrological exchange processes into water management planning, regional agencies will be better prepared to manage the limited water resources among competing priorities that include salmon recovery, flood control, irrigation supply, hydropower production, and recreation. Copyright © 2007 John Wiley & Sons, Ltd.     

An empirical model of event scale cohesive bank profile evolution

The highly stochastic nature of riverbank erosion has driven the need for spatially explicit empirical models. Detailed bank profile surveys along a meander bend of the Brandywine Creek in Pennsylvania, USA, before and after 28 high flow events over a 2·5 year period are used to develop an empirical model of cohesive bank profile erosion. Two hundred and thirty‐six bank erosion observations are classified as hydraulic erosion or subaerial erosion. Threshold conditions required to initiate bank erosion cannot be defined based on field measurements. Using the near‐bank velocity and the number of freeze–thaw cycles as predictors, regression equations are derived for hydraulic erosion that specify the length, thickness, and location on the bank face of eroded blocks. An empirical discriminant function defines the critical geometry of overhang failures, and the volumes removed by overhang failures are computed using another regression equation. All the regression equations are significant, but have low correlation coefficients, suggesting that cohesive bank erosion has a strong stochastic component. Individual events typically remove small masses of soil (average volume 0·084 m³/m) a few centimeters thick (median = 0·057 m) and a few decimeters in length (median = 0·50 m) from the lower third of the bank. Hydraulic erosion is responsible for 87% of all erosion. When applied to three survey sites not used in its development, the profile model predicts the total volume of erosion with errors of 23%, 5% and 1%. Twenty‐four percent of computed erosion volumes for single events are within 50% of observed volumes at these three sites. Extending the approach to decadal timescales and to entire bends will require three‐dimensional observations of bank failure, and spatially and temporally explicit methods to account for the influence of individual large trees on bank failures and near‐bank hydraulic processes. Copyright © 2009 John Wiley & Sons, Ltd.     

Determination of bankfull discharge magnitude and frequency: comparison of methods on 16 gravel‐bed river reaches

Bankfull discharge is identified as an important parameter for studying river morphology, sediment motion, flood dynamics and their ecological impacts. In practice, the determination of this discharge and its hydrological characteristics is not easy, and a choice has to be made between several existing methods. To evaluate the impact of the choice of methods, five bankfull elevation definitions and four hydrological characterizations (determination of duration and frequency of exceedance applied to instantaneous or mean daily data) were compared on 16 gravel‐bed river reaches located in France (the catchment sizes vary from 10 km² to 1700 km²). The consistency of bankfull discharge estimated at reach scale and the hydraulic significance of the five elevation definitions were examined. The morphological definitions (Bank Inflection, Top of Bank) were found more relevant than the definitions based on a geometric criterion. The duration of exceedance was preferred to recurrence intervals (partial duration series approach) because it is not limited by the independency of flood events, especially for low discharges like those associated with the Bank Inflection definition. On average, the impacts of the choice of methods were very important for the bankfull discharge magnitude (factor of 1·6 between Bank Inflection and Top of Bank) and duration of exceedance or frequency (respectively a factor 1·8 and 1·9 between mean daily and instantaneous discharge data). The choice of one combination of methods rather than another can significantly modify the conclusions of a comparative analysis in terms of bankfull discharge magnitude and its hydrological characteristics, so that one must be cautious when comparing results from different studies that use different methods. Copyright © 2006 John Wiley & Sons, Ltd.     

Particle clusters in gravel‐bed rivers: an experimental morphological approach to bed material transport and stability concepts

Structured gravel river beds clearly exert a major influence on bed stability. Indexing structural stability by field measurements of bed strength neglects the processes operating to entrain and transport bed material in different parts of each structure. This study takes a morphological approach to interpreting the critical processes, using particle tracing to determine the movement of individual cluster particles over a range of flood event magnitudes and durations. The experiment was carried out on the River South Tyne, UK; it uses flow hydrographs measured nearby and also benefits from previous studies of historical development, channel morphology and sediment transport at the same site. More than 30 clusters were monitored over a seven‐month period during which clusters occupied 7–16 per cent of the bed. Threshold flows delimiting three apparently contrasting bed sediment process regimes for cluster particles are tentatively set at 100 m³ s^?1 and 183 m³ s^?1; durations of flow at these levels are critical for cluster development, rather than flow peak values. Wake particles are transported most easily. Flow straightening in the wandering channel planform reduces the stability of clusters, since mechanical strength is markedly reduced by this change of direction. The overall area covered by clusters between significant transport events varies little, implying a dynamic equilibrium condition. Copyright © 2005 John Wiley & Sons, Ltd.     

Computational fluid dynamics modelling of boundary roughness in gravel‐bed rivers: an investigation of the effects of random variability in bed elevation

Results from a series of numerical simulations of two‐dimensional open‐channel flow, conducted using the computational fluid dynamics (CFD) code FLUENT, are compared with data quantifying the mean and turbulent characteristics of open‐channel flow over two contrasting gravel beds. Boundary roughness effects are represented using both the conventional wall function approach and a random elevation model that simulates the effects of supra‐grid‐scale roughness elements (e.g. particle clusters and small bedforms). Results obtained using the random elevation model are characterized by a peak in turbulent kinetic energy located well above the bed (typically at y/h = 0·1–0·3). This is consistent with the field data and in contrast to the results obtained using the wall function approach for which maximum turbulent kinetic energy levels occur at the bed. Use of the random elevation model to represent supra‐grid‐scale roughness also allows a reduction in the height of the near‐bed mesh cell and therefore offers some potential to overcome problems experienced by the wall function approach in flows characterized by high relative roughness. Despite these benefits, the results of simulations conducted using the random elevation model are sensitive to the horizontal and vertical mesh resolution. Increasing the horizontal mesh resolution results in an increase in the near‐bed velocity gradient and turbulent kinetic energy, effectively roughening the bed. Varying the vertical resolution of the mesh has little effect on simulated mean velocity profiles, but results in substantial changes to the shape of the turbulent kinetic energy profile. These findings have significant implications for the application of CFD within natural gravel‐bed channels, particularly with regard to issues of topographic data collection, roughness parameterization and the derivation of mesh‐independent solutions. Copyright © 2001 John Wiley & Sons, Ltd.     

Grain size and topographical differences between static and mobile armour layers

A series of laboratory flume experiments under conditions of sediment starvation (zero sediment feeding) and recirculation were conducted in order to identify the temporal evolution and surface properties of static and mobile armour layers. The experiments were carried out in an 8 m long flume using a bimodal grain‐size mixture (D₅₀ = 6·2 mm) and a range of shear stresses ranging from 4·0 to 8·6 N m^–2. The results confirm that a static armour layer is coarser than a mobile one, and that the grain size of a mobile armour layer is rather insensitive to changes in the imposed flow strength. An analysis of laser scan bed surveys revealed the highly structured and imbricated nature of the static armour layer. Under these conditions the vertical roughness length scale of the bed diminished and it became topographically less complex at higher forming discharges. The topography of mobile armour layers created by rising discharges differed. They exhibited a greater roughness length scale and were less organized, despite the fact that the grain size of the surface material maintained an approximately constant value during recirculation. Also, the mobile armour tended to create larger cluster structures than static armour layers when formed by higher discharges. These differences were mainly due to the transport of the coarser fraction of bed sediments, which diminished to zero over the static armour because of being hidden within the bed, whereas in the mobile armour the coarser particles protruded into the flow and were actively transported, increasing the vertical roughness length scale. Overall, the results show that an examination of the grain size characteristics of armour layers cannot be used to infer sediment mobility and bed roughness. Detailed elevation models of exposed surfaces of gravel‐bed rivers are required to provide critical insight on the sediment availability and sedimentation processes. Copyright © 2011 John Wiley & Sons, Ltd.     

Texture‐based image segmentation applied to the quantification of superficial sand in salmonid river gravels

The presence of ?ne sediment in river gravels is widely recognized as being detrimental to salmonid habitat quality. In order to facilitate quanti?cation of sand presence at larger scales, this paper presents an application of image processing allowing for rapid and accurate assessments of super?cial sand presence in dry exposed ?uvial gravels. Images for the process are acquired with a 35 mm SLR ?lm camera and then scanned with a desktop scanner. Texture‐based segmentation is then applied to differentiate between sand and clast areas. Results show that the method is accurate and therefore it offers an alternative to bulk sampling in cases where rapid assessments of sand presence are required. Copyright © 2005 John Wiley & Sons, Ltd.     

Morphodynamic evolution of a lower Mississippi River channel bar after sand mining

In‐channel sand mining by dredge removes large quantities of bed sediment and alters channel morphodynamic processes. While the reach‐scale impacts of dredging are well documented, the effects of the dredged borrow pit on the local flow and sediment transport are poorly understood. These local effects are important because they control the post‐dredge evolution of the borrow pit, setting the pit lifespan and affecting reach‐scale channel morphology. This study documents the observed morphological evolution of a large (1·46 million m³) borrow pit mined on a lateral sandbar in the lower Mississippi River using a time‐series of multibeam bathymetric surveys. During the 2·5 year time‐series, 53% of the initial pit volume infilled with sediment, decreasing pit depth by an average of 0·88 m yr^?1. To explore the controls of the observed infilling, a morphodynamic model (Delft3D) was used to simulate flow and sediment transport within the affected river reach. The model indicated that infilling rates were primarily related to the riverine sediment supply and pit geometry. The pit depth and length influenced the predicted magnitude of the pit bed shear stress relative to its pre‐dredged value, i.e. the bed‐stress reduction ratio (R*), a metric that was correlated with the magnitude and spatial distribution of infilling. A one‐dimensional reduced‐complexity model was derived using predicted sediment supply and R* to simulate patterns of pit infilling. This simplified model of borrow‐pit evolution was able to closely approximate the amount and patterns of sediment deposition during the study period. Additional model experiments indicate that, for a borrow pit of a set volume, creating deep, longitudinally‐shorter borrow pits significantly increased infilling rates relative to elongated pits. Study results provide insight into the resilience of alluvial river channels after a disturbance and the sustainability of sand mining as a sediment source for coastal restoration. Copyright © 2015 John Wiley & Sons, Ltd.     

Interactions between sediment delivery,channel change,climate change and flood risk in a temperate upland environment

This paper uses numerical simulation of flood inundation based on a coupled one‐dimensional–two‐dimensional treatment to explore the impacts upon flood extent of both long‐term climate changes, predicted to the 2050s and 2080s, and short‐term river channel changes in response to sediment delivery, for a temperate upland gravel‐bed river. Results show that 16 months of measured in‐channel sedimentation in an upland gravel‐bed river cause about half of the increase in inundation extent that was simulated to arise from climate change. Consideration of the joint impacts of climate change and sedimentation emphasized the non‐linear nature of system response, and the possibly severe and synergistic effects that come from combined direct effects of climate change and sediment delivery. Such effects are likely to be exacerbated further as a result of the impacts of climate change upon coarse sediment delivery. In generic terms, these processes are commonly overlooked in flood risk mapping exercises and are likely to be important in any river system where there are high rates of sediment delivery and long‐term transfer of sediment to floodplain storage (i.e. alluviation involving active channel aggradation and migration). Similarly, attempts to reduce channel migration through river bank stabilization are likely to exacerbate this process as without bank erosion, channel capacity cannot be maintained. Finally, many flood risk mapping studies rely upon calibration based upon combining contemporary bed surveys with historical flood outlines, and this will lead to underestimation of the magnitude and frequency of floodplain inundation in an aggrading system for a flood of a given magnitude. Copyright © 2006 John Wiley & Sons, Ltd.     

Video‐based gravel transport measurements with a flume mounted light table

The study of bedload transport processes is constrained by an inability to monitor the mass, volume and grain size distribution of sediment in transport at high temporal frequencies. Building upon a previously published design, we have integrated a high‐resolution (1392 × 1024 pixels) video camera with a light table to continuously capture images of 2–181 mm material exiting a flume. The images are continuously recorded at a rate of 15 to 20 frames per second and are post‐processed using LabView^(?) software, yielding continuous grain‐size‐specific transport information on a per second basis. The video capture rate is sufficient to record multiple images of each grain leaving the flume so that particle velocities can be measured automatically. No manual image processing is required. After calibration the method is accurate and precise for sediment in the 2 mm through to 45 mm grain size classes compared with other means of measuring bedload. Based on a set of validation samples, no statistically significant difference existed between the D₁₀, D₁₆, D₂₅, D₅₀, D₇₅, D₈₄, D₉₀ and D₉₅ determined by sieving captured samples and the D_i values determined with the system. On average the system overpredicted transport by 4 per cent (n = 206, SD = 42%). This error can be corrected easily by simply weighing the mass of sediment that leaves the flume. The technology is relatively inexpensive and provides high‐resolution data on coarse sediment transport out of a flume. Copyright © 2008 John Wiley & Sons, Ltd.     

Large grains matter: contrasting bed stability and morphodynamics during two nearly identical experiments

While the stabilizing function of large grains in step‐pool streams has long been recognized, the role they play in gravel‐bed streams is less clear. Most researchers have ignored the role of large grains in gravel‐bed streams, and have assumed that the median bed surface size controls the erodibility of alluvial boundaries. The experiments presented herein challenge this convention. Two experiments were conducted that demonstrate the significant morphodynamic implications of a slight change to the coarse tail of the bed material. The two distributions had the same range of particle sizes, and nearly identical bulk d₅₀ values (1.6 mm); however the d₉₀ of experiment GSD1 was slightly finer (3.7 mm) than that for experiment GSD2 (3.9 mm). Transport rates during GSD1 were nearly four times greater than during GSD2 (even though the dimensionless shear stress was slightly lower), and the channel developed a sinuous pattern with well‐developed riffles, pools and bars. During GSD2 the initial rectangular channel remained virtually unchanged for the duration of the experiment. The relative stability of GSD2 seems to be associated with a slightly larger proportion of stable (large) grains on the bed surface: at the beginning of GSD1, 3.5% of the bed was immobile, while almost twice as much of it (6.1%) was immobile at the beginning of GSD2. The results demonstrate that the largest grains (not the median size) exert first‐order control on channel stability. Copyright © 2017 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.     

The bed morphology of upland single‐thread channels in semi‐arid environments: evidence of repeating bedforms and their wider implications for gravel‐bed rivers

Single‐thread, gravel‐bed streams of moderate slope in the northern Negev are characterized by three channel units: bars exhibit steeper than average slopes and poorly sorted mixtures of small–medium cobbles and coarse–very coarse pebbles; flats are associated with more gentle slopes and well‐sorted medium–fine pebbles and granules; and transitional units have intermediate slopes and grain size. In general, all three units are planar, span the full channel width and have well‐defined boundaries. Bars and flats are more common than the transitional units and alternate downstream for distances of several hundred metres, forming sequences that are reminiscent of the riffle–pool structure commonly observed in humid‐temperate gravel‐bed rivers. A notable contrast is the absence of significant bed relief: bars lack crests and flats lack depressions. The relative lack of bed relief in bar–flat sequences is attributed to the high rate of sediment supply from the sparsely vegetated hillslopes which promotes the infilling of depressions and to the erosion of crests under conditions of intense transport. This reduction of bed relief lowers channel roughness, which in turn increases flow velocity and, therefore, the ability of the channel to transmit the large sediment loads it receives. Although our analyses pertain to a semi‐arid river system, the results have wider implications for understanding the adjustment of channel bedform to high sediment loads in other fluvial environments. Copyright © 2012 John Wiley & Sons, Ltd.     

Life of a bifurcation in a gravel‐bed braided river

Channel bifurcation is a key element in braided rivers, determining the water and sediment distribution and hence controlling the morphological evolution. Recent theoretical and experimental findings, as well as field observations, showed that bifurcations in gravel‐bed braided rivers are often asymmetrical and highly unstable. In this paper field data are presented on a bifurcation in the Tagliamento River, northeast Italy. The planform configuration of the bifurcation and its temporal evolution was monitored by an automatic digital camera during a series of seven floods with different magnitudes. This remote sensing technique allowed a high temporal resolution (pictures were acquired every hour) that was proved to be essential in a highly dynamic system as the one considered here. Digitized maps of the channels provided information on the location of the bifurcation, the width of the anabranches, the angle between them, along with the occurrence and migration of sediment bars. Data were acquired at two different water levels, giving the possibility to compare low and high flow conditions. The monitored bifurcation is largely unstable and shows sudden changes in the water distribution, mainly driven by the bar migrating in the upstream channel and entering the distributaries. A relationship between width asymmetry and flood magnitude was observed, confirming previous analyses. Moreover, recent theoretical findings were applied, in order to test the possibility to estimate general trends in bifurcation evolution. The analysis pointed out the relevance of a correct assessment of the characteristic temporal scales, as the bifurcation evolves on a timescale similar to that of bar migration and flood duration. Understanding the interactions between these processes is therefore crucial in order to increase the ability to model and predict the morphological evolution of a braided network. Copyright © 2012 John Wiley & Sons, Ltd.