Does the permeability of gravel river beds affect near‐bed hydrodynamics?

The permeability of river beds is an important control on hyporheic flow and the movement of fine sediment and solutes into and out of the bed. However, relatively little is known about the effect of bed permeability on overlying near‐bed flow dynamics, and thus on fluid advection at the sediment–water interface. This study provides the first quantification of this effect for water‐worked gravel beds. Laboratory experiments in a recirculating flume revealed that flows over permeable beds exhibit fundamental differences compared with flows over impermeable beds of the same topography. The turbulence over permeable beds is less intense, more organised and more efficient at momentum transfer because eddies are more coherent. Furthermore, turbulent kinetic energy is lower, meaning that less energy is extracted from the mean flow by this turbulence. Consequently, the double‐averaged velocity is higher and the bulk flow resistance is lower over permeable beds, and there is a difference in how momentum is conveyed from the overlying flow to the bed surface. The main implications of these results are three‐fold. First, local pressure gradients, and therefore rates of material transport, across the sediment–water interface are likely to differ between impermeable and permeable beds. Second, near‐bed and hyporheic flows are unlikely to be adequately predicted by numerical models that represent the bed as an impermeable boundary. Third, more sophisticated flow resistance models are required for coarse‐grained rivers that consider not only the bed surface but also the underlying permeable structure. Overall, our results suggest that the effects of bed permeability have critical implications for hyporheic exchange, fluvial sediment dynamics and benthic habitat availability. © 2017 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

On determining the geometric and kinematic characteristics of coherent flow structures over a gravel bed: a new approach using combined PLIF‐PIV

In gravel‐bed rivers, the microtopography of the bed is known to exert a significant influence on the generation of turbulent flow structures that owe their origin to fluid shear generated near the bed. Although field and laboratory measurements have indicated that flows over gravel beds contain a range of coherent flow structures, the kinematic and dynamic properties of these structures are still poorly understood. This paper describes a new experimental methodology to quantify simultaneously both the kinematic and dynamic characteristics of coherent flow structures based upon combined planar laser‐induced fluorescence and particle imaging velocimetry (PLIF‐PIV). The results confirm that the primary generative mechanism of coherent flow structures is at the bed, where merging hairpin vortices form around bed clasts and generate larger‐scale fluid motions that advect downstream. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of turbulence on the transport of individual particles as bedload in a gravel‐bed river

This study investigates the association between mean and turbulent flow variables and the movement of individual particles in a gravel‐bed river. The experimental design implemented in the Eaton‐North River (Québec, Canada) is based on the simultaneous observations at a high temporal resolution of both particle movements as bedload using an underwater video camera and of the streamwise and vertical flow velocity components using a vertical array of three electromagnetic current meters (ECMs). The frequency and distance of displacement of particles larger than 20 mm that were sliding or rolling on the bed were measured from a 10 minutes long film. Mean and turbulent flow properties obtained for periods without sediment transport are compared to those when particles were sliding and rolling. When particles are sliding, weak differences are present for the mean streamwise velocity and normal vertical stresses. Instantaneous Reynolds shear stresses are significantly lower for sliding events which was not expected but could be explained by the important dominance of Quadrant 3 events (inward interactions). When particles are rolling, only the vertical normal stresses show a weak difference from those observed in the absence of transport but they tend to occur when Quadrant 2 (ejections) dominate the flow field. For both sliding and rolling particles, vertical and/or streamwise fluid accelerations show high magnitude values when compared to periods without transport. For sliding particles, streamwise acceleration is mostly negative and combines most of the time with a positive vertical acceleration. For rolling particles, streamwise and vertical acceleration are predominantly of opposite sign. These results suggest that fluid acceleration or deceleration in the streamwise and vertical directions is affecting the pressure field around particles. In this study, fluid acceleration seems to play a more important role than Reynolds shear stress or normal stresses for bedload movements. Copyright © 2010 John Wiley & Sons, Ltd.     

Particle size distribution of bed materials in the sandy river bed of alluvial rivers

The particle size distribution of bed materials in the sandy river bed of alluvial rivers is important in the study of topics such as friction, river bed evolution, erosion, and siltation. It also can reflect the dependency relation between river bed sediment and flow intensity. In this paper, the critical pattern of sediment movement in the near-wall region of a sandy river bed was analyzed. According to the principle of momentum balance, the critical settling-rising condition of bed material in a sandy river bed was found to be instantaneous turbulent velocity equal to 2.7 times the sediment settling velocity in quiescent water. Based on a vertical instantaneous turbulent velocity with a Gaussian distribution, a theoretical relation for calculating the particle size distribution of bed materials in a sandy river bed without pre-known characteristic grain sizes was developed by solving a stochastic equation. The for-mula is verified using measured data, and the results show that the proposed formula was in accordance with the measured data. This study has theoretical significance and practical value for determining the bed material particle size distribution of the sandy bed of alluvial rivers.     

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.     

What is suspended sediment?

Suspended sediment is conventionally regarded as that sediment transported by a fluid that it is fine enough for turbulent eddies to outweigh settling of the particles through the fluid. Early work in the fluvial field attributed suspension to turbulence, and led to the notion of a critical threshold for maintaining sediment in suspension. However, research on both turbulence structures and the interactions between suspended sediment and bedforms in rivers has shown a more complex story and, although there appear to have been no studies of the impact of bedforms on aeolian suspended sediment concentrations, turbulent flow structures and transport rates of saltating particles have been shown to be affected. This research indicates that suspended sediment neither travels with the same velocity as the flow in which it is suspended, nor is it likely to remain in suspension in perpetuity, even under conditions of steady flow or in unsteady flow the where dimensionless critical threshold is permanently exceeded. Rather, like bedload, it travels in a series of hops, and is repeatedly deposited on the bed where it remains until it is re‐entrained. Is there, therefore, a qualitative difference between suspended and saltating sediment, or is it just a quantitative difference in the size of the jump length and the frequency of re‐entrainment? It is our contention that the distinction of suspension as a separate class of sediment transport is both arbitrary and an unhelpful anthropocentric artefact. If we recognize that sediment transport is a continuum and applies to any fluid medium rather than split into different “processes” based on arbitrary thresholds and fluids, then recognizing the continuity will enable development of an holistic approach sediment transport, and thus sediment‐transport models that are likely to be viable across a wider range of conditions than hitherto. Copyright © 2015 John Wiley & Sons, Ltd.     

Sheet flow hydrodynamics over a non-uniform sand bed channel

The current study experimentally investigates the flow characteristics and temporal variations in the sheet flow profile of a non-uniform sand bed channel. Experiments were done to explore turbulent structures in the presence of a sheet flow layer with and without seepage. The turbulent events, such as stream wise velocity, Reynolds shear stresses, and turbulence intensities were found to be increasing and vertical velocity was found decreasing with a sheet layer. The presence of a sheet layer also effects the turbulent energy production and energy dissipation. All the turbulence parameters with and without a sheet layer have also been influenced by the presence of downward seepage. The rate of sheet flow movement is increased with seepage, owing to increased turbulence with seepage. The current study used wavelet analysis on temporally lagged spatial bed elevation profiles obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuation at different spatial scales. A spatial cross correlation analysis at multiple scales, based on the wavelet coefficients, has been done on these bed elevation datasets to observe the effect of downward seepage on the dynamic behavior of sheet flow at different length scales. It is found that seepage increases average bed celerity and also increases the celerity of sheet flow of similar length scales. This increase in the celerity has been hypothesized as the increase of sheet flow movement as well as the increase in turbulent parameters with seepage, which destabilizes the bed particles resulting in a disruption in the continuous propagation pattern of the sheet flow. The increase of sheet flow celerity with seepage is confirmed from the saturation level of the wavelet power spectra of the bed elevation series. The presence of seepage also affects the non-uniformity of collective sheet material.     

Three‐dimensional flow structure around small‐scale bedforms in a simulated gravel‐bed environment

Pebble clusters are common small‐scale morphological features in gravel‐bed rivers, occupying as much as 10 per cent of the bed surface. Important links exist between the presence of pebble clusters and the development of flow structures. These links are poorly understood at the three‐dimensional level. Particularly neglected has been the effect of clusters on the lateral flow characteristics. A laboratory study was conducted using a hydraulic flume, within which simulated pebble clusters were superimposed onto a plane bed of gravel material. High‐resolution three‐dimensional flow data were collected above the bed at two different flow depths using an acoustic Doppler velocimeter. The results present evidence of the importance of lateral flow in the development of turbulent flow structure. Narrow regions of high lateral and downstream turbulence intensity exist to both sides of clusters and in a three‐dimensional separation zone in their lee. This may indicate the presence of horseshoe‐type vortical structures analogous to those identified in less hydraulically rough environments. However, it is likely that these structures are more complicated given the mutual interference of the surrounding medium. The lateral flow was also identified as a key component in the upwelling identified by other authors in the lee of pebble clusters. The results of the vertical flow analysis confirm the hypothesis that six regions with distinct vertical flow characteristics exist above clusters: flow acceleration up the stoss‐side of the cluster; recirculation behind the cluster in the wake region; vortex shedding from the pebble crest and shear layer; flow reattachment downstream of the cluster; upwelling of flow downstream of the point of reattachment; and recovery of flow. Copyright © 2001 John Wiley & Sons, Ltd.     

Gravel‐bed rivers: from particles to patterns

The papers in this special issue reflect several of the major themes and topics from the 7^th International Workshop on Gravel‐Bed Rivers. The papers focus primarily on aspects of bed material transport in gravel‐bed rivers and larger scale morpho‐dynamics. Research in gravel‐bed rivers is increasingly integrating processes over a wide range of temporal and spatial scales by combining field observation, lab experimentation, numerical modeling and theory testing in a range of river types, aided by new technological developments in particle tracking, computational modeling and high resolution spatial data. This is leading to greater understanding of the processes leading to distinctive morpho‐dynamics of river types and a more reliable basis for river management. Copyright © 2012 John Wiley & Sons, Ltd.     

Turbulent flow structures in alluvial channels with curved cross‐sections under conditions of downward seepage

Experimental investigations have been done to analyze turbulent structures in curved sand bed channels with and without seepage. Measures of turbulent statistics such as time‐averaged near‐bed velocities, Reynolds stresses, thickness of roughness sublayer and shear velocities were found to increase with application of downward seepage. Turbulent kinetic energy and Reynolds normal stresses are increased in the streamwise direction under the action of downward seepage, causing bed particles to move rapidly. Analysis of bursting events shows that the relative contributions of all events (ejections, sweeps and interactions) increase throughout the boundary layer, and the thickness of the zone of dominance of sweep events, which are responsible for the bed material movement, increases in the case of downward seepage. The increased sediment transport rate due to downward seepage deforms the cross‐sectional geometry of the channel made of erodible boundaries, which is caused by an increase in flow turbulence and an associated decrease in turbulent kinetic energy dissipation and turbulent diffusion.     

Reduced‐complexity flow routing models for sinuous single‐thread channels: intercomparison with a physically‐based shallow‐water equation model

Reduced‐complexity models of fluvial processes use simple rules that neglect much of the underlying governing physics. This approach is justified by the potential to use these models to investigate long‐term and/or fundamental river behaviour. However, little attention has been given to the validity or realism of reduced‐complexity process parameterizations, despite the fact that the assumptions inherent in these approaches may limit the potential for elucidating the behaviour of natural rivers. This study presents two new reduced‐complexity flow routing schemes developed specifically for application in single‐thread rivers. Output from both schemes is compared with that from a more sophisticated model that solves the depth‐averaged shallow water equations. This comparison provides the first demonstration of the potential for deriving realistic predictions of in‐channel flow depth, unit discharge, energy slope and unit stream power using simple flow routing schemes. It also highlights the inadequacy of modelling unit stream power, shear stress or sediment transport capacity as a function of local bed slope, as has been common practice in a number of previous reduced‐complexity models. Copyright © 2009 John Wiley & Sons, Ltd.     

River ice cover influence on sediment transportation at present and under projected hydroclimatic conditions

A large number of rivers are frozen annually, and the river ice cover has an influence on the geomorphological processes. These processes in cohesive sediment rivers are not fully understood. Therefore, this paper demonstrates the impact of river ice cover on sediment transport, i.e. turbidity, suspended sediment loads and erosion potential, compared with a river with ice‐free flow conditions. The present sediment transportation conditions during the annual cycle are analysed, and the implications of climate change on wintertime geomorphological processes are estimated. A one‐dimensional hydrodynamic model has been applied to the Kokemäenjoki River in Southwest Finland. The shear stress forces directed to the river bed are simulated with present and projected hydroclimatic conditions. The results of shear stress simulations indicate that a thermally formed smooth ice cover diminishes river bed erosion, compared with an ice‐free river with similar discharges. Based on long‐term field data, the river ice cover reduces turbidity statistically significantly. Furthermore, suspended sediment concentrations measured in ice‐free and ice‐covered river water reveal a diminishing effect of ice cover on riverine sediment load. The hydrodynamic simulations suggest that the influence of rippled ice cover on shear stress is varying. Climate change is projected to increase the winter discharges by 27–77% on average by 2070–2099. Thus, the increasing winter discharges and possible diminishing ice cover periods both increase the erosion potential of the river bed. Hence, the wintertime sediment load of the river is expected to become larger in the future. Copyright © 2015 John Wiley & Sons, Ltd.     

Quantifying the Morphological Print of Bedload Transport

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

Sensitivity of interfacial hydraulics to the microtopographic roughness of water‐lain gravels

Flow within the interfacial layer of gravel‐bed rivers is poorly understood, but this zone is important because the hydraulics here transport sediment, generate flow structures and interact with benthic organisms. We hypothesized that different gravel‐bed microtopographies generate measurable differences in hydraulic characteristics within the interfacial layer. This was tested using a high density of spatially and vertically distributed, velocity time series measured in the interfacial layers above three surfaces of contrasting microtopography. These surfaces had natural water‐worked textures, captured in the field using a casting procedure. Analysis was repeated for three discharges, with Reynolds numbers between 165000 and 287000, to evaluate whether discharge affected the impact of microtopography on interfacial flows. Relative submergence varied over a small range (3.5 to 8.1) characteristic of upland gravel‐bed rivers. Between‐surface differences in the median and variance of several time‐averaged and turbulent flow parameters were tested using non‐parametric statistics. Across all discharges, microtopographic differences did not affect spatially averaged (median) values of streamwise velocity, but were associated with significant differences in its spatial variance, and did affect spatially averaged (median) turbulent kinetic energy. Sweep and ejection events dominated the interfacial region above all surfaces at all flows, but there was a microtopographic effect, with Q2 and Q4 events less dominant and structures less persistent above the surface with the widest relief distribution, especially at the highest Reynolds number flow. Results are broadly consistent with earlier work, although this analysis is unique because of the focus on interfacial hydraulics, spatially averaged ‘patch scale’ metrics and a statistical approach to data analysis. An important implication is that observable differences in microtopography do not necessarily produce differences in interfacial hydraulics. An important observation is that appropriate roughness parameterizations for gravel‐bed rivers remain elusive, partly because the relative contributions to flow resistance of different aspects of bed microtopography are poorly constrained. © 2014 The Authors. Earth Surface Processes and Landforms Published by John Wiley & Sons Ltd.     

Effect of seepage on flow and bedforms dynamics

This study, using an experimental approach, focuses on the effect of downward seepage on a threshold alluvial channel morphology and corresponding turbulent flow characteristics. In all the experiments, we observed that the streamwise time‐averaged velocities and Reynolds shear stresses were increased under the influence of downward seepage. Scales of eddy length and eddy turnover time were significantly increased with the application of downward seepage, leading to sediment transport and initiation of bedforms along the channel length. As the amount of seepage discharge increased, eddy length and turnover time were further increased, causing the development of larger bedforms. It was revealed that the geometry of bedforms was linked with the size of eddies. In this work, statistics of bedform dynamics are presented in terms of multi‐scalar bedforms in the presence of seepage. These multi‐scalar ubiquitous bedforms cast a potential impact on flow turbulence as well as stream bed morphology in channels. We used wavelet to analyse temporally lagged spatial bed elevation profiles that were obtained from a set of laboratory experiments and synchronized the wavelet coefficients with bed elevation fluctuations at different length scales. A spatial cross‐correlation analysis, based on the wavelet coefficients, was performed on these bed elevation datasets to observe the effect of downward seepage on the dynamic behaviour of bedforms at different length scales. It was found that celerity of bedforms reduced with increase in seepage percentage. Bedform celerity was best approximated by a probability density function such as Rayleigh distribution under varying downward seepage. Further, statistical analysis of physical parameters of bedforms ascertained that the reduction in bedform celerity was a result of increased bedform size. Copyright © 2017 John Wiley & Sons, Ltd.     

An experimental study on the ripple–dune transition

Flume experiments were conducted on different bed stages across the ripple–dune transition. As flow velocity increases, an initially flat bed surface (made of fairly uniform sandy material) is gradually transformed into a two‐dimensional rippled bed. With further increase in velocity, two‐dimensional ripples are replaced by irregular, linguoid ripples. As the average velocity necessary for the ripple–dune transition to occur is imposed on the bed surface, these non‐equilibrium linguoid ripples are further transformed into larger, two‐dimensional dunes. For each of these stages across the transition, a concrete mould of the bed was created and the flow structure above each fixed bed surface investigated. An acoustic Doppler velocimeter was used to study the flow characteristics above each bed surface. Detailed profiles were used along a transect located in the middle of the channel. Results are presented in the form of spatially averaged profiles of various flow characteristics and of contour maps of flow fields (section view). They clearly illustrate some important distinctions in the flow structure above the different bedform types associated with different stages during the transition. Turbulence intensity and Reynolds stresses gradually increase throughout the transition. Two‐dimensional ripples present a fairly uniform spatial distribution of turbulent flow characteristics above the bed. Linguoid ripples induce three‐dimensional turbulence structure at greater heights above the bed surface and turbulence intensity tends to increase steadily with height above bed surface in the wake region. A very significant increase in turbulence intensity and momentum exchange occurs during the transition from linguoid ripples to dunes. The turbulent flow field properties above dunes are highly dependent on the position along and above the bed surface and these fields present a very high degree of spatial variability (when compared with the rippled beds). Further investigations under natural conditions emphasizing sediment transport mechanisms and rates during the transition should represent the next step of analysis, together with an emphasis on quadrant analysis. Copyright © 2001 John Wiley & Sons, Ltd.     

Fluvial dunes: initiation,characterization, flow structure

Research into fluvial dunes spans disciplines, studies at grain to reach scales, and methodological approaches that include theoretical, experimental, numerical and field investigations. Despite significant research efforts to date, it remains difficult to provide definitive answers to some fundamental questions regarding dunes. This paper reviews three notable challenges that remain regarding fluvial dunes, namely scale‐consistent linking of bed morphologies with turbulent flow fields, the intriguing question of what causes trains of highly‐ordered sediment waves to form in beds of river sediments, and how to define the important characteristics of a dune‐covered bed, including lengths, shapes, and their statistical nature. In each case, the particular challenge is discussed and then recent research and ways forward are presented. Copyright © 2010 John Wiley & Sons, Ltd.     

Effects of submerged flexible vegetation and solid structure bars on channel bed scour

The effects of different submerged obstacle longitudinal bars with different arrangement densities on the flow profile and morphology of a scour hole were investigated under clear water conditions. Acoustic Doppler velocimetry(ADV) data were applied to plot the vertical distributions of three-dimensional velocities and turbulent contours.The experimental results indicate that arrangement density(also can represent porosity),structural material(flexible or solid),and the sidewall effect are the main factors impacting turbulent kinetic energy and the morphology of scour holes.For flexible vegetation,the maximum turbulent kinetic energy near the bed surface increased with the arrangement density.For the same structure,the depth and the magnitude of the lateral expansion of the scour hole also increased with the arrangement density.The flexible vegetation reduced the depth of the scour hole because of deflection and arrangement density.The larger volumes of scour found in the upstream and middle sections of solid structures compare well to those in flexible vegetation.The deflection of porous flexible vegetation transported the turbulent kinetic energy downstream,reduced the turbulent kinetic energy near the sediment bed,and increased the stability of the bars.Flexible vegetation bars are able to protect the bank and the bed of a river under normal conditions,making them a good alternative design in the management and restoration of rivers.     

Detection and reconstruction of large‐scale coherent flow structures in gravel‐bed rivers

Large‐scale flow structures (LSFS) in the streamwise direction are important features of gravel‐bed river flows, because they may contribute to sediment transport and gas exchange. In the present study, these structures are detected using Huang's empirical mode decomposition and reconstructed with phase‐averaging techniques based on a Hilbert transform of the velocity signal. The analysis is based on the fluctuating component of 15 quasi‐instantaneous velocity profiles measured with a three‐dimensional (3D) acoustic Doppler velocity profiler (ADVP) in an armoured gravel‐bed river with a low relative submergence of 2.9 (ratio between flow depth and bed grain diameter). LSFS were identified in most of the measured profiles and consistently showed similar features. We were able to characterize the geometry of these large‐scale coherent structures: the front has a vertical linear shift in the time domain and a vertical profile corresponding to a first quarter moon with the apex situated at z/h ≈ 0.4. In the vertical, the front scales with flow depth h, and in the streamwise direction, LSFS scale with three to seven times the mean flow depth. On the bed, the effect of LSFS is a periodic non‐linear variation of the friction velocity on average between 0.90 and 1.10 times the mean value. A model for the friction velocity cycle resulting from LSFS oscillation is presented. Copyright © 2014 John Wiley & Sons, Ltd.     

Variability of internally generated turbulence in an estuary,from 100 days of continuous observations

We present detailed observations of internally generated turbulence in a sheared, stratified natural flow, as well as an analysis of the external factors leading to its generation and temporal variability. Multi-month time series of vertical profiles of velocity, acoustic backscatter (0.5 Hz), and turbulence parameters were collected with two moored acoustic Doppler current profilers (ADCPs) in the Hudson River estuary, and estuary-long transects of water density were collected 30 times. ADCP backscatter is used for visualization of coherent turbulent structures and evaluation of surface wave biases to the turbulence measurements. Benefits of the continuous long-term turbulence record include our capturing: (1) the seasonality of turbulence due to changing riverflow, (2) hysteresis in stratification and turbulence over the fortnightly cycle of tidal range, and (3) intermittent events such as breaking internal waves. Internal mixing layers (IMLs) are defined as turbulent regions above the logarithmic velocity layer, and the bottom boundary layer (BBL) is defined as the continuously turbulent range of heights above the bed. A cross-correlation analysis reveals how IML and BBL turbulence vary with stratification and external forcing from tidal range, river flow, and winds. Turbulence in both layers is maximal at spring tide and minimal when most stratified, with one exception—IML turbulence at a site with changing channel depth and width is maximal at times of maximum stratification and freshwater input.