Bedload transport processes in a braided gravel-bed river model

A 1:50 scale hydraulic model was designed, based on Froude number similarity and using hydrological and sediment data from a small braided gravel-bed river (the North Branch of the Ashburton River, Canterbury, New Zealand). Eighteen experiments were conducted; seven using steady flows, and eleven using unsteady flows. The experiments were carried out in a 20 m × 3 m tilting flume equipped with a continuous sediment feed and an automated data acquisition and control system. In all experiments water at 30°C was used to reduce viscosity-related scale effects. Analyses of the experimental data revealed that bedload transport rates in braided channels are highly variable, with relative variability being inversely related to mean bedload transport rate. Variability was also found to be cyclic with short-term variations being caused by the migration of bedforms. Bedload transport was found to be more efficient under steady flow than under unsteady flow, and it was postulated that this is caused by a tendency for channel form to evolve towards a condition which maximizes bedload transport for the occurring flow. Average bedload transport rate was found to vary with channel form, although insufficient measurements were made to define a relationship.     

TEMPERATURE EFFECTS ON AIRGUN SIGNATURES1

Experiments in an 850 litre water tank were performed in order to study temperature effects on airgun signatures, and to achieve a better understanding of the physical processes that influence an airgun signature. The source was a bolt airgun with a chamber volume of 1.6 cu.in. The pressure used was 100 bar and the gun depth was 0.5 m. The water temperature in the tank was varied between 5°C and 45°C. Near-field signatures were recorded at different water temperatures. Typical signature characteristics such as the primary-to-bubble ratio and the bubble time period increased with increasing water temperature. For comparison and in order to check whether this is valid for larger guns, computer modelling of airguns with chamber volumes of 1.6 and 40 cu.in. was performed. In the modelling the same behaviour of the signatures with increasing water temperature can be observed. The increase in the primary-to-bubble ratio and the bubble time period with increasing water temperature can be explained by an increased mass transfer across the bubble wall.     

Physical problems of the benthic boundary layer

Since the boundary layer at the sea bed has a number of features in common with boundary layers found in laboratory scale flows and in meteology, a brief review is given first of the properties which may be inferred from experience in these fields or from theoretical studies. Measurements of velocity profiles, turbulence, and shearing stress which have been made near the bottom, in deep water, and on the continental shelf, are described in relation to this background. In particular, the logarithmic form of the velocity profile near the bed and deductions from it appear to be valid in certain conditions, but the occurrence of ripples and other bed forms is a complicating feature. The relation of the dynamical aspects of the flow to the transport of sediment as bed load and in suspension is discussed. The diflusive properties of the layer are then considered, in relation to fluxes near the sea-sediment interface and to the formation of nepheloid layers or layers well mixed in temperature and salinity.At Dept. of Oceanography, Univ. of Washington, Seattle, WA, 98195, U.S.A., Jan. to Aug. 1978.     

Variability of shallow overland flow velocity and soil aggregate transport observed with digital videography

Field experiments at Tiramoana station 30 km north of Christchurch, New Zealand using an erosion plot 16·5 m long, 0·6 m wide, and with a slope of 14–14·5° on rendzina soil aimed to measure the variability of flow velocity and of soil aggregates transport rate in shallow overland flow. Discharge/cross‐section area ratio was used to estimate mean velocity, and high‐speed digital video camera and image analysis provided information about flow and sediment transport variability. Six flow runs with 0·5–3·0 L s^?1 discharges were supercritical with Froude numbers close to or more than 1. Mean flow velocity followed Poiseuille law, float numbers were more than 1·5 and hydraulic resistance was an inverse proportional function of the Reynolds number, which is typical for laminar flows. Hence actual velocity varied through time significantly and the power spectrum was of ‘red‐noise’, which is typical for turbulent flow. Sediment transport rates had even higher variability, and soil aggregates transport was a compound Poisson process. Copyright © 2008 John Wiley & Sons, Ltd.     

Nature and origin of cone-forming volcanic breccias in the Te Herenga Formation, Ruapehu, New Zealand

 Volcanic breccias form large parts of composite volcanoes and are commonly viewed as containing pyroclastic fragments emplaced by pyroclastic processes or redistributed as laharic deposits. Field study of cone-forming breccias of the andesitic middle Pleistocene Te Herenga Formation on Ruapehu volcano, New Zealand, was complemented by paleomagnetic laboratory investigation permitting estimation of emplacement temperatures of constituent breccia clasts. The observations and data collected suggest that most breccias are autoclastic deposits. Five breccia types and subordinate, coherent lava-flow cores constitute nine, unconformity-bounded constructional units. Two types of breccia are gradational with lava-flow cores. Red breccias gradational with irregularly shaped lava-flow cores were emplaced at temperatures in excess of 580  °C and are interpreted as aa flow breccias. Clasts in gray breccia gradational with tabular lava-flow cores, and in some places forming down-slope-dipping avalanche bedding beneath flows, were emplaced at varying temperatures between 200 and 550  °C and are interpreted as forming part of block lava flows. Three textural types of breccia are found in less intimate association with lava-flow cores. Matrix-poor, well-sorted breccia can be traced upslope to lava-flow cores encased in autoclastic breccia. Unsorted boulder breccia comprises constructional units lacking significant exposed lava-flow cores. Clasts in both of these breccia types have paleomagnetic properties generally similar to those of the gray breccias gradational with lava-flow cores; they indicate reorientation after acquisition of some, or all, magnetization and ultimate emplacement over a range of temperatures between 100 and 550  °C. These breccias are interpreted as autoclastic breccias associated with block lava flows. Matrix-poor, well-sorted breccia formed by disintegration of lava flows on steep slopes and unsorted boulder breccia is interpreted to represent channel-floor and levee breccias for block lava flows that continued down slope. Less common, matrix-rich, stratified tuff breccias consisting of angular blocks, minor scoria, and a conspicuously well-sorted ash matrix were generally emplaced at ambient temperature, although some deposits contain clasts possibly emplaced at temperatures as high as 525  °C. These breccias are interpreted as debris-flow and sheetwash deposits with a dominant pyroclastic matrix and containing clasts likely of mixed autoclastic and pyroclastic origin. Pyroclastic deposits have limited preservation potential on the steep, proximal slopes of composite volcanoes. Likewise, these steep slopes are more likely sites of erosion and transport by channeled or unconfined runoff rather than depositional sites for reworked volcaniclastic debris. Autoclastic breccias need not be intimately associated with coherent lava flows in single outcrops, and fine matrix can be of autoclastic rather than pyroclastic origin. In these cases, and likely many other cases, the alternation of coherent lava flows and fragmental deposits defining composite volcanoes is better described as interlayered lava-flow cores and cogenetic autoclastic breccias, rather than as interlayered lava flows and pyroclastic beds. Reworked deposits are probably insignificant components of most proximal cone-forming sequences. Received: 1 October 1998 / Accepted: 28 December 1998     

Barchan ripples under unidirectional water flows in the laboratory: formation and planar morphology

The purpose of the present study is to investigate the formation and planar shape of barchan ripples generated by laboratory unidirectional water currents and to compare the morphology with that of barchans formed by laboratory waves and flows in natural environments. A thin veneer of sand as a sediment source was placed initially at the upstream part of a recirculating flume and the change in bed configurations by the flow was recorded by a video camera. Two types of formative process were observed: one was that barchan ripples grew from barchanoid bedforms and the other was that they developed from a small aggregate of sand particles. The barchan ripples in the present test had the same characteristics in the height–width relation and in the migration speed as previously reported from the desert environment and wave‐flume studies. An examination of the planar shape of the barchans led to the result that the barchan ripples had larger values of body‐length/width, compared with those of barchans formed by water waves in the laboratory and by airflow in natural deserts. The horn‐length/width of the barchan ripples was smaller than that of barchan dunes in deserts but larger than that of barchans in oscillatory laboratory flows. Copyright © 2005 John Wiley & Sons, Ltd.     

On the modelling of sand bedforms using the semivariogram

This study shows the usefulness of the semivariogram for modelling sand ripples created by water flows of varied flow intensity. A combination of two mathematical functions is fitted to each sample semivariogram, that is an exponential (or stochastic) component and a periodic component. The parameters of each of these components have direct physical meaning. A non-dimensional ratio combining the two parameters of the exponential model is interpreted as a regularity index (which increases with the degree of regularity of bedform arrangement). This regularity index is inversely related to the Froude number of the flow. The non-dimensional wavelength, estimated from the dominant periodic function, is also inversely and closely related to the Froude number. The wave height, accurately estimated from properties of the two fitted components, is a direct function of flow velocity and is also proportional to the standard deviation of bed elevations. The bedform shape introduces a considerable discrepancy between the generally assumed normal frequency distribution and the empirical distributions of bed height. The series of bed elevations are generally characterized by a mixture of normal distributions having the same variance but different means. The calculation of a covariance assuming a constant and single mean (as in spectral analysis) can therefore be misleading and the problem may be avoided by using the semivariogram.     

Jet flow over foredunes

Jet flows, which are localized flows exhibiting a high speed maxima, are relatively common in nature, and in many devices. They have only been occasionally observed on dunes, and their dynamics are poorly known. This paper examines computational fluid dynamic (CFD) two‐dimensional (2D) modelling of jet flow over a foredune topography. Flow was simulated in 10° increments from onshore (0°) to highly oblique alongshore (70°) incident wind approach angles. CFD modelling reveals that the formation of a jet is not dependent on a critical wind speed, and an increase in incident wind velocity does not affect the magnitude of jet flow. A jet is first formed at ~1.0 m seawards of the foredune crest on the Prince Edward Island foredune morphology example examined here. A jet is not developed when the incident wind is from an oblique approach angle greater than ~50° because there is significantly less flow acceleration across a much lower slope at this incident angle. The presence of a scarp does influence the structure of the crest jet, in that the jet is more pronounced where a scarp is present. Surface roughness affects the magnitude of jet expansion and jets are better developed on bare surfaces compared to vegetated ones. Copyright © 2016 John Wiley & Sons, Ltd.     

Effect of saltating sediment on flow resistance and bed roughness in overland flow

The acceleration of saltating grains by overland flow causes momentum to be transferred from the flow to the grains, thereby increasing flow resistance and bed roughness. To assess the impact of saltating sediment on overland flow hydraulics, velocity profiles in transitional and turbulent flows on a fixed sand-covered bed were measured using hot-film anemometry. Five discharges were studied. At each discharge, three flows were measured: one free of sediment, one with a relatively low sediment load, and one with a relatively high sediment load. In these flows from 83 to 90 per cent of the sediment was travelling by saltation. As a result, in the sediment-laden flows the near-bed velocities were smaller and the velocity profiles steeper than those in the equivalent sediment-free flows. Sediment loads ranged up to 87·0 per cent of transport capacity and accounted for as much as 20·8 per cent of flow resistance (measured by the friction factor) and 89·7 per cent of bed roughness (measured by the ratio of the roughness length to median grain diameter). It is concluded that saltating sediment has a considerable impact on overland flow hydraulics, at least on fixed granular beds. Saltation is likely to have a relatively smaller effect on overland flow on natural hillslopes and agricultural fields where form and wave resistance dominate. Still, saltation is generally of greater significance in overland flow than in river flow, and for this reason its effect on overland flow hydraulics is deserving of further study. © 1998 John Wiley & Sons, Ltd.     

Characterization of wood-laden flows in rivers

Inorganic sediment is not the only solid-fraction component of river flows; flows may also carry significant amounts of large organic material (i.e. large wood), but the characteristics of these wood-laden flows (WLFs) are not well understood yet. With the aim to shed light on these relatively unexamined phenomena, we collected home videos showing natural flows with wood as the main solid component. Analyses of these videos as well as the watersheds and streams where the videos were recorded allowed us to define for the first time WLFs, describe the main characteristics of these flows and broaden the definition of wood transport regimes (adding a new regime called here hypercongested wood transport). According to our results, WLFs may occur repeatedly, in a large range of catchment sizes, generally in steep, highly confined single thread channels in mountain areas. WLFs are typically highly unsteady and the log motion is non-uniform, as described for other inorganic sediment-laden flows (e.g. debris flows). The conceptual integration of wood into our understanding of flow phenomena is illustrated by a novel classification defining the transition from clear water to hypercongested, wood and sediment-laden flows, according to the composition of the mixture (sediment, wood, and water). We define the relevant metrics for the quantification and modelling of WLFs, including an exhaustive discussion of different modelling approaches (i.e. Voellmy, Bingham and Manning) and provide a first attempt to simulate WLFs. We draw attention to WLF phenomena to encourage further field, theoretical, and experimental investigations that may contribute to a better understanding of flows in river basins, leading to more accurate predictions, and better hazard mitigation and management strategies. © 2019 John Wiley & Sons, Ltd.     

Large wood inhibits debris flow runout in forested southeast Alaska

Due to their potentially long runout, debris flows are a major hazard and an important geomorphic process in mountainous environments. Understanding runout is therefore essential to minimize risk in the near-term and interpret the pace and pattern of debris flow erosion and deposition over geomorphic timescales. Many debris flows occur in forested landscapes where they mobilize large volumes of large woody debris (LWD) in addition to sediment, but few studies have quantitatively documented the effects of LWD on runout. Here, we analyze recent and historic debris flows in southeast Alaska, a mountainous, forested system with minimal human alteration. Sixteen debris flows near Sitka triggered on August 18, 2015 or more recently had volumes of 80 to 25 000 m³ and limited mobility compared to a global compilation of similarly-sized debris flows. Their deposits inundated 31% of the planimetric area, and their runout lengths were 48% of that predicted by the global dataset. Depositional slopes were 6°–26°, and mobility index, defined as the ratio of horizontal runout to vertical elevation change, ranged from 1.2 to 3, further indicating low mobility. In the broader southeast Alaskan region consisting of Chichagof and Baranof Islands, remote sensing-based analysis of 1061 historic debris flows showed that mobility index decreased from 2.3–2.5 to 1.4–1.8 as average forest age increased from 0 to 416 years. We therefore interpret that the presence of LWD within a debris flow and standing trees, stumps, and logs in the deposition zone inhibit runout, primarily through granular phenomena such as jamming due to force chains. Calibration of debris flow runout models should therefore incorporate the ecologic as well as geologic setting, and feedbacks between debris flows and vegetation likely control the transport of sediment and organic material through steep, forested catchments over geomorphic time. © 2020 John Wiley & Sons, Ltd.     

Climate impacts on river flow: projections for the Medway catchment,UK, with UKCP09 and CATCHMOD

The potential impact of climate change on areas of strategic importance for water resources remains a concern. Here, river flow projections for the River Medway, above Teston in southeast England are presented, which is just such an area of strategic importance. The river flow projections use climate inputs from the Hadley Centre Regional Climate Model (HadRM3) for the time period 1960–2080 (a subset of the early release UKCP09 projections). River flow predictions are calculated using CATCHMOD, the main river flow prediction tool of the Environment Agency (EA) of England and Wales. In order to use this tool in the best way for climate change predictions, model setup and performance are analysed using sensitivity and uncertainty analysis. The model's representation of hydrological processes is discussed and the direct percolation and first linear storage constant parameters are found to strongly affect model results in a complex way, with the former more important for low flows and the latter for high flows. The uncertainty in predictions resulting from the hydrological model parameters is demonstrated and the projections of river flow under future climate are analysed. A clear climate change impact signal is evident in the results with a persistent lowering of mean daily river flows for all months and for all projection time slices. Results indicate that a projection of lower flows under future climate is valid even taking into account the uncertainties considered in this modelling chain exercise. The model parameter uncertainty becomes more significant under future climate as the river flows become lower. This has significant implications for those making policy decisions based on such modelling results. Copyright © 2010 John Wiley & Sons, Ltd.     

Observations of bedforms on a dissipative macrotidal beach

Field measurements of wave ripples and megaripples were made with a Sand Ripple Profiler in the surf and shoaling zones of a sandy macrotidal dissipative beach at Perranporth, UK in depths 1–6 m and significant wave heights up to 2.2 m. A frequency domain partitioning approach allowed quantification of height (η), length (λ) and migration rate of ripples and megaripples. Wave ripples with heights up to 2 cm and wavelengths ～20 cm developed in low orbital velocity conditions (u _m?<?0.65 m/s) with mobility number ψ?<?25. Wave ripple heights decreased with increasing orbital velocity and were flattened when mean currents were >0.1 m/s. Wave ripples were superimposed on top of megaripples (η?=?10 cm, λ?=?1 m) and contributed up to 35 % of the total bed roughness. Large megaripples with heights up to 30 cm and lengths 1–1.8 m developed when the orbital velocity was 0.5–0.8 m/s, corresponding to mobility numbers 25–50. Megaripple heights and wavelengths increased with orbital velocity but reduced when mean current strengths were >0.15 m/s. Wave ripple and megaripple migrations were generally onshore directed in the shoaling and surf zones. Onshore ripple migration rates increased with onshore-directed (+ve) incident wave skewness. The onshore migration rate reduced as offshore-directed mean flows (undertow) increased in strength and reached zero when the offshore-directed mean flow was >0.15 m/s. The migration pattern was therefore linked to cross-shore position relative to the surf zone, controlled by competition between onshore-directed velocity skewness and offshore-directed mean flow.     

CORRECTION FACTORS IN THE DETERMINATION OF MEAN VELOCITY OF OVERLAND FLOW

The velocity of overland flow has been conventionally measured using tracers, but it is difficult to measure the mean flow velocity directly because the centroid of the tracer plume is not easily identified. Consequently, previous investigators have measured the velocity of the leading edge of the plume and multiplied it by a correction factor α to obtain an estimate of mean velocity. An alternative method is to measure the velocity of the peak concentration in the tracer plume and multiply this velocity by another correction factor β to estimate mean velocity. To investigate the controls of α and β and develop predictive models for these correction factors, 40 experiments were performed in a flume with a mobile sand bed. Multiple regression analyses reveal that both α and β vary inversely with slope and directly with Reynolds number. The derived regression equations may be used to calculate the mean velocity of other shallow overland flows, at least within the range of slope and Reynolds number for which the equations were developed. In the experiments, slope ranged from 2.7;° to 10° and Reynolds number from 1900 to 12 600.     

Comparing the transitional behaviour of kaolinite and bentonite suspension flows

Past research has demonstrated the dramatic effects that variations in suspended clay can have on the properties of flow by producing a range of transitional flows between turbulent and laminar states, depending on clay concentration and fluid shear. Past studies have been restricted to kaolinite flows, a clay mineral that has relatively weak cohesive properties. This paper extends these studies to suspension flows of bentonite, a clay mineral that attains higher viscosities at far lower volumetric concentrations within a flow. The results show that the types of transitional flow behaviour recognized in past studies can also be found in bentonite suspension flows, but at lower suspended sediment concentrations, thus demonstrating an even more dramatic effect on flow properties, and potentially on sediment transport and resulting bed morphology, than kaolinite flows. The paper proposes new stability diagrams for the phase space of bentonite flows and compares these to past work on kaolinite suspension flows. These new data suggest that the transitional‐flow Reynolds number can be used to delineate the types of transitional flow across different clay types and assess modern and ancient clay‐suspension flows. © 2016 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

Meander hydrodynamics initiated by river restoration deflectors

River restoration projects have installed j‐hook deflectors along the outer bank of meander bends to reduce hydraulic erosion, and in this study we use a computational fluid dynamics (CFD) model to document how these deflectors initiate changes in meander hydrodynamics. We validated the CFD with streamwise and cross‐channel bankfull velocities from a 193° meander bend flume (inlet at 0°) with a fixed point bar and pool equilibrium bed but no j‐hooks, and then used the CFD to simulate changes to flow initiated by bank‐attached boulder j‐hooks (1^st attached at 70°, then a 2^nd at 160°). At bankfull and half bankfull flow the j‐hooks flattened transverse water surface slopes, formed backwater pools upstream of the boulders, and steepened longitudinal water slopes across the boulders and in the conveyance region off the mid‐channel boulder tip. Streamwise velocity and mass transport jets upstream of the j‐hooks were stilled, mid‐channel jets were initiated in the conveyance region, eddies with a cross‐channel axis formed below boulders, and eddies with a vertical axis were shed into wake zones downstream of the point bar and outer bank boulders. At half bankfull depth conveyance region flow cut toward the outer bank downstream of the j‐hook boulders and the secondary circulation cells were reshaped. At bankfull depth the j‐hook at 160° was needed to redirect bank‐impinging flow sent by the upstream j‐hook. The hooked boulder tip of both j‐hooks funneled surface flow into mid‐channel plunging jets, which reversed the secondary circulation cells and initiated 1 to 3 counter rotating cells through the entire meander. The main outer bank collision zone centered at 50° without the j‐hook was moved by the j‐hook to within and just beyond the 70° j‐hook boulder region, which displaced other mass transport zones downstream. J‐hooks re‐organized water surface slopes, streamwise and cross‐channel velocities, and mass transport patterns, to move shear stress from the outer bank and into the conveyance and mid‐channel zones at bankfull flow. At half bankfull flows a patch of high shear re‐attached to the outer bank below the downstream j‐hook. J‐hook geometry and placement within natural meanders can be analyzed with CFD models to help restoration teams reach design goals and understand hydraulic impacts. Copyright © 2011 John Wiley & Sons, Ltd.     

The European Soil Erosion Model (EUROSEM): a dynamic approach for predicting sediment transport from fields and small catchments

The European Soil Erosion Model (EUROSEM) is a dynamic distributed model, able to simulate sediment transport, erosion and deposition over the land surface by rill and interill processes in single storms for both individual fields and small catchments. Model output includes total runoff, total soil loss, the storm hydrograph and storm sediment graph. Compared with other erosion models, EUROSEM has explicit simulation of interill and rill flow; plant cover effects on interception and rainfall energy; rock fragment (stoniness) effects on infiltration, flow velocity and splash erosion; and changes in the shape and size of rill channels as a result of erosion and deposition. The transport capacity of runoff is modelled using relationships based on over 500 experimental observations of shallow surface flows. EUROSEM can be applied to smooth slope planes without rills, rilled surfaces and surfaces with furrows. Examples are given of model output and of the unique capabilities of dynamic erosion modelling in general. © 1998 John Wiley & Sons, Ltd.     

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.     

Macroscale surface roughness and frictional resistance in overland flow

The hydraulics of overland flow on rough granular surfaces can be modelled and evaluated using the inundation ratio rather than the flow Reynolds number, as the primary dimensionless group determining the flow behaviour. The inundation ratio describes the average degree of submergence of the surface roughness and is used to distinguish three flow regimes representing partially inundated, marginally inundated and well-inundated surfaces. A heuristic physical model for the flow hydraulics in each regime demonstrates that the three states of flow are characterized by very different functional dependencies of frictional resistance on the scaled depth of flow. At partial inundation, flow resistance is associated with the drag force derived from individual roughness and therefore increases with depth and percentage cover. At marginal inundation, the size of the roughness elements relative to the depth of flow controls the degree of vertical mixing in the flow so that frictional resistance tends to decrease very rapidly with increasing depth of flow. Well-inundated flows are described using rough turbulent flow hydraulics previously developed for open channel flows. These flows exhibit a much more gradual decrease in frictional resistance with increasing depth than that observed during marginal inundation. A data set compiled from previously published studies of overland flow hydraulics is used to assess the functional dependence of frictional resistance on inundation ratio over a wide range of flow conditions. The data confirm the non-monotonic dependence predicted by the model and support the differentiation of three flow regimes based on the inundation ratio. Although the percentage cover and the surface slope may be of importance in addition to the inundation ratio in the partially and marginally inundated regimes, the Reynolds number appears to be of significance only in describing well-inundated flows at low to moderate Reynolds numbers. As these latter conditions are quite rare in natural environments, the inundation ratio rather than the Reynolds number should be used as the primary dimensionless group when evaluating the hydraulics of overland flow on rough surfaces. © 1997 by John Wiley & Sons, Ltd.     

Evaluating competing hypotheses for the origin and dynamics of river anastomosis

Anastomosing rivers have multiple interconnected channels that enclose flood basins. Various theories potentially explain this pattern, including an increased discharge conveyance and sediment transport capacity of multiple channels, deltaic branching, avulsion forced by base‐level rise, or a tendency to avulse due to upstream sediment overloading. The former two imply a stable anabranching channel pattern, whereas the latter two imply disequilibrium and evolution towards a single‐channel pattern in the absence of avulsion. Our objective is to test these hypotheses on morphodynamic scenario modelling and data of a well‐documented case study: the upper Columbia River. Proportions of channel and floodplain sediments along the river valley were derived from surface mapping. Initial and boundary conditions for the modelling were derived from field data. A 1D network model was built based on gradually varied flow equations, sediment transport prediction, mass conservation, transverse slope and spiral meander flow effects at the bifurcations. The number of channels and crevasse splays decreases in a downstream direction. Also, measured sediment transport is higher at the upstream boundary than downstream. These observations concur with bed sediment overloading from upstream, which can have caused channel aggradation above the surrounding floodplain and subsequent avulsion. The modelling also indicates that avulsion was likely caused by upstream overloading. In the model, multi‐channel systems inevitably evolve towards single‐channel systems within centuries. The reasons are that symmetric channel bifurcations are inherently unstable, while confluenced channels have relatively less friction than two parallel channels, so that more discharge is conveyed through the path with more confluences and less friction. Furthermore, the present longitudinal profile curvature of the valley could only be reproduced in the model by temporary overfeeding. We conclude that this anastomosing pattern is the result of time‐varying sediment overloading and is not an equilibrium pattern feature, and suggest this is valid for many anastomosing rivers. Copyright © 2012 John Wiley & Sons, Ltd.