DMMFc {DMMF} | R Documentation |
Daily based Morgan–Morgan–Finney (DMMF) soil erosion model
Description
Estimating surface runoff and sediment budget using the Daily based Morgan–Morgan–Finney soil erosion model using the algorithm from Choi et al. (2017)
Usage
DMMFc(DEM, R, RI, ET, P_c, P_z, P_s, theta_init, theta_sat, theta_fc, SD, K,
P_I, n_s, CC, GC, IMP, PH, D, NV, d_a, DK_c, DK_z, DK_s, DR_c, DR_z, DR_s,
SoilMap, LULCMap, Breaking, Init_point, Sinks, R_Type = 0, slpMode = 2, ALL = TRUE)
Arguments
DEM |
RasterLayer object of digital elevation model (DEM) |
R |
Numeric vector or RasterBrick object of daily rainfall per unit area (unit: mm) |
RI |
Numeric vector or RasterBrick object of mean hourly rainfall intensity during a day (unit: mm/h) |
ET |
Numeric vector or RasterBrick object of evapotranspiration per unit area (unit: mm) |
P_c |
Numeric vector or RasterLayer object of proportion of clay particles in surface soil (0-1) |
P_z |
Numeric vector or RasterLayer object of proportion of silt particles in surface soil (0-1) |
P_s |
Numeric vector or RasterLayer object of proportion of sand particles in surface soil (0-1) |
theta_init |
Numeric vector or RasterBrick object of initial soil water content of entire soil profile per unit area (unit: vol/vol) |
theta_sat |
Numeric vector or RasterLayer object of saturated soil water content of entire soil profile per unit area (unit: vol/vol) |
theta_fc |
Numeric vector or RasterLayer object of soil water content at field capacity of entire soil profile per unit area (unit: vol/vol) |
SD |
Numeric vector or RasterLayer object of soil depth of entire soil profile (unit: m) |
K |
Numeric vector or RasterLayer object of saturated lateral hydraulic conductivity of entire soil profile (unit: m/d) |
P_I |
Numeric vector or Raster* object of proportion of permanent interception area of rainfall (0-1) |
n_s |
Numeric vector or Raster* object of Manning's roughness coefficient of the soil surface (unit: |
CC |
Numeric vector or Raster* object of proportion of area with canopy cover (0-1) |
GC |
Numeric vector or Raster* object of proportion of area with pervious vegetated ground cover (0-1) |
IMP |
Numeric vector or Raster* object of proportion of area with impervious ground cover (0-1) |
PH |
Numeric vector or Raster* object of average height of vegetation or crop cover where leaf drainage start to fall (unit: m) |
D |
Numeric vector or Raster* object of average diameter of individual plant elements at the surface (unit: m) |
NV |
Numeric vector or Raster* object of number of individual plant elements per unit area (unit: |
d_a |
Numeric vector or Raster* object of typical flow depth of surface runoff (unit: m) |
DK_c |
Numeric vector or RasterLayer object of detachability of clay particles by rainfall (unit: g/J) |
DK_z |
Numeric vector or RasterLayer object of detachability of silt particles by rainfall (unit: g/J) |
DK_s |
Numeric vector or RasterLayer object of detachability of sand particles by rainfall (unit: g/J) |
DR_c |
Numeric vector or RasterLayer object of detachability of clay particles by surface runoff (unit: g/mm) |
DR_z |
Numeric vector or RasterLayer object of detachability of silt particles by surface runoff (unit: g/mm) |
DR_s |
Numeric vector or RasterLayer object of detachability of sand particles by surface runoff (unit: g/mm) |
SoilMap |
Categorical raster map for soil characteristics |
LULCMap |
Categorical raster map for Land use and vegetation types |
Breaking |
Numeric vector of starting points of each simulation period (0: during the period points, 1: breaking points), If |
Init_point |
Numeric vector of starting points of each rainfall event (0: during the period points, 1: breaking points), If |
Sinks |
RasterLayer object of sinks (e.g., streams, reservoirs, and lakes) of sediments and surface runoff (optional). |
R_Type |
Integer object of each rainfall type for estimating kinetic energy of direct throughfall (DT) (0-8) (see details) |
slpMode |
Integer object of slope calculation algorithm (2: second order, 3: third order) |
ALL |
Logical object of selection for the entire output results |
Details
R_Type
can be chosen among options appropriate for each regions according to Morgan (2001) and Morgan (2005). Default option adopts recent universal relationships between kinetic energy density and mean hourly rainfall intensity (RI
) from Shin et al. (2016). Default option is recommended when users have measured RI
. Detailed descriptions are described below.
Option | Region | Kinetic energy density of throughfall [\mathrm{J/m^{2}/mm} ] |
0 | Universal (Default) | 10.3 \cdot \mathrm{RI}^{2/9} |
1 | North America | 11.87 + 8.73 \cdot \log_{10}(\mathrm{RI}) |
2 | North-west Europe | 8.95 + 8.44 \cdot \log_{10}(\mathrm{RI}) |
3 | Mediterranean | 9.81 + 11.25 \cdot \log_{10}(\mathrm{RI}) |
4 | West Mediterranean | 35.9 \cdot (1.0 - 0.56 \cdot \exp(-0.034 \cdot \mathrm{RI})) |
5 | Tropics | 29.8 - ( 127.5 / \mathrm{RI} ) |
6 | East Asia | 9.81 + 10.60 \cdot \log_{10}( \mathrm{RI} ) |
7 | Temperate Southern hemisphere | 29.0 \cdot ( 1.0 - 0.6 \cdot \exp( -0.04 \cdot \mathrm{RI} ) ) |
8 | Universal | 28.3 \cdot ( 1.0 - 0.52 \cdot \exp( -0.042 \cdot \mathrm{RI} ) )
|
Multi-flow direction algorithm (MD\infty
) from Seibert and McGlynn (2007) is built in the DMMF
model.
Value
The output of the function DMMF
is a list of RasterLayer or RasterBrick objects containing the following elements:
-
A
: Surface area size of each element considering slope (unit:\mathrm{m^2}
) -
Rf
: The amount of effective rainfall reaching each element (unit:\mathrm{mm}
) -
SW_c
: Surface water infiltration capacity of each element (unit:\mathrm{mm}
-
theta_r
: Remaining soil water content of each element (unit: vol/vol) -
TC
: Transport capacity of the runoff of each element (unit:\mathrm{kg/m^2}
) -
Q_in
: Volume of surface runoff flowing into each element (unit: L) -
Q_out
: Volume of surface runoff flowing from each element (unit: L) -
IF_in
: Volume of subsurface water flowing into each element (unit: L) -
IF_out
: Volume of subsurface water flowing from each element (unit: L) -
SS_c
: Area density of clay delivered to surface runoff in unit area of each element (unit:\mathrm{kg/m^2}
) -
SS_z
: Area density of silt delivered to surface runoff in unit area of each element (unit:\mathrm{kg/m^2}
) -
SS_s
: Area density of sand delivered to surface runoff of each element (unit:\mathrm{kg/m^2}
) -
G_c
: Area density of available clay for transport by surface runoff of each element (unit:\mathrm{kg/m^2}
) -
G_z
: Area density of available silt for transport by surface runoff of each element (unit:\mathrm{kg/m^2}
) -
G_s
: Area density of available sand for transport by surface runoff of each element (unit:\mathrm{kg/m^2}
) -
SL_c_in
: Mass of clay inputs into each element (unit: kg) -
SL_z_in
: Mass of silt inputs into each element (unit: kg) -
SL_s_in
: Mass of sand inputs into each element (unit: kg) -
SL_in
: Mass of sum of total sediments inputs into each element (unit: kg) -
SL_c_out
: Mass of clay outputs from each element (unit: kg) -
SL_z_out
: Mass of silt outputs from each element (unit: kg) -
SL_s_out
: Mass of sand outputs from each element (unit: kg) -
SL_out
: Mass of sum of total sediments outputs from each element (unit: kg)
Author(s)
Kwanghun Choi and Bjoern Reineking
References
Choi, K., Arnhold, S., Huwe, B., Reineking, B. (2017). Daily based Morgan–Morgan–Finney (DMMF) model: A spatially distributed conceptual soil erosion model to simulate complex soil surface configurations. Water, 9(4), 278.
Shin, S. S., Park, S. D., and Choi, B. K. (2016). Universal power law for relationship between rainfall kinetic energy and rainfall intensity. Advances in Meteorology, Article ID 2494681, 11 pages.
Seibert, J., McGlynn, B. L. (2007). A new triangular multiple flow direction algorithm for computing upslope areas from gridded digital elevation models. Water Resources Research, 43(4):W04501.
Morgan, R. P. C. (2005). Soil erosion and conservation. Blackwell Publishing, Malden, MA, 3rd ed. ISBN 1-4051-1781-8
Morgan, R. P. C. (2001) A simple approach to soil loss prediction: a revised Morgan–Morgan–Finney model. Catena, 44(4):305–322.
See Also
DMMF
for general use of the algorithm.
Examples
## Not run:
## Datasets for DMMFc will be uploaded.
## End(Not run)