Time-averaged flow structure in the central region of a stream confluence: a discussion

This paper is a discussion of Rhoads and Kenworthy (1998) ‘Time-averaged flow structure in the central region of a stream confluence’ Earth Surface Processes and Landforms, 23 , 171–191, that focuses upon the methods used to identify secondary circulation in river channel confluences. It argues that the Rozovskii method that Rhoads and Kenworthy use to rotate their field data to allow identification of secondary circulation cells is flawed, and can result in misleading conclusions about the nature of flow processes in confluences. It recommends that there is a re-emphasis upon helical as opposed to secondary circulation, and that recent developments in both field monitoring and numerical modelling may help significantly in this respect. Copyright © 1999 John Wiley & Sons, Ltd.     

Experimental study on the influence of river flow confluences on the open channel stage–discharge relationship

ABSTRACT

Accurate assessment of stage–discharge relationships in open channel flows is important to the design and management of hydraulic structures and engineering. Flow junctions commonly occur at the confluence of natural rivers or streams. The effect of flow junctions on the stage–discharge relationship at mountain river confluences was found by measuring velocity fields and water levels in experimental models. The results show that the backwater and accumulation–separation at flow junctions affect the flow structures and patterns in the channel; also, flow confluences may induce complex flow characteristics of backwater and flow separation at river junctions, indicating potential submerged flooding disasters within the confluence zone. The impacts of flow junctions on the stage–discharge relationship are investigated for two physical confluence models built from river confluence prototype systems in southwest China. The results show that the presence of tributary river inflows tends to increase the water level of the main river. This is important for flood control, flood-risk evaluation and engineering (e.g. hydropower station construction) in mountain rivers. Finally, a comparative quantitative analysis based on flow motion equations is conducted to evaluate the stage–discharge relationship in both uniform and regular confluence systems. The results indicate that more accurate prediction can be made when taking into account the flow non-uniformity induced by flow separation, backwater and distorted bed in the junction region.     

Assessment of flow paths and confluences for saltwater intrusion in a deltaic river network

Saltwater intrusion is a serious issue in estuarine deltas all over the world due to rapid urban sprawl and water shortage. Therefore, detecting the major flow paths or locations at risk of saltwater intrusion in estuarine ecosystems is important for mitigating saltwater intrusion. In this paper, we introduce a centrality index, the betweenness centrality (BC), to address this problem. Using the BC as the weighted attribute of the river network, we identify the critical confluences for saltwater intrusion and detect the preferential flow paths for saltwater intrusion through the least‐cost‐path algorithm from a graph theory approach. Moreover, we analyse the responses of the BC values of confluences calculated in the river network to salinity. Our results show that the major flow paths and critical confluences for saltwater intrusion in a deltaic river network can be represented by the least cost paths and the BC values of confluences, respectively. In addition, a significant positive correlation between the BC values of confluences and salinity is determined in the Pearl River Delta. Changes in the salinity can produce significant variation in the BC values of confluences. Therefore, freshwater can be diverted into these major flow paths and critical confluences to improve river network management under saltwater intrusion. Copyright © 2015 John Wiley & Sons, Ltd.     

Swimming behavior of juvenile silver carp near the separation zone of a channel confluence

Knowledge of locomotion of fish near river confluences is important for prediction of fish distribution in a river network.The flow separation zone near the confluence of a river network is a favorite habitat and feeding place for silver carp,which is one of the four major species of Chinese carp and usually provides positive rheotaxis to water flow.In the current study,a series of laboratory experiments were done to determine the behavioral responses of juvenile silver carp to the hydrodynamic ...     

Linking interannual river flow river variability across New Zealand to the Southern Annular Mode, 1979–2011

River flow constitutes an important element of the terrestrial branch of the hydrological cycle, yet knowledge regarding the extent to which its variability, at a range of timescales, is linked to a number of modes of atmospheric circulation is meagre. This is especially so in the Southern Hemisphere where strong candidates, such as El Niño Southern Oscillation and the Southern Annular Mode (SAM), for influencing climate and thus river flow variability can be found. This paper presents the results of an analysis of the impact of the SAM on winter and summer river flow variability across New Zealand, purposefully controlling for the influence of El Niño Southern Oscillation and the tendency for the SAM to adopt a positive phase over the last 10–20 years. Study results, based on identifying hydrological regions and applying circulation‐to‐environment and environment‐to‐circulation approaches commonly used in synoptic climatology, reveal a seasonal asymmetry of the response of river flow variability to the SAM; winter flows demonstrate a higher degree of statistical association with the SAM compared to summer flows. Further, because of the complex orography of New Zealand and its general disposition normal to zonal flows of moisture bearing winds, there are intraseasonal spatial variations in river flow SAM associations with clear rain shadow effects playing out in resultant river flow volumes. The complexity of SAM river flow associations found in this study warns against using indices of large scale modes of atmospheric circulation as blunt tools for hydroclimatological prediction at scales beyond hydroclimatological regions or areas with internal hydrological consistency.     

Three-dimensional turbulent structures at a medium-sized confluence with and without an ice cover

River confluences are characterized by a complex mixing zone with three-dimensional (3D) turbulent structures which have been described as both streamwise-oriented structures and Kelvin–Helmholtz (KH) vertical-oriented structures. The latter are visible where there is a turbidity difference between the two tributaries, whereas the former are usually derived from mean velocity measurements or numerical simulations. Few field studies recorded turbulent velocity fluctuations at high frequency to investigate these structures, particularly at medium-sized confluences where logistical constraints make it difficult to use devices such as acoustic doppler velocimeter (ADV). This study uses the ice cover present at the confluence of the Mitis and Neigette Rivers in Quebec (Canada) to obtain long-duration, fixed measurements along the mixing zone. The confluence is also characterized by a marked turbidity difference which allows to investigate the mixing zone dynamics from drone imagery during ice-free conditions. The aim of the study is to characterize and compare the flow structure in the mixing zone at a medium-sized (~40 m) river confluence with and without an ice cover. Detailed 3D turbulent velocity measurements were taken under the ice along the mixing plane with an ADV through eight holes at around 20 positions on the vertical. For ice-free conditions, drone imagery results indicate that large (KH) coherent structures are present, occupying up to 50% of the width of the parent channel. During winter, the ice cover affects velocity profiles by moving the highest velocities towards the centre of the profiles. Large turbulent structures are visible in both the streamwise and lateral velocity components. The strong correlation between these velocity components indicates that KH vortices are the dominating coherent structures in the mixing zone. A spatio-temporal conceptual model is presented to illustrate the main differences on the 3D flow structure at the river confluence with and without the ice cover. © 2019 John Wiley & Sons, Ltd.     

Quantifying the Morphological Print of Bedload Transport

Bedload and river morphology interact in a strong feedback manner. Bedload conditions the development of river morphology along different space and time scales; however, by concentrating the flow in preferential paths, a given morphology controls bedload for a given discharge. As bedload is a non‐linear response of shear stress, local morphology is likely to have a strong impact on bedload prediction when the shear stress is averaged over the section, as is usually done. This was investigated by comparing bedload measured in different bed morphologies (step‐pool, plane bed, riffle‐pool, braiding, and sand beds), with bedload measured in narrow flumes in the absence of any bed form, used here as a reference. The initial methodology consisted of fitting a bedload equation to the flume data. Secondly, the morphological signature of each river was studied as the distance to this referent equation. It was concluded that each morphology affects bedload in a different way. For a given average grain shear stress, the larger the river, the larger the deviation from the flume transport. Narrow streams are those morphologies that behave more like flumes; this is particularly true with flat beds, whereas results deviate from flumes to a greater extent in step‐pools. The riffle‐pool's morphology impacts bedload at different levels depending on the degree of bar development, considered here through the ratio D₈₄/D₅₀ which is used as a proxy for the local bed patchiness and morphology. In braiding rivers morphological effects are important but difficult to assess because width is dependent on transport rate. Bed morphology was found to have negligible effects in sand bed rivers where the Shields stress is usually sufficiently high to minimize the non‐linearity effects when hydraulics is averaged over the section. Copyright © 2015 John Wiley & Sons, Ltd.     

A two‐dimensional morphodynamic model of gravel‐bed river with floodplain vegetation

Artificially straight river channels tend to be unstable, and ultimately develop into river meanders through bank erosion and point‐bar deposition. In this paper account is taken of the effects of riparian and floodplain vegetation on bank strength, floodplain flow resistance, shear stress partitioning, and bedload transport. This is incorporated into an existing 2D hydrodynamic‐morphological model. By applying the new model to an initially straight and single‐threaded channel, the way that its planform and cross‐sectional geometry evolve for different hydraulic and floodplain vegetation conditions is demonstrated. The results show the formation and upstream migration of gravel bars, confluence scouring and the development of meandering and braiding channel patterns. In cases where the channel becomes unstable, the instability grows out of bar formation. The resulting braiding patterns are similar to analytical results. The formation of a transition configuration requires a strong influence from vegetation. Copyright © 2010 John Wiley & Sons, Ltd.     

Extreme flood‐driven fluvial bank erosion and sediment loads: direct process measurements using integrated Mobile Laser Scanning (MLS) and hydro‐acoustic techniques

This methods paper details the first attempt at monitoring bank erosion, flow and suspended sediment at a site during flooding on the Mekong River induced by the passage of tropical cyclones. We deployed integrated mobile laser scanning (MLS) and multibeam echo sounding (MBES), alongside acoustic Doppler current profiling (aDcp), to directly measure changes in river bank and bed at high (~0.05 m) spatial resolution, in conjunction with measurements of flow and suspended sediment dynamics. We outline the methodological steps used to collect and process this complex point cloud data, and detail the procedures used to process and calibrate the aDcp flow and sediment flux data. A comparison with conventional remote sensing methods of estimating bank erosion, using aerial images and Landsat imagery, reveals that traditional techniques are error prone at the high temporal resolutions required to quantify the patterns and volumes of bank erosion induced by the passage of individual flood events. Our analysis reveals the importance of cyclone‐driven flood events in causing high rates of erosion and suspended sediment transport, with a c. twofold increase in bank erosion volumes and a fourfold increase in suspended sediment volumes in the cyclone‐affected wet season. Copyright © 2016 John Wiley & Sons, Ltd.     

Hydro‐morphodynamic evolution in a 90° movable bed discordant confluence with low discharge ratio

Confluences with low discharge and momentum ratios, where narrow steep tributaries with high sediment load join a wide low‐gradient main channel that provides the main discharge, are often observed in high mountain regions such as in the upper‐Rhone river catchment in Switzerland. Few existing studies have examined the hydro‐morphodynamics of this type of river confluence while considering sediment discharge in both confluent channels. This paper presents the evolution of the bed morphology and hydrodynamics as observed in an experimental facility with a movable bed. For that purpose, one experiment was carried out in a laboratory confluence with low discharge and momentum ratios, where constant sediment rates were supplied to both flumes. During the experiment, bed topography and water surface elevations were systematically recorded. When the bed topography reached a steady state (so‐called equilibrium) and the outgoing sediment rate approximated the incoming rate, flow velocity was measured at 12 different points distributed throughout the confluence, and the grain size distribution of the bed surface was analyzed. Typical morphodynamic features of discordant confluences such as a bank‐attached bar and a flow deflection zone are identified in this study. Nevertheless, the presence of a marked scour hole in the discordant confluence and distinct flow regimes for the tributary and main channel, differ from results obtained in previous studies. Strong acceleration of the flow along the outer bank of the main channel is responsible for the scour hole. This erosion is facilitated by the sediment discharge into the confluence from the main channel which inhibits bed armoring in this region. The supercritical flow regime observed in the tributary is the hydrodynamic response to the imposed sediment rate in the tributary. Copyright © 2015 John Wiley & Sons, Ltd.     

Secondary flow in anabranch confluences of a braided,gravel-bed stream

Measurements of the primary and secondary velocity components were máde in two, active, braided river anabranch confluences with a simple Y-shaped plan form, in the gravelly Sunwapta River (D₅₀ of approximately 30 mm). Flow velocity was measured at regularly-spaced intervals using a bidirectional electromagnetic current meter and the measured downstream and cross-stream velocities were converted to primary and secondary velocities to yield the secondary circulation. The primary (downstream) velocity field shows two high velocity streams from the two tributaries which merge (and, in some cases, accelerate) into a single high velocity core over the thalweg. Primary flow velocity declines as the flow expands and diverges at the downstream end of the confluence. The secondary circulation is dominated by two helical cells, back-to-back, plunging over the thalweg and diverging at the bed. This is the first confirmation of this flow structure in confluences, based on field measurements. The strength of the secondary cells declines downstream through each confluence, and laterally away from the thalweg area in cross-section. There is also a tendency for one cell, from the larger of the tributaries, to override the other. The secondary and primary flow structure and strength differs slightly between the two confluences and this is reflected in differences in scour hole form.     

Hydraulic and geomorphic processes in an overbank flood along a meandering,gravel‐bed river: implications for chute formation

Hydraulic interactions between rivers and floodplains produce off‐channel chutes, the presence of which influences the routing of water and sediment and thus the planform evolution of meandering rivers. Detailed studies of the hydrologic exchanges between channels and floodplains are usually conducted in laboratory facilities, and studies documenting chute development are generally limited to qualitative observations. In this study, we use a reconstructed, gravel‐bedded, meandering river as a field laboratory for studying these mechanisms at a realistic scale. Using an integrated field and modeling approach, we quantified the flow exchanges between the river channel and its floodplain during an overbank flood, and identified locations where flow had the capacity to erode floodplain chutes. Hydraulic measurements and modeling indicated high rates of flow exchange between the channel and floodplain, with flow rapidly decelerating as water was decanted from the channel onto the floodplain due to the frictional drag provided by substrate and vegetation. Peak shear stresses were greatest downstream of the maxima in bend curvature, along the concave bank, where terrestrial LiDAR scans indicate initial floodplain chute formation. A second chute has developed across the convex bank of a meander bend, in a location where sediment accretion, point bar development and plant colonization have created divergent flow paths between the main channel and floodplain. In both cases, the off‐channel chutes are evolving slowly during infrequent floods due to the coarse nature of the floodplain, though rapid chute formation would be more likely in finer‐grained floodplains. The controls on chute formation at these locations include the flood magnitude, river curvature, floodplain gradient, erodibility of the floodplain sediment, and the flow resistance provided by riparian vegetation. Copyright © 2015 John Wiley & Sons, Ltd.     

Morphological signatures of normal faulting in low‐gradient alluvial rivers in south‐eastern Louisiana,USA

We explore the fluvial response to faulting in three low‐gradient, sand‐bed rivers in south‐eastern Louisiana, USA, that flow across active normal faults from footwall (upstream) to hangingwall (downstream). We calculate sinuosity, migration rate and migration direction in order to identify anomalies spatially associated with fault scarps. In two of the rivers we model one‐dimensional steady water flow to identify anomalies in surface water slope, width‐to‐depth ratio, and shear stress. In each of these rivers there is one location where flow modeling suggests potential channel incision through the footwall, as indicated by relatively high surface water slopes and shear stress values. In one of these footwall locations, the river straightens and width‐to‐depth ratios decrease, likely contributing to higher surface water slopes and shear stress. This is in contrast to previous studies that have proposed increased sinuosity across fault footwalls and decreased sinuosity across hangingwalls. However, in two hangingwall locations we also observe relatively less sinuous channels. Other planform changes on the hangingwall include topographic steering of channels along and towards the fault and one example of an avulsion. The most notable anomaly in migration rate occurs on the hangingwall of a fault where a river has cut off a meander loop. Although fluvial response to faulting varies here, comparatively large and small channels exhibit similar responses. Further, Pleistocene fault slip rates are orders of magnitude lower than the channel migration rates, suggesting that faulting should not be a major influence on the fluvial evolution. Nonetheless, notable channel anomalies exist near faults, suggesting that recent fault slip rates are higher than Pleistocene rates, and/or that low‐gradient alluvial channels are more sensitive to faulting than previous studies have suggested. Rivers appear to be influenced by faulting in this setting, however background rates of meander loop cutoff may be just as influential as faulting. Copyright © 2015 John Wiley & Sons, Ltd.     

Network concepts to describe channel importance and change in multichannel systems: test results for the Jamuna River,Bangladesh

Most of the largest rivers on Earth have multiple active channels connected at bifurcations and confluences. At present a method to describe a channel network pattern and changes in the network beyond the simplistic braiding index is unavailable. Our objectives are to test a network approach to understand the character, stability and evolution of a multi‐channel river pattern under natural discharge conditions. We developed a semi‐automatic method to derive a chain‐like directional network from images that represent the multi‐channel river and to connect individual network elements through time. The Jamuna River was taken as an example with a series of Landsat TM and ETM+ images taken at irregular intervals between 1999 and 2004. We quantified the overall importance of individual channels in the entire network using a centrality property. Centrality showed that three reaches can be distinguished along the Jamuna with a different network character: the middle reach has dominantly one important channel, while upstream and downstream there are about two important channels. Temporally, relatively few channels changed dramatically in both low‐flow and high‐flow periods despite the increase of braiding index during a flood. Based on the centrality we calculated a weighted braiding index that represents the number of important channels in the network, which is about two in the Jamuna River and which is larger immediately after floods. We conclude that the network measure centrality provides a novel characterization of river channel networks, highlighting properties and tendencies that have morphological significance. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of pit removal methods on river network position

Digital elevation models often contain depressions that result in areas described as having no drainage, referred to as sinks or pits. These depressions disrupt the drainage surface, which disrupt routing of flow over the surface. Most of the attributes that can be extracted from a digital elevation model rely on flow‐routing algorithms to calculate the upslope contributing area. There is little information on the influence of the various algorithms on the position and on the connectivity of the extracted networks. The aim of this study was to assess the effects of pit removal methods, data sources and flow‐routing algorithms on the position of river networks. The results show that all factors and methods have an impact on the position of the extracted networks. The pit removal method combining filling and carving extracted river networks closer to the reference, as well the elevation models with higher resolution. Single‐flow direction methods provided more accurate positioning of river network, in this test area where the drainage is generally well defined. Copyright © 2011 John Wiley & Sons, Ltd.     

Observations of bedload movement,bar development and sediment supply in the braided Ohau river

Relative bedload transport rate and hydraulic parameters were measured on two occasions in a reach of the braided, gravel-bedded Ohau River. Each reach contained a deep, fast-flowing chute leading to an area of diverging, shallow flow which contained a middle bar. The measurements are self-consistent, and indicate that, where flow is concentrated in a deep chute, shear stress is high, but where flow diverges, depth, slope and shear stress decrease. In the first survey, the bed was scoured in the chute and sediment transport rates were high, but where flow diverged sediment transport rate decreased. It appears that deposition leads to bar growth, bar growth in turn reduces slope and depth, deposition is encouraged, and the bar continues to grow vertically, laterally and upstream. In the second survey no sediment transport was observed, despite hydraulic conditions very similar to the first survey. The absence of sediment transport is attributed to the cessation of sediment supply to the river channel.     

On the use of horizontal acoustic Doppler profilers for continuous bed shear stress monitoring

Continuous monitoring of bed shear stress in large river systems may serve to better estimate alluvial sediment transport to the coastal ocean.Here we explore the possibility of using a horizontally deployed acoustic Doppler current profiler(ADCP) to monitor bed shear stress,applying a prescribed boundary layer model,previously used for discharge estimation.The model parameters include the local roughness length and a dip correction factor to account for sidewall effects.Both these parameters depend on river stage and on the position in the cross-section, and were estimated from shipborne ADCP data.We applied the calibrated boundary layer model to obtain bed shear stress estimates over the measuring range of the HADCP.To validate the results,co-located coupled ADCPs were used to infer bed shear stress,both from Reynolds stress profiles and from mean velocity profiles. From HADCP data collected over a period of 1.5 years,a time series of width profiles of bed shear stress was obtained for a tidal reach of the Mahakam River,East Kalimantan,Indonesia.A smaller dataset covering 25 hours was used for comparison with results from the coupled ADCPs.The bed shear stress estimates derived from Reynolds stress profiles appeared to be strongly affected by local effects causing upflow and downflow,which are not included in the boundary layer model used to derive bed shear stress with the horizontal ADCP.Bed shear stresses from the coupled ADCP are representative of a much more localized flow,while those derived with the horizontal ADCP resemble the net effect of the flow over larger scales.Bed shear stresses obtained from mean velocity profiles from the coupled ADCPs show a good agreement between the two methods,and highlight the robustness of the method to uncertainty in the estimates of the roughness length.     

Influence of junction angle on three‐dimensional flow structure and bed morphology at confluent meander bends during different hydrological conditions

Recent field and modeling investigations have examined the fluvial dynamics of confluent meander bends where a straight tributary channel enters a meandering river at the apex of a bend with a 90° junction angle. Past work on confluences with asymmetrical and symmetrical planforms has shown that the angle of tributary entry has a strong influence on mutual deflection of confluent flows and the spatial extent of confluence hydrodynamic and morphodynamic features. This paper examines three‐dimensional flow structure and bed morphology for incoming flows with high and low momentum‐flux ratios at two large, natural confluent meander bends that have different tributary entry angles. At the high‐angle (90°) confluent meander bend, mutual deflection of converging flows abruptly turns fluid from the lateral tributary into the downstream channel and flow in the main river is deflected away from the outer bank of the bend by a bar that extends downstream of the junction corner along the inner bank of the tributary. Two counter‐rotating helical cells inherited from upstream flow curvature flank the mixing interface, which overlies a central pool. A large influx of sediment to the confluence from a meander cutoff immediately upstream has produced substantial morphologic change during large, tributary‐dominant discharge events, resulting in displacement of the pool inward and substantial erosion of the point bar in the main channel. In contrast, flow deflection is less pronounced at the low‐angle (36°) confluent meander bend, where the converging flows are nearly parallel to one another upon entering the confluence. A large helical cell imparted from upstream flow curvature in the main river occupies most of the downstream channel for prevailing low momentum‐flux ratio conditions and a weak counter‐rotating cell forms during infrequent tributary‐dominant flow events. Bed morphology remains relatively stable and does not exhibit extensive scour that often occurs at confluences with concordant beds. Copyright © 2014 John Wiley & Sons, Ltd.     

Factors influencing the calculation of periodic secondary circulation in a tidal river: numerical modelling of the lower Sacramento River,USA

Secondary circulation is the component of three‐dimensional (3D) flow in river channels perpendicular to the primary flow direction. Secondary circulation calculated from acoustic Doppler current profiler (ADCP) transects is sensitive to the calculation method and is affected by the transect angle relative to the mean flow direction and variations in the flow direction along a transect. To quantify bounds on transect alignment relative to river flow for field data collection and examine tidal time‐scale variability in secondary circulation, the 3D hydrodynamic model UnTRIM was applied to simulate the hydrodynamics in the lower reach of the Sacramento River (CA, USA). Secondary circulation was calculated using the Rozovskii and the zero net discharge methods on repeated transects extracted from the model results in regions of both relatively uniform and complex flows. When the depth‐averaged flow direction along a transect varied by more than about 5 °, occurring when the transect was as little as 10 to 20 ° out of normal to the mean flow direction, the Rozovskii method produced more realistic secondary circulation than the zero net discharge method. Analysis indicated that ADCP transects should be within 20 ° of perpendicular to the mean flow direction when calculating secondary circulation. Secondary circulation strength around two tidally influenced bends generally increased with increasing flow and broke down near slack water. However, the strength of the secondary circulation was not only a function of the flow magnitude, but also depended on the direction of the water flow and the transect location relative to the river curvature, which varied with the tidal flow direction. Copyright © 2015 John Wiley & Sons, Ltd.