| MUSLE | R Documentation |
Simulation of soil erosion on a daily timestep.
MUSLE(Q, qp, A, K, LS, C = 0.085, P = 0.40, a = 11.8, b = 0.56)
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. |
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.
Sediment yield (t/day).
Sarah Goslee
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.
MUSLE.K, MUSLE.LS, peak
# 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)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.