Use of the Universal Soil Loss Equation
The USLE computes upland erosion from small watersheds on an average annual basis.
It includes the detachment and transport components, but it does not account for the deposition component.
Therefore, the USLE cannot be used to compute sediment yield.
For example, in a 1000-mi2 drainage basin, only 5 percent of the soil loss computed by the USLE may appear as sediment yield at the basin outlet.
The remaining 9S percent is redistributed on uplands or flood plains and does not constitute a net soil loss from the drainage basin.
Example 15-3
Assume a 600-ac watershed above a proposed floodwater-retarding structure in Fountain County, Indiana.
Compute the average annual soil loss by the universal soil loss equation for the following conditions: (1) cropland, 280 ac, contour strip-cropped, soil is Fayette silt loam, slopes are 8% and 200 ft long; (2) pasture, 170 ac, 50% canopy cover, 80% ground cover with grass, soil is Fayette silt loam, slopes are 8% and 200 ft long; and (3) forest, 150 ac, soil is Marshall silt loam, 30% tree canopy cover, slopes are 12% and 100 ft long.
- From Fig. 15-2, R = 185.
From Table 15-4, K = 0.38.
From Fig. 15-3, LS = 1.4.
The value of C for cropland is obtained from local sources; assume C = 0.12 for this example.
From Table 15-7, P = 0.25.
Using Eq. 15-6: A = 185 X 0.38 X 1.4 X 0.12 X 0.25 = 2.95 tons/ac/y.
- R = 185; K = 0.38: LS = 1.4.
From Table 15-5, C = 0.012.
No value of P has been established for pasture; therefore, P = 1.
Using Eq. 15-6: A = 185 X 0.38 X 1.4 X 0.012 X 1.0 = 1.18 tons/ac/y.
- R = 185.
From Table 15-4, K = 0.33.
From Fig. 15-3, LS = 1.8.
From Table 15-6. C = 0.006.
No value of P has been established for forest.
Using Eq. 15-6: A = 185 X 0.33 X 1.5 X 0.006 X 1.0 = 0.66 tons/ac/y.
The total sheet and rill erosion from the 600-ac watershed is (280 X 2.95) + (170 X 1.18) + (150 X 0.66) = 1126 tons/y.
Channel Erosion
Channel erosion includes gully erosion, streambank erosion, streambed degradation, floodplain scour, and other sources of sediment, excluding upland erosion. Gullies are incipient channels in process of development.
Gully growth is usually accelerated by severe climatic events, improper land use, or changes in stream base levels.
Most of the significant gully activity, in terms of the quantities of sediment produced and delivered to downstream locations, is found in regions of moderate to steep topography having thick soil mantles.
The total sediment outflow from gullies is usually less than sheet and rill erosion [31].
Streambank erosion and streambed degradation can be significant in certain cases.
Changes in channel alignment and/ or removal of natural vegetation from stream banks may cause increased bank erosion.
Streambed degradation, typically downstream of reservoirs, can also constitute an additional source of sediment.
Methods for determining soil loss due to the various types of channel erosion include the following: (1) comparing aerial photographs taken at different times to assess the growth rate of channels, (2) performing river cross-sectional surveys to determine changes in cross-sectional area, (3) assembling historical data to determine the average age and growth rate of channels, and (4) performing field studies to evaluate the annual growth rate of channels.
Field surveys can often provide sufficient data to estimate streambank erosion as follows [41]:
S = HLR (15-7)
in which S = annual volume of streambank erosion; H = average height of bank; L = length of eroded bank, each side of channel if both sides are eroding; and R = annual rate of bank recession (net rate if one side is eroding while the other is depositing).
Streambed degradation can be estimated as follows [41]:
S = WLD (15-8)
in which S = annual volume of streambed degradation; W = average bottom width of degrading channel reach; L = length of degrading channel reach; and D = annual rate of streambed degradation.
Accelerated Erosion Due to Strip-mining and Construction Activities
Strip-mining and construction activities greatly accelerate erosion rates.
For instance, Collier et al. [10] found that a watershed with 10.4 percent of its area strip-mined eroded 76 times more sediment than a similar undisturbed watershed.
Wolman and Shick [45] found that sediment concentrations in streams draining construction areas ranged from 3000 to 150,000 ppm, compared to concentrations of 2000 ppm in comparable natural settings.
These studies indicate that human-induced land disturbances have a substantial impact on sediment production.
With a careful choice of factors, the USLE can be used to compute soil loss from disturbed lands.
Sediment Yield
In engineering applications, the quantity of sediment eroded at the sources is not as important as the quantity of sediment delivered to a downstream point, i.e., the sediment yield.
The sediment yield is calculated by multiplying the gross sediment production, which includes all types of erosion (sheet, rill , gully, and channel erosion) by a sediment- delivery ratio that varies in the range 0 to 1 (it can also be expressed as a percentage).
Therefore, a calculation of sediment yield hinges upon an estimate of gross sediment production (from the various sources) and an appropriate sediment-delivery ratio.
Sediment-delivery Ratio
The sediment-delivery ratio (SDR) is largely a function of (1) sediment source, (2) proximity of sediment source to the fluvial transport system, (3) density and condition of the fluvial transport system, (4) sediment size and texture, and (5) catchment characteristics.
The sediment source has an influence on the delivery ratio.
Not all sediments originating in sheet and rill erosion are likely to enter the fluvial transport system.
However, sediments produced by channel erosion are generally closer to the transport system and are more likely to be delivered to downstream points.
The proximity of the sediment source to the transport system is also an important variable in the estimation of SD R.
The amounts of sediment delivered to downstream points will depend to a large extent on the ability of the fluvial transport system to entrain and hold on to the sediment particles.
Silt and clay particles can be transported much more readily than sand particles; therefore, the delivery of silts and clays is more likely to occur than that of sands.
Catchment characteristics also affect sediment-delivery ratios.
High relief often indicates both a high erosion rate and a high SDR.
High channel density is usually an indication of an efficient transport system and, consequently, of a high SDR.
Estimation of Sediment-delivery Ratios.
The SDR is the ratio of sediment yield to gross sediment production. Sediment yield can be evaluated by one of several methods. At reservoir locations, estimates of sediment yield can be obtained by reservoir sedimentation surveys. Alternatively, sediment yield can be evaluated by direct measurement of sediment load at the point of interest. Estimates of gross sediment production from upland sources can be obtained using either the USLE formula or a regionally derived formula for sheet and rill erosion. When warranted, this estimate can be augmented by field estimates of gully and channel erosion.
In the absence of actual measurements, statistical analysis can be used to develop regional regression equations to predict SDRs.
The simplest SDR prediction equation is that based solely on drainage area, as shown in Fig. 15-4.
This figure shows that SDR varies approximately in inverse proportion to the 1/5 power of the drainage area.
Other sources, however. have quoted values of this power as low as 1/8 [2].
The fact remains that the greater the drainage area, the smaller the catchment relief and the greater the chances for sediment deposition within the catchment.
Consequently, the smaller the catchment's SDR.
Rough estimates of SDR can be obtained from Fig. 15-4. but caution is recommended for more refined studies.
An example of the use of statistical analysis for the estimation of sediment-delivery ratios is given by Roehl [34].
Using data from the southeast Piedmont region of the United States, he developed the following predictive equation:
(15-9)
in which SDR = sediment-delivery ratio in percentage; A = drainage area in square miles; L/R = dimensionless ratio of catchment length-to-relief (length measured parallel to main drainageway, relief measured as the elevation difference between drainage divide and outlet); and B = weighted mean bifurcation ratio, defined as the the ratio of the number of streams in a given order to the number of streams in the next higher order.
Values of SDR measured by Roehl in the Piedmont area were in the range 3.7 to 59.4 percent.
Empirical Formulas for Sediment Yield
As with the SDRs, statistical analysis can be used to develop regional equations for the prediction of sediment yield.
A study by Dendy and Bolton [13] showed that sediment yield can be related to catchment area and mean annual runoff.
Sediment Yield versus Drainage Area.
Dendy and Bolton studied sedimentation data from about 1500 reservoirs, ponds, and sediment detention basins.
In developing their formulas, they used data from about 800 of these reservoirs with drainage areas greater than or equal to 1 mi2.
The smaller watersheds-those of drainage area less than 1 mi2 - were excluded because of their large variability of sediment yield, reflecting the diverse effects of soils, local terrain, vegetation, land use, and agricultural practices.
For drainage areas between 1 and 30,000 mi2, Dendy and Bolton found that the annual sediment yield per unit area was inversely related to the 0.16 power of the drainage area:
(15-10)
in which S = sediment yield in tons per square mile per year; SR = reference sediment yield corresponding to a 1-mi2 drainage area. equal to 1645 tons per year; A = drainage area in square miles; and AR = reference drainage area (1 mi2).
Sediment Yield versus Mean Annual Runoff.
Dendy and Bolton studied sedimentation data from 50S reservoirs having mean annual runoff data.
Annual sediment yield per unit area was shown to increase sharply as mean annual runoff Q increased from 0 to 2 in.
Thereafter, for mean annual runoff from 2 to 50 in. annual sediment yield per unit area decreased exponentially.
This led to the following equations.
(15-11a)
(15-11b)
in which QR = reference mean annual runoff, QR = 2 in.
Dendy and Bolton combined Eqs. 15-10 and 15-11 into a set of equations to vexpress sediment yield in terms of drainage area and mean annual runoff.
(15-12a)
(15-12b)
For SR = 1645 tons/mi2/y, QR = 2 in., and AR = 1 mi2, Eq. 15-12 reduces to the following:
For Q < 2 in.:
S = 1280Q¡.46(1.43 - 0.26 log A)
For Q ~ 2 in.:
S = 1965e-o.osSQ(l.43 - 0.26 log A)
(15-13a)
(l5-13b)
Equations 5-12 and 5-13 are based on average values of grouped data; therefore, they should be used with caution.
In certain cases, local factors such as soils, geology, topography, land use, and vegetation may have a greater influence on sediment yield than either mean annual runoff or drainage area.
Nevertheless, these equations provide a first approximation to the regional assessment of sediment yield for watershed planning purposes.
Example 15-4.
Calculate the sediment yield by the Dendy and Bolton formula for a 150-mi2 watershed with 3.5 in. of mean annual runoff.
The application of Eq. l5-13b leads to:
S = 1965 X e(-o.osS XJ,S) [1.43 - 0.26 log (150)] = 1400 ton/ mi2/y
Therefore, the sediment yield is 210,000 ton/y.
Other widely used sediment yield prediction formulas are the modified universal soil loss equation (MUSLE), developed by Williams [42], and the Flaxman formula [17].
Unlike the USLE, which is based on annual values, the MUSLE is intended for use with individual storms.
The Flaxman formula was developed using data from the western United States and is therefore particularly applicable to that region.
15.3 SEDIMENT TRANSPORT
Sediment transport refers to the entraInment and movement of sediment by flowing water.
An understanding of the principles of sediment transport is essential for the interpretation and solution of many hydraulic, hydrologic, and water resources engineering problems.
The study of sediment transport can be divided into (1) sediment transport mechanics, (2) sediment transport prediction, and (3) sediment routing.
Sediment transport mechanics refers to the fundamental processes by which sediment is entrained and transported by flowing water.
Sediment transport prediction refers to the methods and techniques to predict the equilibrium or steady rate of sediment transport in streams and rivers.
The prediction of sediment transport is accomplished by means of a sediment transport formula.
Sediment routing refers to the nonequilibrium or unsteady sediment transport processes, the net result of which is either the aggradation or degradation of stream and river beds.
The description of sediment transport is based on principles of fluid mechanics, river mechanics, and fluvial geomorphology.
The energy and turbulence of the flow gives streams and rivers the capacity to entrain and transport sediment.
The sediment being transported can originate in either (a) upland sources or (b) channel sources.
A significant feature of sediment transport is the entrainment and transport by the flow of the material constituting the channel bed.
Thus, sediment transport serves not only as the means for the movement of sediment from upstream to downstream but also as the mechanism by which streams and rivers determine their own crosssectional shapes and boundary roughness.
While the transport of sediment is a fluid mechanics subject, the interaction between flowing stream and its boundary is a river mechanics subject.
Sediment Transport Mechanics
Sediment load or sediment discharge is the total amount of sediment transported by a stream or river past a given point, expressed in terms of weight per unit time. Based on the predominant mode of transport, sediment load can be classified into (a) bed load and (b) suspended load.
Bed load is the fraction of sediment load that moves by saltation and rolling along the channel bed, primarily by action of bottom shear stresses caused by vertical velocity gradients.
Suspended load is the fraction of sediment load that moves in suspension by the action of turbulence.
Particles transported as bed load are coarser than particles transported as suspended load.
However, the distinction between bed load and suspended load is not all-exclusive; some particles may move as bed load at one point, as suspended load at another, and vice versa.
Based on whether the particle sizes are represented in the channel bed, sediment load can be classified into (a) bed-material load and (b) fine-material load.
Bed-material load is the fraction of sediment load whose particle sizes are significantly represented in the channel bed.
Conversely, fine-material load - commonly referred to as wash load-is the fraction of sediment load whose particle sizes are not significantly represented in the channel bed.
Stated in other terms, bed-material load is the coarser fraction of sediment load that may have originated in the channel bed and that may be subject to deposition under certain flow conditions.
Wash load is the finer fraction of sediment load that has not originated in the channel bed and that is not likely to deposit.
Wash load is then, washed through the reach, largely unaffected by the hydraulics of the flow.
The relationship between the two classifications of sediment load is shown in Fig. 15-5 [11].
This figure shows that the concepts of bed load and wash load are mutually exclusive.
The middle overlap is the suspended bed-material load, i.e., the fraction of sediment load that moves in suspended mode and is composed of particle sizes that are represented in the channel bed.
Initiation of Motion.
Water flowing over a streambed has a marked vertical velocity gradient near the streambed.
This velocity gradient exerts a shear stress on the particles lying on the streambed, i.e., a bottom shear stress.
For wide channels, the bottom shear stress can be approximated by the following formula [5]:
To = 'YdSo (15-14)
in which τ0 = bottom shear stress; γ = specific weight of water; d = flow depth; and So = equilibrium or energy slope of the channel flow.