Three‐dimensional flow structure and patterns of bed shear stress in an evolving compound meander bend

Compound meander bends with multiple lobes of maximum curvature are common in actively evolving lowland rivers. Interaction among spatial patterns of mean flow, turbulence, bed morphology, bank failures and channel migration in compound bends is poorly understood. In this paper, acoustic Doppler current profiler (ADCP) measurements of the three‐dimensional (3D) flow velocities in a compound bend are examined to evaluate the influence of channel curvature and hydrologic variability on the structure of flow within the bend. Flow structure at various flow stages is related to changes in bed morphology over the study timeframe. Increases in local curvature within the upstream lobe of the bend reduce outer bank velocities at morphologically significant flows, creating a region that protects the bank from high momentum flow and high bed shear stresses. The dimensionless radius of curvature in the upstream lobe is one‐third less than that of the downstream lobe, with average bank erosion rates less than half of the erosion rates for the downstream lobe. Higher bank erosion rates within the downstream lobe correspond to the shift in a core of high velocity and bed shear stresses toward the outer bank as flow moves through the two lobes. These erosion patterns provide a mechanism for continued migration of the downstream lobe in the near future. Bed material size distributions within the bend correspond to spatial patterns of bed shear stress magnitudes, indicating that bed material sorting within the bend is governed by bed shear stress. Results suggest that patterns of flow, sediment entrainment, and planform evolution in compound meander bends are more complex than in simple meander bends. Moreover, interactions among local influences on the flow, such as woody debris, local topographic steering, and locally high curvature, tend to cause compound bends to evolve toward increasing planform complexity over time rather than stable configurations. Copyright © 2016 John Wiley & Sons, Ltd.     

Modulation of the flow structure by progressive bedforms in the Kinoshita meandering channel

An in‐house fully three‐dimensional general‐purpose finite element model is applied to solve the hydrodynamic structure in a periodic Kinoshita‐generated meandering channel. The numerical model solves the incompressible Reynolds‐averaged Navier–Stokes equations for mass and momentum, while solving the k ? ε equations for turbulence. The free surface is described by the rigid‐lid approximation (using measured water surface data) for flat (smooth‐bed) and self‐formed (rough‐bed) conditions. The model results are compared against experimental measurements in the ‘Kinoshita channel’, where three‐dimensional flow velocities and turbulence parameters were measured. This validation was carried out for the upstream‐valley meander bend orientation under smooth (flat bed) conditions. After validation, several simulations were carried out to predict the hydrodynamics in conditions where either it was not possible to perform measurements (e.g. applicability of the laboratory acoustic instruments) and to extrapolate the model to other planform configurations. For the flat smooth‐bed case, a symmetric (no skewness) planform configuration was modeled and compared to the upstream‐skewed case. For the self‐formed rough‐bed case, prediction of the hydrodynamics during the progression of bedforms was performed. It appears that the presence of bedforms on a bend has the following effects: (i) the natural secondary flow of the bend is disrupted by the presence of the bedforms, thus depending on the location of the dune, secondary flows might differ completely from the traditional orientation; (ii) an increment on both the bed and bank shear stresses is induced, having as much as 50% more fluvial erosion, and thus a potential increment on the migration rate of the bend. Implications on sediment transport and bend morphodynamics are also discussed in the paper. Copyright © 2013 John Wiley & Sons, Ltd.     

Chute channel dynamics in large,sand‐bed meandering rivers

Meander bends of many large, sand‐bed meandering rivers are partitioned by chute channels that convey permanent flow, and co‐exist with the mainstem for decades. As a first step toward understanding the dynamics and morphodynamic implications of these ‘bifurcate meander bends’, this study applied binary logistic regression analysis to determine whether it is possible to predict chute initiation based on attributes of meander bend character and dynamics. Regression models developed for the Strickland River, Papua New Guinea, the lower Paraguay River, Paraguay/Argentina, and the Beni River, Bolivia, revealed that the probability of chute initiation at a meander bend is a function of the bend extension rate (the rate at which a bend elongates in a direction perpendicular to the valley axis trend). Image analyses of all rivers and field observations from the Strickland suggest that the majority of chute channels form during scroll–slough development. Rapid extension is shown to favour chute initiation by breaking the continuity of point bar deposition and vegetation encroachment at the inner bank, resulting in widely‐spaced scrolls with intervening sloughs that are positively aligned with primary over‐bar flow. The rivers plot in order of increasing chute activity on an empirical meandering‐braided pattern continuum defined by potential specific stream power (ω_pv) and bedload calibre (D₅₀). Increasing stream power is considered to result in higher bend extension rates, with implications for chute initiation. In addition, chute stability is shown to depend on river sediment load relative to flow discharge (Qs/Q), such that while the Beni may plot in the region of highly braided rivers by virtue of a high potential specific stream power, the formation of stable chute channels is suppressed by the high sediment load. This tendency is consistent with previous experimental studies, and results in a planform that is transitional between single‐thread meandering and braided. Copyright © 2011 John Wiley & Sons, Ltd.     

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

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

Spatial variability of three‐dimensional Reynolds stresses in a developing channel bend

Experimental results of the mean flow field and turbulence characteristics for flow in a model channel bend with a mobile sand bed are presented. Acoustic Doppler velocimeters (ADVs) were used to measure the three components of instantaneous velocities at multiple cross sections in a 135° channel bend for two separate experiments at different stages of clear water scour conditions. With measurements at multiple cross sections through the bend it was possible to map the changes in both the spatial distribution of the mean velocity field and the three Reynolds shear stresses. Turbulent stresses are known to contribute to sediment transport and the three‐dimensionality inherent to flow in open channel bends presents a useful case for determining specific relations between three‐dimensional turbulence and sediment entrainment and transport. These measurements will also provide the necessary data for validating numerical simulations of turbulent flow and sediment transport. The results show that the magnitude and distribution of three‐dimensional Reynolds stresses increase through the bend, with streamwise‐cross stream and cross stream‐vertical components exceeding the maximum principal Reynolds stress through the bend. The most intriguing observation is that near‐bed maximum positive streamwise‐cross stream Reynolds stress coincides with the leading edge of the outer bank scour hole (or thalweg), while maximum cross stream‐vertical Reynolds stress (in combination with high negative streamwise‐cross stream Reynolds stress near the bend apex) coincides with the leading edge of the inner bank bar. Maximum Reynolds stress and average turbulent kinetic energy appear to be greater and more localized over the scour hole before final equilibrium scour is reached. This suggests that the turbulent energy in the flow is higher while the channel bed is developing, and both lower turbulent energy and a broader distribution of turbulent stresses near the bed are required for cessation of particle mobilization and transport. Copyright © 2010 John Wiley & Sons, Ltd.     

Spatial variation of flow characteristics in a subarctic meandering river in ice-covered and open-channel conditions: A 2D hydrodynamic modelling approach

To be able to understand year-round river channel evolution both at present and in the future, the spatial variation of the flow characteristics and their sediment transport capabilities under ice cover need to be detected. As the measurements done through cross-sectional drill holes cover only a small portion of the river channel area, the numerical simulations give insight into the wider spatial horizontal variation of the flow characteristics. Therefore, we simulate the ice-covered flow with a hydrodynamic two-dimensional (2D) model in a meandering subarctic river (Pulmanki River, Finland) in mid-winter conditions and compare them to the pre-winter open-channel low flow situation. Based on the simulations, which are calibrated with reference measurements, we aim to detect (1) how ice-covered mid-winter flow characteristics vary spatially and (2) the erosion and sedimentation potential of the ice-covered flow compared to open-channel conditions. The 2D hydrodynamic model replicated the observed flow characteristics in both open-channel and ice-covered conditions. During both seasons, the greatest erosional forces locate in the shallow sections. The narrow, freely flowing channel area found in mid-winter cause the main differences in the spatial flow variation between seasons. Despite the causes of the horizontal recirculating flow structures being similar in both seasons, the structures formed in different locations depended on whether the river was open or ice covered. The critical thresholds for particle entrainment are exceeded more often in open-channel conditions than during ice-covered flow. The results indicate spatially extensive sediment transport in open-channel conditions, but that the spatial variability and differences in depositional and erosional locations increase in ice-covered conditions. Asymmetrical bends and straight reaches erode throughout the year, whereas symmetrical, smaller bends mainly erode in open-channel conditions and are prone to deposition in winter. The long ice-covered season can greatly affect the annual morphology of the submerged channel. © 2019 John Wiley & Sons, Ltd.     

Influence of experimental removal of large woody debris on spatial patterns of three‐dimensional flow in a meander bend

This study reports the results of a large woody debris (LWD) removal experiment in a meander bend along a low‐energy stream in the Midwestern United States. The LWD obstacle was located in the center of the channel at the bend exit and consisted of a mature tree with an intact soil‐covered root wad and a large accumulation of logs, branches and pieces of lumber on top of and adjacent to the main tree. The results indicate that the LWD obstruction influenced 3D flow structure in this bend at all flow stages. The main effect of LWD is to dramatically decelerate flow throughout the majority of the bend, while locally accelerating flow where it passes through the narrow chute at the downstream end of the LWD obstruction. Results from the LWD removal experiment indicate that patterns of three‐dimensional flow structure in meander bends are sensitive to complete removal of LWD. After the removal of LWD from the bend, both downstream and secondary velocities increased and, though still weak, secondary flow intensified. Large, relatively stable, obstructions that span a significant portion of the channel may act as natural dams, effectively ponding water upstream of the LWD, thereby producing substantial convective deceleration of the flow. This research is the first to document three‐dimensional flow structure before and after a controlled removal of LWD from a meander bend. Studies of the type reported here represent a first step toward determining the ensemble of process interactions between LWD and bend dynamics. Copyright © 2006 John Wiley & Sons, Ltd.     

Distorted Froude‐scaled flume analysis of large woody debris

This paper presents the results of a movable‐boundary, distorted, Froude‐scaled hydraulic model based on Abiaca Creek, a sand‐bedded channel in northern Mississippi. The model was used to examine the geomorphic and hydraulic impact of simplified large woody debris (LWD) elements. The theory of physical scale models is discussed and the method used to construct the LWD test channel is developed. The channel model had bed and banks moulded from 0·8 mm sand, and flow conditions were just below the threshold of motion so that any sediment transport and channel adjustment were the result of the debris element. Dimensions and positions of LWD elements were determined using a debris jam classification model. Elements were attached to a dynamometer to measure element drag forces, and channel adjustment was determined through detailed topographic surveys. The fluid drag force on the elements decreased asymptotically over time as the channel boundary eroded around the elements due to locally increased boundary shear stress. Total time for geomorphic adjustment computed for the prototype channel at the Q₂ discharge (discharge occurring once every two years on average) was as short as 45 hours. The size, depth and position of scour holes, bank erosion and bars created by flow acceleration past the elements were found to be related to element length and position within the channel cross‐section. Morphologies created by each debris element in the model channel were comparable with similar jams observed in the prototype channel. Published in 2001 John Wiley & Sons, Ltd.     

Sensitivity of simulated flow fields and bathymetries in meandering channels to the choice of a morphodynamic model

Morphodynamic models are used by river practitioners and scientists to simulate geomorphic change in natural and artificial river channels. It has long been recognized that these models are sensitive to the choice of parameter values, and proper calibration is now common practice. This paper investigates the less recognized impact of the choice of the model itself. All morphodynamic models purport to simulate the same flow and sediment dynamics, often relying on the same governing equations. Yet in solving these equations, the models have different underlying assumptions, for example regarding spatial discretization, turbulence, sediment inflow, lateral friction, and bed load transport. These differences are not always considered by the average model user, who might expect similar predictions from calibrated models. Here, a series of numerical simulations in meandering channels was undertaken to test whether six morphodynamic codes (BASEMENT, CCHE‐2D, NAYS, SSIIM‐1, TELEMAC‐2D and TELEMAC‐3D) would yield significantly different equilibrium bathymetries if subjected to identical, initial flow conditions. We found that, despite producing moderately similar velocity patterns on a fixed‐flat bed (regression coefficient r of 0.77 ± 0.20), the codes disagree substantially with respect to simulated bathymetries (r = 0.49 ± 0.31). We relate these discrepancies to differences in the codes' assumptions. Results were configuration specific, i.e. codes that perform well for a given channel configuration do not necessarily perform well with higher or lower sinuosity configurations. Finally, limited correlation is found between accuracy and code complexity; the inclusion of algorithms that explicitly account for the effects of local bed slope and channel curvature effects on transport magnitude and direction does not guarantee accuracy. The range of solutions obtained from the evaluated codes emphasizes the need for carefully considering the choice of code. We recommend the creation of a central repository providing universal validation cases and documentation of recognized fluvial codes in commonly studied fluvial settings. Copyright © 2016 John Wiley & Sons, Ltd.     

Two‐dimensional and three‐dimensional computational models in hydrodynamic and morphodynamic reconstructions of a river bend: sensitivity and functionality

This study assesses hydrodynamic and morphodynamic model sensitivity and functionality in a curved channel. The sensitivity of a depth‐averaged model to user‐defined parameters (grain size, roughness, transverse bed slope effect, transport relations and secondary flow) is tested. According to the sensitivity analysis, grain size, transverse bed slope effect and sediment transport relations are critical to simulated meander bend morphodynamics. The parametrization of grain size has the most remarkable effect: field‐based grain size parametrization is necessary in a successful morphodynamic reconstruction of a meander bend. The roughness parametrization method affects the distribution of flow velocities and therefore also morphodynamics. The combined effect of various parameters needs further research. Two‐dimensional (2D) and three‐dimensional (3D) reconstructions of a natural meander bend during a flood event are assessed against field measurements of acoustic Doppler current profiler and multi‐temporal mobile laser scanning data. The depth‐averaged velocities are simulated satisfactorily (differences from acoustic Doppler current profiler velocities 5–14%) in both 2D and 3D simulations, but the advantage of the 3D hydrodynamic model is unquestionable because of its ability to model vertical and near‐bed flows. The measured and modelled near‐bed flow, however, differed notably from each other's, the reason of which was left open for future research. It was challenging to model flow direction beyond the apex. The 3D flow features, which also affected the distribution of the bed shear stress, seem not to have much effect on the predicted morphodynamics: the 2D and 3D morphodynamic reconstructions over the point bar resembled each other closely. Although common features between the modelled and measured morphological changes were also found, some specific changes that occurred were not evident in the simulation results. Our results show that short‐term, sub‐bend scale morphodynamic processes of a natural meander bend are challenging to model, which implies that they are affected by factors that have been neglected in the simulations. The modelling of short‐term morphodynamics in natural curved channel is a challenge that requires further study. Copyright © 2014 John Wiley & Sons, Ltd.     

Experimental study on transient and steady‐state dynamics of bedforms in supply limited configuration

The purpose of the present study is to investigate experimentally the development of bedforms in a configuration where the sediment supply is limited. The experimental setup is a rectangular closed duct combining an innovative system to control the rate of sediment supply Q_in , and a digitizing system to measure in real time the 3D bedform topography. We carried out different sets of experiments with two sediment sizes (100 µm and 500 µm) varying both the sediment supply and the water flow rate to obtain a total of 46 different configurations. After a transient phase, steady sub‐centimeter bedforms of various shapes have been observed: barchans dunes, straight transverse dune, linguoid transverse dunes and bedload sheets. Height, spacing, migration speed, and mean bed elevation of the equilibrium bedforms were measured. For a given flow rate, two regimes were identified with fine sediment: (i) a monotonic increasing regime where the equilibrium bedform height and velocity increase with the sediment supply rate Q_in and (ii) an invariant regime for which both parameters are almost independent of Q_in. For coarse sediment, only the first regime is observed. We interpret the saturation of height and velocity for fine sediment bedforms as the transition from a supply‐limited regime to a transport‐limited regime in which the bedload flux has reached its maximum value under the prevailing flow conditions. We also demonstrate that all experiments can be rescaled if the migration speed and height of the bedforms are, respectively, divided and multiplied by the cube of the shear velocity. This normalization is independent of grain size and of bedform morphology. These experimental results provide a new quantification of the factors controlling equilibrium height and migration speed of bedforms in supply‐limited conditions against which theoretical and numerical models can be tested.     

Shear velocity estimates in rough‐bed open‐channel flow

Shear velocity u_* is an important parameter in geophysical flows, in particular with respect to sediment transport dynamics. In this study, we investigate the feasibility of applying five standard methods [the logarithmic mean velocity profile, the Reynolds stress profile, the turbulent kinetic energy (TKE) profile, the wall similarity and spectral methods] that were initially developed to estimate shear velocity in smooth bed flow to turbulent flow over a loose bed of coarse gravel (D₅₀ = 1·5 cm) under sub‐threshold conditions. The analysis is based on quasi‐instantaneous three‐dimensional (3D) full depth velocity profiles with high spatial and temporal resolution that were measured with an Acoustic Doppler Velocity Profiler (ADVP) in an open channel. The results of the analysis confirm the importance of detailed velocity profile measurements for the determination of shear velocity in rough‐bed flows. Results from all methods fall into a range of ± 20% variability and no systematic trend between methods was observed. Local and temporal variation in the loose bed roughness may contribute to the variability of the logarithmic profile method results. Estimates obtained from the TKE and Reynolds stress methods reasonably agree. Most results from the wall similarity method are within 10% of those obtained by the TKE and Reynolds stress methods. The spectral method was difficult to use since the spectral energy of the vertical velocity component strongly increased with distance from the bed in the inner layer. This made the choice of the reference level problematic. Mean shear stress for all experiments follows a quadratic relationship with the mean velocity in the flow. The wall similarity method appears to be a promising tool for estimating shear velocity under rough‐bed flow conditions and in field studies where other methods may be difficult to apply. This method allows for the determination of u_* from a single point measurement at one level in the intermediate range (0·3 < h < 0·6). Copyright © 2013 John Wiley & Sons, Ltd.     

Effect of non‐uniformity of flow on velocity and turbulence intensities over a cobble‐bed

Non‐uniform flows encompassing both accelerating and decelerating flows over a cobble‐bed flume have been experimentally investigated in a flume at a scale of intermediate relative submergence. Measurements of mean longitudinal flow velocity u, and determinations of turbulence intensities u′, v′, w′, and Reynolds shear stress ?u_fw_f have been made. The longitudinal velocity distribution was divided into the inner zone close to the bed and the outer zone far from the bed. In the inner zone of the boundary layer (near the bed) the velocity profile closely followed the ‘Log Law’; however, in the outer zone the velocity distribution deviated from the Log Law consistently for both accelerating and decelerating flows and the changes in bed slopes ranging from ?2% to + 2% had no considerable effect on the outer zone. For a constant bed slope (S = ±0·015), the larger the flow rate, the smaller the turbulence intensities. However, no detectable pattern has been observed for u′, v′ and w′ distributions near the bed. Likewise, for a constant flow rate (Q = 0·040 m³/s), with variation in bed slope the longitudinal turbulent intensity profile in the longitudinal direction remained concave for both accelerating and decelerating flows; whereas vertical turbulent intensity (w′) profile presented no specific form. The results reveal that the positions of maximum values of turbulence intensities and the Reynolds shear stress depend not only on the flow structure (accelerating or decelerating) but also on the intermediate relative submergence scale. Copyright © 2009 John Wiley & Sons, Ltd.     

Predictors of urban variable source area: a cross‐sectional analysis of urbanized catchments in the United States

Many studies have empirically confirmed the relationship between urbanization and changes to the hydrologic cycle and degraded aquatic habitats. While much of the literature focuses on extent and configuration of impervious area as a causal determinant of degradation, in this article, I do not attribute causes of decreased watershed storage on impervious area a priori. Rather, adapting the concept of variable source area (VSA) and its relationship to incremental storage to the particular conditions of urbanized catchments, I develop a statistically robust linear regression‐based methodology to detect evidence of VSA‐dominant response. Using the physical and meteorological characteristics of the catchments as explanatory variables, I then use logistic regression to statistically analyze significant predictors of the VSA classification. I find that the strongest predictor of VSA‐type response is the percent of undeveloped area in the catchment. Characteristics of developed areas, including total impervious area, percent‐developed open space and the type of drainage infrastructure, do not add to the explanatory power of undeveloped land in predicting VSA‐type response. Within only developed areas, I find that total impervious area and percent‐developed open space both decrease the odds of a catchment exhibiting evidence of VSA‐type response and the effect of developed open space is more similar to that of total impervious area than undeveloped land in predicting VSA response. Different types of stormwater management infrastructure, including combined sewer systems and infiltration, retention and detention infrastructure are not found to have strong statistically significant effects on probability of VSA‐type response. VSA‐type response is also found to be stronger during the growing season than the dormant season. These findings are consistent across a national cross‐section of urbanized watersheds, a higher resolution dataset of Baltimore Metropolitan Area watersheds and a subsample of watersheds confirmed not to be served by (combined sewer systems). Copyright © 2016 John Wiley & Sons, Ltd.     

Experimental study of the ground motion on a large scale topographic hill at Kitherion (Greece)

We present the results of an experiment about the effects on the seismic ground motion of a hill of large dimensions (6km long, 3km large and 700m high) near Corinth (Greece). We installed 7 seismometers across the hill, and analyzed the ground motion with different methods: the classical spectral ratios (CSR) and the horizontal to vertical spectral ratios calculated both on noise (HVNR) and earthquake data (RF). The amplification at the resonant frequency (around 0.7Hz) is low, as expected, but one of the two stations installed at the top exhibits large amplification (up to a factor of 10) at 3Hz. The comparison between the H/V results and the CSR ones suggests that the former method is able to point out the fundamental frequencies of a hill. Amplifications and spatial localization predicted by numerical 3D modeling are consistent with the experimental data only at frequencies lower than 1Hz.     

Mean and turbulent flow fields in a simulated ice‐covered channel with a gravel bed: some laboratory observations

Northern rivers experience freeze‐up over the winter, creating asymmetric under‐ice flows. Field and laboratory measurements of under‐ice flows typically exhibit flow asymmetry and its characteristics depend on the presence of roughness elements on the ice cover underside. In this study, flume experiments of flows under a simulated ice cover are presented. Open water conditions and simulated rough ice‐covered flows are discussed. Mean flow and turbulent flow statistics were obtained from an Acoustic Doppler Velocimeter (ADV) above a gravel‐bed surface. A central region of faster flow develops in the middle portion of the flow with the addition of a rough cover. The turbulent flow characteristics are unambiguously different when simulated ice covered conditions are used. Two distinct boundary layers (near the bed and in the vicinity of the ice cover, near the water surface) are clearly identified, each being characterized by high turbulent intensity levels. Detailed profile measurements of Reynolds stresses and turbulent kinetic energy indicate that the turbulence structure is strongly influenced by the presence of an ice cover and its roughness characteristics. In general, for y/d > 0·4 (where y is height above bed and d is local flow depth), the addition of cover and its roughening tends to generate higher turbulent kinetic energy values in comparison to open water flows and Reynolds stresses become increasingly negative due to increased turbulence levels in the vicinity of the rough ice cover. The high negative Reynolds stresses not only indicate high turbulence levels created by the rough ice cover but also coherent flow structures where quadrants one and three dominate. Copyright © 2012 John Wiley & Sons, Ltd.     

Seasonal dynamics and exports of elements from a first‐order stream to a large inland lake in Michigan

Headwater streams are critical components of drainage systems, directly connecting terrestrial and downstream aquatic ecosystems. The amount of water in a stream can alter hydrologic connectivity between the stream and surrounding landscape and is ultimately an important driver of what constituents headwater streams transport. There is a shortage of studies that explore concentration–discharge (C‐Q) relationships in headwater systems, especially forested watersheds, where the hydrological and ecological processes that control the processing and export of solutes can be directly investigated. We sought to identify the temporal dynamics and spatial patterns of stream chemistry at three points along a forested headwater stream in Northern Michigan and utilize C‐Q relationships to explore transport dynamics and potential sources of solutes in the stream. Along the stream, surface flow was seasonal in the main stem, and perennial flow was spatially discontinuous for all but the lowest reaches. Spring snowmelt was the dominant hydrological event in the year with peak flows an order of magnitude larger at the mouth and upper reaches than annual mean discharge. All three C‐Q shapes (positive, negative, and flat) were observed at all locations along the stream, with a higher proportion of the analytes showing significant relationships at the mouth than at the mid or upper flumes. At the mouth, positive (flushing) C‐Q shapes were observed for dissolved organic carbon and total suspended solids, whereas negative (dilution) C‐Q shapes were observed for most cations (Na⁺, Mg²⁺, Ca²⁺) and biologically cycled anions (NO₃^?, PO₄^3?, SO₄^2?). Most analytes displayed significant C‐Q relationships at the mouth, indicating that discharge is a significant driving factor controlling stream chemistry. However, the importance of discharge appeared to decrease moving upstream to the headwaters where more localized or temporally dynamic factors may become more important controls on stream solute patterns.     

Three‐dimensional flow structure and channel change in an asymmetrical compound meander loop,Embarras River,Illinois

The planform dynamics of meandering rivers produce a complex array of meander forms, including elongated meander loops. Thus far, few studies have examined in detail the flow structure within meander loops and the relation of flow structure to patterns of planform change. This field‐based investigation examines relations between three‐dimensional fluid motion and channel change within an elongated, asymmetrical meander loop containing multiple pool–riffle structures. The downstream velocity field is characterized by a high‐velocity core that shifts slightly outward as flow moves through individual lobes of the loop. For some of the measured flows this core becomes submerged below the water surface downstream of the lobe apexes. Vectors of cross‐stream/vertical velocities indicate that skew‐induced helical motion develops within the pools near lobe apexes and decays over riffles where channel curvature is less pronounced. Maximum rates of bank retreat generally occur near lobe apexes where impingement of the flow on the outer channel bank is greatest. However, maximum rates and loci of bank retreat differ for upstream and downstream lobes of the loop, leading to increasing asymmetry of loop geometry over time—a finding consistent with experimental investigations of loop evolution. Copyright © 2003 John Wiley & Sons, Ltd.     

The effects of ice cover on flow characteristics in a subarctic meandering river

The effects of ice cover on flow characteristics in meandering rivers are still not completely understood. Here, we quantify the effects of ice cover on flow velocity, the vertical and spatial flow distribution, and helical flow structure. Comparison with open‐channel low flow conditions is performed. An acoustic doppler current profiler (ADCP) is used to measure flow from up to three meander bends, depending on the year, in a small sandy meandering subarctic river (Pulmanki River) during two consecutive ice‐covered winters (2014 and 2015). Under ice, flow velocities and discharges were predominantly slower than during the preceding autumn open‐channel conditions. Velocity distribution was almost opposite to theoretical expectations. Under ice, velocities reduced when entering deeper water downstream of the apex in each meander bend. When entering the next bend, velocities increased again together with the shallower depths. The surface velocities were predominantly greater than bottom/riverbed velocities during open‐channel flow. The situation was the opposite in ice‐covered conditions, and the maximum velocities occurred in the middle layers of the water columns. High‐velocity core (HVC) locations varied under ice between consecutive cross‐sections. Whereas in ice‐free conditions the HVC was located next to the inner bank at the upstream cross‐sections, the HVC moved towards the outer bank around the apex and again followed the thalweg in the downstream cross‐sections. Two stacked counter‐rotating helical flow cells occurred under ice around the apex of symmetric and asymmetric bends: next to the outer bank, top‐ and bottom‐layer flows were towards the opposite direction to the middle layer flow. In the following winter, no clear counter‐rotating helical flow cells occurred due to the shallower depths and frictional disturbance by the ice cover. Most probably the flow depth was a limiting factor for the ice‐covered helical flow circulation, similarly, the shallow depths hinder secondary flow in open‐channel conditions. Copyright © 2016 John Wiley & Sons, Ltd.     

Some early Russian studies of subsurface storm‐flow processes

Subsurface flow is an important component of the catchment hydrological cycle. Various mechanisms of this process and their role in storm‐flow generation attracted the attention of many researchers throughout the twentieth century. The results of most of these studies are well documented. However, similar studies conducted in the past by many Russian hydrologists have never been made available to the English speaking hydrological community. This paper attempts to fill this gap and briefly review some of these investigations, focusing on their main results. It starts with the review of the early experimental studies (after 1930s–1950s), which allow the main characteristics of subsurface storm flow to be established. This is followed by a review of the research conducted in 1960s, which resulted in some conceptualization of the subsurface flow mechanisms. The paper also draws some parallels between this, mostly unknown, subsurface flow research in Russia and the better‐known contemporary studies of this process. Copyright © 2002 John Wiley & Sons, Ltd.