5.4 Modified USLE (MUSLE) & Single-Storm Sediment Yield

5.4 Modified USLE (MUSLE) & Single-Storm Sediment Yield

Quick Reference: While USLE and RUSLE predict average annual gross hillslope detachment ($A$, tons/ac/yr), civil and environmental engineers must size sediment basins, traps, and impoundments for discrete design storm events (e.g., the 2-year or 10-year, 24-hour storm). In 1975, Jimmy R. Williams developed the Modified Universal Soil Loss Equation (MUSLE) by replacing the rainfall energy factor ($R$) with a runoff energy factor: $Y = 95 \times (Q \times q_p)^{0.56} \times K \times LS \times C \times P$. Here, $Y$ is the single-storm delivered sediment yield in tons, $Q$ is runoff volume in acre-feet, and $q_p$ is peak discharge in cfs. Because MUSLE is driven directly by runoff energy, it models delivered sediment at a watershed outfall without requiring a separate Sediment Delivery Ratio ($SDR$).

---

Fundamental Limitations of Annual USLE / RUSLE

While USLE and RUSLE serve as indispensable tools for long-term conservation planning and SWPPP erosion risk modeling, they exhibit three critical limitations when applied to civil construction engineering:

1. Average Annual Time Horizon: USLE/RUSLE yields an annualized erosion rate based on multi-decade meteorological averages. However, construction permits (NPDES Construction General Permits) and municipal engineering codes require designing temporary sediment control structures (sediment basins, sediment traps, sediment barriers) to withstand specific regulatory design storms (e.g., the 10-year, 24-hour storm event).
2. Rainfall Energy vs. Runoff Realities: In USLE, erosion is driven by rainfall kinetic energy ($R = EI_{30}$). However, rainfall does not equal runoff. A massive, high-intensity rainstorm falling upon highly permeable, dry sandy soil or an unpaved gravel pad may generate substantial raindrop splash detachment, but virtually zero surface runoff ($Q = 0$). In such cases, detached particles simply settle in place a few inches away. Conversely, a moderate rainstorm falling upon saturated clay subsoils or an extensive paved commercial pad generates massive runoff volume and high peak velocity, scouring and transporting enormous sediment volumes off-site.
3. Gross Detachment vs. Delivered Sediment Yield: USLE predicts gross detachment across hillslope surfaces. It cannot determine how much of that detached soil actually travels through the drainage network to reach a property boundary, a sediment basin inlet, or a receiving stream outfall.

---

The Williams (1975) Modified USLE (MUSLE) Formulation

To overcome these limitations for hydrological engineering, Jimmy R. Williams (USDA-ARS, Temple, Texas) developed the Modified Universal Soil Loss Equation (MUSLE) in 1975. Williams replaced the rainfall erosivity factor ($R$) with an empirical runoff energy term that combines total runoff volume and peak runoff rate.

The MUSLE Equation (U.S. Customary Units)

$Y = 95 \times (Q \times q_p)^{0.56} \times K \times LS \times C \times P$

Where:

• $Y$ = Single-Storm Sediment Yield: Delivered to the watershed or catchment outlet for the specific design storm event, expressed in tons (U.S. short tons = $2,000\text{ lbs}$). (Note: $Y$ represents total mass of delivered sediment, not tons per acre per year).
• $Q$ = Total Runoff Volume: Generated by the contributing drainage basin for the design storm, expressed in acre-feet ($ ext{ac-ft}$). One acre-foot equals $43,560\text{ cubic feet}$.
• $q_p$ = Peak Discharge Rate: The crest runoff flow rate exiting the watershed outfall for the design storm, expressed in cubic feet per second ($ ext{cfs}$). Typically calculated using the NRCS TR-55 Graphical Peak Discharge method or HEC-HMS hydrologic modeling.
• $K$ = Soil Erodibility Factor: Area-weighted average soil erodibility factor for the contributing watershed (same standard USLE units).
• $LS$ = Topographic Factor: Representative topographic length and steepness factor for the contributing drainage area.
• $C$ = Cover Management Factor: Area-weighted cover factor representing site surface stabilization conditions during the design storm.
• $P$ = Support Practice Factor: Area-weighted conservation/support practice factor.
• $95$ = Empirical Calibration Coefficient: Dimensional conversion constant for U.S. customary units.

Metric (SI) Formulation

In international and metric engineering applications, the MUSLE equation is written as:

$Y_{\text{metric}} = 11.8 \times (Q_{\text{m}^3} \times q_{\text{m}^3/\text{s}})^{0.56} \times K \times LS \times C \times P$

Where $Y$ is delivered sediment yield in metric tonnes, $Q$ is runoff volume in cubic meters ($ ext{m}^3$), and $q_p$ is peak discharge in cubic meters per second ($ ext{m}^3/\text{s}$). The coefficient $11.8$ reflects the dimensional unit conversion from the English coefficient $95$.

---

Physical Meaning of the Runoff Energy Term $(Q \times q_p)^{0.56}$

The runoff energy factor fundamentally improves sediment prediction because flowing surface water performs two distinct physical tasks:

1. Runoff Volume ($Q$, acre-feet): Quantifies the total volumetric mass of water traversing the drainage network. A large runoff volume sustains hydraulic flow depth and flow duration, preventing suspended silt and fine sand particles from settling out along the flow path.
2. Peak Discharge Rate ($q_p$, cfs): Dictates maximum flow depth, flow velocity, and peak tractive shear stress ($\tau = \gamma R_h S$) exerted on the soil surface at the crest of the runoff hydrograph. Peak velocity determines the maximum particle size that runoff can detach and transport without deposition.
3. The Interactive Product $(Q \times q_p)$: Combines total flow mass and peak hydraulic power. The exponent $0.56$ reflects the non-linear hydraulic scaling between flow energy and sediment transport capacity.
4. Zero Runoff Boundary Condition: If a storm produces no surface runoff ($Q = 0$), the runoff energy term evaluates to zero ($95 \times (0)^{0.56} = 0$), yielding $Y = 0\text{ tons}$. Unlike USLE, MUSLE correctly predicts zero sediment yield when no runoff occurs, regardless of rainfall kinetic energy.

---

Sediment Delivery Ratio (SDR) & Deposition Mechanics

A central concept in watershed geomorphology and CPESC exam practice is the Sediment Delivery Ratio ($SDR$).

Gross Hillslope Erosion (A) ──► [ Upland Deposition in Swales, Buffers, Flats ] ──► Delivered Sediment Yield (Y)
                                    SDR = Y / A  (Typically 0.05 to 0.85)

Mathematical Definition of SDR

The Sediment Delivery Ratio ($SDR$) is defined as the dimensionless ratio of delivered sediment yield ($Y$) measured at a designated watershed outlet or property line to the total gross upland hillslope erosion ($A$) occurring across the entire contributing watershed:

$SDR = \frac{Y}{A}$

$SDR$ is expressed as a decimal ($0.0\text{ to } 1.0$) or as a percentage ($0%\text{ to } 100%$).

Deposition Along Upland Flow Paths

Gross erosion measures every grain of soil detached on a hillside. However, as sediment-laden sheet flow moves downslope, it encounters changes in terrain. Runoff velocity decelerates whenever:

• Slope gradient flattens (e.g., at the toe of a cut slope or in broad concave swales);
• Surface hydraulic roughness ($n$) increases (e.g., entering a dense vegetated buffer strip);
• Runoff spreads out across wide, shallow overbank areas.

When flow velocity drops below the settling velocity ($v_s$) of a soil particle (governed by Stokes' Law), the particle drops out of suspension and deposits on the land surface. Consequently, only a fraction of detached gross erosion ever reaches the watershed outlet.

Watershed Drainage Area Relationship

Empirical research by the USDA-SCS and Vito A. Vanoni demonstrated that $SDR$ decreases systematically as total drainage basin area ($DA$) increases:

$SDR = \alpha \times (DA)^{-\beta} \quad \text{commonly expressed as} \quad SDR = 0.41 \times (DA)^{-0.28}$

Where $DA$ is contributing drainage area in square miles (or acres).

Other Factors Governing SDR

• Relief-to-Length Ratio ($R/L$): Watersheds with steep vertical relief and short mainstream channel lengths maintain high flow velocities and exhibit significantly higher $SDR$ values.
• Drainage Density: Watersheds with highly incised, dense networks of drainage gullies and storm pipes deliver runoff directly to channels without overland filtration, increasing $SDR$.
• Sediment Particle Size Distribution: Coarse sands and gravels deposit almost immediately ($SDR < 0.10$), whereas colloidal clays and fine silts remain in suspension for days, exhibiting an $SDR$ approaching $1.00$.

Core Distinction for the Exam: When using USLE/RUSLE, you must manually multiply gross soil loss ($A$) by an empirical $SDR$ to estimate delivered sediment: $Y = A \times SDR$. When using MUSLE, you do NOT multiply by an $SDR$ because the runoff energy factor ($95(Q \cdot q_p)^{0.56}$) already accounts for watershed transport efficiency and upland deposition!

---

Sediment Basin Sizing & Trap Efficiency (\eta)

Under NPDES Phase I and Phase II regulations, temporary sediment basins are mandated for construction drainage areas disturbing $10\text{ or more acres}$. Sizing these structures requires matching basin volume to incoming design storm sediment yield.

Basin Trap Efficiency (\eta)

The Trap Efficiency (\eta) of a sediment basin is the percentage of total incoming sediment yield ($Y_{ ext{in}}$) that settles out and is retained within the basin pool:

$\eta = \left( \frac{Y_{\text{trapped}}}{Y_{\text{in}}} \right) \times 100$

Trap efficiency is governed by the Camp-Hazen sedimentation model, which depends on the ratio of particle settling velocity ($v_s$) to basin surface overflow rate ($q / A_s$), the basin length-to-width ratio (minimum $2:1$, recommended $4:1$ or greater utilizing porous baffles), and the dead storage volume.

Net sediment escaping the basin outfall and entering receiving waters is:

$Y_{\text{discharged}} = Y_{\text{in}} \times (1 - \eta)$

Converting Sediment Mass ($Y$, tons) to Storage Volume ($V_{\text{sed}}$, cubic yards)

To size the physical sediment storage zone (cleanout volume) of an engineered basin, the CPESC practitioner must convert sediment yield mass ($Y$, tons) into physical settled volume:

$V_{\text{sed}} (\text{cu ft}) = \frac{Y (\text{tons}) \times 2,000\text{ lbs/ton}}{\gamma_{\text{dry}} (\text{lbs/cu ft})}$

$V_{\text{sed}} (\text{cu yds}) = \frac{V_{\text{sed}} (\text{cu ft})}{27\text{ cu ft/cu yd}}$

Where $\gamma_{\text{dry}}$ is the settled dry bulk density of deposited sediment, which typically ranges from $65\text{ to } 95\text{ lbs/cu ft}$ (commonly assumed as $75\text{ lbs/cu ft}$ for mixed construction subsoils, or roughly $1.01\text{ tons per cubic yard}$).

---

Step-by-Step Worked Calculation: MUSLE Single-Storm Sediment Yield & Basin Sizing

Design Problem Scenario

A civil grading contractor is developing a $12.0\text{-acre}$ commercial logistics park in suburban Atlanta, Georgia. Local drainage ordinances require the on-site temporary sediment basin to fully contain the sediment volume generated by a 10-year, 24-hour regulatory design storm without exceeding the designated sediment cleanout elevation.

• Site Hydrologic Modeling (NRCS TR-55 Method):
  • Design Storm: 10-year, 24-hour rainfall depth = $5.20\text{ inches}$
  • Total Runoff Volume: $Q = 2.50\text{ acre-feet}$
  • Peak Runoff Discharge Rate: $q_p = 22.0\text{ cfs}$
• Physical Catchment Parameters:
  • Area-Weighted Soil Erodibility: $K = 0.28$ (sandy clay loam subsoil)
  • Representative Topographic Factor: $LS = 2.10$ ($150\text{ ft}$ average slope length at $4:1$ [$25%$] gradient)
  • Cover Management Factor: $C = 0.35$ (active grading with temporary hydromulch applied to $40%$ of perimeter areas)
  • Support Practice Factor: $P = 0.90$ (contour surface roughening / track-walking)
• Settled Sediment Dry Density: $\gamma_{\text{dry}} = 75.0\text{ lbs/cu ft}$

Step 1: Calculate the Runoff Energy Term

Compute the interactive product of runoff volume and peak discharge:

$(Q \times q_p) = 2.50\text{ ac-ft} \times 22.0\text{ cfs} = 55.0$

Raise the product to the $0.56$ power:

$(55.0)^{0.56} = 9.4074$

Multiply by the empirical coefficient $95$:

$\text{Runoff Energy Factor} = 95 \times 9.4074 = 893.70$

Step 2: Compute the Product of Site Specific Factors ($K \times LS \times C \times P$)

$K \times LS = 0.28 \times 2.10 = 0.588$ $C \times P = 0.35 \times 0.90 = 0.315$ $\text{Site Factor Product} = 0.588 \times 0.315 = 0.18522$

Step 3: Compute Delivered Single-Storm Sediment Yield ($Y$)

$Y = 893.70 \times 0.18522 = 165.53\text{ tons of sediment}$

The 10-year, 24-hour storm will deliver approximately $165.5\text{ tons}$ of eroded sediment into the temporary sediment basin.

Step 4: Convert Sediment Mass to Basin Storage Volume ($V_{\text{sed}}$)

Convert sediment tonnage into pounds:

$\text{Total Mass} = 165.53\text{ tons} \times 2,000\text{ lbs/ton} = 331,060\text{ lbs}$

Compute volume in cubic feet using the dry bulk density of $75.0\text{ lbs/cu ft}$:

$V_{\text{sed}} = \frac{331,060\text{ lbs}}{75.0\text{ lbs/cu ft}} = 4,414.13\text{ cu ft}$

Convert cubic feet to cubic yards:

$V_{\text{sed}} = \frac{4,414.13\text{ cu ft}}{27\text{ cu ft/cu yd}} = 163.49\text{ cu yds} \approx 164\text{ cu yds}$

Engineering Conclusion & Design Specification

To ensure compliance with local sediment retention criteria and prevent premature basin failure during a 10-year storm, the CPESC engineer must specify:

1. A dedicated sediment storage dead pool volume of at least $4,415\text{ cubic feet}$ ($164\text{ cubic yards}$) below the principal spillway crest or skimmer invert;
2. A clear sediment cleanout marker stake positioned at the $50%$ storage volume elevation ($82\text{ cubic yards}$ / $2,208\text{ cu ft}$), mandating mechanical dredging whenever accumulated sediment reaches this marker.