MUSLE {VFS} | R Documentation |
Modified Universal Soil Loss Equation
Description
Simulation of soil erosion on a daily timestep.
Usage
MUSLE(Q, qp, A, K, LS, C = 0.085, P = 0.40, a = 11.8, b = 0.56)
Arguments
Q |
Runoff volume (m^3/d). |
qp |
Runoff peak discharge (m^3/s). |
A |
Field area (ha). |
K |
Soil erodibility factor. |
LS |
Landscape factor. |
C |
Crop management factor. Default is for a corn field. |
P |
Erosion control practice factor. |
a |
Location coefficient. Default value from Williams 1975. |
b |
Location coefficient. Default value from Williams 1975. |
Details
Uses the Modified Universal Soil Loss Equation to estimate daily sediment yield. If K and LS are not known, they can be estimated from soil or field properties using MUSLE.K
and MUSLE.LS
.
If the location coefficients are known from measured sedimentation data, more accurate estimates can be made.
Value
Sediment yield (t/day).
Author(s)
Sarah Goslee
References
Williams, J. R. (1975) Sediment-yield prediction with universal equation using runoff energy factor. Pp. 244–251 in: Present and prospective technology for predicting sediment yield and sources. ARS.S-40, US Gov. Print. Office, Washington, DC. 244-252.
Wischmeier, W. H. and Smith, D. D. (1978) Predicting rainfall erosion losses-a guide to conservation planning. U.S. Department of Agriculture, Agriculture Handbook No. 537.
See Also
Examples
# Approximate erodibility factor from soil texture.
Kf <- MUSLE.K(.3, .5, .2)
# Calculate landscape factor from field size and shape.
# 100-m field length with 2% slope
# note that MUSLE.LS takes feet
LS <- MUSLE.LS(100 * 3.28, .02)
# assume 0.4 ha cornfield with known rainfall intensity
peakd <- peak(intensity = 55, area = 0.4)
SedYield <- MUSLE(85, qp = peakd, A = .4, K = Kf, LS = LS)