Sediment transport

Last updated
Dust blowing from the Sahara Desert over the Atlantic Ocean towards the Canary Islands Wafrica amo 2007209 lrg.jpg
Dust blowing from the Sahara Desert over the Atlantic Ocean towards the Canary Islands

Sediment transport is the movement of solid particles (sediment), typically due to a combination of gravity acting on the sediment, and the movement of the fluid in which the sediment is entrained. Sediment transport occurs in natural systems where the particles are clastic rocks (sand, gravel, boulders, etc.), mud, or clay; the fluid is air, water, or ice; and the force of gravity acts to move the particles along the sloping surface on which they are resting. Sediment transport due to fluid motion occurs in rivers, oceans, lakes, seas, and other bodies of water due to currents and tides. Transport is also caused by glaciers as they flow, and on terrestrial surfaces under the influence of wind. Sediment transport due only to gravity can occur on sloping surfaces in general, including hillslopes, scarps, cliffs, and the continental shelf—continental slope boundary.

Contents

Sediment transport is important in the fields of sedimentary geology, geomorphology, civil engineering, hydraulic engineering and environmental engineering (see applications, below). Knowledge of sediment transport is most often used to determine whether erosion or deposition will occur, the magnitude of this erosion or deposition, and the time and distance over which it will occur.

Environments

Aeolian

Sand blowing off a crest in the Kelso Dunes of the Mojave Desert, California KelsoSand.JPG
Sand blowing off a crest in the Kelso Dunes of the Mojave Desert, California

Aeolian or eolian (depending on the parsing of æ) is the term for sediment transport by wind. This process results in the formation of ripples and sand dunes. Typically, the size of the transported sediment is fine sand (<1 mm) and smaller, because air is a fluid with low density and viscosity, and can therefore not exert very much shear on its bed.

Bedforms are generated by aeolian sediment transport in the terrestrial near-surface environment. Ripples [1] and dunes [2] form as a natural self-organizing response to sediment transport.

Aeolian sediment transport is common on beaches and in the arid regions of the world, because it is in these environments that vegetation does not prevent the presence and motion of fields of sand.

Wind-blown very fine-grained dust is capable of entering the upper atmosphere and moving across the globe. Dust from the Sahara deposits on the Canary Islands and islands in the Caribbean, [3] and dust from the Gobi desert has deposited on the western United States. [4] This sediment is important to the soil budget and ecology of several islands.

Deposits of fine-grained wind-blown glacial sediment are called loess.

Fluvial

Deep, eroding glaciofluvial deposits alongside the Matanuska River, Alaska Matanuska River 8727.JPG
Deep, eroding glaciofluvial deposits alongside the Matanuska River, Alaska

In geography and geology, fluvial sediment processes or fluvial sediment transport are associated with rivers and streams and the deposits and landforms created by sediments. It can result in the formation of ripples and dunes, in fractal-shaped patterns of erosion, in complex patterns of natural river systems, and in the development of floodplains and the occurrence of flash floods. Sediment moved by water can be larger than sediment moved by air because water has both a higher density and viscosity. In typical rivers the largest carried sediment is of sand and gravel size, but larger floods can carry cobbles and even boulders.

When the stream or rivers are associated with glaciers, ice sheets, or ice caps, the term glaciofluvial or fluvioglacial is used, as in periglacial flows and glacial lake outburst floods. [5] [6] Fluvial sediment processes include the motion of sediment and erosion or deposition on the river bed. [7] [8]

Coastal

Sand ripples, Laysan Beach, Hawaii. Coastal sediment transport results in these evenly spaced ripples along the shore. Monk seal for scale. Laysan beach.jpg
Sand ripples, Laysan Beach, Hawaii. Coastal sediment transport results in these evenly spaced ripples along the shore. Monk seal for scale.

Coastal sediment transport takes place in near-shore environments due to the motions of waves and currents. At the mouths of rivers, coastal sediment and fluvial sediment transport processes mesh to create river deltas.

Coastal sediment transport results in the formation of characteristic coastal landforms such as beaches, barrier islands, and capes. [9]

A glacier joining the Gorner Glacier, Zermatt, Switzerland. These glaciers transport sediment and leave behind lateral moraines. Glacier.zermatt.arp.750pix.jpg
A glacier joining the Gorner Glacier, Zermatt, Switzerland. These glaciers transport sediment and leave behind lateral moraines.

Glacial

As glaciers move over their beds, they entrain and move material of all sizes. Glaciers can carry the largest sediment, and areas of glacial deposition often contain a large number of glacial erratics, many of which are several metres in diameter. Glaciers also pulverize rock into "glacial flour", which is so fine that it is often carried away by winds to create loess deposits thousands of kilometres afield. Sediment entrained in glaciers often moves approximately along the glacial flowlines, causing it to appear at the surface in the ablation zone.

Hillslope

In hillslope sediment transport, a variety of processes move regolith downslope. These include:

These processes generally combine to give the hillslope a profile that looks like a solution to the diffusion equation, where the diffusivity is a parameter that relates to the ease of sediment transport on the particular hillslope. For this reason, the tops of hills generally have a parabolic concave-up profile, which grades into a convex-up profile around valleys.

As hillslopes steepen, however, they become more prone to episodic landslides and other mass wasting events. Therefore, hillslope processes are better described by a nonlinear diffusion equation in which classic diffusion dominates for shallow slopes and erosion rates go to infinity as the hillslope reaches a critical angle of repose. [10]

Debris flow

Large masses of material are moved in debris flows, hyperconcentrated mixtures of mud, clasts that range up to boulder-size, and water. Debris flows move as granular flows down steep mountain valleys and washes. Because they transport sediment as a granular mixture, their transport mechanisms and capacities scale differently from those of fluvial systems.

Applications

Suspended sediment from a stream emptying into a fjord (Isfjorden, Svalbard, Norway) IsfjordenSediment.JPG
Suspended sediment from a stream emptying into a fjord (Isfjorden, Svalbard, Norway)

Sediment transport is applied to solve many environmental, geotechnical, and geological problems. Measuring or quantifying sediment transport or erosion is therefore important for coastal engineering. Several sediment erosion devices have been designed in order to quantify sediment erosion (e.g., Particle Erosion Simulator (PES)). One such device, also referred to as the BEAST (Benthic Environmental Assessment Sediment Tool) has been calibrated in order to quantify rates of sediment erosion. [11]

Movement of sediment is important in providing habitat for fish and other organisms in rivers. Therefore, managers of highly regulated rivers, which are often sediment-starved due to dams, are often advised to stage short floods to refresh the bed material and rebuild bars. This is also important, for example, in the Grand Canyon of the Colorado River, to rebuild shoreline habitats also used as campsites.

Sediment discharge into a reservoir formed by a dam forms a reservoir delta. This delta will fill the basin, and eventually, either the reservoir will need to be dredged or the dam will need to be removed. Knowledge of sediment transport can be used to properly plan to extend the life of a dam.

Geologists can use inverse solutions of transport relationships to understand flow depth, velocity, and direction, from sedimentary rocks and young deposits of alluvial materials.

Flow in culverts, over dams, and around bridge piers can cause erosion of the bed. This erosion can damage the environment and expose or unsettle the foundations of the structure. Therefore, good knowledge of the mechanics of sediment transport in a built environment are important for civil and hydraulic engineers.

When suspended sediment transport is increased due to human activities, causing environmental problems including the filling of channels, it is called siltation after the grain-size fraction dominating the process.

Initiation of motion

Mechanisms of movement in sediment transport. (a) Roll: sediment particle is rotated by shear stress along the ground (b) Lift: Sediment particle is lifted by the shear stress into the bulk (c) Pluck: Sediment particle is extracted from crevices in the ground by shear stress Sediment1.png
Mechanisms of movement in sediment transport. (a) Roll: sediment particle is rotated by shear stress along the ground (b) Lift: Sediment particle is lifted by the shear stress into the bulk (c) Pluck: Sediment particle is extracted from crevices in the ground by shear stress

Stress balance

For a fluid to begin transporting sediment that is currently at rest on a surface, the boundary (or bed) shear stress exerted by the fluid must exceed the critical shear stress for the initiation of motion of grains at the bed. This basic criterion for the initiation of motion can be written as:

.

This is typically represented by a comparison between a dimensionless shear stress and a dimensionless critical shear stress . The nondimensionalization is in order to compare the driving forces of particle motion (shear stress) to the resisting forces that would make it stationary (particle density and size). This dimensionless shear stress, , is called the Shields parameter and is defined as: [12]

.

And the new equation to solve becomes:

.

The equations included here describe sediment transport for clastic, or granular sediment. They do not work for clays and muds because these types of floccular sediments do not fit the geometric simplifications in these equations, and also interact thorough electrostatic forces. The equations were also designed for fluvial sediment transport of particles carried along in a liquid flow, such as that in a river, canal, or other open channel.

Only one size of particle is considered in this equation. However, river beds are often formed by a mixture of sediment of various sizes. In case of partial motion where only a part of the sediment mixture moves, the river bed becomes enriched in large gravel as the smaller sediments are washed away. The smaller sediments present under this layer of large gravel have a lower possibility of movement and total sediment transport decreases. This is called armouring effect. [13] Other forms of armouring of sediment or decreasing rates of sediment erosion can be caused by carpets of microbial mats, under conditions of high organic loading. [14]

Critical shear stress

Original Shields diagram, 1936 Shields diagram.jpg
Original Shields diagram, 1936

The Shields diagram empirically shows how the dimensionless critical shear stress (i.e. the dimensionless shear stress required for the initiation of motion) is a function of a particular form of the particle Reynolds number, or Reynolds number related to the particle. This allows the criterion for the initiation of motion to be rewritten in terms of a solution for a specific version of the particle Reynolds number, called .

This can then be solved by using the empirically derived Shields curve to find as a function of a specific form of the particle Reynolds number called the boundary Reynolds number. The mathematical solution of the equation was given by Dey. [15]

Particle Reynolds number

In general, a particle Reynolds number has the form:

Where is a characteristic particle velocity, is the grain diameter (a characteristic particle size), and is the kinematic viscosity, which is given by the dynamic viscosity, , divided by the fluid density, .

The specific particle Reynolds number of interest is called the boundary Reynolds number, and it is formed by replacing the velocity term in the particle Reynolds number by the shear velocity, , which is a way of rewriting shear stress in terms of velocity.

where is the bed shear stress (described below), and is the von Kármán constant, where

.

The particle Reynolds number is therefore given by:

Bed shear stress

The boundary Reynolds number can be used with the Shields diagram to empirically solve the equation

,

which solves the right-hand side of the equation

.

In order to solve the left-hand side, expanded as

,

the bed shear stress needs to be found, . There are several ways to solve for the bed shear stress. The simplest approach is to assume the flow is steady and uniform, using the reach-averaged depth and slope. because it is difficult to measure shear stress in situ, this method is also one of the most-commonly used. The method is known as the depth-slope product.

Depth-slope product

For a river undergoing approximately steady, uniform equilibrium flow, of approximately constant depth h and slope angle θ over the reach of interest, and whose width is much greater than its depth, the bed shear stress is given by some momentum considerations stating that the gravity force component in the flow direction equals exactly the friction force. [16] For a wide channel, it yields:

For shallow slope angles, which are found in almost all natural lowland streams, the small-angle formula shows that is approximately equal to , which is given by , the slope. Rewritten with this:

Shear velocity, velocity, and friction factor

For the steady case, by extrapolating the depth-slope product and the equation for shear velocity:

,

The depth-slope product can be rewritten as:

.

is related to the mean flow velocity, , through the generalized Darcy–Weisbach friction factor, , which is equal to the Darcy-Weisbach friction factor divided by 8 (for mathematical convenience). [17] Inserting this friction factor,

.

Unsteady flow

For all flows that cannot be simplified as a single-slope infinite channel (as in the depth-slope product, above), the bed shear stress can be locally found by applying the Saint-Venant equations for continuity, which consider accelerations within the flow.

Example

Set-up

The criterion for the initiation of motion, established earlier, states that

.

In this equation,

, and therefore
.
is a function of boundary Reynolds number, a specific type of particle Reynolds number.
.

For a particular particle Reynolds number, will be an empirical constant given by the Shields Curve or by another set of empirical data (depending on whether or not the grain size is uniform).

Therefore, the final equation to solve is:

.

Solution

Some assumptions allow the solution of the above equation.

The first assumption is that a good approximation of reach-averaged shear stress is given by the depth-slope product. The equation then can be rewritten as:

.

Moving and re-combining the terms produces:

where R is the submerged specific gravity of the sediment.

The second assumption is that the particle Reynolds number is high. This typically applies to particles of gravel-size or larger in a stream, and means the critical shear stress is constant. The Shields curve shows that for a bed with a uniform grain size,

.

Later researchers [18] have shown this value is closer to

for more uniformly sorted beds. Therefore the replacement

is used to insert both values at the end.

The equation now reads:

This final expression shows the product of the channel depth and slope is equal to the Shield's criterion times the submerged specific gravity of the particles times the particle diameter.

For a typical situation, such as quartz-rich sediment in water , the submerged specific gravity is equal to 1.65.

Plugging this into the equation above,

.

For the Shield's criterion of . 0.06 * 1.65 = 0.099, which is well within standard margins of error of 0.1. Therefore, for a uniform bed,

.

For these situations, the product of the depth and slope of the flow should be 10% of the diameter of the median grain diameter.

The mixed-grain-size bed value is , which is supported by more recent research as being more broadly applicable because most natural streams have mixed grain sizes. [18] If this value is used, and D is changed to D_50 ("50" for the 50th percentile, or the median grain size, as an appropriate value for a mixed-grain-size bed), the equation becomes:

Which means that the depth times the slope should be about 5% of the median grain diameter in the case of a mixed-grain-size bed.

Modes of entrainment

The sediments entrained in a flow can be transported along the bed as bed load in the form of sliding and rolling grains, or in suspension as suspended load advected by the main flow. [16] Some sediment materials may also come from the upstream reaches and be carried downstream in the form of wash load.

Rouse number

The location in the flow in which a particle is entrained is determined by the Rouse number, which is determined by the density ρs and diameter d of the sediment particle, and the density ρ and kinematic viscosity ν of the fluid, determine in which part of the flow the sediment particle will be carried. [19]

Here, the Rouse number is given by P. The term in the numerator is the (downwards) sediment the sediment settling velocity ws, which is discussed below. The upwards velocity on the grain is given as a product of the von Kármán constant, κ = 0.4, and the shear velocity, u.

The following table gives the approximate required Rouse numbers for transport as bed load, suspended load, and wash load. [19] [20]

Mode of TransportRouse Number
Initiation of motion>7.5
Bed load >2.5, <7.5
Suspended load: 50% Suspended>1.2, <2.5
Suspended load: 100% Suspended>0.8, <1.2
Wash load <0.8

Settling velocity

Streamlines around a sphere falling through a fluid. This illustration is accurate for laminar flow, in which the particle Reynolds number is small. This is typical for small particles falling through a viscous fluid; larger particles would result in the creation of a turbulent wake. Stokes sphere.svg
Streamlines around a sphere falling through a fluid. This illustration is accurate for laminar flow, in which the particle Reynolds number is small. This is typical for small particles falling through a viscous fluid; larger particles would result in the creation of a turbulent wake.
For particles with a small settling velocity, diffusion will increase the complexity of the particle's path to the bottom and the time it takes to settle compared to particles with high settling velocities. Modeled Sediment Trajectories (1).png
For particles with a small settling velocity, diffusion will increase the complexity of the particle's path to the bottom and the time it takes to settle compared to particles with high settling velocities.

The settling velocity (also called the "fall velocity" or "terminal velocity") is a function of the particle Reynolds number. Generally, for small particles (laminar approximation), it can be calculated with Stokes' Law. For larger particles (turbulent particle Reynolds numbers), fall velocity is calculated with the turbulent drag law. Dietrich (1982) compiled a large amount of published data to which he empirically fit settling velocity curves. [21] Ferguson and Church (2006) analytically combined the expressions for Stokes flow and a turbulent drag law into a single equation that works for all sizes of sediment, and successfully tested it against the data of Dietrich. [22] Their equation is

.

In this equation ws is the sediment settling velocity, g is acceleration due to gravity, and D is mean sediment diameter. is the kinematic viscosity of water, which is approximately 1.0 x 10−6 m2/s for water at 20 °C.

and are constants related to the shape and smoothness of the grains.

ConstantSmooth SpheresNatural Grains: Sieve DiametersNatural Grains: Nominal DiametersLimit for Ultra-Angular Grains
18182024
0.41.01.11.2

The expression for fall velocity can be simplified so that it can be solved only in terms of D. We use the sieve diameters for natural grains, , and values given above for and . From these parameters, the fall velocity is given by the expression:

Alternatively, settling velocity for a particle of sediment can be derived using Stokes Law assuming quiescent (or still) fluid in steady state. The resulting formulation for settling velocity is,

A plot showing the relationship between sediment particle diameter and the Stokes settling velocity Settling velocity.png
A plot showing the relationship between sediment particle diameter and the Stokes settling velocity

where  is the gravitational constant; is the density of the sediment;  is the density of water;  is the sediment particle diameter (commonly assumed to be the median particle diameter, often referred to as  in field studies); and  is the molecular viscosity of water. The Stokes settling velocity can be thought of as the terminal velocity resulting from balancing a particle's buoyant force (proportional to the cross-sectional area) with the gravitational force (proportional to the mass). Small particles will have a slower settling velocity than heavier particles, as seen in the figure. This has implications for many aspects of sediment transport, for example, how far downstream a particle might be advected in a river.

Hjulström–Sundborg diagram

The logarithmic Hjulstrom curve Hjulstroms diagram en.PNG
The logarithmic Hjulström curve

In 1935, Filip Hjulström created the Hjulström curve, a graph which shows the relationship between the size of sediment and the velocity required to erode (lift it), transport it, or deposit it. [23] The graph is logarithmic.

Åke Sundborg later modified the Hjulström curve to show separate curves for the movement threshold corresponding to several water depths, as is necessary if the flow velocity rather than the boundary shear stress (as in the Shields diagram) is used for the flow strength. [24]

This curve has no more than a historical value nowadays, although its simplicity is still attractive. Among the drawbacks of this curve are that it does not take the water depth into account and more importantly, that it does not show that sedimentation is caused by flow velocity deceleration and erosion is caused by flow acceleration. The dimensionless Shields diagram is now unanimously accepted for initiation of sediment motion in rivers.

Transport rate

A schematic diagram of where the different types of sediment load are carried in the flow. Dissolved load is not sediment: it is composed of disassociated ions moving along with the flow. It may, however, constitute a significant proportion (often several percent, but occasionally greater than half) of the total amount of material being transported by the stream. Stream Load.gif
A schematic diagram of where the different types of sediment load are carried in the flow. Dissolved load is not sediment: it is composed of disassociated ions moving along with the flow. It may, however, constitute a significant proportion (often several percent, but occasionally greater than half) of the total amount of material being transported by the stream.

Formulas to calculate sediment transport rate exist for sediment moving in several different parts of the flow. These formulas are often segregated into bed load, suspended load, and wash load. They may sometimes also be segregated into bed material load and wash load.

Bed load

Bed load moves by rolling, sliding, and hopping (or saltating) over the bed, and moves at a small fraction of the fluid flow velocity. Bed load is generally thought to constitute 5–10% of the total sediment load in a stream, making it less important in terms of mass balance. However, the bed material load (the bed load plus the portion of the suspended load which comprises material derived from the bed) is often dominated by bed load, especially in gravel-bed rivers. This bed material load is the only part of the sediment load that actively interacts with the bed. As the bed load is an important component of that, it plays a major role in controlling the morphology of the channel.

Bed load transport rates are usually expressed as being related to excess dimensionless shear stress raised to some power. Excess dimensionless shear stress is a nondimensional measure of bed shear stress about the threshold for motion.

,

Bed load transport rates may also be given by a ratio of bed shear stress to critical shear stress, which is equivalent in both the dimensional and nondimensional cases. This ratio is called the "transport stage" and is an important in that it shows bed shear stress as a multiple of the value of the criterion for the initiation of motion.

When used for sediment transport formulae, this ratio is typically raised to a power.

The majority of the published relations for bedload transport are given in dry sediment weight per unit channel width, ("breadth"):

.

Due to the difficulty of estimating bed load transport rates, these equations are typically only suitable for the situations for which they were designed.

Notable bed load transport formulae

Meyer-Peter Müller and derivatives

The transport formula of Meyer-Peter and Müller, originally developed in 1948, [25] was designed for well-sorted fine gravel at a transport stage of about 8. [19] The formula uses the above nondimensionalization for shear stress, [19]

,

and Hans Einstein's nondimensionalization for sediment volumetric discharge per unit width [19]

.

Their formula reads:

. [19]

Their experimentally determined value for is 0.047, and is the third commonly used value for this (in addition to Parker's 0.03 and Shields' 0.06).

Because of its broad use, some revisions to the formula have taken place over the years that show that the coefficient on the left ("8" above) is a function of the transport stage: [19] [26] [27] [28]

[26]
[27] [28]

The variations in the coefficient were later generalized as a function of dimensionless shear stress: [19] [29]

[29]
Wilcock and Crowe

In 2003, Peter Wilcock and Joanna Crowe (now Joanna Curran) published a sediment transport formula that works with multiple grain sizes across the sand and gravel range. [30] Their formula works with surface grain size distributions, as opposed to older models which use subsurface grain size distributions (and thereby implicitly infer a surface grain sorting).

Their expression is more complicated than the basic sediment transport rules (such as that of Meyer-Peter and Müller) because it takes into account multiple grain sizes: this requires consideration of reference shear stresses for each grain size, the fraction of the total sediment supply that falls into each grain size class, and a "hiding function".

The "hiding function" takes into account the fact that, while small grains are inherently more mobile than large grains, on a mixed-grain-size bed, they may be trapped in deep pockets between large grains. Likewise, a large grain on a bed of small particles will be stuck in a much smaller pocket than if it were on a bed of grains of the same size. In gravel-bed rivers, this can cause "equal mobility", in which small grains can move just as easily as large ones. [31] As sand is added to the system, it moves away from the "equal mobility" portion of the hiding function to one in which grain size again matters. [30]

Their model is based on the transport stage, or ratio of bed shear stress to critical shear stress for the initiation of grain motion. Because their formula works with several grain sizes simultaneously, they define the critical shear stress for each grain size class, , to be equal to a "reference shear stress", . [30]

They express their equations in terms of a dimensionless transport parameter, (with the "" indicating nondimensionality and the "" indicating that it is a function of grain size):

is the volumetric bed load transport rate of size class per unit channel width . is the proportion of size class that is present on the bed.

They came up with two equations, depending on the transport stage, . For :

and for :

.

This equation asymptotically reaches a constant value of as becomes large.

Wilcock and Kenworthy

In 2002, Peter Wilcock and T. A. Kenworthy, following Peter Wilcock (1998), [32] published a sediment bed-load transport formula that works with only two sediments fractions, i.e. sand and gravel fractions. [33] A mixed-sized sediment bed-load transport model using only two fractions offers practical advantages in terms of both computational and conceptual modeling by taking into account the nonlinear effects of sand presence in gravel beds on bed-load transport rate of both fractions. In fact, in the two-fraction bed load formula appears a new ingredient with respect to that of Meyer-Peter and Müller that is the proportion of fraction on the bed surface where the subscript represents either the sand (s) or gravel (g) fraction. The proportion , as a function of sand content , physically represents the relative influence of the mechanisms controlling sand and gravel transport, associated with the change from a clast-supported to matrix-supported gravel bed. Moreover, since spans between 0 and 1, phenomena that vary with include the relative size effects producing "hiding" of fine grains and "exposure" of coarse grains. The "hiding" effect takes into account the fact that, while small grains are inherently more mobile than large grains, on a mixed-grain-size bed, they may be trapped in deep pockets between large grains. Likewise, a large grain on a bed of small particles will be stuck in a much smaller pocket than if it were on a bed of grains of the same size, which the Meyer-Peter and Müller formula refers to. In gravel-bed rivers, this can cause "equal mobility", in which small grains can move just as easily as large ones. [31] As sand is added to the system, it moves away from the "equal mobility" portion of the hiding function to one in which grain size again matters. [33]

Their model is based on the transport stage, i.e., or ratio of bed shear stress to critical shear stress for the initiation of grain motion. Because their formula works with only two fractions simultaneously, they define the critical shear stress for each of the two grain size classes, , where represents either the sand (s) or gravel (g) fraction. The critical shear stress that represents the incipient motion for each of the two fractions is consistent with established values in the limit of pure sand and gravel beds and shows a sharp change with increasing sand content over the transition from a clast- to matrix-supported bed. [33]

They express their equations in terms of a dimensionless transport parameter, (with the "" indicating nondimensionality and the "" indicating that it is a function of grain size):

is the volumetric bed load transport rate of size class per unit channel width . is the proportion of size class that is present on the bed.

They came up with two equations, depending on the transport stage, . For :

and for :

.

This equation asymptotically reaches a constant value of as becomes large and the symbols have the following values:

In order to apply the above formulation, it is necessary to specify the characteristic grain sizes for the sand portion and for the gravel portion of the surface layer, the fractions and of sand and gravel, respectively in the surface layer, the submerged specific gravity of the sediment R and shear velocity associated with skin friction .

Kuhnle et al.

For the case in which sand fraction is transported by the current over and through an immobile gravel bed, Kuhnle et al.(2013), [34] following the theoretical analysis done by Pellachini (2011), [35] provides a new relationship for the bed load transport of the sand fraction when gravel particles remain at rest. It is worth mentioning that Kuhnle et al. (2013) [34] applied the Wilcock and Kenworthy (2002) [33] formula to their experimental data and found out that predicted bed load rates of sand fraction were about 10 times greater than measured and approached 1 as the sand elevation became near the top of the gravel layer. [34] They, also, hypothesized that the mismatch between predicted and measured sand bed load rates is due to the fact that the bed shear stress used for the Wilcock and Kenworthy (2002) [33] formula was larger than that available for transport within the gravel bed because of the sheltering effect of the gravel particles. [34] To overcome this mismatch, following Pellachini (2011), [35] they assumed that the variability of the bed shear stress available for the sand to be transported by the current would be some function of the so-called "Roughness Geometry Function" (RGF), [36] which represents the gravel bed elevations distribution. Therefore, the sand bed load formula follows as: [34]

where

the subscript refers to the sand fraction, s represents the ratio where is the sand fraction density, is the RGF as a function of the sand level within the gravel bed, is the bed shear stress available for sand transport and is the critical shear stress for incipient motion of the sand fraction, which was calculated graphically using the updated Shields-type relation of Miller et al.(1977). [37]

Suspended load

Suspended load is carried in the lower to middle parts of the flow, and moves at a large fraction of the mean flow velocity in the stream.

A common characterization of suspended sediment concentration in a flow is given by the Rouse Profile. This characterization works for the situation in which sediment concentration at one particular elevation above the bed can be quantified. It is given by the expression:

Here, is the elevation above the bed, is the concentration of suspended sediment at that elevation, is the flow depth, is the Rouse number, and relates the eddy viscosity for momentum to the eddy diffusivity for sediment, which is approximately equal to one. [38]

Experimental work has shown that ranges from 0.93 to 1.10 for sands and silts. [39]

The Rouse profile characterizes sediment concentrations because the Rouse number includes both turbulent mixing and settling under the weight of the particles. Turbulent mixing results in the net motion of particles from regions of high concentrations to low concentrations. Because particles settle downward, for all cases where the particles are not neutrally buoyant or sufficiently light that this settling velocity is negligible, there is a net negative concentration gradient as one goes upward in the flow. The Rouse Profile therefore gives the concentration profile that provides a balance between turbulent mixing (net upwards) of sediment and the downwards settling velocity of each particle.

Bed material load

Bed material load comprises the bed load and the portion of the suspended load that is sourced from the bed.

Three common bed material transport relations are the "Ackers-White", [40] "Engelund-Hansen", "Yang" formulae. The first is for sand to granule-size gravel, and the second and third are for sand [41] though Yang later expanded his formula to include fine gravel. That all of these formulae cover the sand-size range and two of them are exclusively for sand is that the sediment in sand-bed rivers is commonly moved simultaneously as bed and suspended load.

Engelund–Hansen

The bed material load formula of Engelund and Hansen is the only one to not include some kind of critical value for the initiation of sediment transport. It reads:

where is the Einstein nondimensionalization for sediment volumetric discharge per unit width, is a friction factor, and is the Shields stress. The Engelund–Hansen formula is one of the few sediment transport formulae in which a threshold "critical shear stress" is absent.

Wash load

Wash load is carried within the water column as part of the flow, and therefore moves with the mean velocity of main stream. Wash load concentrations are approximately uniform in the water column. This is described by the endmember case in which the Rouse number is equal to 0 (i.e. the settling velocity is far less than the turbulent mixing velocity), which leads to a prediction of a perfectly uniform vertical concentration profile of material.

Total load

Some authors have attempted formulations for the total sediment load carried in water. [42] [43] These formulas are designed largely for sand, as (depending on flow conditions) sand often can be carried as both bed load and suspended load in the same stream or shoreface.

Bed load sediment mitigation at intake structures

Riverside intake structures used in water supply, canal diversions, and water cooling can experience entrainment of bed load (sand-size) sediments. These entrained sediments produce multiple deleterious effects such as reduction or blockage of intake capacity, feedwater pump impeller damage or vibration, and result in sediment deposition in downstream pipelines and canals. Structures that modify local near-field secondary currents are useful to mitigate these effects and limit or prevent bed load sediment entry. [44]

See also

Related Research Articles

<span class="mw-page-title-main">Bingham plastic</span> Material which is solid at low stress but becomes viscous at high stress

In materials science, a Bingham plastic is a viscoplastic material that behaves as a rigid body at low stresses but flows as a viscous fluid at high stress. It is named after Eugene C. Bingham who proposed its mathematical form.

<span class="mw-page-title-main">Bed load</span> Particles in a flowing fluid that are transported along the bed

The term bed load or bedload describes particles in a flowing fluid that are transported along the stream bed. Bed load is complementary to suspended load and wash load.

<span class="mw-page-title-main">Soil mechanics</span> Branch of soil physics and applied mechanics that describes the behavior of soils

Soil mechanics is a branch of soil physics and applied mechanics that describes the behavior of soils. It differs from fluid mechanics and solid mechanics in the sense that soils consist of a heterogeneous mixture of fluids and particles but soil may also contain organic solids and other matter. Along with rock mechanics, soil mechanics provides the theoretical basis for analysis in geotechnical engineering, a subdiscipline of civil engineering, and engineering geology, a subdiscipline of geology. Soil mechanics is used to analyze the deformations of and flow of fluids within natural and man-made structures that are supported on or made of soil, or structures that are buried in soils. Example applications are building and bridge foundations, retaining walls, dams, and buried pipeline systems. Principles of soil mechanics are also used in related disciplines such as geophysical engineering, coastal engineering, agricultural engineering, hydrology and soil physics.

<span class="mw-page-title-main">Ekman spiral</span> Velocity profile of wind driven current with depth

The Ekman spiral is an arrangement of ocean currents: the directions of horizontal current appear to twist as the depth changes. The oceanic wind driven Ekman spiral is the result of a force balance created by a shear stress force, Coriolis force and the water drag. This force balance gives a resulting current of the water different from the winds. In the ocean, there are two places where the Ekman spiral can be observed. At the surface of the ocean, the shear stress force corresponds with the wind stress force. At the bottom of the ocean, the shear stress force is created by friction with the ocean floor. This phenomenon was first observed at the surface by the Norwegian oceanographer Fridtjof Nansen during his Fram expedition. He noticed that icebergs did not drift in the same direction as the wind. His student, the Swedish oceanographer Vagn Walfrid Ekman, was the first person to physically explain this process.

<span class="mw-page-title-main">Plug flow</span> Simple model of fluid flow in a pipe

In fluid mechanics, plug flow is a simple model of the velocity profile of a fluid flowing in a pipe. In plug flow, the velocity of the fluid is assumed to be constant across any cross-section of the pipe perpendicular to the axis of the pipe. The plug flow model assumes there is no boundary layer adjacent to the inner wall of the pipe.

<span class="mw-page-title-main">Love wave</span> Horizontally polarized surface waves

In elastodynamics, Love waves, named after Augustus Edward Hough Love, are horizontally polarized surface waves. The Love wave is a result of the interference of many shear waves (S-waves) guided by an elastic layer, which is welded to an elastic half space on one side while bordering a vacuum on the other side. In seismology, Love waves (also known as Q waves (Quer: German for lateral)) are surface seismic waves that cause horizontal shifting of the Earth during an earthquake. Augustus Edward Hough Love predicted the existence of Love waves mathematically in 1911. They form a distinct class, different from other types of seismic waves, such as P-waves and S-waves (both body waves), or Rayleigh waves (another type of surface wave). Love waves travel with a lower velocity than P- or S- waves, but faster than Rayleigh waves. These waves are observed only when there is a low velocity layer overlying a high velocity layer/ sub–layers.

<span class="mw-page-title-main">Yield (engineering)</span> Phenomenon of deformation due to structural stress

In materials science and engineering, the yield point is the point on a stress-strain curve that indicates the limit of elastic behavior and the beginning of plastic behavior. Below the yield point, a material will deform elastically and will return to its original shape when the applied stress is removed. Once the yield point is passed, some fraction of the deformation will be permanent and non-reversible and is known as plastic deformation.

The suspended load of a flow of fluid, such as a river, is the portion of its sediment uplifted by the fluid's flow in the process of sediment transportation. It is kept suspended by the fluid's turbulence. The suspended load generally consists of smaller particles, like clay, silt, and fine sands.

The derivation of the Navier–Stokes equations as well as their application and formulation for different families of fluids, is an important exercise in fluid dynamics with applications in mechanical engineering, physics, chemistry, heat transfer, and electrical engineering. A proof explaining the properties and bounds of the equations, such as Navier–Stokes existence and smoothness, is one of the important unsolved problems in mathematics.

<span class="mw-page-title-main">Critical state soil mechanics</span>

Critical state soil mechanics is the area of soil mechanics that encompasses the conceptual models representing the mechanical behavior of saturated remoulded soils based on the critical state concept. At the critical state, the relationship between forces applied in the soil (stress), and the resulting deformation resulting from this stress (strain) becomes constant. The soil will continue to deform, but the stress will no longer increase.

The Cauchy momentum equation is a vector partial differential equation put forth by Cauchy that describes the non-relativistic momentum transport in any continuum.

Shear velocity, also called friction velocity, is a form by which a shear stress may be re-written in units of velocity. It is useful as a method in fluid mechanics to compare true velocities, such as the velocity of a flow in a stream, to a velocity that relates shear between layers of flow.

The depth–slope product is used to calculate the shear stress at the bed of an open channel containing fluid that is undergoing steady, uniform flow. It is widely used in river engineering, stream restoration, sedimentology, and fluvial geomorphology. It is the product of the water depth and the mean bed slope, along with the acceleration due to gravity and density of the fluid.

Ocean dynamics define and describe the flow of water within the oceans. Ocean temperature and motion fields can be separated into three distinct layers: mixed (surface) layer, upper ocean, and deep ocean.

<span class="mw-page-title-main">Stream power</span>

Stream power, originally derived by R. A. Bagnold in the 1960s, is the amount of energy the water in a river or stream is exerting on the sides and bottom of the river. Stream power is the result of multiplying the density of the water, the acceleration of the water due to gravity, the volume of water flowing through the river, and the slope of that water. There are many forms of the stream power formula with varying utilities, such as comparing rivers of various widths or quantifying the energy required to move sediment of a certain size. Stream power is closely related to other criteria such as stream competency and shear stress. Stream power is a valuable measurement for hydrologists and geomorphologists tackling sediment transport issues as well as for civil engineers, who use it in the planning and construction of roads, bridges, dams, and culverts.

The Shields parameter, also called the Shields criterion or Shields number, is a nondimensional number used to calculate the initiation of motion of sediment in a fluid flow. It is a nondimensionalization of a shear stress, and is typically denoted or . This parameter has been developed by Albert F. Shields, and is called later Shields parameter. The Shields parameter is the main parameter of the Shields formula. It is given by:

<span class="mw-page-title-main">Stream competency</span> Concept in hydrology

In hydrology stream competency, also known as stream competence, is a measure of the maximum size of particles a stream can transport. The particles are made up of grain sizes ranging from large to small and include boulders, rocks, pebbles, sand, silt, and clay. These particles make up the bed load of the stream. Stream competence was originally simplified by the “sixth-power-law,” which states the mass of a particle that can be moved is proportional to the velocity of the river raised to the sixth power. This refers to the stream bed velocity which is difficult to measure or estimate due to the many factors that cause slight variances in stream velocities.

<span class="mw-page-title-main">Radiation stress</span> Term in physical oceanography

In fluid dynamics, the radiation stress is the depth-integrated – and thereafter phase-averaged – excess momentum flux caused by the presence of the surface gravity waves, which is exerted on the mean flow. The radiation stresses behave as a second-order tensor.

The multiphase particle-in-cell method (MP-PIC) is a numerical method for modeling particle-fluid and particle-particle interactions in a computational fluid dynamics (CFD) calculation. The MP-PIC method achieves greater stability than its particle-in-cell predecessor by simultaneously treating the solid particles as computational particles and as a continuum. In the MP-PIC approach, the particle properties are mapped from the Lagrangian coordinates to an Eulerian grid through the use of interpolation functions. After evaluation of the continuum derivative terms, the particle properties are mapped back to the individual particles. This method has proven to be stable in dense particle flows, computationally efficient, and physically accurate. This has allowed the MP-PIC method to be used as particle-flow solver for the simulation of industrial-scale chemical processes involving particle-fluid flows.

<span class="mw-page-title-main">Shields formula</span> Parameter (and formula) to describe stability of grains in flowing water

The Shields formula is a formula for the stability calculation of granular material in running water.

References

  1. Anderson, R (1990). "Eolian ripples as examples of self-organization in geomorphological systems". Earth-Science Reviews. 29 (1–4): 77. doi:10.1016/0012-8252(0)90029-U.
  2. Kocurek, Gary; Ewing, Ryan C. (2005). "Aeolian dune field self-organization – implications for the formation of simple versus complex dune-field patterns". Geomorphology. 72 (1–4): 94. Bibcode:2005Geomo..72...94K. doi:10.1016/j.geomorph.2005.05.005.
  3. Goudie, A; Middleton, N.J. (2001). "Saharan dust storms: nature and consequences". Earth-Science Reviews. 56 (1–4): 179. Bibcode:2001ESRv...56..179G. doi:10.1016/S0012-8252(01)00067-8.
  4. "Dust Storm Spreads Out of Gobi Desert". Earthobservatory.nasa.gov. 13 April 2006. Retrieved 2022-05-08.
  5. Neuendorf, Klaus K.E.; Mehl, James P. Jr; Jackson, Julia A., eds. (2011). Glossary of Geology (5th revised ed.). Alexandria, Virginia: American Geological Institute. p. 800. ISBN   978-3-642-06621-4. OCLC   751527782.
  6. Wilson, W.E. & Moore, J.E. 2003. Glossary of Hydrology, American Geological Institute, Springer, 248pp.
  7. Charlton, Ro (2008). Fundamentals of fluvial geomorphology . London: Rutledge. p.  234. ISBN   978-0-415-33454-9.
  8. Wohl, Ellen (2014). Rivers in the Landscape: Science and Management. Wiley-Blackwell. p. 330. ISBN   978-1-118-41489-7.
  9. Ashton, Andrew; Murray, A. Brad; Arnault, Olivier (2001). "Formation of coastline features by large-scale instabilities induced by high-angle waves". Nature. 414 (6861): 296–300. Bibcode:2001Natur.414..296A. doi:10.1038/35104541. PMID   11713526. S2CID   205023325.
  10. Roering, Joshua J.; Kirchner, James W.; Dietrich, William E. (1999). "Evidence for nonlinear, diffusive sediment transport on hillslopes and implications for landscape morphology". Water Resources Research. 35 (3): 853. Bibcode:1999WRR....35..853R. doi: 10.1029/1998WR900090 .
  11. Grant, J.; Walker, T.R.; Hill, P.S.; Lintern, D.G. (2013). "BEAST-A portable device for quantification of erosion in intact sediment cores". Methods in Oceanography. 5: 39–55. doi:10.1016/j.mio.2013.03.001.
  12. Shields, A. (1936) Anwendung der Ähnlichkeitsmechanik und der Turbulenzforschung auf die Geschiebebewegung; In Mitteilungen der Preussischen Versuchsanstalt für Wasserbau und Schiffbau, Heft 26 (Online Archived 2011-07-18 at the Wayback Machine ; PDF; 3,8 MB)
  13. Sharmeen, Saniya; Willgoose, Garry R. (2006). "The interaction between armouring and particle weathering for eroding landscapes". Earth Surface Processes and Landforms . 31 (10): 1195–1210. Bibcode:2006ESPL...31.1195S. doi: 10.1002/esp.1397 . S2CID   91175516.
  14. Walker, T.R.; Grant, J. (2009). "Quantifying erosion rates and stability of bottom sediments at mussel aquaculture sites in Prince Edward Island, Canada". Journal of Marine Systems. 75 (1–2): 46–55. Bibcode:2009JMS....75...46W. doi:10.1016/j.jmarsys.2008.07.009.
  15. Dey S. (1999) Sediment threshold. Applied Mathematical Modelling, Elsevier, Vol. 23, No. 5, 399-417.
  16. 1 2 Hubert Chanson (2004). The Hydraulics of Open Channel Flow: An Introduction. Butterworth-Heinemann, 2nd edition, Oxford, UK, 630 pages. ISBN   978-0-7506-5978-9.
  17. Whipple, Kelin (2004). "Hydraulic Roughness" (PDF). 12.163: Surface processes and landscape evolution. MIT OCW. Retrieved 2009-03-27.
  18. 1 2 Parker, G (1990). "Surface-based bedload transport relation for gravel rivers". Journal of Hydraulic Research. 28 (4): 417–436. Bibcode:1990JHydR..28..417P. doi:10.1080/00221689009499058.
  19. 1 2 3 4 5 6 7 8 Whipple, Kelin (September 2004). "IV. Essentials of Sediment Transport" (PDF). 12.163/12.463 Surface Processes and Landscape Evolution: Course Notes. MIT OpenCourseWare . Retrieved 2009-10-11.
  20. Moore, Andrew. "Lecture 20—Some Loose Ends" (PDF). Lecture Notes: Fluvial Sediment Transport. Kent State. Retrieved 23 December 2009.
  21. Dietrich, W. E. (1982). "Settling Velocity of Natural Particles" (PDF). Water Resources Research. 18 (6): 1615–1626. Bibcode:1982WRR....18.1615D. doi:10.1029/WR018i006p01615. S2CID   128625889.
  22. Ferguson, R. I.; Church, M. (2006). "A Simple Universal Equation for Grain Settling Velocity". Journal of Sedimentary Research. 74 (6): 933–937. doi:10.1306/051204740933.
  23. The long profile – changing processes: types of erosion, transportation and deposition, types of load; the Hjulstrom curve. coolgeography.co.uk. Last accessed 26 Dec 2011.
  24. Special Topics: An Introduction to Fluid Motions, Sediment Transport, and Current-generated Sedimentary Structures; As taught in: Fall 2006. Massachusetts Institute of Technology. 2006. Last accessed 26 Dec 2011.
  25. Meyer-Peter, E; Müller, R. (1948). Formulas for bed-load transport. Proceedings of the 2nd Meeting of the International Association for Hydraulic Structures Research. pp. 39–64.
  26. 1 2 Fernandez-Luque, R; van Beek, R (1976). "Erosion and transport of bedload sediment". J. Hydrol. Res. 14 (2): 127–144. Bibcode:1976JHydR..14..127F. doi:10.1080/00221687609499677.
  27. 1 2 Cheng, Nian-Sheng (2002). "Exponential Formula for Bedload Transport". Journal of Hydraulic Engineering. 128 (10): 942. doi:10.1061/(ASCE)0733-9429(2002)128:10(942). hdl: 10356/83917 .
  28. 1 2 Wilson, K. C. (1966). "Bed-load transport at high shear stress". J. Hydraul. Div. 92 (6). ASCE: 49–59. doi:10.1061/JYCEAJ.0001562.
  29. 1 2 Wiberg, Patricia L.; Dungan Smith, J. (1989). "Model for Calculating Bed Load Transport of Sediment". Journal of Hydraulic Engineering. 115: 101. doi:10.1061/(ASCE)0733-9429(1989)115:1(101).
  30. 1 2 3 Wilcock, Peter R.; Crowe, Joanna C. (2003). "Surface-based Transport Model for Mixed-Size Sediment". Journal of Hydraulic Engineering. 129 (2): 120. doi:10.1061/(ASCE)0733-9429(2003)129:2(120). S2CID   40651180.
  31. 1 2 Parker, G.; Klingeman, P. C.; McLean, D. G. (1982). "Bedload and Size Distribution in Paved Gravel-Bed Streams". Journal of the Hydraulics Division. 108 (4). ASCE: 544–571. doi:10.1061/JYCEAJ.0005854.
  32. Wilcock, P. R. (1998). "Two-fraction model of initial sediment motion in gravel-bed rivers". Science. 280 (5362): 410–412. Bibcode:1998Sci...280..410W. doi:10.1126/science.280.5362.410. PMID   9545213.
  33. 1 2 3 4 5 Wilcock, Peter R.; Kenworthy, T. (2002). "A two-fraction model for the transport of sand/gravel mixtures". Water Resour. Res. 38 (10): 1194. Bibcode:2002WRR....38.1194W. doi: 10.1029/2001WR000684 .
  34. 1 2 3 4 5 Kuhnle, R. A.; Wren, D. G.; Langendoen, E. J.; Rigby, J. R. (2013). "Sand Transport over an Immobile Gravel Substrate". Journal of Hydraulic Engineering. 139 (2): 167–176. doi:10.1061/(ASCE)HY.1943-7900.0000615.
  35. 1 2 Pellachini, Corrado (2011). Modelling fine sediment transport over an immobile gravel bed (phd). Trento: Unitn-eprints.
  36. Nikora, V; Goring, D; McEwan, I; Griffiths, G (2001). "Spatially averaged open-channel flow over rough bed". J. Hydraul. Eng. 127 (2): 123–133. doi:10.1061/(ASCE)0733-9429(2001)127:2(123).
  37. Miller, M.C.; McCave, I.N.; Komar, P.D. (1977). "Threshold of sediment motion under unidirectional currents". Sedimentology. 24 (4): 507–527. Bibcode:1977Sedim..24..507M. doi:10.1111/j.1365-3091.1977.tb00136.x.
  38. Harris, Courtney K. (March 18, 2003). "Lecture 9: Suspended Sediment Transport II" (PDF). Sediment transport processes in coastal environments. Virginia Institute of Marine Science. Archived from the original (PDF) on 28 May 2010. Retrieved 23 December 2009.
  39. Moore, Andrew. "Lecture 21—Suspended Sediment Transport" (PDF). Lecture Notes: Fluvial Sediment Transport. Kent State. Retrieved 25 December 2009.
  40. Ackers, P.; White, W.R. (1973). "Sediment Transport: New Approach and Analysis". Journal of the Hydraulics Division . 99 (11). ASCE: 2041–2060. doi:10.1061/JYCEAJ.0003791.
  41. Ariffin, J.; A.A. Ghani; N.A. Zakaira; A.H. Yahya (14–16 October 2002). "Evaluation of equations on total bed material load" (PDF). International Conference on Urban Hydrology for the 21st Century. Kuala Lumpur.
  42. Yang, C (1979). "Unit stream power equations for total load". Journal of Hydrology. 40 (1–2): 123. Bibcode:1979JHyd...40..123Y. doi:10.1016/0022-1694(79)90092-1.
  43. Bailard, James A. (1981). "An Energetics Total Load Sediment Transport Model For a Plane Sloping Beach". Journal of Geophysical Research. 86 (C11): 10938. Bibcode:1981JGR....8610938B. doi:10.1029/JC086iC11p10938.
  44. Natato, T.; Ogden, F.L. (1998). "Sediment control at water intakes along sand-bed rivers". Journal of Hydraulic Engineering. 126 (6): 589–596. doi:10.1061/(ASCE)0733-9429(1998)124:6(589).