Northern Prairie Wildlife Research Center

Civil Engineering and Environmental Science, University of

Oklahoma, Norman, OK 73019;

U.S. Army Corps of Engineers Waterways Experiment Station,

Vicksburg, MS 39180

Data have been collected to assess the different rates of sediment deposition resulting from drainage areas consisting of native prairie, Conservation Reserve Program land (CRP), buffered wetlands, and tilled watersheds. Characteristics of the sediment, statistical analysis of the sediment for 1993, and other data reported by the U.S. Environmental Protection Agency (EPA) (1) indicate that wetlands with drainage areas consisting of tilled cropland have a greater sediment deposition than CRP cropland.

**Rainfall and Runoff**

Rainfall and runoff provide the basic driving energy for the erosion and
transport processes. In stream channels where the bulk of sediment transport
relationships have been developed, particles are continuously exchanged between
the bed and the flow. In upland erosion, this does not usually occur because
the shear required to detach the particle is much more than that required
for transport. Foster and Huggins (4) reported that the critical shear stress
required to detach soil particles is around 2.9 N/m^{2} whereas transport
only requires around 0.5 N/m^{2} for a sand of 342 µm. Rainfall impact
on the soil surface is necessary for the detachment of particles, and, therefore,
the supply of sediment. Modeling of upland erosion requires the simulation
of supply rates and transport to obtain an estimate of sediment transport
at a particular location in the landscape.

**Topography**

Topography affects the rate of runoff through slope and the length of slope. As slope increases so does the transport capacity. As length increases at a given slope, more drainage area contributes more runoff volume and rates providing more transport capacity. As soil particles are detached in the upper reaches of the watershed and transported down slope to the wetlands, the relative magnitudes of supply and transport rates govern the amount of sediment transported. Slope of the landsurface plays a key role in this sub-process. As slope decreases closer to the wetlands, it is expected that the transport capacity will fall below the supply rate. In the near-wetland part of the landscape, deposition has been observed and is explained by the decrease in transport capacity due to flatter slopes.

**Soils**

Soils have different susceptibilities to erosion forces. Thus, soil type and the management of the soil influences the rate at which soil may be detached by rainfall impact (5). Soil transportability is the ease with which the detached soil particles are transported by the water. Soil is eroded as primary particles (e.g., sand, silt, and clay) or as aggregates of different particle sizes and organic matter. Continual, intensive tillage breaks down the soil structure and depletes organic matter enhancing erodibility. Sediment produced on uplands as opposed to sediment flowing in rivers is a mixture of primary particles and aggregates (6). The proportion of primary particles and aggregates depends on upslope sorting, soil properties, and the amounts coming from rill and interrill areas. Soil particle sizes affect transport and deposition rates. Clays, if they are primary particles, may never be deposited due to resuspension caused by turbulence; whereas silts and sands are the first to be deposited as transport capacity decreases below the supply rate.

**Cover**

Cover, including plant canopy, mulches, plant residue, or densely growing plants, has a greater influence on erosion than any other factor. Cover, in contact with the soil, greatly reduces rill erosion by cushioning raindrop impact and reducing the shear stress exerted by flowing water. Strips of dense mulches or grasses can induce deposition and filter sediment from the runoff. Tillage type, frequency, and elapsed time between tillage and a runoff event affect the detachment of soil particles.

Given these factors that affect the detachment of soil particles (supply) and transport of sediment, the overall goal of simulation is not only to adequately represent the processes but to provide realistic estimates of the impacts of anthropogenic effects on wetlands. We were interested in both the capacity and the supply-limited rates because as runoff from the watershed divide traverses the topography, either control may dominate the transport of sediment. From preliminary analysis of sediment taken from the wetlands, CRP treated watersheds yield less sediment than those that were tilled (1). This indicates that the rate of sediment deposition is dependent on the supply which in turn is affected by tillage and land management practices.

where, qt is the mass transport rate in metric tons per meter of width and
seconds; So is the land surface slope; and q is the unit width flow rate in
m^{2}/s. This equation can be further modified (8) to account for
the soil erodibility, cropping and tillage effects, and conservation practices
as follows

where, K is the soil erodibility, C is the crop management factor, and P is the conservation practice factor. The factors K, C, and P are readily available for specific soils, crop management, and conservation practices. Further, K, C, and P represent the supply-limiting characteristics of the watershed.

Foster and Meyer (17) compared the sediment transport predicted by the Yalin equation to laboratory and erosion plot studies and found excellent agreement. Erosion of natural agricultural soils has been observed to occur by detachment and saltation of soil aggregates due to lifting forces induced by the flowing water. The Yalin equation assumes that sediment motion begins when the lift force of flow exceeds a critical lift force. The particle is transported downstream until the particle weight forces it out of the flow and back in the bed. The number of particles in motion at a given time is a linear function of an excess shear parameter, . The required hydrodynamic parameters are energy slope, and the hydraulic radius, R. The transportability of the soil is determined by the particle (aggregate or primary) specific gravity, the diameter, and the critical lift force, Ycr, given by the Shields diagram which defines incipient motion based on particle and fluid characteristics. The Yalin equation defines a dimensionless transport rate, (equation 3), in which Ws is transport rate in Newtons per second per unit width; SG is the specific gravity; w is the mass density of the water; d is the diameter in meters; g is the acceleration of gravity in m/s2; the shear velocity; V* is as in equation 4;

The excess shear parameter, in equation (3) is

where the lift force, Y is

The critical lift force Ycr is defined as a function of a dimensionless particle diameter, D*

The Shields diagram, usually presented in graphical form, may be expressed as a piece-wise continuous set of functions:

In equation (3) the value of = A * where A is

This defines each of the terms in equation 3. The resulting value of is dimensionless, whereas, the value of Ws computed with parameters in the Meters-Kilogram system, yields the capacity-limited transport rate in units of Newtons per second per unit width [N/(m s)].

To illustrate the transport computations for the W13E1 watershed, assumed
conditions with respect to the supply-limited conditions (parameters K, C,
and P) are made. Two conditions were modeled for CRP and tilled treatments
with the following parameters. The resulting sediment transport is supply-limited
and was the maximum amount that could be transported from the drainage area.
So= 0.02 K= 0.42 Silt Loam with 2% OM C_{CRP}= 0.24 No canopy and
20% cover by undecayed residues. C_{tilled}= 0.45 No canopy and 0%
cover by undecayed residues. P_{CRP}= 0.5 Farming on contour P_{tilled}=
1 No conservation practices.

The parameters chosen are typical values but should be modified for each soil mapping unit, land use/management, and cropping plan in each drainage area.

Constants:

g= 9.80 m/s2 v_{m} 8.94e-07 m^{2}/s P_{w}= 1000.00
kg/m3 Y= 9800 N/m3

Sediment Characteristics:

d_50= 3.30e-05 m Sg= 2.65 D*= 0.90 Ycr= 0.27 A= 0.86 The time-dependent supply-limited transport rate is shown in Figure 2 with the time-dependent capacity-limited transport rate. Note that the capacity-limited transport rate is less than both the CRP and tilled treatment transport rates. This indicates that the capacity-limited transport rate governs. The estimated transport rate is, therefore, the lesser of the two (capacity versus supply); the capacity-limited transport rate. Meyer and Wischmeier (3) found from erosion plot studies that for slopes of less than 3%, transport capacity was the limiting rate. The segment modeled in this example has an average slope of 2%.

Under the conditions analyzed for the W13E1 pothole, both supply-limited and capacity-limited transport may dominate the sediment transport rate. Particularly on flat slopes encountered near the emergent vegetation of the wetland, sediment begins to drop out due to reduced transport capacity. The supply of sediment is also affected by the land treatment. CRP treated watersheds show less sediment transport than tillage without conservation treatment. Supply-limited transport rates for CRP conditions are 33% of the tilled conditions shown in Figure 2. The capacity-limited rate using the Yalin equation shows that it is controlling at all times during the storm for the flow rates encountered at 75 meters downslope from the watershed divide. Further analysis is needed to determine if this is true along the land surface profile from the watershed divide to the edge of the wetland.

Return to Contents