Reach‐scale channel geometry of a mountain river

Mountain rivers can be subject to strong constraints imposed by changes in gradient and grain size supplied by processes such as glaciation and rockfall. Nonetheless, adjustments in the channel geometry and hydraulics of mountain rivers at the reach scale can produce discernible patterns analogous to those in fully alluvial rivers. Mountain rivers can differ in that imposed reach‐scale gradient is an especially important control on reach‐scale channel characteristics, as indicated by examination of North St Vrain Creek in Colorado. North St Vrain Creek drains 250 km² of the Rocky Mountains. We used 25 study reaches within the basin to examine controls on reach‐scale channel geometry. Variables measured included channel geometry, large woody debris, grain size, and mean velocity. Drainage area at the study reaches ranged from 2·2 to 245 km², and gradient from 0·013 to 0·147 m m^?1. We examined correlations among (1) potential reach‐scale response variables describing channel bankfull dimension and shape, hydraulics, bedform wavelength and amplitude, grain size, ?ow resistance, standard deviation of hydraulic radius, and volume of large woody debris, and (2) potential control variables that change progressively downstream (drainage area, discharge) or that are likely to re?ect a reach‐speci?c control (bed gradient). We tested the hypothesis that response variables correlate most strongly with local bed gradient because of the segmented nature of mountain channels. Results from simple linear regression analyses indicate that most response variables correlate best with gradient, although channel width and width/depth ratio correlate best with discharge. Multiple regression analyses using Mallow's C_p selection criterion and log‐transformation of all variables produced similar results in that most response variables correlate strongly with gradient. These results suggest that the hypothesis is partially supported: channel bed gradient is likely to be a good predictor for many reach‐scale response variables along mountain rivers, but discharge is also an important predictor for some response variables. Copyright © 2004 John Wiley & Sons, Ltd.     

Improvement of streams hydro‐geomorphological assessment using LiDAR DEMs

Hydro‐geomorphological assessments are an essential component for riverine management plans. They usually require costly and time‐consuming field surveys to characterize the spatial variability of key variables such as flow depth, width, discharge, water surface slope, grain size and unit stream power throughout the river corridor. The objective of this research is to develop automated tools for hydro‐geomorphological assessments using high‐resolution LiDAR digital elevation models (DEMs). More specifically, this paper aims at developing geographic information system (GIS) tools to extract channel slope, width and discharge from 1 m‐resolution LiDAR DEMs to estimate the spatial distribution of unit stream power in two contrasted watersheds in Quebec: a small agricultural stream (Des Fèves River) and a large gravel‐bed river (Matane River). For slope, the centreline extracted from the raw LiDAR DEM was resampled at a coarser resolution using the minimum elevation value. The channel width extraction algorithm progressively increased the centerline from the raw DEM until thresholds of elevation differences and slopes were reached. Based on the comparison with over 4000 differential global positioning system (GPS) measurements of the water surface collected in a 50 km reach of the Matane River, the longitudinal profile and slope estimates extracted from the raw and resampled LiDAR DEMs were in very good agreement with the field measurements (correlation coefficients ranging from 0 · 83 to 0 · 87) and can thus be used to compute stream power. The extracted width also corresponded very well to the channel as seen from ortho‐photos, although the presence of bars in the Matane River increased the level of error in width estimates. The estimated maximum unit stream power spatial patterns corresponded well with field evidence of bank erosion, indicating that LiDAR DEMs can be used with confidence for initial hydro‐geomorphological assessments. Copyright © 2013 John Wiley & Sons, Ltd.     

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

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

Geometric characteristics and evolution of a tidal channel network in experimental setting

This paper reports on a laboratory study that aims to reproduce a tidal channel network, in order to enhance the understanding of the morphodynamic evolution of the channel characteristics as the network expands and finally reaches equilibrium. A high‐resolution laser system scanned the bed topography at different time steps creating multiple digital elevation models of the channel network. Two hundred and seventy individual channel segments are analyzed and cross‐correlated in terms of their width, depth and length. The laboratory results show positive linear correlations between depth and width as well as between length and width of channel segments of the network configuration at final equilibrium. In a downstream direction, channels appear to widen more than they deepen, indirectly a sign that discharge has a stronger control on channel width than on depth. In contrast to fluvial drainage networks that commonly display fractal and scale‐invariant behavior, the geometric properties of the experimental tidal creek network shows scale dependence. Channel attributes exhibit consistent patterns of exponentially decreasing abundance, with increasing creek length, depth and width. The nature of the observed exponential distributions within creek attributes (width, depth, length) allows for statistical predictability of relative creek attribute dimensions downstream and through time. Copyright © 2010 John Wiley & Sons, Ltd.     

Predicting river width,depth and velocity at ungauged sites in England and Wales using multilevel models

Using a dataset of gauged river discharges taken from sites in England and Wales, linear multilevel models (also known as mixed effects models) were applied to quantify the variability in discharge and the discharge‐hydraulic geometry relationships across three nested spatial scales. A jackknifing procedure was used to test the ability of the multilevel models to predict hydraulic geometry, and therefore width, mean depth and mean velocity, at ungauged stations. These models provide a framework for making predictions of hydraulic geometry parameters, with associated levels of uncertainty, using different levels of data availability. Results indicate that as one travels downstream along a river there is greater variability in hydraulic geometry than is the case between rivers of similar sizes. This indicates that hydraulic geometry (and therefore hydrology) is driven by catchment area, to a greater extent than by natural geomorphological variations in the streamwise direction at the mesoscale, but these geomorphological variations can still have a major impact on channel structure. Copyright © 2008 John Wiley & Sons, Ltd.     

Wood distribution in neotropical forested headwater streams of La Selva,Costa Rica

Surveys of wood along 30 forested headwater stream reaches in La Selva Biological Station in north‐eastern Costa Rica represent the first systematic data reported on wood loads in neotropical streams. For streams with drainage areas of 0·1–8·5 km² and gradients of 0·2–8%, wood load ranged from 3 to 34·7 m³ wood/100 m channel and 41–612 m³ wood/ha channel. These values are within the range reported for temperate streams. The variables wood diameter/flow depth, stream power, the presence of backflooding, and channel width/depth are consistently selected as significant predictors by statistical models for wood load. These variables explain half to two‐thirds of the variability in wood load. These results, along with the spatial distribution of wood with respect to the thalweg, suggest that transport processes exert a greater influence on wood loads than recruitment processes. Wood appears to be more geomorphically effective in altering bed elevations in gravel‐bed reaches than in reaches with coarser or finer substrate. Copyright © 2009 John Wiley & Sons, Ltd.     

Bed scour by debris flows: experimental investigation of effects of debris‐flow composition

Debris flows can grow greatly in size by entrainment of bed material, enhancing their runout and hazardous impact. Here, we experimentally investigate the effects of debris‐flow composition on the amount and spatial patterns of bed scour and erosion downstream of a fixed to erodible bed transition. The experimental debris flows were observed to entrain bed particles both grain by grain and en masse, and the majority of entrainment was observed to occur during passage of the flow front. The spatial bed scour patterns are highly variable, but large‐scale patterns are largely similar over 22.5–35° channel slopes for debris flows of similar composition. Scour depth is generally largest slightly downstream of the fixed to erodible bed transition, except for clay‐rich debris flows, which cause a relatively uniform scour pattern. The spatial variability in the scour depth decreases with increasing water, gravel (= grain size) and clay fraction. Basal scour depth increases with channel slope, flow velocity, flow depth, discharge and shear stress in our experiments, whereas there is no correlation with grain collisional stress. The strongest correlation is between basal scour and shear stress and discharge. There are substantial differences in the scour caused by different types of debris flows. In general, mean and maximum scour depths become larger with increasing water fraction and grain size, and decrease with increasing clay content. However, the erodibility of coarse‐grained experimental debris flows (gravel fraction = 0.64) is similar on a wide range of channel slopes, flow depths, flow velocities, discharges and shear stresses. This probably relates to the relatively large influence of grain‐collisional stress to the total bed stress in these flows (30–50%). The relative effect of grain‐collisional stress is low in the other experimental debris flows (<5%), causing erosion to be largely controlled by basal shear stress. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical simulation of three‐dimensional flow hydraulics in a braided channel

Results are presented from a numerical simulation of three‐dimensional flow hydraulics around a mid‐channel bar carried out using the FLUENT/UNS computational fluid dynamics (CFD) software package. FLUENT/UNS solves the three‐dimensional Reynolds‐averaged form of the Navier–Stokes equations. Turbulence closure is achieved using a RNG k–ϵ model. Simulated flow velocities are compared with measured two‐dimensional velocities (downstream and cross‐stream) obtained using an electromagnetic current meter (ECM). The results of the simulation are qualitatively consistent with the flow structures observed in the field. Quantitative comparison of the simulated and measured velocity magnitudes indicates a strong positive correlation between the two (r=0·88) and a mean difference of 0·09 m s⁻¹. Deviations between simulated and measured velocities may be identified that are both random and systematic. The former may reflect a number of factors including subgrid‐scale natural spatial variability in flow velocities associated with local bed structures and measurement uncertainty resulting from problems of ECM orientation. Model mesh configuration, roughness parameterization and inlet boundary condition uncertainty may each contribute to systematic differences between simulated and measured flow velocities. These results illustrate the potential for using CFD software to simulate flow hydraulics in natural channels with complex configurations. They also highlight the need for detailed spatially distributed datasets of three‐dimensional flow variables to establish the accuracy and applicability of CFD software. Copyright © 1999 John Wiley & Sons, Ltd.     

Evolution of a bifurcation in a meandering river with adjustable channel widths,Rhine delta apex,The Netherlands

Rivers may dramatically change course on a fluvial plain. Such an avulsion temporarily leads to two active channels connected at a bifurcation. Here we study the effect of dynamic meandering at the bifurcation and the effect of channel width adjustment to changing discharge in both downstream branches on the evolution of a bifurcation and coexisting channels. As an example, we reconstructed the last major avulsion at the Rhine delta apex. We combined historical and geological data to reconstruct a slowly developing avulsion process spanning 2000 years and involving channel width adjustment and meandering at the bifurcation. Based on earlier idealised models, we developed a one‐dimensional model for long‐term morphodynamic prediction of upstream channel and bifurcates connected at the bifurcation node. The model predicts flow and sediment partitioning at the node, including the effect of migrating meanders at the bifurcation and channel width adjustment. Bifurcate channel width adaptation to changing discharge partitioning dramatically slows the pacing of bifurcation evolution because the sediment balance for width adjustment and bed evolution are coupled. The model further shows that meandering at the bifurcation modulates channel abandonment or enlargement periodically. This explains hitherto unrecognised reactivation signals in the sedimentary record of the studied bifurcation meander belts, newly identified in our geological reconstruction. Historical maps show that bifurcation migration due to meander bend dynamics increases the bifurcation angle, which increases the rate of closure of one bifurcate. The combination of model and reconstruction identifies the relevant timescales for bifurcation evolution and avulsion duration. These are the time required to fill one downstream channel over one backwater length, the time to translate one meander wavelength downstream and, for strong river banks, the adaptation timescale to adjust channel width. The findings have relevance for all avulsions where channel width can adjust to changing discharge and where meandering occurs. Copyright © 2011 John Wiley & Sons, Ltd.     

Channel adjustments of the lower Trinity River,Texas, downstream of Livingston Dam

Channel cross‐sectional changes since construction of Livingston Dam and Lake Livingston in 1968 were studied in the lower Trinity River, Texas, to test theoretical models of channel adjustment, and to determine controls on the spatial extent of channel response. High and average flows were not significantly modified by the dam, but sediment transport is greatly reduced. The study is treated as an opportunistic experiment to examine the effects of a reduction in sediment supply when discharge regime is unchanged. Channel scour is evident for about 60 km downstream, and the general phenomena of incision, widening, coarsening of channel sediment and a decrease in channel slope are successfully predicted, in a qualitative sense, by standard models of channel response. However, there is no consistent channel response within this reach, as various qualitatively different combinations of increases, decreases or no change in width, depth, slope and roughness occur. These multiple modes of adjustment are predicted by the unstable hydraulic geometry model. Between about 60 km and the Trinity delta 175 km downstream of the dam, no morphological response to the dam is observed. Rather than a diminution of the dam's effects on fluvial processes, this is due to a fundamental change in controls of the fluvial system. The downstream end of the scour zone corresponds to the upstream extent of channel response to Holocene sea level rise. Beyond 60 km downstream, the Trinity River is characterized by extensive sediment storage and reduced conveyance capacity, so that even after dam construction sediment supply still exceeds transport capacity. The channel bed of much of this reach is near or below sea level, so that sea level rise and backwater effects from the estuary are more important controls on the fluvial system than upstream inputs. Copyright © 2005 John Wiley & Sons, Ltd.     

An experimental test of whether bar instability contributes to the formation,periodicity and maintenance of pool–riffle sequences

Pool–riffle sequences (PRSs) are periodic river‐bed morphologies with wavelengths several times the channel width. Causes of PRS formation and maintenance are not clearly understood, which may limit the effectiveness of protection and rehabilitation measures. Some confusion has existed about whether the PRS morphology is the same as or distinct from alternate bars. In this paper we investigate whether the bar instability forming alternate bars also contributes to PRS formation, periodicity and maintenance. This was unclear because bar instability occurs only when the ratio of channel width/depth exceeds a critical value, generally understood to be approximately 10, which is larger than the width‐to‐depth ratio of many PRSs. A mobile‐bed physical model is used to test whether bar instability occurs in channels characteristic of PRS morphology, with low width‐to‐depth ratio, and high relative roughness. The physical model was scaled from a prototype PRS in a gravel and cobble bed river. Alternate bars formed in the model at channel width‐to‐depth ratios as low as 3·8. The wavelength of the alternate bars formed was generally 2·2–5 times channel width, which was similar to the prototype PRS. Therefore, bar instability can occur in virtually all PRSs, and it contributes to the widespread formation of periodic PRS morphology. The model showed that maintenance of the bar height in the prototype PRS also depends on variations in channel width. It is concluded that periodic PRSs are formed and maintained by the interaction between bar instability, and flow deflections associated with variations in channel geometry such as width variation. Resonance between bar instability and three‐dimensional bed forms such as alternate bars and variations in channel geometry. Variations in channel geometry are also important in determining the location and dimensions of individual pools and bars. Copyright © 2008 John Wiley & Sons, Ltd.     

Comparison of hydraulic geometry between sand‐ and gravel‐bed rivers in relation to channel pattern discrimination

A comparison has been made between the hydraulic geometry of sand‐ and gravel‐bed rivers, based on data from alluvial rivers around the world. The results indicate a signi?cant difference in hydraulic geometry among sand‐ and gravel‐bed rivers with different channel patterns. On this basis, some diagrams for discrimination of meandering and braided channel patterns have been established. The relationships between channel width and water discharge, between channel depth and water discharge, between width–depth ratio and water discharge and between channel slope and water discharge can all be used for channel pattern discrimination. The relationship between channel width and channel depth can also be used for channel pattern discrimination. However, the accuracy of these relationships for channel pattern discrimination varies, and the depth–discharge relationship is a better discriminator of pattern type than the classic slope–discharge function. The cause for this difference has been explained qualitatively. To predict the development of channel patterns under different natural conditions, the pattern discriminator should be searched on the basis of independent or at least semi‐independent variables. The relationship between stream power and bed material grain size can be used to discriminate channel patterns, which shows a better result than the discriminator using the slope–discharge relationship. Copyright © 2004 John Wiley & Sons, Ltd.     

Flow behaviour,suspended sediment transport and transmission losses in a small (sub‐bank‐full) flow event in an Australian desert stream

The behaviour of a discrete sub‐bank‐full flow event in a small desert stream in western NSW, Australia, is analysed from direct observation and sediment sampling during the flow event and from later channel surveys. The flow event, the result of an isolated afternoon thunderstorm, had a peak discharge of 9 m³/s at an upstream station. Transmission loss totally consumed the flow over the following 7·6 km. Suspended sediment concentration was highest at the flow front (not the discharge peak) and declined linearly with the log of time since passage of the flow front, regardless of discharge variation. The transmission loss responsible for the waning and eventual cessation of flow occurred at a mean rate of 13.2% per km. This is quite rapid, and is more than twice the corresponding figure for bank‐full flows estimated by Dunkerley (1992) on the same stream system. It is proposed that transmission losses in ephemeral streams of the kind studied may be minimized in flows near bank‐full stage, and be higher in both sub‐bank‐full and overbank flows. Factors contributing to enhanced flow loss in the sub‐bank‐full flow studied included abstractions of flow to pools, scour holes and other low points along the channel, and overflow abstractions into channel filaments that did not rejoin the main flow. On the other hand, losses were curtailed by the shallow depth of banks wetted and by extensive mud drapes that were set down over sand bars and other porous channel materials during the flow. Thus, in contrast with the relatively regular pattern of transmission loss inferred from large floods, losses from low flows exhibit marked spatial variability and depend to a considerable extent on streamwise variations in channel geometry, in addition to the depth and porosity of channel perimeter sediments. Copyright © 1999 John Wiley & Sons, Ltd.     

Controls on the transition from meandering to straight channels in the wetlands of the Okavango Delta,Botswana

Previous studies of alluvial rivers have shown that channel patterns form a continuum controlled by interactions among factors such as gradient, discharge, sediment size, and bank strength. Data from channels in the permanent wetlands of the Okavango Delta add to these ?ndings by focusing on pattern transitions in channels with banks formed by sedges and grasses that are rooted in peat and underlain by unconsolidated sediment. Channels are well de?ned, and transport ?ne–medium sand as bedload between the vertical, vegetation‐lined banks. Water depths, velocities, grain sizes, and bankline vegetation do not vary signi?cantly or systematically downstream, but the permeable banks allow water to leak from the channels, contributing to an overall downstream decrease in discharge and width. In addition, as the Okavango River ?ows from the <12 km wide ‘Panhandle’ and splits into distributaries in the broader ‘Fan’, valley gradient steepens by c. 60 per cent. These downstream changes result in channel pattern adjustments. In the Panhandle, the Okavango River is a relatively wide (c. 30–100 m), actively meandering, sinuous channel (P > 2·0), but further downstream in the Fan, the narrower (<40 m) distributaries follow laterally stable, less sinuous (‘straight’) courses (P < 1·75). Some channel pattern discrimination diagrams based on simple indices of gradient, discharge, sediment size or stream power are inadequate for analysing the meandering–straight transition in the Okavango but Parker's (1976) approach, based on ratios of depth–width and slope–Froude number, accurately characterizes the transition. Our ?eld observations, combined with the results from previous experimental studies, suggest that in relatively wide channels (w/d > 10), thalweg meandering results in scour of the unconsolidated sediment at the bank base, leading to undermining and collapse of the vegetation, and to slow meander migration. However, as channels narrow downstream (w/d < 10) with discharge losses, proportionally increasing sidewall drag exerted by bankline vegetation suppresses thalweg meandering and bank scour, and channels follow stable, less sinuous courses. Copyright © 2004 John Wiley & Sons, Ltd.     

Planform dynamics of braided streams

The high dynamism and complexity of braided networks poses a series of open questions, significant for river restoration and management. The present work is aimed at the characterization of the morphology of braided streams, in order to assess whether the system reaches a steady state under constant flow conditions and, in that case, to determine how it can be described and on which parameters it depends. A series of 14 experimental runs were performed in a laboratory physical model with uniform sand, varying the discharge and the longitudinal slope. Planimetric and altimetric configurations were monitored in order to assess the occurrence of a steady state. A set of parameters was considered, such as the braid‐plain width and the number and typology of branches and nodes. Results point out that a relationship exists between braiding morphology and two dimensionless parameters, related to total water discharge and stream power. We found that network complexity increases at higher values of water discharge and a larger portion of branches exhibits morphological activity. Results are then compared to the outputs of a simple one‐dimensional model, that allows to easily predict the average network complexity, once the bed topography is known. Model computations permit also the investigation of the effect of water discharge variations and to compare different width definitions. The at‐a‐station variability of planimetric parameters shows a peculiar behaviour, both regarding number of branches and wetted width. In particular, the analysis of the relationship between width and discharge highlighted relevant differences in comparison to single thread channel. Copyright © 2009 John Wiley & Sons, Ltd.     

Investigation of spatially distributed braided river flows using a two‐dimensional hydraulic model

This paper reports the application of a two‐dimensional hydraulic model to a braided reach of the Avoca River, New Zealand. Field measurements of water surface elevation, depth and velocity obtained at low flow were used to validate the model and to optimize the parameterization of bed friction. The main systematic trends in the measured flow variables are reproduced by the model. However, field data are characterized by greater spatial variability than model output reflecting differences in the scale of processes measured in the field and represented by the model. Additional model runs were conducted to simulate flow patterns within the study reach at five higher discharges. The purpose of these simulations was to evaluate the potential for using two‐dimensional hydraulic models to quantify the reach‐scale hydraulic characteristics of braided rivers and their dependence on discharge. Changes in flow depth and velocity with increasing discharge exhibit trends that are consistent with the results of previous field investigations, although the tendency for the wetted area of the braidplain within particular depth and velocity categories to remain fixed as discharge rises, as has been noted for several braided rivers in New Zealand, was not observed. Modelled shear stress frequency distributions fit gamma functions that incorporate a distribution shape parameter, the value of which follows clear systematic trends with rising discharge. These results illustrate both the problems of, and potential for, using two‐dimensional hydraulic models in braided river applications. This leads to something of a paradox in that while such models provide a means of generating hydraulic information that would be difficult to obtain in the field at an equivalent spatial resolution, they are, due to the problems inherent to data collection, difficult to validate conclusively. Despite this limitation, the application of spatially distributed models to investigate relationships between discharge and reach‐scale form and process variables appears to have considerable potential. Copyright © 2003 John Wiley & Sons, Ltd.     

At‐a‐station hydraulic geometry relations, 1: theoretical development

This paper, the first of two, hypothesizes that: (1) the temporal variation of stream power of a river channel at a given station with varying discharge is accomplished by the temporal variation in channel form (flow depth and channel width) and hydraulic variables, including energy slope, flow velocity and friction; (2) the change in stream power is distributed among the changes in flow depth, channel width, flow velocity, slope, and friction, depending on the boundary conditions that the channels has to satisfy. The second hypothesis is a result of the principle of maximum entropy and the theory of minimum energy dissipation or its simplified minimum stream power. These two hypotheses lead to families of at‐a‐station hydraulic geometry relations. The conditions under which these families of relations can occur in the field are discussed. Copyright © 2007 John Wiley & Sons, Ltd.     

Annual bed‐elevation regime in the alluvial channel of Squamish River,Southwestern British Columbia,Canada

The aim of this study is to examine the annual regime of channel scour and fill by monitoring bed‐elevation changes in a reach of Squamish River in southwestern British Columbia, Canada. Sonar surveys of 13 river cross‐sections in a sandy gravel‐bed single‐channel study reach were repeated biweekly over a full hydrologic year (1995/6). The survey results show that bedload movement occurs as waves or pulses forming bedwaves that appear to maintain an overall coherence with movement downstream. These bedwaves propagate downstream by a mode here termed pulse scour and pulse fill, a process distinguished from the conventional mode of scour and fill commonly associated with flood events (here termed local scour and local fill). Bedwave celerity was estimated to be about 15·5 m d⁻¹ corresponding to a bedwave residence time in the study reach of almost one hydrologic year. The total amount of local bed‐elevation change ranged between 0·22 m and 2·41 m during the period of study. Analysis of the bed‐elevation and flow data reveals that, because of the bedwave phenomenon, there is no simple relation between the mean bed‐elevation and discharge nor any strong linear correlation among cross‐sectional behaviour. The bed‐elevation data also suggest that complex changes to the bed within a cross‐section are masked when the bed is viewed in one dimension, although no definitive trends in bed behaviour were found in the two‐dimensional analysis. Although a weak seasonal effect is evident in this study, the bed‐elevation regime is dominated by sediment supply‐driven fluctuations in bedload transport occurring at timescales shorter than the seasonal fluctuation in discharge. The study also indicates that bed‐elevation monitoring on Squamish River, and others like it, for purposes of detecting and measuring aggradation/degradation must take into account very considerable and normal channel‐bed variability operating at timescales from hours to months. Copyright © 2000 John Wiley & Sons, Ltd.     

Channel change,bankfull and effective discharges on a vertically accreting,meandering, gravel‐bed river

Along the lower reaches of the Waipaoa River, New Zealand, cross‐section survey data indicate there was a 23 per cent decrease in bankfull width and a 22 per cent reduction in channel cross‐section area between 1948 and 2000, as the channel responded to increased inputs of fine (suspended) sediment following deforestation of the headwaters in late C19 and early C20. We determined the bankfull discharge within a ～39 km long reach by routing known discharges through the one‐dimensional MIKE 11 flow model. The model runs suggest that the bankfull discharge varies between ～800 and ～2300 m³ s^?1 and that the average recurrence interval is 4 ± 2 years on the annual maximum series; by contrast, the effective flow (360 m³ s^?1) is equaled or exceeded three times a year. The variability in bankfull discharge arises because the banks tend to be lower in places where flood flows are constricted than in reaches where overbank flow is dispersed over a wide area, and because scour has counteracted aggradation in some locations. There is no downstream variation in Shields stress, or in relative shear stress, within the study reach. Bankfull shear stress is, on average, five times greater than the shear stress required to initiate motion. At the effective discharge it is more than twice the threshold value. The effective discharge probably has more relevance than the bankfull discharge to the overall picture of sediment movement in the lower reaches of the Waipaoa River but, because width is constrained by the stability and resistance of the bank material to erosion during high flows that also scour the bed, the overall channel geometry is likely determined by discharges at or near bankfull. Copyright © 2006 John Wiley & Sons, Ltd.     

Organic matter and fine grains as possible determinants of spatial and seasonal variability in bed sediment fauna: A case study from a Hercynian gravel stream

The spatial and seasonal distributions of organic matter and fine grains were tested as possible determinants of fauna distribution in bed sediment of a Hercynian gravel stream. Invertebrate densities and the amounts of fine grains and organic carbon were assessed in freeze-core samples taken along 70 cm depth profiles at three different positions in the stream channel. Sampling was conducted on five occasions of low discharge over two years. The variability in invertebrate community composition was analysed using Detrended Correspondence Analysis with posterior projection of explanatory variables; Variation Partitioning was used to estimate the independent and shared effects of the explanatory variables. We found that the best predictors of the invertebrate community were spatial variables (depth, position in the channel) and then variables influenced by seasonal patterns (surface water temperature and discharge). The influence of organic matter and fine grain content was significant only after eliminating spatial autocorrelation. High amounts of organic matter, randomly accumulated in the sediment, improved the model by explaining high fauna densities. The fine grain content was not a limiting factor to fauna at our study site. It is possible that the large amount of mica flakes in the sediment has caused the arrangement of grains with a pore space sufficient for fauna even when fine grain content was high.