Local‐scale variability of seepage and hydraulic conductivity in a shallow gravel‐bed river

Seepage rate and direction measured with a seepage metre modified for use in flowing water were greatly variable along a 300‐m reach of a shallow, gravel‐bed river and depended primarily on the local‐scale bed topography. The median value of seepage measured at 24 locations was 24 cm/day, but seepage measured at specific sites ranged from ?340 to +237 cm/day. Seepage also varied substantially over periods of hours to days and occasionally reversed direction in response to evolution of the sediment bed. Vertical hydraulic conductivity was related to seepage direction and was larger during upward seepage than during downward seepage; with differences ranging from 4 to 40% in areas of active sediment transport to more than an order of magnitude in areas where current was too slow to mobilize bed sediment. Seepage was poorly related to hydraulic gradient measured over vertical distances of 0·3 m and appeared to be opposite the hydraulic gradient at 18% of the locations where both parameters were measured. Results demonstrate the scale dependence of these measurements in coarse‐grained hyporheic settings and indicate that hydraulic gradients should be determined over a much shorter vertical increment if used to indicate exchange across the sediment–water interface. Published in 2009 by 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.     

Spatial and temporal intensification of lateral hyporheic flux in narrowing intra‐meander zones

Meander bends in alluvial rivers morphologically evolve towards meander cut‐off with narrowing intra‐meander necks, and this should steepen hydraulic gradients and intensify intra‐meander hyporheic flux. This research used dye tracking and head loss measurements in a 1 : 500 planimetrically scaled laboratory river table to quantify the spatial and temporal intensification of intra‐meander flux rates at two evolution ages. The younger meander bend, M1, had a sinuosity of 2.3, a river neck width of 0.39 cm, and 0.6% river slope, and the older meander bend, M3, had a sinuosity of 5.2, a river neck width of 0.12 cm, and 0.5% river slope. Flux into and out of the meander bend was estimated along the normalized curvilinear distance s*, with the meander neck at s* = 0.1 and s* = 0.9, the meander centroid at s* = 0.37 and s* = 0.63, and the apex at s* = 0.5. Between the meander centroid and neck, we documented a 60% spatial intensification for M1 and a 90% spatial intensification for M3. Between M1 and M3, we documented a 135% temporal intensification at the neck and a 100% intensification at the centroid. Our empirical spatial and temporal intensification rates involving the M1 and the M3 scenario were one to three times lower than theoretical rates derived from a river evolution model with equivalent M1 and M3 planimetry. Overestimation by the theoretical model was attributed to exaggerated head loss caused by the model neglecting groundwater contributions to river stage. Hyporheic exchange provides critical ecosystem services, and its spatial and temporal variation with meander evolution should be considered in river management. Copyright © 2012 John Wiley & Sons, Ltd.     

Prospective modelling of 3D hyporheic exchange based on high‐resolution topography and stream elevation

River managers and scientists interested in hyporheic processes need adequate tools for characterizing hyporheic exchange flow (HEF) at local sites where only poor information on subsurface properties are available. This study evaluates a three‐dimensional modelling approach, on the basis of detailed surface parameterization and a simplified subsurface structure, for comparison of potential HEF characteristics at three experimental reaches at the channel‐unit scale. First, calibration is conducted to determine the best fit‐of‐heads given the model simplification, then the structure of residuals are used to evaluate the origin of the misfit, and finally, a sensitivity analysis is conducted to identify inter‐site differences in HEF. Results show that such an approach can highlight potential magnitude differences in HEF characteristics between reaches. The sensitivity analysis is successful in delineating the small area of exchange that remains under conditions of high groundwater discharge. In this case, however, the calibrated model performs poorly in representing the exchange pattern at the sediment–water interface, thus suggesting that the approach is less adequate for a deterministic simulation of observed heads. The summary statistics are in the range of similar published models, for which the reported indicator is the root mean square error on heads normalized by the head drop over the reach. We recommend, however, that modellers use a more comparable indicator, such as a measure of the residuals normalized by a measure of observed vertical head differences. Overall, when subsurface data are unavailable or sparse, a three‐dimensional groundwater model based on high‐resolution topographic data combined with a sensitivity analysis appears as a useful tool for a preliminary characterization of HEF. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of the hyporheic zone on the phosphorus dynamics of a large gravel‐bed river,Garonne River,France

Phosphorus (P) concentrations in sediments and in surface and interstitial water from three gravel bars in a large river (Garonne River, southern France) were measured daily, downstream of a wastewater treatment plant for a city of 740 000 inhabitants (Toulouse). Measurements were made of vertical hydraulic gradient (VHG), total dissolved phosphorus (TDP), soluble reactive phosphorus (SRP) and total phosphorus (TP) in water and of three extractable forms of phosphorus (water extractable, NaOH extractable and H₂SO₄ extractable) in hyporheic sediments from the gravel bars. Dissolved phosphorus was the major contributor to TP (74–79%) in both interstitial and surface waters on all sampling dates, and in most cases surface water P concentrations were significantly higher than interstitial concentrations. Hyporheic sediment TP concentrations ranged between 269 and 465 µg g^?1 and were highest in fine sediment fractions. Acid‐extractable P, a non‐bioavailable form, represented at least 95% of sediment TP. A positive relationship was observed between VHG and TP in two of the gravel bars, with wells that were strongly downwelling having lower TP concentrations. These results suggest that in downwelling zones, hyporheic sediments can trap surface‐derived dissolved P, and that much of this P becomes stored in refractory particulate forms. Bioavailable P is mainly present in dissolved form and only occupies a small fraction of total P, with particulate P comprising the majority of total P. Copyright © 2009 John Wiley & Sons, Ltd.     

The effects of discharge and slope on hyporheic flow in step‐pool morphologies

This paper focuses on surface–subsurface water exchange in a steep coarse‐bedded stream with a step‐pool morphology. We use both flume experiments and numerical modelling to investigate the influence of stream discharge, channel slope and sediment hydraulic conductivity on hyporheic exchange. The model step‐pool reach, whose topography is scaled from a natural river, consists of three step‐pool units with 0.1‐m step heights, discharges ranging between base and over‐bankfull flows (scaled values of 0.3–4.5 l/s) and slopes of 4% and 8%. Results indicate that the deepest hyporheic flow occurs with the steeper slope and at moderate discharges and that downwelling fluxes at the base of steps are highest at the largest stream discharges. In contrast to findings in a pool‐riffle morphology, those in this study show that steep slopes cause deeper surface–subsurface exchanges than gentle slopes. Numerical simulation results show that the portion of the hyporheic zone influenced by surface water temperature increases with sediment hydraulic conductivity. These experiments and numerical simulations emphasize the importance of topography, sediment permeability and roughness elements along the channel surface in governing the locations and magnitude of downwelling fluxes and hyporheic exchange. Our results show that hyporheic zones in these steep streams are thicker than previously expected by extending the results from streams with pool‐riffle bed forms. Copyright © 2014 John Wiley & Sons, Ltd.     

Active width of gravel‐bed braided rivers

Previous analyses have identified the active width of braided rivers, the bed area over which bed load flux and short‐term morphological change occurs, as an important element of braiding dynamics and predictions of bed load flux. Here we compare theoretical predictions of active width in gravel‐bed braided rivers with observations from Sunwapta River, and from a generic physical model of gravel braided rivers, to provide general observations of the variation in active width, and to develop an understanding of the causes of variation. Bed topography was surveyed daily along a 150 m reach of the pro‐glacial Sunwapta River for a total of four weeks during summer when flow was above threshold for morphological activity. In the laboratory, detailed digital elevation models (DEMs) were derived from photogrammetric survey at regular intervals during a constant discharge run. From the field and flume observations there is considerable local and circumstantial variation in active width, but also a general trend in average active width with increasing discharge. There is also a clear relationship of active width with active braiding index (number of active branches in the braided channel network), and with dimensionless stream power, which appears to be consistent across the range of data from field and physical models. Thus there is a link between active width and the river morphology and dynamics, and the possibility of a general relationship for estimating active width from channel pattern properties or reach‐scale stream power values, from which approximate bedload flux calculations may be made. The analysis also raises questions about differences between hydraulically‐based numerical model computations of instantaneous active width and observation of time‐integrated morphological active width. Understanding these differences can give insight into the nature of bedload transport in braided rivers and the relationship to morphological processes of braiding. Copyright © 2011 John Wiley & Sons, Ltd.     

Determining hydraulic resistance in gravel‐bed rivers from the dynamics of their water surfaces

Traditionally, approaches to account for the effect of the boundary roughness of a gravel‐bed river have used a grain‐size index of the bed surface as a surrogate for hydraulic resistance. The use of a single grain‐size does not take into account the spatial heterogeneity in the bed surface and how this heterogeneity imparts resistance on the flow, nor the way in which this relationship changes with variables such as flow stage. A new technique to remotely quantify hydraulic resistance is proposed. It is based on measuring the dynamics of a river's water surface and relating this to the actual hydraulic resistance created by a rough sediment boundary. The water surface dynamics are measured using a new acoustic technique, grazing angle sound propagation (GRASP). This proposed method to measure hydraulic resistance is based on a greater degree of physical reasoning, and this is discussed in the letter. By measuring acoustically the temporal dynamics of turbulent water surfaces over a water‐worked gravel bed in a laboratory flume, a dependency is demonstrated between the temporal variation in the reflected acoustic pressure and measured hydraulic resistance. It is shown that the standard deviation in acoustic pressure decreases with increasing hydraulic resistance. This is shown to apply for a range of relative submergences and bed slopes that are typical of gravel‐bed rivers. This remote sensing technique is both rapid and inexpensive, and has the potential to be applied to natural river channels and to other environmental turbulent flows, such as overland flows. A whole new class of low‐cost, remote and non‐intrusive instruments could be developed as a result and used in a wide range of hydraulic and hydrological applications. Copyright © 2006 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.     

Seasonal variation in cascade‐driven hyporheic exchange,northern Honduras

A characterization of hyporheic exchange for dry and wet season baseflow, as well as partially dewatered discharge, was done in Prieta Creek, a first‐order cascade in northern Honduras. The cascade had discharges from 1 to 15 l s^?1, had average slopes of 12%, pool spacing of 3 m, and shallow substrate of sand and gravel. Tracer tests were conducted in a 15‐m sub‐reach, a length considered to be adequate for the experiment based on the DaI test, a ratio of exchange and transport processes. In the three tests, between 9 and 18% of tracer was not recovered, possibly due to entrainment in flowpaths passing beneath the downstream monitoring location. Tracer data were analysed by the one‐dimensional transport with inflow and storage (OTIS) transient storage model (TSM) to derive standard exchange parameters, and by the solute transport in rivers (STIR) model to examine hyporheic residence time distributions (RTDs). The best fit of the observed tracer breakthrough curves was obtained by using the STIR model with a combination of two exponential RTDs to represent hyporheic retention. With increasing discharge, the OTIS model predicted increasing storage exchange fluxes and exchange coefficients and decreasing storage zone areas and transient storage times, which are trends supported by riparian and streambed piezometric head data. Riparian water levels rose during the transition from the dry to wet season, which could constrict the hyporheic storage zone. Thirteen of the 19 streambed piezometers recorded seasonal changes in hydraulic gradients and flux direction, with fewer yet stronger upwelling zones during higher discharges. The MODFLOW model missed the observed seasonal changes, possibly due to subtle changes in the seasonal change in water surface profiles. We conclude that partially dewatered dry season exchange, compared to wet season exchange, was initiated and terminated with smaller pressure gradients and, in different streambed locations, was smaller in volume, had longer residence times, and may connect with deeper and longer flow paths. Copyright © 2010 John Wiley & Sons, Ltd.     

Regime theories in gravel‐bed rivers: models,controlling variables,and applications in disturbed Italian rivers

Downstream hydraulic geometry relationships describe the shape of alluvial channels in terms of bankfull width, flow depth, flow velocity, and channel slope. Recent investigations have stressed the difference in spatial scales associated with these variables and thus the time span required for their adjustment after a disturbance. The aim of this study is to explore the consequences in regime models considering the hypothesis that while channel width and depth adjust quickly to changes in water and sediment supply, reach slope requires a longer time span. Three theoretical models were applied. One model incorporates an extremal hypothesis (Millar RG. 2005. Theoretical regime equations for mobile gravel‐bed rivers with stable banks. Geomorphology 64 : 207–220), and the other two are fully physically based (Ikeda S, Parker G, Kimura Y. 1988. Stable width and depth of straight gravel rivers with heterogeneous bed materials. Water Resources Research 24 : 713–722; Parker G, Wilcock PR, Paola C, Dietrich W, Pitlick J. 2007. Physical basis for quasi universal relations describing bankfull hydraulic geometry of single‐thread gravel‐bed rivers. Journal of Geophysical Research 112 , DOI: 10.1029/2006JF000549). In order to evaluate the performance of models introducing the slope as an independent variable, we propose two modifications to previous models. The performance of regime models was tested against published data from 142 river reaches and new hydraulic geometry data from gravel‐bed rivers in Patagonia (Argentina) and north‐eastern Italy. Models that assume slope as a control (Ikeda et al., 1988; or Millar, 2005) predict channel depth and width reasonably well. Parker et al.'s (2007) model improved predictions because it filters the scatter in slope data with a relation slope–discharge. The extremal hypothesis model of Millar (2005) predicts comparably to the other physically based models. Millar's model was chosen to describe the recent changes in the Piave and Brenta rivers due to human intervention – mainly in‐channel gravel mining. The change in sediment supply and recovery was estimated for these rivers. This study supports the interpretation that sediment supply is the key factor guiding morphological changes in these rivers. Copyright © 2013 John Wiley & Sons, Ltd.     

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.     

The spatial and temporal patterns of aggradation in a temperate,upland, gravel‐bed river

Intensive field monitoring of a reach of upland gravel‐bed river illustrates the temporal and spatial variability of in‐channel sedimentation. Over the six‐year monitoring period, the mean bed level in the channel has risen by 0·17 m with a maximum bed level rise of 0·5 m noted at one location over a five month winter period. These rapid levels of aggradation have a profound impact on the number and duration of overbank flows with flood frequency increasing on average 2·6 times and overbank flow time increasing by 12·8 hours. This work raises the profile of coarse sediment transfer in the design and operation of river management, specifically engineering schemes. It emphasizes the need for the implementation of strategic monitoring programmes before engineering work occurs to identify zones where aggradation is likely to be problematic. Exploration of the sediment supply and transfer system can explain patterns of channel sedimentation. The complex spatial, seasonal and annual variability in sediment supply and transfer raise uncertainties into the system's response to potential changes in climate and land‐use. Thus, there is a demand for schemes that monitor coarse sediment transfer and channel response. Copyright © 2009 John Wiley & Sons, Ltd.     

The variability in the morphological active width: Results from physical models of gravel‐bed braided rivers

The morphological active width, defined as the lateral extent of bed material displacement over time, is a fundamental parameter in multi‐threaded gravel‐bed rivers, linking complex channel dynamics to bedload transport. Here, results are presented from five constant discharge experiments, and three event hydrographs, covering a range of flow strengths and channel configurations for which morphological change, bedload transport rates, and stream power were measured in a physical model. Changes in channel morphology were determined via differencing of photogrammetrically‐derived digital elevation models (DEMs) of the model surface generated at regular intervals over the course of ~115 h of experimental runs. Independent measures of total bedload output were made using downstream sediment baskets. Results indicate that the morphological active width increases with total and dimensionless stream power and is strongly and positively correlated with bulk change (total volume of bed material displaced over time) and active braiding intensity (ABI). Although there is considerable scatter due to the inherent variability in braided river morphodynamics, the active width is positively correlated with independent measurements of bedload transport rate. Active width, bulk change, and bedload transport rates were all negligible below a dimensionless stream power threshold value of ~ 0.09, above which all increase with flow strength. Therefore, the active width could be used as a general predictor of bulk change and bedload transport rates, which in turn could be approximated from total and dimensionless stream power or ABI in gravel‐bed braided rivers. Furthermore, results highlight the importance of the active width, rather than the morphological active depth, in predicting volumes of change and bedload transport rates. The results contribute to the larger goals of better understanding braided river morphodynamics, creating large high‐resolution datasets of channel change for model calibration and validation, and developing morphological methods for predicting bedload transport rates in braiding river systems. Copyright © 2018 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.     

Heat tracing of heterogeneous hyporheic exchange adjacent to in‐stream geomorphic features

There are many field techniques used to quantify rates of hyporheic exchange, which can vary in magnitude and direction spatially over distances of only a few metres, both within and between morphological features. We used in‐stream mini‐piezometers and heat transport modelling of stream and streambed temperatures to quantify the rates and directions of water flux across the streambed interface upstream and downstream of three types of in‐stream geomorphic features: a permanent dam, a beaver dam remnant and a stream meander. We derived hyporheic flux estimates at three different depths at six different sites for a month and then paired those flux rates with measurements of gradient to derive hydraulic conductivity (K) of the streambed sediments. Heat transport modelling provided consistent daily flux estimates that were in agreement directionally with hydraulic gradient measurements and also identified vertical heterogeneities in hydraulic conductivity that led to variable hyporheic exchange. Streambed K varied over an order of magnitude (1·9 × 10^?6 to 5·7 × 10^?5 m/s). Average rates of hyporheic flux ranged from static (q < ±0·02 m/day) to 0·42 m/day. Heat transport modelling results suggest three kinds of flow around the dams and the meander. Copyright © 2010 John Wiley & Sons, Ltd.     

Distribution of clay-sized sediments in streambeds and influence of fine sediment clogging on hyporheic exchange

In this work, the deposition of clay-sized fine particles (d₅₀ = 0.006 mm) and its subsequent influence on the dune-induced hyporheic exchange are investigated. Fine sand (D₅₀ = 0.28 mm), coarse sand (D₅₀ = 1.7 mm), and gravel (D₅₀ = 5.5 mm) grains were used to form homogenous model streambeds; one control - no clay input, and two treatments - increasing clay inputs for each grain type. The results indicate that the clogging profiles of clay-sized sediments may not be predicted accurately using the previously proposed metric based on the relative sizes of infiltrating and substrate sediments. Further, the depositional patterns vary with the initial concentration of clay particles in the surface water. The assessment of clogging profiles in coarse-grained model streambeds also reveals a preferential infiltration of the clay particles in the hyporheic downwelling regions. The results from the dye tracer test suggest that the accumulation of clay particles altered the exchange characteristics in the treatment flumes. For each grain size, the treatment flumes exhibit lower hyporheic flux and higher median residence times compared to their respective control flumes. The dye penetration depths were lower in treatment flumes with fine and coarse sand compared to their respective control flumes. Interestingly, higher penetration depths were observed in treatment flumes with gravel compared to their respective control flume potentially due to the generation of preferential flow paths in the partially clogged gravel beds. The clogging altered the hyporheic fluxes and residence times in the coarse-grained model beds to a greater degree in comparison to the fine sand beds. Overall, our findings indicate that the properties of both fine and substrate sediments influence the clogging patterns in streambeds, and the subsequent influence of fine sediment clogging on hyporheic exchange and associated processes may vary across stream ecosystems.     

Effects of sand content on initial gravel motion in gravel‐bed rivers

When fine sediments are present in gravel streambeds (gravel‐framework beds), the gravel can be more easily removed from its original position, compared with gravel in a streambed without fine sediment but otherwise under the same hydraulic conditions. In this study, the effect of the presence of sand on the initiation of gravel motion in gravel riverbeds was investigated using flume experiments. The relationship between the critical Shields stress for gravel motion initiation and the fraction of sand in the bed was determined experimentally. The results can be summarized as follows. (1) When the fraction of sand in the bed is smaller than about 0.4, the critical Shields stress for the initiation of gravel motion decreases with increasing fraction of sand. The critical Shields stress increases, however, with increasing fraction of sand when it is larger than about 0.4. (2) The difference between the value of the critical Shields stress predicted by the Egiazaroff equation and the value obtained from the experimental data becomes maximum at about 0.4 of the fraction of sand. Here an empirical relation between the critical Shields stress and the fraction of sand is proposed so as to consider the effects of the ratio of the characteristic gravel size to the mean size of the bed material on the critical Shields stress. (3) Gravel in armored beds can be more easily mobilized by supplying sand as part of a sediment augmentation scheme. The sand fraction in the subsurface layer of the bed appears to reduce the friction angle of exposed particles. Sediment augmentation using sand has been recently demonstrated to be a viable alternative for mobilizing gravel for the restoration of gravel‐bed rivers downstream of dams. The quantitative evaluation obtained through the experiments reported here may be useful for the design of augmentation schemes. Copyright © 2017 John Wiley & Sons, Ltd.     

A coupled sediment routing and lateral migration model for gravel‐bed rivers

We describe additions made to a multi‐size sediment routing model enabling it to simulate width adjustment simultaneously alongside bed aggradation/incision and fining/coarsening. The model is intended for use in single thread gravel‐bed rivers over annual to decadal timescales and for reach lengths of 1–10 km. It uses a split‐channel approach with separate calculations of flow and sediment transport in the left and right sides of the channel. Bank erosion is treated as a function of excess shear stress with bank accretion occurring when shear stress falls below a second, low, threshold. A curvature function redistributes shear stress to either side of the channel. We illustrate the model through applications to a 5·6‐km reach of the upper River Wharfe in northern England. The sediment routing component with default parameter values gives excellent agreement with field data on downstream fining and down‐reach reduction in bedload flux, and the width‐adjustment components with approximate calibration to match maximum observed rates of bank shifting give plausible patterns of local change. The approach may be useful for exploring interactions between sediment delivery, river management and channel change in upland settings. Copyright © 2011 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.