Burst-like sediment suspension events in a sand bed river

In a much quoted paper, Jackson (1976) hypothesized that turbulent [bursting] motions such as those documented in laboratory boundary layers play a major role in alluvial sediment suspension. To date, the hypothesis remained largely untested, due to difficulties in monitoring turbulent suspension in rivers. This study provides field data documenting burst-like turbulent motions over a sandy bed channel and quantifying the role of these motions in sand suspension. The data were collected in a 10 m deep channel of the Fraser River near Mission, British Columbia, Canada. Turbulent fluctuations of both flow components, downstream and normal to the bed, along with the output of an optical suspended sediment sensor, were monitored 1 m above the river bed. Typical flow velocities averaged 0·9 ms⁻¹ at the sensors, where mean suspended sediment concentrations were 500 mgl⁻¹; decimetre height small dunes on the backs of larger, half-metre amplitude dunes covered the channel bed in the area. Brief but intense, burst-like [ejection and inrush] events were identified in the flow records, where they are responsible for a high degree of [intermittency] in shear stress over the dunes: 80 per cent of the turbulent momentum exchange across the 1 m level can be ascribed to such brief (3-8 s duration) events, active under 12 per cent of the time. In addition, the record of fluctuating sediment concentrations reveals these burst-like motions to be highly effective in vertically mixing suspended sediment and thus, ultimately, in maintaining suspended sediment transport above the dune bed. The bulk (60 and 90 per cent in two deployments) of the vertical sediment mixing was accomplished by intense events active some 10 per cent of the time. No discrete recurrence timescale for these ‘burst-like’ mixing events is evident, however. Rather, a continuous variation of return periods was observed as a function of the magnitude of vertical mixing event considered. To that extent, conceptual models of sediment transport in terms of burst events with a predictable recurrence such as proposed by Jackson (1976) may be misleading.     

Turbulent unsteady flow profiles over an adverse slope

When an unsteady free surface flow encounters an adverse slope, it results in a decelerating flow up the adverse slope. The time dependent turbulent flow is treated here by appropriately reducing the two-dimensional Reynolds averaged Navier-Stokes equation along with the equation of continuity considering turbulence closure. With suitable choice of parameters, the resulting differential equations are numerically solved to compute free surface and streamwise velocity profiles with time. It is found that initially the advancing free surface is convex upwards for a short time, followed by a jump of the free surface with a negative streamwise velocity that is a backwater rolling breaker due to deceleration of flow. At later time, however, the velocity becomes positive, that is, the breakers roll forward. This dual feature of motion, that is a surge followed by rolling breakers, is repeated for sometime before the jumps stop. The theoretical analysis presented here is motivated by tidal bores propagating upstream in an estuarine river.     

Pyroclastic flows from the 1991 eruption of Unzen volcano,Japan

Pyroclastic flows from Unzen were generated by gravitational collapse of the growing lava dome. As soon as the parental lobe failed at the edge of the dome, spontaneous shattering of lava occurred and induced a gravity flow of blocks and finer debris. The flows had a overhanging, tongue-like head and cone- or rollershaped vortices expanding outward and upward. Most of the flows traveled from 1 to 3 km, but some flows reached more than 4 km, burning houses and killing people in the evacuated zone of Kita-kamikoba on the eastern foot of the volcano. The velocities of the flows ranged from 15 to 25 m/s on the gentle middle flank. Observations of the flows and their deposits suggest that they consisted of a dense basal avalanche and an overlying turbulent ash cloud. The basal avalanche swept down a topographic low and formed to tongue-like lobe having well-defined levees; it is presumed to have moved as a non-Newtonian fluid. The measured velocities and runout distances of the flows can be matched to a Bingham model for the basal avalanche by the addition of turbulent resistance. The rheologic model parameters for the 29 May flow are as follows: the density is 1300 kg/m³, the yield strength is 850 Pa, the viscosity is 90 Pa s, and the thickness of the avalanche is 2 m. The ash cloud is interpreted as a turbulent mixing layer above the basal avalanche. The buoyant portions of the cloud produced ash-fall deposits, whereas the dense portions moved as a surge separated from the parental avalanche. The ash-cloud surges formed a wide devastated zone covered by very thin debris. The initial velocities of the 3 June surges, when they detached from avalanches, are determined by the runout distance and the angle of the energy-line slope. A comparison between the estimated velocities of the 3 June avalanches and the surges indicates that the surges that extended steep slopes along the avalanche path, detached directly from the turbulent heads of the avalanches. The over-running surge that reached Kita-Kamikoba had an estimated velocity higher than that of the avalanche; this farther-travelled surge is presumed to have been generated by collapse of a rising ash-cloud plume.     

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.     

Fluid flow and sediment entrainment in the Garonne River bore and tidal bore collision

A detailed field study was carried out on a tidal bore to document the turbulent processes and sediment entrainment which occurred. The measured bore, within the Arcins Channel of the Garonne River (France), was undular in nature and was followed by well‐defined secondary wave motion. Due to the local river geometry a collision between the Arcins channel tidal bore and the bore which formed within the main Garonne River channel was observed about 800 m upstream of the sampling site. This bore collision generated a transient standing wave with a black water mixing zone. Following this collision the bore from the main Garonne River channel propagated ‘backward’ to the downstream end of the Arcins channel. Velocity measurements with a fine temporal resolution were complemented by measurements of the sediment concentration and river level. The instantaneous velocity data indicated large and rapid fluctuations of all velocity components during the tidal bore. Large Reynolds shear stresses were observed during and after the tidal bore passage, including during the 'backward' bore propagation. Large suspended sediment concentration estimates were recorded and the suspended sediment flux data showed some substantial sediment motion, consistent with the murky appearance of the flood tide waters. Copyright © 2015 John Wiley & Sons, Ltd.     

Experimental investigation on local shear stress and turbulence intensities over a rough non-uniform bed with and without sediment using 2D particle image velocimetry

The current work focuses on locally resolving velocities,turbulence,and shear stresses over a rough bed with locally non-uniform character.A nonporous subsurface layer and fixed interfacial sublayer of gravel and sand were water-worked to a nature-like bed form and additionally sealed in a hydraulic flume.Two-dimensional Particle Image Velocimetry(2 D-PIV) was applied in the vertical plane of the experimental flume axis.Runs with clear water and weak sediment transport were done under slightly supercritical flow to ensure sediment transport conditions without formation of considerable sediment deposits or dunes.The study design included analyzing the double-averaged flow parameters of the entire measurement domain and investigating the flow development at 14 consecutive vertical subsections.Local geometrical variabilities as well the presence of sediment were mainly reflected in the vertical velocity component.Whereas the vertical velocity decreased over the entire depth in presence of sediment transport,the streamwise velocity profile was reduced only within the interfacial sublayer.In the region with decelerating flow conditions,however,the streamwise velocity profile systematically increased along the entire depth extent.The increase in the main velocity(reduction of flow resistance)correlated with a decrease of the turbulent shear and main normal stresses.Therefore,effects of rough bed smoothening and drag force reduction were experimentally documented within the interfacial sublayer due to mobile sediment.Moreover,the current study leads to the conclusion that in nonuniform flows the maximum Reynolds stress values are a better predictor for the bed shear stress than the linearly extrapolated Reynolds stress profile.This is an important finding because,in natural flows,uniform conditions are rare.     

Submerged hydraulic jump with sediment-laden flow

In this experimental study,the effect of suspended sediment concentration on the characteristics of a submerged hydraulic jump in a rectangular channel has been investigated,based on the sediment conce...     

Linear and angular momentum conservation in hydraulic jump in diverging channels

This paper addresses the integral conservation of linear and angular momentum in the steady hydraulic jump in a linearly diverging channel.The flow is considered to be divided into a mainstream that conveys the total liquid discharge, and a roller where no average mass transport occurs. It is assumed that no macroscopic rheological relationship holds, so mass, momentum and angular momentum integral balances are independent relationships. Normal stresses are assumed to be hydrostatic on vertical, cylindrical surfaces. Viscous stresses are assumed to be negligible with respect to turbulent stresses. Assuming that the horizontal velocity distribution in the mainstream is uniform and that the horizontal momentum and angular momentum in the roller are negligible with respect to their mainstream counterparts, an analytical solution is obtained for the free surface profile of the flow. This solution is fundamental for finding the sequent depths and their positions. Consequently, it permits solving for the length of the jump, which is assumed to be equal to the length of the roller. Mainstream and roller thicknesses can also be derived from the present solution. This model may also be theoretically used to derive the average shear stresses exerted by the roller on the mainstream and the power losses per unit weight. This second relationship, which returns the well-known classical expression for total power loss in the jump, demonstrates that the strongly idealized mechanical model proposed here is internally consistent.     

Reconciling pyroclastic flow and surge: the multiphase physics of pyroclastic density currents

Two end-member types of pyroclastic density current are commonly recognized: pyroclastic surges are dilute currents in which particles are carried in turbulent suspension and pyroclastic flows are highly concentrated flows. We provide scaling relations that unify these end-members and derive a segregation mechanism into basal concentrated flow and overriding dilute cloud based on the Stokes number (S_T), the stability factor (Σ_T) and the dense-dilute condition (D_D). We recognize five types of particle behaviors within a fluid eddy as a function of S_T and Σ_T: (1) particles sediment from the eddy, (2) particles are preferentially settled out during the downward motion of the eddy, but can be carried during its upward motion, (3) particles concentrate on the periphery of the eddy, (4) particles settling can be delayed or ‘fast-tracked’ as a function of the eddy spatial distribution, and (5) particles remain homogeneously distributed within the eddy. We extend these concepts to a fully turbulent flow by using a prototype of kinetic energy distribution within a full eddy spectrum and demonstrate that the presence of different particle sizes leads to the density stratification of the current. This stratification may favor particle interactions in the basal part of the flow and D_D determines whether the flow is dense or dilute. Using only intrinsic characteristics of the current, our model explains the discontinuous features between pyroclastic flows and surges while conserving the concept of a continuous spectrum of density currents.     

Spatial and temporal analysis of hillslope–channel coupling and implications for the longitudinal profile in a dryland basin

The long‐term evolution of channel longitudinal profiles within drainage basins is partly determined by the relative balance of hillslope sediment supply to channels and the evacuation of channel sediment. However, the lack of theoretical understanding of the physical processes of hillslope–channel coupling makes it challenging to determine whether hillslope sediment supply or channel sediment evacuation dominates over different timescales and how this balance affects bed elevation locally along the longitudinal profile. In this paper, we develop a framework for inferring the relative dominance of hillslope sediment supply to the channel versus channel sediment evacuation, over a range of temporal and spatial scales. The framework combines distinct local flow distributions on hillslopes and in the channel with surface grain‐size distributions. We use these to compute local hydraulic stresses at various hillslope‐channel coupling locations within the Walnut Gulch Experimental Watershed (WGEW) in southeast Arizona, USA. These stresses are then assessed as a local net balance of geomorphic work between hillslopes and channel for a range of flow conditions generalizing decadal historical records. Our analysis reveals that, although the magnitude of hydraulic stress in the channel is consistently higher than that on hillslopes, the product of stress magnitude and frequency results in a close balance between hillslope supply and channel evacuation for the most frequent flows. Only at less frequent, high‐magnitude flows do channel hydraulic stresses exceed those on hillslopes, and channel evacuation dominates the net balance. This result suggests that WGEW exists mostly (~50% of the time) in an equilibrium condition of sediment balance between hillslopes and channels, which helps to explain the observed straight longitudinal profile. We illustrate how this balance can be upset by climate changes that differentially affect relative flow regimes on slopes and in channels. Such changes can push the long profile into a convex or concave condition. © 2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

An analysis of the characteristics of rough bed turbulent shear stresses in an open channel

Entrainment of sediment particles from channel beds into the channel flow is influenced by the characteristics of the flow turbulence which produces stochastic shear stress fluctuations at the bed. Recent studies of the structure of turbulent flow has recognized the importance of bursting processes as important mechanisms for the transfer of momentum into the laminar boundary layer. Of these processes, the sweep event has been recognized as the most important bursting event for entrainment of sediment particles as it imposes forces in the direction of the flow resulting in movement of particles by rolling, sliding and occasionally saltating. Similarly, the ejection event has been recognized as important for sediment transport since these events maintain the sediment particles in suspension. In this study, the characteristics of bursting processes and, in particular, the sweep event were investigated in a flume with a rough bed. The instantaneous velocity fluctuations of the flow were measured in two-dimensions using a small electromagnetic velocity meter and the turbulent shear stresses were determined from these velocity fluctuations. It was found that the shear stress applied to the sediment particles on the bed resulting from sweep events depends on the magnitude of the turbulent shear stress and its probability distribution. A statistical analysis of the experimental data was undertaken and it was found necessary to apply a Box-Cox transformation to transform the data into a normally distributed sample. This enabled determination of the mean shear stress, angle of action and standard error of estimate for sweep and ejection events. These instantaneous shear stresses were found to be greater than the mean flow shear stress and for the sweep event to be approximately 40 percent greater near the channel bed. Results from this analysis suggest that the critical shear stress determined from Shield's diagram is not sufficient to predict the initiation of motion due to its use of the temporal mean shear stress. It is suggested that initiation of particle motion, but not continuous motion, can occur earlier than suggested by Shield's diagram due to the higher shear stresses imposed on the particles by the stochastic shear stresses resulting from turbulence within the flow.     

The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.     

DISCONTINUOUS FLOW OF TURBID DENSITY CURRENTS II.INTERNAL HYDRAULIC JUMP

1 INTERNAL HYDRAULIC JUMPS AND LOCAL WATER ENTRAINMENT Traveling and stationary internal hydraulic jumps in two-layer flow can occur in the following situations (Fig. 1). When a stratified fluid flows over an obstacle created by a sudden blockage, the layer of increased depth within the heavier fluid will propagate upstream as a traveling hydraulic jump. When the heavier fluid overtops the obstacle it will flow along the downstream face of the obstacle as a supercritical flow…     

DISCONTINUOUS FLOW OF TURBID DENSITY CURRENTS Ⅱ. INTERNAL HYDRAULIC JUMP

Traveling and stationary internal hydraulic jumps in density currents with positive or negative entrainment coefficients were analyzed based on simple assumptions. An expression of internal hydraulic jumps with entrainment coefficients was derived. Experimental data, published in literature, of stationary internal hydraulic jumps in turbid, thermal and saline density currents including measured values of water entrainment were used to compare with theory. Comparison was also made of traveling internal hydraulic jumps between measured data and theory.     

Intermittent surging movements of a coastal landslide

Since the late 1960s there has been a large mudslide in coastal cliffs of Permo-Triassic strata (conglomerate overlying mudstone) at West Down Beacon, 2 km west of Budleigh Salterton, Devon. A total displacement of 100 m was achieved between 1981 and 1985 as the lobate toe of the mudslide pushed forward across the beach in eleven surges of movement. Each surge involved a displacement of between 5 and 15 m and was completed within a few hours. The toe moved by planar sliding, possibly on more than one seaward-dipping shear surface or zone. Occasionally there were additional relatively minor displacements (less than 1 m), but normally the mudslide was stationary between the major surges. After each surge high-oblique aerial photographs were obtained for stereoscopic interpretation. Debris falls from the cliff appear to have triggered some of the mudslide surges. A displacement recording of one surge has indicated that the rate of movement of the mudslide may have been partly controlled by variations in sea level during the tidal cycle.     

Decadal behaviour of a washover fan,Skallingen Denmark

In this study, the decadal evolution of a washover fan on the west coast of Denmark is examined from its initial generation in 1990 until 2015. Since its inception, the bare and flat washover fan surface has recovered and accreted slowly due to re-activation by overwash during surges and due to aeolian activity and dune formation, stimulated by vegetation growth. The volume of sand on the washover has increased steadily at an average rate of about 23 m³/yr per unit length of shoreline, and a total of 175,000 m³ of sand is now deposited on the fan, while at the same time the shoreline has receded by some 250 m. The evolution can be divided into three stages: 1) An initiation phase when storm surge levels and energetic wave conditions caused a breach in the foredunes and overwash processes formed a washover fan with a relatively low elevation above mean sea level; 2) An initial recovery phase during which waves supplied sand to the fan during frequent overwash activity and winds transported this sand into marginal dunes surrounding the fan; and 3) A later recovery phase when the surface of the fan had accreted to a level where vegetation could survive and trap sediment into new foredune growth across the fan. The rate of accretion has been overall linear but scales with neither annual overwash frequency, nor with aeolian transport potential. Instead, the linear accretion is more closely related to the steady onshore migration of nearshore bars that weld to the beach and provide a sand supply for transfer to the fan. The fan evolution demonstrates the importance of washover fans in preserving barrier resilience during transgressional phases caused by increasing mean sea level. © 2019 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.     

A climatic dataset of ocean vertical turbulent mixing coefficient based on real energy sources

Using data on wind stress, significant height of combined wind waves and swell, potential temperature, salinity and seawater velocity, as well as objectively-analyzed in situ temperature and salinity, we established a global ocean dataset of calculated wind- and tide-induced vertical turbulent mixing coefficients. We then examined energy conservation of ocean vertical mixing from the point of view of ocean wind energy inputs, gravitational potential energy change due to mixing (with and without artificially limiting themixing coefficient), and K-theory vertical turbulent parameterization schemes regardless of energy inputs. Our research showed that calculating the mixing coefficient with average data and artificial limiting the mixing coefficient can cause a remarkable lack of energy conservation, with energy losses of up to 90% and changes in the energy oscillation period. The data also show that wind can introduce a huge amount of energy into the upper layers of the Southern Ocean, and that tidesdo so in regions around underwater mountains. We argue that it is necessary to take wind and tidal energy inputs into account forlong-term ocean climate numerical simulations. We believe that using this ocean vertical turbulent mixing coefficient climatic dataset is a fast and efficient method to maintain the ocean energy balance in ocean modeling research.     

The dependence of UK extreme sea levels and storm surges on the North Atlantic Oscillation

The role of the North Atlantic Oscillation (NAO) in effecting changes in winter extreme high and low waters and storm surges in UK waters has been investigated with the use of a depth-averaged tide+surge numerical model. Spatial patterns of correlation of extreme high and low waters (extreme still water sea levels) with the NAO index are similar to those of median or mean sea level studied previously. Explanations for the similarities, and for differences where they occur, are proposed. Spatial patterns of correlations of extreme high and low and median surge with the NAO index are similar to the corresponding extreme sea-level patterns. Suggestions are made as to which properties of surges (frequency, duration, magnitude) are linked most closely to NAO variability. Several climate models suggest higher (more positive) average values of NAO index during the next 100 years. However, the impact on the UK coastline in terms of increased flood risk should be low (aside from other consequences of climate change such as a global sea-level rise) if the existing relationships between extreme high waters and NAO index are maintained.     

Mt. Pelée, martinique: may 8 and 20, 1902, pyroclastic flows and surges

The distribution, stratigraphic relationships and fragmental components of the May 8 and 20, 1902, pyroclastic flows from Mt. Pelée, Martinique, together with eyewitness accounts, suggest the following explanation for those eruptions. The eruptions were vertically directed magmatic (perhaps initiated phreatically), and contained abundant juvenile lithics from congealed magma of the dome and neck. This resulted in a two-part eruption column having (1) a dense, lithic-charged part which collapsed into the crater and flowed out of a pre-existing notch in its side, giving rise to pyrochlastic flows, and (2) a magmatically derived column containing gases, juvenile vitric material and crystals which largely by-passed the neck and dome and escaped into the atmosphere. All of the energy of the flows was apparently focused through the notch. They emerged fully turbulent and flowed down Rivière Blanche. Gravity segregation of large and abundant fragments soon resulted in a dense, high-concentration, poorly fluidized block-and-ash flow confined to the valley, leaving above a fully turbulent, high-energy ash-cloud surge. As the ash-cloud surge moved down the mountain, it continued to expand outward. The process of gravity segregation continued as the ash-cloud surge expanded, resulting in secondary block-and-ash underflows. Toward St. Pierre, the secondary block-and-ash flows developed on a gently sloping upland surface 100 m or more above the valley of Rivière Blanche. The turbulent, fragment-depleted surges above the secondary block-and-ash flows maintained sufficient energy to devastate the landscape outward to about 3000 m, including St. Pierre. The surges refracted around obstacles and in one place, moved up a small valley in a direction opposite to the main flows.