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.     

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.     

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.     

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.     

Turbulent flow structure in a gravel-bed river: Markov chain analysis of the fluctuating velocity profile

Three electromagnetic current meter probes were deployed in a Canadian gravel-bed river to obtain simultaneous records at 10 Hz of streamwise (u) and vertical (v) velocity components at three heights above the bed. By looking at the positive and negative signs of the instantaneous fluctuations from the time-average values of each velocity component at each height, the fluctuating velocity profile of u or v can be treated as a Markov chain with eight states and its statistical properties can be tested against null hypotheses based on the absence of spatial structure. We report results of this novel approach. The most common states of the u profile were those with either higher-than-average or lower-than-average velocities at all heights; these ‘high speed’ and ‘low speed’ states persisted for up to 3 s. The most common v profiles were all-upwards or all-downwards, but these persisted for shorter times than the high speed and low speed u profiles. Analysis of transition probabilities shows statistically significant tendencies for acceleration from the low speed u profile, and change from all-upwards to all-downwards v profile, to take place progressively from the uppermost probe downwards, in a sweep-like way. Deceleration from the high speed to low speed u profile and change from all-downwards to all-upwards v profile (burst-like behaviour) do not show such clear patterns. The results are interpreted in terms of the advection of inverted wedges of relatively high-momentum fluid, followed by more chaotic structures. A separate set of flow visualization experiments over a mixed gravel bed in a flume supports the presence of advected wedge structures, the decelerating part of the sequence corresponding to irregular ejections of near-bed fluid.     

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.     

Depth‐dependent hydraulic conductivity distribution patterns of a streambed

Characterization of streambed hydraulic conductivity from the channel surface to a great depth below the channel surface can provide needed information for the determination of stream‐aquifer hydrologic connectedness, and it is also important to river restoration. However, knowledge on the streambed hydraulic conductivity for sediments 1 m below the channel surface is scarce. This study describes a method that was used to determine the distribution patterns of streambed hydraulic conductivity for sediments from channel surface to a depth of 15 m below. The method includes Geoprobe's direct‐push techniques and Permeameter tests. Direct‐push techniques were used to generate the electrical conductivity (EC) logs and to collect sequences of continuous sediment cores from river channels, as well as from the alluvial aquifer connected to the river. Permeameter tests on these sediment cores give the profiles of vertical hydraulic conductivity (K_v) of the channel sediments and the aquifer materials. This method was applied to produce K_v profiles for a streambed and an alluvial aquifer in the Platte River Valley of Nebraska, USA. Comparison and statistical analysis of the K_v profiles from the river channel and from the proximate alluvial aquifer indicates a special pattern of K_v in the channel sediments. This depth‐dependent pattern of K_v distribution for the channel sediments is considered to be produced by hyporheic processes. This K_v‐distribution pattern implied that the effect of hyporheic processes on streambed hydraulic conductivity can reach the sediments about 9 m below the channel surface. Copyright © 2010 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.     

Sediment transport in high‐speed flows over a fixed bed: 2. Particle impacts and abrasion prediction

Single bed load particle impacts were experimentally investigated in supercritical open channel flow over a fixed planar bed of low relative roughness height simulating high‐gradient non‐alluvial mountain streams as well as hydraulic structures. Particle impact characteristics (impact velocity, impact angle, Stokes number, restitution and dynamic friction coefficients) were determined for a wide range of hydraulic parameters and particle properties. Particle impact velocity scaled with the particle velocity, and the vertical particle impact velocity increased with excess transport stage. Particle impact and rebound angles were low and decreased with transport stage. Analysis of the particle impacts with the bed revealed almost no viscous damping effects with high normal restitution coefficients exceeding unity. The normal and resultant Stokes numbers were high and above critical thresholds for viscous damping. These results are attributed to the coherent turbulent structures near the wall region, i.e. bursting motion with ejection and sweep events responsible for turbulence generation and particle transport. The tangential restitution coefficients were slightly below unity and the dynamic friction coefficients were lower than for alluvial bed data, revealing that only a small amount of horizontal energy was transferred to the bed. The abrasion prediction model formed by Sklar and Dietrich in 2004 was revised based on the new equations on vertical impact velocity and hop length covering various bed configurations. The abrasion coefficient k_v was found to be vary around k_v ~ 10⁵ for hard materials (tensile strength f_t > 1 MPa), one order of magnitude lower than the value assumed so far for Sklar and Dietrich's model. Copyright © 2017 John Wiley & Sons, Ltd.     

Field studies on factors affecting very fine and ultra fine particulate organic matter deposition in low‐gradient sand‐bed streams

The knowledge on particle deposition in streams is mainly based on investigations in mountain streams. No data exist from low‐gradient sand‐bed streams that largely differ in the morphological and hydraulic factors proposed to affect deposition. To identify physical control on particle deposition in low‐gradient streams, we assessed deposition of very fine and ultra fine organic particulate matter in 18 sand‐bed stream reaches. We added particles derived from lake sediment and assessed the mean transport distance S_P and the deposition velocity v_dep. Additionally, reach hydraulics were estimated by injections of a conservative solute tracer (NaCl). Among the low‐gradient streams, particle deposition kinetics were variable but similar to deposition in mountain streams. S_P was solely related to the flow velocity. This relation was confirmed when comprising published data on deposition of fine organic particles. An association between particle deposition and transient storage factors was insignificant. We found significance of the transient storage to S_P only for repeated measures within a single reach, when flow velocity and benthic conditions were nearly constant. Measured v_dep/v_fall ratios were much larger than unity in most reaches. Evidence from this relation suggests that the vertical transport of very fine and ultra fine organic particulate matter through the water column was caused mainly by vertical mixing. Copyright © 2006 John Wiley & Sons, Ltd.     

Velocity distribution in a gradually accelerating free surface flow

This study investigates the interaction of the vertical velocity v and the streamwise velocity u in a gradually accelerating flow. The analytical result shows that the momentum of uv driven by the mean velocities in a non-uniform flow is not negligible. This additional momentum directly results in the concave profiles of Reynolds shear stress in gradually accelerating flows, a departure from the expected linear profile. Consequently, this momentum causes the maximum velocity to be located below the free surface, i.e., the dip-phenomenon. This paper investigated the interactions of the Reynolds shear stress, non-zero vertical velocity and dip-phenomenon, it is found that the non-zero vertical velocity causes the dip-phenomenon. The approach is tested using the experimental data of Song and others, and good agreements between the predicted and measured velocity profiles have been achieved.     

Ensemble modelling of ice cover for a reservoir affected by pumped‐storage operation and climate change

Ensemble modelling was used to assess the robustness of projected impacts of pumped‐storage (PS) operation and climate change on reservoir ice cover. To this end, three one‐dimensional and a two‐dimensional laterally averaged hydrodynamic model were set up. For the latter, the strength of the impacts with increasing distance from the dam was also investigated. Climate change effects were simulated by forcing the models with 150 years of synthetic meteorological time series created with a weather generator based on available air temperature scenarios for Switzerland. Future climate by the end of the 21^st century was projected to shorten the ice‐covered period by ~2 months and decrease ice thicknesses by ~13 cm. Under current climate conditions, the ice cover would already be affected by extended PS operation. For example, the average probability of ice coverage on a specific day was projected to decrease by ~13% for current climate and could further be reduced from ~45% to ~10% for future climate. Overall, the results of all models were consistent. Although the number of winters without ice cover was projected to increase for all one‐dimensional models, studying individual segments of the two‐dimensional model showed that the impact was pronounced for segments close to the PS intake/outlet. In summary, the reservoir's ice cover is expected to partially vanish with higher probability of open water conditions closer to the PS intake/outlet.     

Estimating the mean speed of laminar overland flow using dye injection‐uncertainty on rough surfaces

A common method for estimating mean flow speeds in studies of surface runoff is to time the travel of a dye cloud across a measured flow path. Motion of the dye front reflects the surface flow speed, and a correction must be employed to derive a value for the profile mean speed, which is always lower. Whilst laminar flow conditions are widespread in the interrill zone, few data are available with which to establish the relationship linking surface and profile mean speeds, and there are virtually none for the flow range 100 < Re < 500 (Re = Reynolds number) which is studied here. In laboratory experiments on a glued sand board, mean flow speeds were estimated from both dye speeds and the volumetric flow relation v = Q/ wd with d measured using a computer‐controlled needle gauge at 64 points. In order to simulate conditions applicable to many dryland soils, the board was also roughened with plant litter and with ceramic tiles (to simulate surface stone cover). Results demonstrate that in the range 100 < Re < 500, there is no consistent relation between surface flow speeds and the profile mean. The mean relationship is v = 0·56 v _surf, which departs significantly from the theoretical smooth‐surface relation v = 0·67 v _surf, and exhibits a considerable scatter of values that show a dependence on flow depth. Given the inapplicability of any fixed conversion factor, and the dependence on flow depth, it is suggested that the use of dye timing as a method for estimating v be abandoned in favour of precision depth measurement and the use of the relation v = Q/ wd , at least within the laminar flow range tested. Copyright © 2001 John Wiley & Sons, Ltd.     

Three-dimensional turbulent structures at a medium-sized confluence with and without an ice cover

River confluences are characterized by a complex mixing zone with three-dimensional (3D) turbulent structures which have been described as both streamwise-oriented structures and Kelvin–Helmholtz (KH) vertical-oriented structures. The latter are visible where there is a turbidity difference between the two tributaries, whereas the former are usually derived from mean velocity measurements or numerical simulations. Few field studies recorded turbulent velocity fluctuations at high frequency to investigate these structures, particularly at medium-sized confluences where logistical constraints make it difficult to use devices such as acoustic doppler velocimeter (ADV). This study uses the ice cover present at the confluence of the Mitis and Neigette Rivers in Quebec (Canada) to obtain long-duration, fixed measurements along the mixing zone. The confluence is also characterized by a marked turbidity difference which allows to investigate the mixing zone dynamics from drone imagery during ice-free conditions. The aim of the study is to characterize and compare the flow structure in the mixing zone at a medium-sized (~40 m) river confluence with and without an ice cover. Detailed 3D turbulent velocity measurements were taken under the ice along the mixing plane with an ADV through eight holes at around 20 positions on the vertical. For ice-free conditions, drone imagery results indicate that large (KH) coherent structures are present, occupying up to 50% of the width of the parent channel. During winter, the ice cover affects velocity profiles by moving the highest velocities towards the centre of the profiles. Large turbulent structures are visible in both the streamwise and lateral velocity components. The strong correlation between these velocity components indicates that KH vortices are the dominating coherent structures in the mixing zone. A spatio-temporal conceptual model is presented to illustrate the main differences on the 3D flow structure at the river confluence with and without the ice cover. © 2019 John Wiley & Sons, Ltd.     

Bedload transport resistance in rough open‐channel flows

During bedload movement by saltation, streamwise momentum is transferred from the ?ow to the saltating grains. When the grains collide with other grains on the bed or in the ?ow, streamwise momentum is reduced, and there is a decrease in streamwise ?ow velocity and an increase in ?ow resistance, herein termed bedload transport resistance f_bt. Based on experiments in two ?umes with ?xed and mobile plane beds and previously published data, an equation is developed that may be used to predict f_bt for both capacity and non‐capacity ?ows. The variables in this equation are identi?ed by dimensional analysis and the coef?cients are determined by non‐linear regression. This equation applies to rough turbulent open‐channel ?ows, where the relative submergence is between 1 and 20 and the entire sediment load moves by saltation. An investigation of the relative magnitudes of f_bt and grain resistance f_c suggests that where dimensionless shear stress θ is less than 1 and saltation is the dominant mode of bedload transport, f_bt/f_c increases with θ but never exceeds 1. Copyright © 2004 John Wiley & Sons, Ltd.     

Characteristics of turbulence over the semi-fixed desert area north of Xinjiang,China

As the largest fixed and semi-fixed desert in China, the Gurbantünggüt Desert undergoes a long period of snow cover in the winter and the rapid growth of ephemeral plants in the spring, presenting obvious seasonal changes in the underlying desert surface type, which can lead to variation in the turbulence of the near-surface boundary layer turbulence over the desert. In this study, gradient tower data and eddy covariance data from 2017 were analysed to investigate the turbulence characteristics of the different surface boundary layers in the hinterland of the Gurbantünggüt Desert. The results indicate that stable atmospheric conditions in the desert occur exclusively during the early morning and at night in the desert, and the onset and duration of this stable state varies seasonally. Two regimes of intermittent turbulence occur during the night, a weak turbulent regime that occurs when the wind speed is less than the threshold and a strong turbulent regime when the wind speed exceeds the threshold, and different wind speed thresholds were observed at each level. These parameters follow a seasonal pattern of summer (July) > spring (April) > autumn (October) > winter (January) in terms of magnitude. The mean turbulence intensities of the along-wind, cross-wind and vertical wind are 0.5, 0.47 and 0.14, respectively, with I_u > I_v > I_w. The normalized standard deviation of the wind velocity components (σ_u, σ_v and σ_w) generally satisfies a 1/3 power-law relation. Our results show that the night-time turbulence regime classification for the Gurbantünggüt Desert strongly depends on meteorological and orographic features, and the intermittent turbulent events have the non-stationarity of the flow in common. The results can contribute to the study of land surface processes, climate change and desertification in inland arid desert areas.     

Scaling and self‐similarity in one‐dimensional unsteady open channel flow

The conditions under which the Saint Venant equations system for unsteady open channel flow, as an initial–boundary value problem, becomes self‐similar are investigated by utilizing one‐parameter Lie group of point scaling transformations. One of the advantages of this methodology is that the self‐similarity conditions due to the initial and boundary conditions can also be investigated thoroughly in addition to the conditions due to the governing equation. The obtained self‐similarity conditions are compared with the scaling relationships that are derived through the Froude similitude. It is shown that the initial–boundary value problem of a one‐dimensional unsteady open channel flow process in a prototype domain can be self‐similar with that of several different scaled domains. However, the values of all the flow variables (at specified time and space) under different scaled domains can be upscaled to the same values in the prototype domain (at the corresponding time and space), as shown in this study. Distortion in scales of different space dimensions has been implemented extensively in physical hydraulic modelling, mainly because of cost, space and time limitations. Unlike the traditional approach, the distinction is made between the longitudinal–horizontal and transverse–horizontal length scales in this study. The scaled domain obtained by the proposed approach, when scaling ratios of channel width and water depth are equal, is particularly important for the similarity of flow characteristics in a cross‐section because the width‐to‐depth ratio and the inclination angles of the banks are conserved in a cross‐section. It is also shown that the scaling ratio of the roughness coefficient under distorted channel conditions depends on that of hydraulic radius and longitudinal length. The proposed scaling relations obtained by the Lie group scaling approach may provide additional spatial, temporal and economical flexibility in setting up physical hydraulic models. Copyright © 2013 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.     

Pressure‐driven,shoreline currents in a perennially ice‐covered,pro‐glacial lake in Antarctica,identified from a LiCl tracer injected into a pro‐glacial stream

The distribution of streamwater within ice‐covered lakes influences sub‐ice currents, biological activity and shoreline morphology. Perennially ice‐covered lakes in the McMurdo Dry Valleys, Antarctica, provide an excellent natural laboratory to study hydrologic–limnologic interactions under ice cover. For a 2 h period on 17 December 2012, we injected a lithium chloride tracer into Andersen Creek, a pro‐glacial stream flowing into Lake Hoare. Over 4 h, we collected 182 water samples from five stream sites and 15 ice boreholes. Geochemical data showed that interflow travelled West of the stream mouth along the shoreline and did not flow towards the lake interior. The chemistry of water from Andersen Creek was similar to the chemistry of water below shoreline ice. Additional evidence for Westward flow included the morphology of channels on the ice surface, the orientation of ripple marks in lake sediments at the stream mouth and equivalent temperatures between Andersen Creek and water below shoreline ice. Streamwater deflected to the right of the mouth of the stream, in the opposite direction predicted by the Coriolis force. Deflection of interflow was probably caused by the diurnal addition of glacial runoff and stream discharge to the Eastern edge of the lake, which created a strong pressure gradient sloping to the West. This flow directed stream momentum away from the lake interior, minimizing the impact of stream momentum on sub‐ice currents. It also transported dissolved nutrients and suspended sediments to the shoreline region instead of the lake interior, potentially affecting biological productivity and bedform development. Copyright © 2014 John Wiley & Sons, Ltd.     

Modeling of terracette‐hillslope soil moisture as a function of aspect,slope and vegetation in a semi‐arid environment

In the semi‐arid western United States, water availability plays a defining role in land use. Soil moisture, vegetation, and microtopography are key variables in the hydrologic function of these ecosystems. Previous research has not addressed the influence of site‐specific aspect, vegetation, or slope gradient on terracette soil moisture patterns in semi‐arid rangelands. Therefore, the objectives of this study were to: (1) assess the influence of terracette site aspect, vegetation cover, and slope on soil moisture; (2) conceptualize conditions at the hillslope scale given terracette morphology; and (3) estimate the extent of terracettes at a regional scale. The Simultaneous Heat and Water (SHAW) model was used to simulate soil water dynamics of terracettes given variations in site conditions. These results were coupled with time‐of‐flight laser scans to quantify terracette bench and riser percent‐area, and statewide assessments of terracette extent using digital orthoimagery and a geographical information system (GIS). Modeling results indicated site aspect had minimal influence (±0.005 m³ m^?3) on terracette soil moisture. Vegetation, represented as leaf area index (LAI), had the single‐most influential effect on terracette volumetric water content (θ _v) demonstrated by an inverse relationship of LAI to mean terracette hillslope θ _v; and slope increases of ≥15% on northern azimuths increased mean θ _v which contrasted with southern azimuths for similar slope increases. Laser scanning results indicated bench width and riser length could be estimated from mean site slope (R ² = 0.82 risers and R ² = 0.93 benches). Aerial orthoimagery/GIS assessments estimated >159 000 ha of terracettes throughout the State of Idaho, with >41 000 ha (~26%) occurring on lands managed as grazing allotments. These findings provide an increased understanding of rangeland hydrologic processes as influenced by cattle density, vegetation, and terracettes which can aide land managers in their selection and application of best management practices on these lands. Copyright © 2017 John Wiley & Sons, Ltd.