Riffle‐pool morphometry and stage‐dependant morphodynamics of a large floodplain river (Vereinigte Mulde,Sachsen‐Anhalt,Germany)

Riffle‐pool sequences are a common feature of gravel‐bed rivers. However, mechanisms of their generation and maintenance are still not fully understood. In this study a monitoring approach is employed that focuses on analysing cross‐sectional and longitudinal channel geometry of a large floodplain river (Vereinigte Mulde, Sachsen‐Anhalt, Germany) with a high temporal and spatial resolution, in order to conclude from stage‐dependant morphometric changes to riffle and pool maintaining processes. In accordance with previous authors, pool cross‐sections of the Mulde River are narrow and riffle cross‐sections are wide suggesting that they should rather be addressed as two general types of channel cross‐sections than solely as bedforms. At high flows, riffles and pools in the study reaches changed in length and height but not in position. Pools were scoured and riffles aggraded, a development which was reversed during receding flows below the threshold of 0·4Q_bf (40% bankfull discharge). An index for the longitudinal amplitude of riffle‐pool sequences, the bed undulation intensity or bedform amplitude, is introduced and proved to be highly significant as a form parameter, its first derivative as a process parameter. The process of pool scour and riffle fill is addressed as bedform maintenance or bedform accentuation. It is indicated by increasing longitudinal bed amplitudes. According to the observed dynamics of bed amplitudes, maintenance of riffle‐pool sequences lags behind discharge peaks. Maximum bed amplitudes may be reached with a delay of several days after peak discharges. Increasing bed undulation intensity is interpreted to indicate bed mobility. Post‐flood decrease of the bed undulation intensity indicates a retrograde phase when transport from pools to riffles has ceased and bed mobility is restricted to riffle tails and heads of pools. This type of transport behaviour is referred to as disconnected mobility. The comparison of two river reaches, one with undisturbed sediment supply, the other with sediment deficit, suggests that high bed undulation intensity values at low flows indicate sediment deficit and potentially channel degrading conditions. It is more generally hypothesized that channel bed undulations constitute a major component of form roughness and that increased bed amplitudes are an important feature of channel bed adjustment to sediment deficit be it temporally during late floods or permanently due to a supply limitation of bedload. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

Recent variation in reach‐scale bankfull discharge in the Lower Yellow River

Bankfull discharge is a key parameter in the context of river engineering and geomorphology, as an indicator of flood discharge capacity in alluvial rivers, and varying in response to the incoming flow and sediment regimes. Bankfull channel dimensions have significantly adjusted along the Lower Yellow River (LYR) due to recent channel degradation, caused by the operation of the Xiaolangdi Reservoir, which has led to longitudinal variability in cross‐sectional bankfull discharges. Therefore, it is more representative to describe the flood discharge capacity of the LYR, using the concept of reach‐averaged bankfull discharge. Previous simple mean methods to estimate reach‐scale bankfull discharge cannot meet the condition of flow continuity or account for the effect of different spacing between two sections. In this study, a general method to calculate cross‐sectional bankfull discharge using the simulated stage‐discharge relation is outlined briefly, and an integrated method is then proposed for estimating reach‐scale bankfull discharge. The proposed method integrates a geometric mean based on the log‐transformation with a weighted average based on the spacing between two consecutive sections, which avoids the shortcomings of previous methods. The post‐flood reach‐scale bankfull discharges in three different channel‐pattern reaches of the LYR were estimated annually during the period from 1999 to 2010 using the proposed method, based on surveyed post‐flood profiles at 91 sedimentation sections and the measured hydrological data at seven hydrometric sections. The calculated results indicate that: (i) the estimated reach‐scale bankfull discharges can effectively represent the flood discharge capacity of different reaches, with their ranges of variation being less than those of typical cross‐sectional bankfull discharges; and (ii) the magnitude of the reach‐scale bankfull discharge in each reach can respond well to the accumulative effect of incoming flow and sediment conditions. Finally, empirical relationships for different reaches in the LYR were developed between the reach‐scale bankfull discharge and the previous four‐year average discharge and incoming sediment coefficient during flood seasons, with relatively high correlation coefficients between them being obtained, and the reach‐scale bankfull discharges in different reaches predicted by the delayed response model were also presented for a comparison. These relations for the prediction of reach‐scale bankfull discharges were validated using the cross‐sectional profiles and hydrological data measured in 2011. Copyright © 2013 John Wiley & Sons, Ltd.     

Sediment routing hypothesis for pool‐riffle maintenance

This paper provides comprehensive evidence that sediment routing around pools is a key mechanism for pool‐riffle maintenance in sinuous upland gravel‐bed streams. The findings suggest that pools do not require a reversal in energy for them to scour out any accumulated sediments, if little or no sediments are fed into them. A combination of clast tracing using passive integrated transponder (PIT) tagging and bedload traps (positioned along the thalweg on the upstream riffle, pool entrance, pool exit and downstream riffle) are used to provide information on clast pathways and sediment sorting through a single pool‐riffle unit. Computational fluid dynamics (CFD) is also used to explore hydraulic variability and flow pathways. Clast tracing results provide a strong indication that clasts are not fed through pools, rather they are transported across point bar surfaces, or around bar edges (depending upon previous clast position, clast size, and event magnitude). Spatial variations in bedload transport were found throughout the pool‐riffle unit. The pool entrance bedload trap was often found to be empty, when the others had filled, further supporting the notion that little or no sediment was fed into the pool. The pool exit slope trap would occasionally fill with sediment, thought to be sourced from the eroding outer bank. CFD results demonstrate higher pool shear stresses (τ ≈ 140 N m^–2) in a localized zone adjacent to an eroding outer bank, compared to the upstream and downstream riffles (τ ≈ 60 N m^–2) at flows of 6 · 2 m³ s^–1 (≈ 60% of the bankfull discharge) and above. There was marginal evidence for near‐bed velocity reversal. Near‐bed streamlines, produced from velocity vectors indicate that flow paths are diverted over the bar top rather than being fed through the thalweg. Some streamlines appear to brush the outer edge of the pool for the 4 · 9 m³ s^–1 to 7 · 8 m³ s^–1 (between 50 and 80% of the bankfull discharge) simulations, however complete avoidance was found for discharges greater than this. Copyright © 2013 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.     

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.     

Spatial variations in surface sediment structure in riffle–pool sequences: a preliminary test of the Differential Sediment Entrainment Hypothesis (DSEH)

Riffle–pool sequences are maintained through the preferential entrainment of sediment grains from pools rather than riffles. This preferential entrainment has been attributed to a reversal in the magnitude of velocity and shear stress under high flows; however the Differential Sediment Entrainment Hypothesis (DSEH) postulates that differential entrainment can instead result from spatial sedimentological contrasts. Here we use a novel suite of in situ grain‐scale field measurements from a riffle–pool sequence to parameterize a physically‐based model of grain entrainment. Field measurements include pivoting angles, lift forces and high resolution digital elevation models (DEMs) acquired using terrestrial laser scanning, from which particle exposure, protrusion and surface roughness were derived. The entrainment model results show that grains in pools have a lower critical entrainment shear stress than grains in either pool exits or riffles. This is because pool grains have looser packing, hence greater exposure and lower pivoting angles. Conversely, riffle and pool exit grains have denser packing, lower exposure and higher pivoting angles. A cohesive matrix further stabilizes pool exit grains. The resulting predictions of critical entrainment shear stress for grains in different subunits are compared with spatial patterns of bed shear stress derived from a two‐dimensional computational fluid dynamics (CFD) model of the reach. The CFD model predicts that, under bankfull conditions, pools experience lower shear stresses than riffles and pool exits. However, the difference in sediment entrainment shear stress is sufficiently large that sediment in pools is still more likely to be entrained than sediment in pool exits or riffles, resulting in differential entrainment under bankfull flows. Significantly, this differential entrainment does not require a reversal in flow velocities or shear stress, suggesting that sedimentological contrasts alone may be sufficient for the maintenance of riffle–pool sequences. This finding has implications for the prediction of sediment transport and the morphological evolution of gravel‐bed rivers. Copyright © 2012 John Wiley & Sons, Ltd.     

Morphological effects of different channel‐forming discharges in a gravel‐bed river

The study analyses the morphological response of a gravel‐bed river to discharges of different magnitude (from moderate events that occur several times a year to a 12‐year flood) and so defines the range of formative discharges for single morphological units (channels, bars, islands) and a range of magnitude of morphological activity from the threshold discharges for gravel transport and minor bar modification up to flows causing major morphological changes. The study was conducted on the Tagliamento River, a large gravel‐bed river in north‐eastern Italy, using two different methods, analysis of aerial photographs and field observation of painted gravel particles. The available photographs (five flights from August 1997 to November 2002) and the two commissioned flights (June 2006 and April 2007) do not define periods with a single flood event, but the intervals are short enough (11 to 22 months) to have a limited number of flood events in each case. The fieldwork, which involved cross‐section survey, grain‐size analysis and observation of painted sediments, complemented the aerial surveys by allowing analysis of channel response to single flood events. Substantial morphological changes (e.g. bank erosion of several tens of metres up to more than 100 m) associated with flood events with a recurrence interval between 1·1 year and 12 years have been documented. Multiple forming discharges were defined based on the activity of different morphological units. Discharges equal to 20–50% of the bankfull discharge are formative for the channels, whereas the bankfull discharge (1·1 year flood in this case of the Tagliamento River) is formative for low bars. Larger floods, but still relatively frequent (with a recurrence interval less than five years), are required for full gravel transport on high bars and significant morphological changes of islands. Copyright © 2009 John Wiley & Sons, Ltd.     

Inbank and overbank velocity conditions in an arid zone anastomosing river

Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q_2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.     

Random controls on semi‐rhythmic spacing of pools and riffles in constriction‐dominated rivers

Average pool spacing between five and seven bankfull widths has been documented in environments throughout the world, but has limited theoretical justification in coarse‐bedded and bedrock environments. Pool formation in coarse‐bedded and bedrock channels has been attributed to bedrock and boulder constrictions. Because the spacing of these constrictions may be irregular in nature, it is difficult to reconcile pool‐formation processes with the supposedly rhythmic spacing of pools and riffles. To address these issues, a simulation model for pool and riffle formation is used to demonstrate that semi‐rhythmic spacing of pools with an approximate spacing of five to seven bankfull widths can be recreated from a random distribution of obstructions and minimum pool‐ and riffle‐length criteria. It is assumed that a pool–riffle couplet will achieve a minimum length based on dominant‐discharge conditions. Values for the minimum‐length assumption are based on field data collected in New England and California, while the theoretical basis relies on the demonstrated hydraulic response of individual pools to elongation. Results from the simulations show that the location of pools can be primarily random in character, but still assume an average spacing between four and eight bankfull widths for a variety of conditions. Field verification data generally support the model but highlight a highly skewed distribution of pool‐forming elements and pool spacing. The relation between pool spacing and bankfull widths is attributed to the common geometric response of these features to dominant‐discharge conditions. Copyright © 2001 John Wiley & Sons, Ltd.     

Flood‐driven topographic changes in a gravel‐cobble river over segment,reach, and morphological unit scales

Regulated rivers generally incise below dams that cut off sediment supply, but how that happens and what the consequences are at different spatial scales is poorly understood. Modern topographic mapping at meter‐scale resolution now enables investigation of the details of spatial processes. In this study, spatial segregation was applied to a meter‐scale raster map of topographic change from 1999 to 2008 on the gravel‐cobble, regulated lower Yuba River in California to answer specific scientific questions about how a decadal hydrograph that included a flood peak of 22 times bankfull discharge affected the river at segment, reach, and morphological unit scales. The results show that the river preferentially eroded sediment from floodplains compared to the channel, and this not only promoted valley‐wide sediment evacuation, but also facilitated the renewal and differentiation of morphological units, especially in the channel. At the reach scale, area of fill and mean net rate of elevational change were directly correlated with better connectivity between the channel and floodplain, while the mean rate of scour in scour areas was influenced by the ratio of slope to bankfull Froude number, a ratio indicative of lateral migration versus vertical downcutting. Hierarchical segregation of topographic change rasters proved useful for understanding multi‐scalar geomorphic dynamics. Copyright © 2016 John Wiley & Sons, Ltd.     

One‐dimensional and two‐dimensional hydrodynamic modeling derived flow properties: impacts on aquatic habitat quality predictions

Studies of the effects of hydrodynamic model dimensionality on simulated flow properties and derived quantities such as aquatic habitat quality are limited. It is important to close this knowledge gap especially now that entire river networks can be mapped at the microhabitat scale due to the advent of point‐cloud techniques. This study compares flow properties, such as depth and velocity, and aquatic habitat quality predicted from pseudo‐2D and fully 2D hydrodynamic modeling. The models are supported by high‐resolution, point‐cloud derived bathymetries, from which close‐spaced cross‐sections were extracted for the 1D modeling, of three morphologically and hydraulically different river systems. These systems range from small low‐gradient meandering pool–riffle to large steep confined plane‐bed rivers. We test the effects of 1D and 2D models on predicted hydraulic variables at cross‐sections and over the full bathymetry to quantify the differences due to model dimensionality and those from interpolation. Results show that streambed features, whose size is smaller than cross‐sectional spacing, chiefly determine the different results of 1D and 2D modeling whereas flow discharge, stream size, morphological complexity and model grid sizes have secondary effects on flow properties and habitat quality for a given species and life stage predicted from 1D and 2D modeling. In general, the differences in hydraulic variables are larger in the bathymetric than in the cross‐sectional analysis, which suggests that some errors are introduced from interpolation of spatially disaggregated simulated variables with a 1D model, instead of model dimensionality 1D or 2D. Flow property differences are larger for velocity than for water surface elevation and depth. Differences in weighted usable area (WUA) derived from 1D and 2D modeling are relatively small for low‐gradient meandering pool–riffle systems, but the differences in the spatial distribution of microhabitats can be considerable although clusters of same habitat quality are spatially comparable. Copyright © 2014 John Wiley & Sons, Ltd.     

Bankfull discharge in mountain streams: theory and practice

The results are presented of an investigation of bankfull discharge in two Polish Carpathian streams: Skawica and Krzyworzeka. Existing definitions of river bankfull were reviewed and applied in tests carried out on selected cross‐sections of the streams. The Woodyer method was given special attention, with a correspondingly detailed survey of plants characterizing river benches. Riley's bench index method and the methods of Williams, Wolman, Schumm and Brown, and Woloszyn were tested. The report concludes that bankfull discharge value for a mountain stream should not be reported as a single number, but rather as a range of discharges within which one could expect the bankfull value to lie. Copyright © 2002 John Wiley & Sons, Ltd.     

Morphology of riffle–pool sequences in the River Severn,England

Despite the occurrence of riffle–pool sequences in many rivers there are few data concerning riffle–pool unit morphology. Of many criteria proposed to identify riffle–pool units, only three methods can be regarded as objective and robust. These are the ‘zero‐crossing’, the ‘spectral analysis’ and the ‘control‐point’ methods. In this paper statistics are developed using the first two of these methods to describe the streamwise morphology of 275 riffles and 285 pools which form a continuous 32·1 km reach of the bed of the River Severn in Shropshire, England. Yalin's theoretical relationship between the average riffle:pool unit length ( λ _p) and channel width ( W ), λ _p = 3 W , applies to the River Severn. Reach‐average riffle height ( H ) is a constant proportion of bankfull depth ( h ); typically H ≅ 0·16 h . Riffle height is a positive function of riffle length. Pool depth is a positive function of pool length. However, both riffle length and pool length increase more rapidly than the bed‐level amplitude, such that long riffles or pools are relatively ‘flat’. As channel gradient reduces, bedforms flatten and become more asymmetric as riffle stoss sides and the proximal slope of pools lengthen at the expense of riffle lee sides and pool distal slopes. The statistical relationships between riffle steepness (H/L) and water depth are similar to those for equilibrium subaqueous dunes. The Severn data are consistent with Yalin's theoretical analysis relating riffle bedform length (L _r) to water depth, i.e. L _r = α2π h, wherein α ≅ 1 for steep near‐equilibrium bedforms but α ≅ 2 to 3 as the relative depth decreases and riffles become long, low features. Theoretical consideration and turbulence data indicate that the frequency of coherent turbulent‐flow structures associated with the riffle–pool mixing length in the Severn should be of the order of 50 to 100 s. The morphological similarity of the steepest River Severn riffles with dunes raises intriguing questions with respect to self‐similar, convergent organization of periodic alluvial bedforms and to bedform dynamic classification particularly. Copyright © 2000 John Wiley & Sons, Ltd.     

Velocity characteristics along a small step–pool channel

This paper summarizes measurements of velocity along three reaches of a small mountain channel with step–pool bedforms. A one‐dimensional electromagnetic current meter was used to record velocity fluctuations at 37 fixed measurement points during five measurement intervals spanning the peak of the annual snowmelt hydrograph. Measurement cross‐sections were located upstream from a bed‐step, at the step lip, downstream from the step, and in a uniform‐gradient run. Data analyses focused on characteristics of velocity profiles, and on correlations between velocity characteristics and the potential control variables bedform type, reach gradient and flow depth. To test the hypothesis that velocity characteristics are related to channel bedform types, ANOVA and ANCOVA tests were performed for the average velocity and coefficient of variation of point velocity data. Results indicate that high frequency velocity variations correlate to some degree with both channel characteristics and discharge. Velocity became more variable as stage increased, particularly at low‐gradient reaches with less variable bed roughness. Velocity profiles suggest that locations immediately downstream from bed‐steps are dominated by wake turbulence from mid‐profile shear layers. Locations immediately upstream from steps, at step lips, and in runs are dominated by bed‐generated turbulence. Adverse pressure gradients upstream and downstream from steps may be enhancing turbulence generation, whereas favourable pressure gradients at steps are suppressing turbulence. The bed‐generated turbulence and skin friction of runs appear to be less effective energy dissipators than the wake‐generated turbulence and form drag of step–pool bedforms. Copyright © 2000 John Wiley & Sons, Ltd.     

The linkage between velocity patterns and sediment entrainment in a forced‐pool and riffle unit

A field‐based project was initiated in order to characterize velocities and sediment entrainment in a forced‐pool and riffle sequence. Three‐dimensional velocities and turbulence intensities were measured with an acoustic Doppler velocimeter at 222 different points at three similar flows that averaged approximately 4·35 m³ s⁻¹ within a large pool–riffle unit on North Saint Vrain Creek, Colorado. Sediment‐sorting patterns were observed with the introduction of 500 tracer particles painted according to initial seeding location. Tracer particles moved sporadically during a 113 day period in response to the annual snowmelt peak flow, which reached a maximum level of 14·8 m³ s⁻¹. Velocity data indicate high instantaneous velocities and turbulence levels in the centre of pools. Patterns of sediment deposition support the notion that stream competence is higher in the pool than the downstream riffle. Flow convergence around a large channel constriction appears to play a major role in multiple processes that include helical flow development and sediment routing, and backwater development with low velocities and turbulence levels above the constriction that may locally limit sediment supply. Jet flow, flow separation, vortex scour and turbulence generation enhance scour in the centre of pools. Ultimately, multiple processes appear to play some role in maintenance of this forced pool and the associated riffle. Copyright © 2008 John Wiley & Sons, Ltd.     

Temporal and spatial response of hyporheic zone geochemistry to a storm event

Although there has been recent focus on understanding spatial variability in hyporheic zone geochemistry across different morphological units under baseflow conditions, less attention has been paid to temporal responses of hyporheic zone geochemistry to non‐steady‐state conditions. We documented spatial and temporal variability of hyporheic zone geochemistry in response to a large‐scale storm event, Tropical Storm Irene (August 2011), across a pool–riffle–pool sequence along Chittenango Creek in Chittenango, NY, USA. We sampled stream water as well as pore water at 15 cm depth in the streambed at 14 locations across a 30 m reach. Sampling occurred seven times at daily intervals: once during baseflow conditions, once during the rising limb of the storm hydrograph, and five times during the receding limb. Principal component analysis was used to interpret temporal and spatial changes and dominant drivers in stream and pore water geochemistry (n = 111). Results show the majority of spatial variance in hyporheic geochemistry (62%) is driven by differential mixing of stream and ground water in the hyporheic zone. The second largest driver (17%) of hyporheic geochemistry was temporal dilution and enrichment of infiltrating stream water during the storm. Hyporheic sites minimally influenced by discharging groundwater (‘connected’ sites) showed temporal changes in water chemistry in response to the storm event. Connected sites within and upstream of the riffle reflected stream geochemistry throughout the storm, whereas downstream sites showed temporally lagged responses in some conservative and biogeochemically reactive solutes. This suggests temporal changes in hyporheic geochemistry at these locations reflect a combination of changes in infiltrating stream chemistry and hyporheic flowpath length and residence time. The portion of the study area strongly influenced by groundwater discharge increased in size throughout the storm, producing elevated Ca²⁺ and concentrations in the streambed, suggesting zones of localized groundwater inputs expand in response to storms. Copyright © 2013 John Wiley & Sons, Ltd.     

Hyporheic exchange flow around constructed in‐channel structures and implications for restoration design

In‐channel rock vane structures are widely used in stream restoration as a way to reduce stream channel erosion and create pool or riffle features. When these structures change hydraulic gradients they may affect ecological stream functions, such as hyporheic exchange flow (HEF) patterns. A study of constructed in‐channel structure controls on HEF was conducted in the third‐order Batavia Kill, New York using stream and hyporheic temperature amplitude analysis and computational fluid dynamics (CFD) hydraulic simulations. Temperature monitors were installed in the water column and channel bed at six locations around each of seven in‐channel restoration structures (three cross‐vanes and four J‐hooks) at baseflow in 2007. Elevation surveys of the structures were then used to simulate HEF using CFD. The results indicate a pattern of pronounced upwelling in the run section just below the structure, upwelling transitioning to downwelling within the pool, and pronounced downwelling in the glide out of the pool. This pattern is consistent with natural riffle pool sequences. The direction of HEF inferred from the temperature amplitude analysis agreed with the direction of flow simulated with CFD at 80% of the locations, and the few disagreements were expected due to model limitations. CFD simulation demonstrated that increasing stream flows result in changes in HEF spatial patterns and magnitude at each structure. This work illustrates how CFD simulations can guide design of in‐channel restoration structures for HEF function. Copyright © 2009 John Wiley & Sons, Ltd.