Flow modelling in gravel‐bed rivers: rethinking the bottom boundary condition

A key problem in computational fluid dynamics (CFD) modelling of gravel‐bed rivers is the representation of multi‐scale roughness, which spans the range from grain size, through bedforms, to channel topography. These different elements of roughness do not clearly map onto a model mesh and use of simple grain‐scale roughness parameters may create numerical problems. This paper presents CFD simulations for three cases: a plane bed of fine gravel, a plane bed of fine gravel including large, widely‐spaced pebble clusters, and a plane gravel bed with smaller, more frequent, protruding elements. The plane bed of fine gravel is modelled using the conventional wall function approach. The plane bed of fine gravel including large, widely‐spaced pebble clusters is modelled using the wall function coupled with an explicit high‐resolution topographic representation of the pebble clusters. In these cases, the three‐dimensional Reynolds‐averaged continuity and Navier–Stokes equations are solved using the standard k ? ε turbulence model, and model performance is assessed by comparing predicted results with experimental data. For gravel‐bed rivers in the field, it is generally impractical to map the bed topography in sufficient detail to enable the use of an explicit high‐resolution topography. Accordingly, an alternative model based on double‐averaging is developed. Here, the flow calculations are performed by solving the three‐dimensional double‐averaged continuity and Navier‐Stokes equations with the spatially‐averaged 〈k ? ε〉 turbulence model. For the plane bed of fine gravel including large, widely‐spaced pebble clusters, the model performance is assessed by comparing the spatially‐averaged velocity with the experimental data. The case of a plane gravel bed with smaller, more frequent, protruding elements is represented by a series of idealized hypothetical cases. Here, the spatially‐averaged velocity and eddy viscosity are used to investigate the applicability of the model, compared with using the explicit high‐resolution topography. The results show the ability of the model to capture the spatially‐averaged flow field and, thus, illustrate its potential for representing flow processes in natural gravel‐bed rivers. Finally, practical data requirements for implementing such a model for a field example are given. Copyright © 2011 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.     

Gravel‐bed river evolution in earthquake‐prone regions subject to cycled hydrographs and repeated sediment pulses

Sediment often enters rivers in the form of sediment pulses associated with landslides and debris flows. This is particularly so in gravel‐bed rivers in earthquake‐prone mountain regions, such as Southwest China. Under such circumstances, sediment pulses can rapidly change river topography and leave the river in repeated states of gradual recovery. In this paper, we implement a one‐dimensional morphodynamic model of river response to pulsed sediment supply. The model is validated using data from flume experiments, so demonstrating that it can successfully reproduce the overall morphodynamics of experimental pulses. The model is then used to explore the evolution of a gravel‐bed river subject to cycled hydrographs and repeated sediment pulses. These pulses are fed into the channel in a fixed region centered at a point halfway down the calculational domain. The pulsed sediment supply is in addition to a constant sediment supply at the upstream end. Results indicate that the river can reach a mobile‐bed equilibrium in which two regions exist within which bed elevation and surface grain size distribution vary periodically in time. One of these is at the upstream end, where a periodic discharge hydrograph and constant sediment supply are imposed, and the other is in a region about halfway down the channel where periodic sediment pulses are introduced. Outside these two regions, bed elevation and surface grain size distribution reach a mobile‐bed equilibrium that is invariant in time. The zone of fluctuation‐free mobile‐bed equilibrium upstream of the pulse region is not affected by repeated sediment pulses under the scenarios tested, but downstream of the pulse region, the channel reaches different fluctuation‐free mobile‐bed equilibriums under different sediment pulse scenarios. The vertical bed structure predicted by the simulations indicates that the cyclic variation associated with the hydrograph and sediment pulses can affect the substrate stratigraphy to some depth. Copyright © 2017 John Wiley & Sons, Ltd.     

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

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

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

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.     

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.     

Evolution of gravel surfaces in a sediment‐recirculating flume

The assumption of the equilibrium state of gravel surfaces in flume experiments under feeding or recirculating conditions is generally justified by three equilibrium criteria based on sediment transport, slopes, and bed features. When these parameters become stable, an experiment is expected to reach equilibrium. This equilibrium state, however, is based on a one‐dimensional model, the Exner equation, which may not truly reflect the equilibrium state of the system considering the complex flow and sediment processes. In this paper, the evolutionary process of a gravel surface is investigated based on a large‐scale recirculating flume experiment. The performances of the three equilibrium criteria are evaluated first, and then the evolution of the bed morphology is studied. The key findings include the following: (1) the sediment transport rate, slopes of water and bed surfaces, and one‐dimensional morphological features reach equilibrium roughly simultaneously; (2) two‐dimensional morphology continually evolves after these characteristics reach equilibrium, which is confirmed by the characteristics of the sediment transport process; and (3) the results from a numerical simulation suggest that a much longer time is required to reach an equilibrium state. Our results suggest that sufficient experimental time is required to investigate the equilibrium morphological characteristics of gravel surfaces, which is much longer than the equilibrium time reflected by the one‐dimensional equilibrium criteria. Copyright © 2017 John Wiley & Sons, Ltd.     

Observations of bedload transport in a gravel bed river during high flow using fiber‐optic DTS methods

The question: ‘how does a streambed change over a minor flood?’ does not have a clear answer due to lack of measurement methods during high flows. We investigate bedload transport and disentrainment during a 1.5‐year flood by linking field measurements using fiber optic distributed temperature sensing (DTS) cable with sediment transport theory and an existing explicit analytical solution to predict depth of sediment deposition from amplitude and phase changes of the diurnal near‐bed pore‐water temperature. The method facilitates the study of gravel transport by using near‐bed temperature time series to estimate rates of sediment deposition continuously over the duration of a high flow event coinciding with bar formation. The observations indicate that all gravel and cobble particles present were transported along the riffle at a relatively low Shields Number for the median particle size, and were re‐deposited on the lee side of the bar at rates that varied over time during a constant flow. Approximately 1–6% of the bed was predicted to be mobile during the 1.5‐year flood, indicating that large inactive regions of the bed, particularly between riffles, persist between years despite field observations of narrow zones of local transport and bar growth on the order ~3–5 times the median particle size. In contrast, during a seven‐year flood approximately 8–55% of the bed was predicted to become mobile, indicating that the continuous along‐stream mobility required to mobilize coarse gravel through long pools and downstream to the next riffle is infrequent. Copyright © 2017 John Wiley & Sons, Ltd.     

Form and growth of bars in a wandering gravel‐bed river

The changing form of developing alluvial river bars has rarely been studied in the field, especially in the context of the fixed, compound, mainly alternate gravel bars that are the major morphological feature of the wandering style. Century scale patterns of three‐dimensional growth and development, and the consequent scaling relations of such bars, are examined along the gravel‐bed reach of lower Fraser River, British Columbia, Canada. A retrospective view based on maps and aerial photographs obtained through the twentieth century shows that individual bars have a life history of about 100 years, except in certain, protected positions. A newly formed gravel bar quickly assumes its ultimate thickness and relatively quickly approaches its equilibrium length. Growth continues mainly by lateral accretion of unit bars, consistent with the lateral style of instability of the river. Bar growth is therefore allometric. Mature bars approach equilibrium dimensions and volume that scale with the overall size of the channel. Accordingly, the bars conform with several published criteria for the ultimate dimensions of alternate barforms. Sand bars, observed farther downstream, have notably different morphology. Fraser River presents a typical wandering channel planform, exhibiting elements of both meandered and low‐order braided channels. Hydraulic criteria to which the Fraser bars conform illustrate why this planform develops and persists. The modest rate of bed material transfer along the channel – typical of the wandering type – determines a century‐length time scale for bar development. This time scale is consistent with estimates that have been made for change of the macroform elements that determine the overall geometry of alluvial channels. Copyright © 2009 John Wiley & Sons, Ltd.     

Network‐scale dynamics of grain‐size sorting: implications for downstream fining,stream‐profile concavity,and drainage basin morphology

We explore the link between channel‐bed texture and river basin concavity in equilibrium catchments using a numerical landscape evolution model. Theory from homogeneous sediment transport predicts that river basin concavity directly increases with bed sediment size. If the effective grain size on a river bed governs its concavity, then natural phenomena such as grain‐size sorting and channel armouring should be linked to concavity. We examine this hypothesis by allowing the bed sediment texture to evolve in a transport‐limited regime using a two grain‐size mixture of sand and gravel. Downstream ?ning through selective particle erosion is produced in equilibrium. As the channel‐bed texture adjusts downstream so does the local slope. Our model predicts that it is not the texture of the original sediment mixture that governs basin concavity. Rather, concavity is linked to the texture of the sorted surface layer. Two different textural regimes are produced in the experiments: a transitional regime where the mobility of sand and gravel changes with channel‐bed texture, and a sand‐dominated region where the mobility of sand and gravel is constant. The concavity of these regions varies depending on the median gravel‐ or sand‐grain size, erosion rate, and precipitation rate. The results highlight the importance of adjustments in both surface texture and slope in natural rivers in response to changes in ?uvial and sediment inputs throughout a drainage network. This adjustment can only be captured numerically using multiple grain sizes or empirical downstream ?ning rules. Copyright © 2004 John Wiley & Sons, Ltd.     

Suspended sediment transport and deposition over a dune: Río Paraná, Argentina

Field data from the Rio Paraná, Argentina, are used to examine patterns of suspended sediment transport over a sand dune. Measurements of three‐dimensional velocity are made with an acoustic Doppler current profiler whilst suspended sediment concentration and particle size have been quantified using a laser in situ sediment scattering transmissometer. Suspended sediment concentration and streamwise and vertical sediment flux are highest close to the bed, with an upward vertical flux over the stoss side of the dune and downward flux over the lee side. Suspended sediment concentrations are higher over the crest compared with the trough and suspended sediment is coarsest near the bed. About 17% of the suspended‐load transported over the crest is deposited in the lee side before it reaches the trough. Most of this deposited sand is coarser sediment that originates close to the bed over the crest, a result consistent with simulations based on the model of Mohrig and Smith (Water Resources Research 1996; 32: 3207–3217) for the excursion lengths of sediment dispersed in the lee side of a dune. Copyright © 2009 John Wiley & Sons, Ltd.     

Hydraulic Tomography: Continuity and Discontinuity of High‐K and Low‐K Zones

Hydraulic tomography is an emerging field and modeling method that provides a continuous hydraulic conductivity (K) distribution for an investigated region. Characterization approaches that rely on interpolation between one‐dimensional (1D) profiles have limited ability to accurately identify high‐K channels, juxtapositions of lenses with high K contrast, and breaches in layers or channels between such profiles. However, locating these features is especially important for groundwater flow and transport modeling, and for design and operation of in situ remediation in complex hydrogeologic environments. We use transient hydraulic tomography to estimate 3D K in a volume of 15‐m diameter by 20‐m saturated thickness in a highly heterogeneous unconfined alluvial (clay to sand‐and‐gravel) aquifer with a K range of approximately seven orders of magnitude at an active industrial site in Assemini, Sardinia, Italy. A modified Levenberg‐Marquardt algorithm was used for geostatistical inversion to deal with the nonlinear nature of the highly heterogeneous system. The imaging results are validated with pumping tests not used in the tomographic inversion. These tests were conducted from three of five clusters of continuous multichannel tubing (CMTs) installed for observation in the tomographic testing. Locations of high‐K continuity and discontinuity, juxtaposition of very high‐K and very low‐K lenses, and low‐K “plugs” are evident in regions of the investigated volume where they likely would not have been identified with interpolation from 1D profiles at the positions of the pumping well and five CMT clusters. Quality assessment methods identified a suspect high‐K feature between the tested volume and a lateral boundary of the model.     

Secondary flow deflection in the lee of transverse dunes with implications for dune morphodynamics and migration

Measurements of lee‐side airflow response from an extensive array of meteorological instruments combined with smoke and flow streamer visualization is used to examine the development and morphodynamic significance of the lee‐side separation vortex over closely spaced transverse dune ridges. A differential deflection mechanism is presented that explains the three‐dimensional pattern of lee‐side airflow structure for a variety of incident flow angles. These flow patterns produce reversed, along‐dune, and deflected surface flow vectors in the lee that are inferred to result in net ‘lateral diversion’ of sand transport over one dune wavelength for incident angles as small as 10° from crest‐transverse (i.e. 80° from the crest line). This lateral displacement increases markedly with incident flow angle when expressed as the absolute value of the total deflection in degrees. Reversed and multi‐directional flow occurs for incident angles between 90° and 50°. These results document the three‐dimensional nature of flow and sand transport over transverse dunes and provide empirical evidence for an oblique migration model. Copyright © 2013 John Wiley & Sons, Ltd.     

Geomorphic response to gravel augmentation and high‐flow dam release in the Trinity River,California

The geomorphic effect of introducing a gravel augmentation totaling 520 m³ into a gravel‐bed stream during a dam‐controlled flood in May of 2015 was monitored with bedload transport measurements, an array of seismometers, and repeated topographic surveys. Half of the augmented gravel was injected into the flow with front‐end loaders on the rising limb of the flood and the other half was injected on the first day of the peak. Virtually all of the gravel transported past the injection point was deposited within about 7 to 10 channel widths of the injection point. Most of the injected gravel deposited along the left bank of the river whereas the right half of the channel bed was dominated by scour. The downstream third of the depositional area consisted of a small dune field that developed prior to the second gravel injection and subsequently migrated about one channel width downstream. A second depositional front was observed upstream from the gravel injection point, where a delta‐like wedge of bed material developed in the first hours of the flow release and changed little over the remainder of the release. These two depositional areas represent small‐scale bed‐material storage reservoirs with the potential to accumulate and periodically release packets of bed material. Interactions with such storage reservoirs are hypothesized to cause large bed‐material pulses to disperse by fragmenting into multiple smaller pulses. As a refinement to the conceptual model that views sediment pulse evolution in terms of dispersion and translation, the concept of pulse fragmentation has practical implications for gravel management. It implies that gravel augmentations can produce morphologic changes at locations that are separated from the augmentation point by arbitrarily long reaches, and it highlights the dependence of pulse propagation rates on the nature and distribution of the bed‐material storage reservoirs in the channel system. Copyright © 2017 John Wiley & Sons, Ltd.     

Field application of close‐range digital photogrammetry (CRDP) for grain‐scale fluvial morphology studies

In situ measurement of grain‐scale fluvial morphology is important for studies on grain roughness, sediment transport and the interactions between animals and the geomorphology, topics relevant to many river practitioners. Close‐range digital photogrammetry (CRDP) and terrestrial laser scanning (TLS) are the two most common techniques to obtain high‐resolution digital elevation models (DEMs) from fluvial surfaces. However, field application of topography remote sensing at the grain scale is presently hindered mainly by the tedious workflow challenges that one needs to overcome to obtain high‐accuracy elevation data. A recommended approach for CRDP to collect high‐resolution and high‐accuracy DEMs has been developed for gravel‐bed flume studies. The present paper investigates the deployment of the laboratory technique on three exposed gravel bars in a natural river environment. In contrast to other approaches, having the calibration carried out in the laboratory removes the need for independently surveyed ground‐control targets, and makes for an efficient and effective data collection in the field. Optimization of the gravel‐bed imagery helps DEM collection, without being impacted by variable lighting conditions. The benefit of a light‐weight three‐dimensional printed gravel‐bed model for DEM quality assessment is shown, and confirms the reliability of grain roughness data measured with CRDP. Imagery and DEM analysis evidences sedimentological contrasts between gravel bars within the reach. The analysis of the surface elevations shows the effect variable grain‐size and sediment sorting have on the surface roughness. By plotting the two‐dimensional structure functions and surface slopes and aspects we identify different grain arrangements and surface structures. The calculation of the inclination index allows determining the surface‐forming flow direction(s). We show that progress in topography remote sensing is important to extend our knowledge on fluvial morphology processes at the grain scale, and how a technique customized for use by fluvial geomorphologists in the field benefits this progress. Copyright © 2016 John Wiley & Sons, Ltd.     

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Examining the physical components of boundary shear stress for water‐worked gravel deposits

It is argued in this commentary that, in order to understand better the physical mechanisms that generate boundary shear stress over water‐worked gravel beds, flow velocity data should be re‐evaluated by spatial averaging the Reynolds equations to produce time‐ and space‐averaged (double‐averaged) momentum equations. A series of laboratory experiments were conducted in which the flow velocities were measured using a PIV system over two water‐worked gravel deposits. Combined with detailed data on the bed surface topography and vertical porosity, the physical components of shear stress were obtained. This enabled the various momentum transfer mechanisms present above, within and at the interface of a porous, fluvial deposit, to be quantified. This included the examination of the relevant contributions of temporal and spatial fluctuations in velocity and surface drag to the overall momentum transfer. It is demonstrated that double‐averaging represents a logical framework for assessing the fluid forces responsible for sediment entrainment and for investigating intragravel flow and sediment–water interface exchange mechanisms within the roughness layer in water‐worked gravel deposits. By considering the physical components of shear stress and their relative sizes it was possible to provide a physically based explanation for existing observations of enhanced mobility of gravel–sand mixtures and the transfer of solutes into porous, gravel deposits. This analysis reveals the importance of obtaining co‐located, high quality spatial data on the flow field and bed surface topography in order to gain a physical understanding of the mechanisms which generate boundary shear stress. Copyright © 2010 John Wiley & Sons, Ltd.     

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

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