Relation between flow,surface‐layer armoring and sediment transport in gravel‐bed rivers

This study investigates trends in bed surface and substrate grain sizes in relation to reach‐scale hydraulics using data from more than 100 gravel‐bed stream reaches in Colorado and Utah. Collocated measurements of surface and substrate sediment, bankfull channel geometry and channel slope are used to examine relations between reach‐average shear stress and bed sediment grain size. Slopes at the study sites range from 0·0003 to 0·07; bankfull depths range from 0·2 to 5 m and bankfull widths range from 2 to 200 m. The data show that there is much less variation in the median grain size of the substrate, D_50s, than there is in the median grain size of the surface, D₅₀; the ratio of D₅₀ to D_50s thus decreases from about four in headwater reaches with high shear stress to less than two in downstream reaches with low shear stress. Similar trends are observed in an independent data set obtained from measurements in gravel‐bed streams in Idaho. A conceptual quantitative model is developed on the basis of these observations to track differences in bed load transport through an idealized stream system. The results of the transport model suggest that downstream trends in total bed load flux may vary appreciably, depending on the assumed relation between surface and substrate grain sizes. Copyright © 2007 John Wiley & Sons, Ltd.     

Coarse bed load transport in an alluvial gravel bed stream,Dupuyer Creek,Montana

Coarse bed load was sampled in a gravel/cobble bed stream during two major floods in the snowmelt runoff season. The channel is characterized by high rates of bank erosion and, therefore, high rates of sediment supply and bed load flux. Peak discharge reached four times bank‐full, and bed load was sampled at flows 0·7–1·7 times bank‐full. A large aperture bed load sampler (1 m by 0·45 m) captured the largest particles in motion, and specifically targeted the coarse bed load size distribution by using a relatively large mesh (32 mm or D₂₅ of streambed surface size distribution). Bed load flux was highly variable, with a peak value of 0·85 kg/s/m for the coarse fraction above 38 mm. Bed load size distribution and maximum particle size was related to flow strength. Entrainment was size selective for particles D₇₀ and larger (88–155 mm), while particles in the range D₃₀–D₇₀ (35–88 mm) ceased to move at essentially the same flow. Bed load flux was size selective in that coarse fractions of the streambed surface were under‐represented in or absent from the bed load. Painted tracer particles revealed that the streambed surface in the riffles could remain stable even during high rates of bed load transport. These observations suggest that a large proportion of bed load sediments was sourced from outside the riffles. Repeat surveys confirmed major scour and fill in pools (up to 0·75 m), and bank erosion (>2 m), which together contributed large volumes of sediment to the bed load. Copyright © 2007 John Wiley & Sons, Ltd.     

In situ quantification of spatial and temporal variability of hyporheic exchange in static and mobile gravel‐bed rivers

Seepage meters modified for use in flowing water were used to directly measure rates of exchange between surface and subsurface water in a gravel‐ and cobble bed river in western Pennsylvania, USA (Allegheny River, Q_mean = 190 m³/s) and a sand‐ and gravel‐bed river in Colorado, USA (South Platte River, Q_mean = 9·7 m³/s). Study reaches at the Allegheny River were located downstream from a dam. The bed was stable with moss, algae, and river grass present in many locations. Median seepage was + 0·28 m/d and seepage was highly variable among measurement locations. Upward and downward seepage greatly exceeded the median seepage rate, ranging from + 2·26 (upward) to ? 3·76 (downward) m/d. At the South Platte River site, substantial local‐scale bed topography as well as mobile bedforms resulted in spatial and temporal variability in seepage greatly in exceedence of the median groundwater discharge rate of 0·24 m/d. Both upward and downward seepage were recorded along every transect across the river with rates ranging from + 2·37 to ? 3·40 m/d. Despite a stable bed, which commonly facilitates clogging by fine‐grained or organic sediments, seepage rates at the Allegheny River were not reduced relative to those at the South Platte River. Seepage rate and direction depended primarily on measurement position relative to local‐ and meso‐scale bed topography at both rivers. Hydraulic gradients were small at nearly all seepage‐measurement locations and commonly were not a good indicator of seepage rate or direction. Therefore, measuring hydraulic gradient and hydraulic conductivity at in‐stream piezometers may be misleading if used to determine seepage flux across the sediment‐water interface. Such a method assumes that flow between the well screen and sediment‐water interface is vertical, which appears to be a poor assumption in coarse‐grained hyporheic settings. Copyright © 2011 John Wiley & Sons, Ltd.     

Dimensionless critical shear stress evaluation from flume experiments using different gravel beds

Flume experiments were conducted using four different gravel beds (D₅₀ + 12–39 mm) and a range of marked particles (10–65 mm). The shear stresses were evaluated from friction velocities, when initial movement of marked particles occurred. Two kinds of equations were produced: first for the threshold of initial movement, and second for generalized movement. Equations of the type 0_c + a(D_i/D₅₀)^b, as proposed by Andrews (1983) are applicable even if the material is relatively well sorted. However, the values of a and b are lower (respectively 0·050 and -0·70) for initial movement. Generalized movement requires a higher shear stress (a + 0·068 and b + -0·80). D₉₀ of the bed material and y₀ (the bed roughness parameter) were also used as reference values in place of D₅₀. They produced lower values than in natural streams, mainly owing to the fact that the material used in the flume is better sorted: clusters are less well developed and the bed roughness is lower.     

High frequency stream bed mobility of a low‐gradient agricultural stream with implications on the hyporheic zone

Little Kickapoo Creek (LKC), a low‐gradient stream, mobilizes its streambed–fundamentally altering its near‐surface hyporheic zone–more frequently than do higher‐gradient mountain and karst streams. LKC streambed mobility was assessed through streambed surveys, sediment sampling, and theoretical calculations comparing basal shear stress (τ_b) with critical shear stress (τ_c). Baseflow τ_b is capable of entraining a d₅₀ particle; bankfull flow could entrain a 51·2 mm particle. No particle that large occurs in the top 30 cm of the substrate, suggesting that the top 30 cm of the substrate is mobilized and redistributed during bankfull events. Bankfull events occur on average every 7·6 months; flows capable of entraining d₅₀ and d₈₅ particles occur on average every 0·85 and 2·1 months, respectively. Streambed surveys verify streambed mobility at conditions below bankfull. While higher gradient streams have higher potential energy than LKC, they achieve streambed‐mobilization thresholds less frequently. Heterogeneous sediment redistribution creates an environment where substrate hydraulic conductivity (K) varies over four orders of magnitude. The frequency and magnitude of the substrate entrainment has implications on hyporheic zone function in fluid, solute and thermal transport models, interpretations of hyporheic zone stability, and understanding of LKC's aquatic ecosystem. Copyright © 2008 John Wiley & Sons, Ltd.     

ADV measurements of velocity distributions in a gravel‐bed flume

Velocity measurements carried out by an acoustic doppler velocimeter (ADV) in a rectangular laboratory ?ume having a gravel bed are presented. The velocity pro?les are measured in six verticals of the channel cross‐section having an increasing distance (from 4 to 38·5 cm) from the ?ume wall. The experimental runs are carried out for ?ve different bed arrangements, characterized by different concentrations of coarser elements, and for the two conditions of small‐ and large‐scale roughness. For both hydraulic conditions, the velocity measurements are ?rst used to test the applicability of the Dean pro?le and of the logarithmic pro?le corrected by a divergence function proposed in this paper. Then, for each value of the depth sediment ratio h/d₈₄, the non‐dimensional friction factor parameter is calculated by integration of the measured velocity distributions in the different verticals of the cross‐section. Finally a semi‐logarithmic ?ow resistance equation is empirically deduced. Copyright © 2003 John Wiley & Sons, Ltd.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

A casting procedure for reproducing coarse‐grained sedimentary surfaces

In the ?eld, the measurement of near‐bed hydraulics remains problematic. Greater precision is possible in the laboratory, but, in the case of gravels, it is dif?cult to create a water‐worked channel‐bed that is realistic enough to replicate faithfully the conditions found in nature. In this paper, a technique to reproduce coarse‐grained sedimentary fabrics of large areal extent is described. It involves moulding natural river‐bed surfaces from which facsimiles are cast. Remarkably realistic casts with dimensions of 1 m by 2 m have been produced and their quality assessed using spatial data derived using automated digital photogrammetry. The casts reproduce the prototype surfaces with errors at millimetre scale (0·5 per cent of the microrelief). The technique has facilitated the introduction of sedimentary surfaces that incorporate natural, complex structures of grains up to cobble size into experimental channels where detailed studies of near‐bed hydraulics can be carried out. Copyright © 2003 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.     

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.     

Effect of bed surface roughness on longitudinal dispersion in artificial open channels

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.     

Shear velocity estimates in rough‐bed open‐channel flow

Shear velocity u_* is an important parameter in geophysical flows, in particular with respect to sediment transport dynamics. In this study, we investigate the feasibility of applying five standard methods [the logarithmic mean velocity profile, the Reynolds stress profile, the turbulent kinetic energy (TKE) profile, the wall similarity and spectral methods] that were initially developed to estimate shear velocity in smooth bed flow to turbulent flow over a loose bed of coarse gravel (D₅₀ = 1·5 cm) under sub‐threshold conditions. The analysis is based on quasi‐instantaneous three‐dimensional (3D) full depth velocity profiles with high spatial and temporal resolution that were measured with an Acoustic Doppler Velocity Profiler (ADVP) in an open channel. The results of the analysis confirm the importance of detailed velocity profile measurements for the determination of shear velocity in rough‐bed flows. Results from all methods fall into a range of ± 20% variability and no systematic trend between methods was observed. Local and temporal variation in the loose bed roughness may contribute to the variability of the logarithmic profile method results. Estimates obtained from the TKE and Reynolds stress methods reasonably agree. Most results from the wall similarity method are within 10% of those obtained by the TKE and Reynolds stress methods. The spectral method was difficult to use since the spectral energy of the vertical velocity component strongly increased with distance from the bed in the inner layer. This made the choice of the reference level problematic. Mean shear stress for all experiments follows a quadratic relationship with the mean velocity in the flow. The wall similarity method appears to be a promising tool for estimating shear velocity under rough‐bed flow conditions and in field studies where other methods may be difficult to apply. This method allows for the determination of u_* from a single point measurement at one level in the intermediate range (0·3 < h < 0·6). Copyright © 2013 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.     

Persistence of the surface texture of a gravel‐bed river during a large flood

We present herein clear field evidence for the persistence of a coarse surface layer in a gravel‐bed river during flows capable of transporting all grain sizes present on the channel bed. Detailed field measurements of channel topography and bed surface grain size were made in a gravel‐bed reach of the Colorado River prior to a flood in 2003. Runoff produced during the 2003 snowmelt was far above average, resulting in a sustained period of high flow with a peak discharge of 27 m³/s (170% of normal peak flow); all available grain sizes within the study reach were mobilized in this period of time. During the 2003 peak flow, the river avulsed immediately upstream of the study reach, thereby abandoning approximately one half kilometer of the former channel. The abandonment was rapid (probably within a few hours), leaving the bed texture essentially frozen in place at the peak of the flood. All locations sampled prior to the flood were resampled following the stream abandonment. In response to the high flow, the surface median grain size (D_50s) coarsened slightly in the outer part of the bend while remaining nearly constant along the inner part of the bend, resulting in an overall increase from 18 to 21 mm for the study reach. Thus, the coarse bed surface texture persisted despite shear stresses throughout the bend that were well above the critical entrainment value. This may be explained because the response of the bed texture to increases in flow strength depends primarily upon the continued availability of the various grain size percentiles in the supply, which in this case was essentially unlimited for all sizes present in the channel. Copyright © 2007 John Wiley & Sons, Ltd.     

The impact of stress history on bed structure

Recent research has started to focus on how prolonged periods of sub‐threshold flows may be capable of imparting structural changes that contribute to increased bed stability. To date, this effect (termed ‘stress history’) has been found to be significant in acting to increase a bed's critical shear stress at entrainment threshold. However, it is supported by only limited, qualitative and often speculative information on the mechanisms of this stabilization process in grade‐specific studies. As such, this paper uses high resolution laser scanning to quantitatively ascertain the granular mechanics underpinning the relationship between stress history and entrainment threshold for beds of a range of grain size distributions. Employing a bed slope of 1/200, three grain size distributions with median grain sizes (D₅₀) of 4·8 mm [uniform (σ_g = (D₈₄/D₁₆)^0.5 = 1·13; bimodal (σ_g = 2·08); and, unimodal (σ_g = 1·63)] were exposed to antecedent stress histories of 60 and 960 minutes duration. Antecedent shear stress magnitude was set at 50% of the critical shear stress for the D₅₀ when no stress history period was employed. Two laser displacement scans of the bed surface (approximate area 100 mm × 117 mm) were taken, one prior to the antecedent period and one after this period, so that changes to surface topography could be quantified (resolution of x = 0·10 mm, y = 0·13 mm and z = 0·24 mm). Rearrangement of bed surface structure is described using statistical analysis and two‐dimensional (2D) semi‐variograms to analyse scaling behaviour. Results reveal vertical settlement, changes to bed roughness and particle repositioning. However, the bed grain size distribution influences the relative importance of each mechanism in determining stress history induced bed stability; this is the focus of discussion in this paper. Copyright © 2012 John Wiley & Sons, Ltd.     

A wind tunnel investigation of the influences of fetch length on the flux profile of a sand cloud blowing over a gravel surface

Wind tunnel tests were conducted to examine the fetch effect of a gravel surface on the ?ux pro?le of the sand cloud blowing over it using typical dune sand. The results suggest that the ?ux pro?le of blown sand over a gravel surface differs from that over a sandy surface and is characterized by a peak ?ux at a height above the surface while that over a sandy surface decreases exponentially with height. The ?ux pro?le of a sand cloud over a gravel surface can be expressed by a Gaussian peak function: q = a + b exp (?0·5((h ? c)/d)²), where q is the sand transport rate at height h, and a, b, c and d are regression coef?cients. The signi?cance of the coef?cients in the function could be de?ned in accordance with the fetch length of the gravel surface and wind velocity. Coef?cient c represents the peak ?ux height and increases with both wind velocity and fetch length, implying that the peak ?ux height is related to the bounce height of the particles in the blowing sand cloud. Coef?cient d shows a tendency to increase with both wind velocity and fetch length. The sum of a and b, representing the peak ?ux, increases with wind velocity but decreases with fetch length. The average saltation height derived from the cumulative percentage curve shows a tendency to increase with both the fetch length and wind velocity. For any fetch length of a gravel surface the sand transport equation is expressed as Q = C(1 ? U_t/U)(ρ/g)U³, where Q is the sand transport rate, U is the wind velocity, U_t is the threshold velocity measured at the same height as U, g is the gravitational acceleration, ρ is the air density, C is a proportionality coef?cient that decreases with the fetch length of the gravel surface. At a given wind velocity, the sand transport rate over a gravel surface is only 52–68 per cent of that over a sandy surface. The ?ux rate in true creep over a gravel surface increases with wind velocity but decreases with the fetch length, whereas the creep proportion (the ratio of creep ?ux to the sand transport rate) decreases with both the wind velocity and fetch length. Two‐variable (including fetch length and wind velocity) equations were developed to predict the peak ?ux height, average saltation height and transport rate. Copyright © 2004 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.     

Displacement and transport of marked pebbles,cobbles and boulders during floods in a steep mountain stream

In a small experimental catchment of the Dolomites (Rio Cordon, 5 km²) field observations have been carried out on the movement of various sized bed material particles. Displacement length of 860 marked pebbles, cobbles and boulders (0·032 < D < 0·512 m) has been measured along the river bed during individual snowmelt and flood events in the periods 1993–1994 and 1996–1998. Floods were grouped into two categories. The first includes ‘ordinary’ events, which are characterized by peak discharges with a return period of 1–5 years and by an hourly bedload rate not exceeding 20 m³ h^?1. The second refers to ‘exceptional’ events with a return period of 50–60 years. A flood of this latter type occurred on 14 September 1994, with a peak discharge of 10·4 m³ s^?1 and average hourly bedload rate of 324 m³ h^?1. The variation according to grain size of total displacement length L_i depends on the degree of mobilization of the individual fractions of the bed surface: L_i is independent of D_i for smaller, fully mobile grain sizes and decreases rapidly with D_i for larger fractions in a state of partial transport. Sustained selective transport without a supply of sediment from upstream leads to the development of a stable coarse armoured surface through progressive winnowing of finer material from the bed surface. With supply unlimited conditions for transport, both the occurrence of extreme events and the duration of a sequences of ‘ordinary’ floods play an important role in the degree of mobilization of the individual fractions of the bed. Copyright © 2004 John Wiley & Sons, Ltd.     

Effect of changes in fine‐grained matrix on bedload sediment transport in a gravel‐bed river

While clay and silt matrices of gravel‐bed rivers have received attention from ecologists concerned variously with the deteriorating environments of benthic and hyporheic organisms, their impact on sediment entrainment and transport has been explored less. A recent increase of such a matrix in the bed of Nahal Eshtemoa, an ephemeral river of the northern Negev, has more than doubled the boundary shear stress needed to initiate bedload, from 7 N m^‐2 (τ^* = 0.027) during the flash floods of 1991–2001 to 15 N m^‐2 (τ^* = 0.059) during those of 2008–2009. The relation between bedload flux and boundary shear stress continues to be well‐defined, but it is displaced. The matrix now contains a significant amount of silt and clay size material. The reasons for the increased entrainment threshold of bedload are explored. Large‐scale laser scanning of the dry bed reveals a reduction in grain‐scale morphological roughness, while artificial in situ tests of matrix integrity indicate considerable cohesion. The implications for adopting bed material sampling strategies that account for matrix development are assessed. Copyright © 2012 John Wiley & Sons, Ltd.