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.     

Influence of bed clusters and size gradation on operational time distribution for non‐uniform bed‐load transport

The operational time distribution (OTD) defines the time for bed‐load sediment spent in motion, which is needed to characterize the random nature of sediment transport. This study explores the influence of bed clusters and size gradation on OTD for non‐uniform bed‐loads. First, both static and mobile bed armouring experiments were conducted in laboratorial flumes to monitor the transport of mixed sand/gravel sediments. Only in the mobile armouring experiment did apparent bed clusters develop, because of stable feeding and a longer transport period. Second, a generalized subordinated advection (GSA) model was applied to quantify the observed dynamics of tracer particles. Results show that for the static armour layer (without sediment feed), the best‐fit OTD assigns more weight to the large displacement of small particles, likely because of the size‐selective entrainment process. The capacity coefficient in the GSA model, which affects the width of the OTD, is space dependent only for small particles whose dynamics can be significantly affected by larger particles and whose distribution is more likely to be space dependent in a mixed sand and gravel system. However, the OTD for the mobile armour layer (with sediment recirculation) exhibited longer tails for larger particles. This is because the trailing edge of larger particles is more resistant to erosion, and their leading front may not be easily trapped by self‐organized bed clusters. The strong interaction between particle–bed may cause the capacity coefficient to be space‐dependent for bed‐load transport along mobile armour layers. Therefore, the combined laboratory experiments and stochastic model analysis show that the OTD may be affected more by particle–bed interactions (such as clusters) than by particle–particle interactions (e.g. hiding and exposing), and that the GSA model can quantify mixed‐size sand/gravel transport along river beds within either static or mobile armour layers. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of bedload transport characteristics of Idaho streams and rivers

Field data are essential in evaluating the adequacy of predictive equations for sediment transport. Each dataset based on the sediment transport rates and other relevant information gives an increased understanding and improved quantification of different factors influencing the sediment transport regime in the specific environment. Data collected for 33 sites on 31 mountain streams and rivers in Central Idaho have enabled the analysis of sediment transport characteristics in streams and rivers with different geological, topographic, morphological, hydrological, hydraulic, and sedimentological characteristics. All of these streams and rivers have armored, poorly sorted bed material with the median particle size of surface layer coarser than the subsurface layer. The fact that the largest particles in the bedload samples did not exceed the median particle size of the bed surface material indicates that the armor layer is stable for the observed flow discharges (generally bankfull or less, and in some cases two times higher than bankfull discharge). The bedload transport is size‐selective. The transport rates are generally low, since sediment supply is less than the ability of flow to move the sediment for one range of flow discharges, or, the hydraulic ability of the stream is insufficient for entrainment of the coarse bed material. Detailed analyses of bedload transport rates, bedload and bed material characteristics were performed for each site. The obtained results and conclusions are used to identify different influences on bedload transport rates in analyzed gravel‐bed rivers. Copyright © 2008 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.     

Evaluating fine sediment mobilization and storage in a gravel‐bed river using controlled reservoir releases

Two controlled flow events were generated by releasing water from a reservoir into the Olewiger Bach, located near Trier, Germany. This controlled release of near bank‐full flows allowed an investigation of the fine sediment (<63 µm) mobilized from channel storage. Both a winter (November) and a summer (June) release event were generated, each having very different antecedent flow conditions. The characteristics of the release hydrographs and the associated sediment transport indicated a reverse hysteresis with more mass, but smaller grain sizes, moving on the falling limb. Fine sediment stored to a depth of 10 cm in the gravels decreased following the release events, indicating the dynamic nature and importance of channel‐stored sediments as source materials during high flow events. Sediment traps, filled with clean natural gravel, were buried in riffles before the release of the reservoir water and the total mass of fine sediment collected by the traps was measured following the events. Twice the mass of fine sediment was retained by the gravel traps compared with the natural gravels, which may be due to their altered porosity. Although the amount of fine sediment collected by the traps was not significantly related to measures of gravel structure, it was found to be significantly correlated to measures of local flow velocity and Froude number. A portion of the traps were fitted with lids to restrict surface exchange of water and sediment. These collected the highest amounts of event‐mobilized sediments, indicating that inter‐gravel lateral flows, not just surface infiltration of sediments, are important in replenishing and redistributing the channel‐stored fines. These findings regarding the magnitude and direction of fine sediment movement in gravel beds are significant in both a geomorphic and a biological context. Copyright © 2006 John Wiley & Sons, Ltd.     

Water‐worked gravel beds in laboratory flumes – a natural analogue?

An investigation has been conducted to identify the key parameters that are likely to scale laboratory sediment deposits to the field scale. Two types of bed formation were examined: one where sediment is manually placed and screeded and the second where sediment is fed into a running flume. This later technique created deposits through sequential cycles of sediment transport and deposition. Detailed bed surface topography measurements have been made over a screeded bed and three fed beds. In addition, bulk subsurface porosity and hydraulic conductivity have been measured. By comparing the four beds, results revealed that certain physical properties of the screeded bed were clearly different from those of the fed beds. The screeded bed had a random organization of grains on both the surface and within the subsurface. The fed beds exhibited greater surface and subsurface organization and complexity, and had a number of properties that closely resembled those found for water‐worked gravel beds. The surfaces were water‐worked and armoured and there was preferential particle orientation and direction of imbrication in the subsurface. This suggested that fed beds are able to simulate, in a simplified manner, both the surface and subsurface properties of established gravel‐bed river deposits. The near‐bed flow properties were also compared. It revealed that the use of a screeded bed will typically cause an underestimation in the degree of temporal variability in the flow. Furthermore, time‐averaged streamwise velocities were found to be randomly organized over the screeded bed but were organized into long streamwise flow structures over the fed beds. It clearly showed that caution should be taken when comparing velocity measurements over screeded beds with water‐worked beds, and that the formation of fed beds offers an improved way of investigating intragravel flow and sediment–water interface exchange processes in gravel‐bed rivers at a laboratory scale. Copyright © 2008 John Wiley & Sons, Ltd.     

Modeling the effects of pulsed versus chronic sand inputs on salmonid spawning habitat in a low‐gradient gravel‐bed river

It is widely recognized that high supplies of fine sediment, largely sand, can negatively impact the aquatic habitat quality of gravel‐bed rivers, but effects of the style of input (chronic vs. pulsed) have not been examined quantitatively. We hypothesize that a continuous (i.e. chronic) supply of sand will be more detrimental to the quality of aquatic habitat than an instantaneous sand pulse equal to the integrated volume of the chronic supply. We investigate this issue by applying a two‐dimensional numerical model to a 1 km long reach of prime salmonid spawning habitat in central Idaho. Results show that in both supply scenarios, sand moves through the study reach as bed load, and that both the movement and depth of sand on the streambed mirrors the hydrograph of this snowmelt‐dominated river. Predictions indicate greater and more persistent mortality of salmonid embryos under chronic supplies than pulse inputs, supporting our hypothesis. However, predicted mortality varies both with salmonid species and location of spawning. We found that the greatest impacts occur closer to the location of the sand input under both supply scenarios. Results also suggest that reach‐scale morphology may modulate the impact of sand loads, and that under conditions of high sand loading climate‐related increases in flow magnitude could increase embryo mortality through sand deposition, rather than streambed scour. Copyright © 2013 John Wiley & Sons, Ltd.     

Visual characterization technique for gravel‐cobble river bed surface sediments; validation and environmental applications Contribution to the programme of CIRSA (Centre Interuniversitaire de Recherche sur le Saumon Atlantique)

This paper presents an evaluation of the feasibility and the reliability of a visual characterization technique for gravel–cobble river bed surface substrate. Based on principal axis regressions, using phi scale (ϕ), comparisons of visual estimation and grid sampling techniques show that useful predictive relations (R² = 0·78–0·88) exist between visual estimates of the surface d₁₆, d₅₀ and d₈₄ and estimates obtained for the same percentiles with the grid sampling technique. Comparisons of visual estimation and the surface‐bulk sampling technique also indicate a predictive relation (R² = 0·70) between the d₅₀ of the two methods. Trained operators can visually estimate gravel–cobble bed surface d₁₆ to uncertainties of 41 per cent, d₅₀ to 15 per cent and d₈₄ to 11 per cent (for example, there is a 5·5 mm error on a d₈₄ size of 50 mm). Furthermore, evidence shows that if operators are properly trained, a calibration relation for each percentile can be applied independently of operators. This visual characterization allows effective detailed mapping of spatial patterns in substrate size distribution along extensive reaches of gravel‐bed rivers. The technique can be very useful in creating terrain models for various geomorphological, hydrological and biological applications such as the determination of entrainment thresholds, hydraulic roughness and substrate suitability for benthic insects or salmonid habitat. Copyright © 2001 John Wiley & Sons, Ltd.     

The particle size characteristics of fine‐grained channel deposits in the River Exe Basin,Devon, UK

Deposition and storage of fine‐grained (<62·5 μm) sediment in the hyporheic zone of gravel bed rivers frequently represents an important cause of aquatic habitat degradation. The particle size characteristics of such fine‐grained bed sediment (FGBS) exert an important control on its hydrodynamic properties and environmental impact. Traditionally, particle size analysis of FGBS in gravel bed rivers has focused on the absolute size distribution of the chemically dispersed mineral fraction. However, recent work has indicated that in common with fluvial suspended sediment, significant differences may exist between the absolute and the in situ, or effective, particle size composition of FGBS, as a result of the existence of aggregates, or composite particles. In the investigation reported in this paper, sealable bed traps that could be remotely opened to sample sediment deposited during specific storm runoff events and a laser back‐scatter probe were used to quantify the temporal and spatial variability of both the absolute and effective particle size composition of FGBS, and the associated suspended sediment from four gravel bed rivers in the Exe Basin, Devon, UK. The absolute particle size distributions of both the FGBS and suspended sediment evidenced c. >95%<62·5 μm sized primary particles and displayed a seasonal winter–summer fining, while the opposite trend was displayed by the effective particle size distribution of the FGBS and suspended sediment. The effective particle size distributions of both were typically highly aggregated, comprising up to 68%>62·5 μm sized particles. Spatial variation in the effective particle size and aggregation parameters was of secondary importance relative to temporal variation. The effective particle size distribution of the FGBS was consistently coarser and more aggregated than the associated suspended sediment and there was evidence of aggregate break‐up in samples of resuspended bed sediment. The implications of these findings for sediment transport modelling are considered. Copyright © 1999 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.     

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.     

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.     

A hybrid method for size characterization of coarse subsurface fluvial sediments

This paper introduces and evaluates a novel method for ascertaining the grain‐size distribution of subsurface sediments that involves profoundly less sampling effort than standard methods. It is based on hybrid sampling principles previously applied to the construction of synthetic surface grain‐size distributions. The method is developed from an empirical demonstration of the approximate similarity of surface and subsurface grain‐size distributions when compared over a common range of sizes. Subsurface hybrid models are found to provide good facsimiles of grain‐size distributions de?ned using standard criteria and to yield distribution percentiles with millimetre accuracy. The technique is presented as an expedient alternative to standard methods for large, perennial gravel‐bed rivers. As this is a new technique, prudent application is advised in lieu of further investigation. Copyright © 2004 John Wiley & Sons, Ltd.     

In situ characterization of grain‐scale fluvial morphology using Terrestrial Laser Scanning

The grain‐scale morphology of fluvial sediments is an important control on the character and dynamics of river systems; however current understanding of its role is limited by the difficulties of robustly quantifying field surface morphology. Terrestrial Laser Scanning (TLS) offers a new methodology for the rapid acquisition of high‐resolution and high‐precision surface elevation data from in situ sediments. To date, most environmental and fluvial applications of TLS have focused on large‐scale systems, capturing macroscale morphologies. Application of this new technology at scales necessary to characterize the complexity of grain‐scale fluvial sediments therefore requires a robust assessment of the quality and sources of errors in close‐range TLS data. This paper describes both laboratory and field experiments designed to evaluate close‐range TLS for sedimentological applications and to develop protocols for data acquisition. In the former, controlled experiments comprising high‐resolution scans of white, grey and black planes and a sphere were used to quantify the magnitude and source of three‐dimensional (3D) point errors resulting from a combination of surface geometry, reflectivity effects and inherent instrument precision. Subsequently, a methodology for the collection and processing of grain‐scale TLS data is described through an application to a coarse grained gravel system, the River Feshie (D₅₀ 32 to 63 mm). This stepwise strategy incorporates averaging repeat scans and filtering scan artefact and non‐surface points using local 3D search algorithms. The sensitivity of the results to the filter parameter values are assessed by careful internal validation of Digital Terrain Models (DTMs) created from the resulting point cloud data. The transferability of this methodology is assessed through application to a second river, Bury Green Brook, dominated by finer gravel (D₅₀ 18 to 33 mm). The factor limiting the resolution of DTMs created from this second dataset was found to be the relative sizes of the laser footprint and smallest grains. Copyright © 2009 John Wiley & Sons, Ltd.     

The response of a gravel‐bed river planform configuration to flow variations and bed reworking: a modelling study

A 2D depth‐averaged model has been developed for simulating water flow, sediment transport and morphological changes in gravel‐bed rivers. The model was validated with a series of laboratory experiments and then applied to the Nove reach of the Brenta River (Northern Italy) to assess its bed material transport, interpret channel response to a series of intensive flood events (R.I. ≈ 10 years) and provide a possible evolutionary scenario for the medium term. The study reach is 1400 m long with a mean slope of 0.0039 m m^?1. High‐resolution digital terrain models were produced combining LiDAR data with colour bathymetry techniques. Extensive field sedimentological surveys were also conducted for surface and subsurface material. Data were uploaded in the model and the passage of two consecutive high intensity floods was simulated. The model was run under several hypotheses of sediment supply: one considering substantial equilibrium between sediment input and transport capacity, and the others reducing the sediment supply. The sediment supply was then calibrated comparing channel morphological changes as observed in the field and calculated by the model. Annual bed material transport was assessed and compared with other techniques. Low‐frequency floods (R.I. ≈ 1.5 years) are expected to produce negligible changes in the channel while high floods may erode banks rather than further incising the channel bed. Location and distribution of erosion and deposition areas within the Nove reach were predicted with acceptable biases stemming from imperfections of the model and the specified initial, boundary and forcing conditions. A medium‐term evolutionary scenario simulation underlined the different response to and impact of a consecutive sequence of floods. Copyright © 2015 John Wiley & Sons, Ltd.     

Detection and reconstruction of large‐scale coherent flow structures in gravel‐bed rivers

Large‐scale flow structures (LSFS) in the streamwise direction are important features of gravel‐bed river flows, because they may contribute to sediment transport and gas exchange. In the present study, these structures are detected using Huang's empirical mode decomposition and reconstructed with phase‐averaging techniques based on a Hilbert transform of the velocity signal. The analysis is based on the fluctuating component of 15 quasi‐instantaneous velocity profiles measured with a three‐dimensional (3D) acoustic Doppler velocity profiler (ADVP) in an armoured gravel‐bed river with a low relative submergence of 2.9 (ratio between flow depth and bed grain diameter). LSFS were identified in most of the measured profiles and consistently showed similar features. We were able to characterize the geometry of these large‐scale coherent structures: the front has a vertical linear shift in the time domain and a vertical profile corresponding to a first quarter moon with the apex situated at z/h ≈ 0.4. In the vertical, the front scales with flow depth h, and in the streamwise direction, LSFS scale with three to seven times the mean flow depth. On the bed, the effect of LSFS is a periodic non‐linear variation of the friction velocity on average between 0.90 and 1.10 times the mean value. A model for the friction velocity cycle resulting from LSFS oscillation is presented. Copyright © 2014 John Wiley & Sons, Ltd.     

Biases of submerged bulk and freeze‐core samples

Freeze‐coring and bulk sampling are routine methods used to sample subsurface spawning gravel under shallow water. Both methods have limitations. Freeze‐coring is not believed to representatively sample coarse grain sizes and the sample volumes are relatively small. Conversely, when bulk sampling, even within an enclosure, some fine sediment is suspended and washed away from the sample. This paper assesses the biases in sampling performance between the two methods and determines whether the loss of fines that occurs when bulk sampling could be predicted and thus corrected for. At six riffles the spawning substrate was sampled under approximately 50 cm of water with a bulk sample and three adjacent freeze‐cores. For each riffle, data from the two samples were combined using the method of Fripp and Diplas (1993) and the resultant composite sample was compared with the original freeze‐core and bulk samples to assess the relative precision and biases of the two techniques. On average, the D₅₀ of the bulk samples was 4 mm larger and a one‐third loss of the <2 mm fraction occurred compared with the composite samples. In contrast, freeze‐core samples contain on average 32% more sediment >16 mm compared with composite samples. Based on six samples, taken from six riffles, the amount of sediment finer than 0·5 mm lost using our bulk sampling technique with an enclosure appears to be predictable and correctable. Copyright © 2005 John Wiley & Sons, Ltd.     

Bridging the gap between turbulence and larger scales of flow motions in rivers

Although flow turbulence in rivers is of critical importance to earth scientists, ecologists and engineers, its relations with larger flow scales are not well understood, thus leaving a fundamental gap in our knowledge. From an analysis of a long time series of the streamwise and vertical flow velocity fluctuations measured in a gravel‐bed river, we show that the signature of the fundamental turbulent flow structures (e.g. ejections and sweeps) is embedded within increasingly larger flow scales in a self‐similar manner. The imbrication of turbulent structures into large flow pulsations of flow acceleration and deceleration covers more than two‐orders of magnitude from a few seconds to nearly 10 minutes. This property is explained by the clustering of turbulent events creating an emergent pattern at larger scales. The size of the larger flow pulsations scales with the spacing of the pools and riffles in the river. This implies a mutual adjustment between turbulence generation mechanisms and long pulsations of flow acceleration and deceleration controlled by the bed morphology. These results bridge a gap in our understanding of flows in rivers and offer a new perspective on the interactions between the turbulent flow with larger scales of flow motion that are critical for sediment transport, habitat selection and fish behaviour. Copyright © 2011 John Wiley & Sons, Ltd.