MMF-SAGA Soil Erosion Model
| Author(s) | (c) 2009-2012 V. Wichmann, M. Setiawan |
| Library ID | sim_erosion |
| Tool ID | 0 |
| Version | 1.0 |
| Menu | Simulation | Erosion |
Description
Soil erosion modelling with a modified MMF (Morgan-Morgan-Finney) model (Morgan & Duzant 2008).
This tool is called MMF-SAGA because some things have been implemented differently compared to the original publication. The most important are:
- the extension of the model to support spatially distributed modelling
- the introduction of a "Channel Network" layer
- the introduction of a "rainfall duration" (time span) parameter
- the exposure of the flow depth parameter
A more detailed description of the model, its modifications, and model application is provided by Setiawan (2012), chapter 6.
Currently, a number of additional grid datasets are outputted to facilitate model evaluation. This can be easily changed within the source code.
References
Morgan, R.P.C. (2001): A simple approach to soil loss prediction: a revised Morgan-Morgan-Finney model. Catena, 44, 305-322.
Morgan, R.P.C., Duzant, J.H. (2008): Modified MMF (Morgan-Morgan-Finney) model for evaluating effects of crops and vegetation cover on soil erosion. Earth Surf. Process. Landforms, 32, 90-106.
Setiawan, M. A. (2012): Integrated Soil Erosion Management in the upper Serayu Watershed, Wonosobo District, Central Java Province, Indonesia. Dissertation at the Faculty of Geo- and Atmospheric Sciences of the University of Innsbruck, Austria. PDF
Parameters
| | Name | Type | Identifier | Description | Constraints |
| Input | Digital Terrain Model | grid, input | DTM | DTM, digital terrain model [m] | - |
| Slope | grid, input | S | S, slope [rad] | - |
| Channel Network | grid, input, optional | CHANNEL | Channel network, all other cells NoData | - |
| Permanent Interception | grid, input | PI | PI, permanent interception expressed as the proportion [between 0-1] of rainfall | - |
| Canopy Cover | grid, input | CC | CC, canopy cover expressed as a portion [between 0-1] of the soil surface protected by vegetation or crop | - |
| Plant Height | grid, input | PH | PH, plant height [m], representing the effective height from which raindrops fall from the crop or vegetation | - |
| Ratio Evapotranspiration | grid, input | EtEo | Et/Eo, ratio of actual to potential evapotranspiration | - |
| Ground cover | grid, input | GC | GC, Ground cover expressed as a portion [between 0-1] of the soil surface protected by vegetation or crop cover on the ground | - |
| Diameter plant elements | grid, input | D | D, Average diameter [m] of the individual plants elements (stem, leaves) at the ground surface | - |
| Number plant elements | grid, input | NV | NV, Number of plant elements per unit area [number/unit area] at the ground surface | - |
| Soil moisture (at FC) | grid, input | MS | MS, Soil moisture at field capacity [% w/w] | - |
| Bulk density top layer | grid, input | BD | BD, Bulk density of the top layer [Mg/m3] | - |
| Effective hydrological depth | grid, input | EHD | EHD, Effective hydrological depth of the soil [m] | - |
| Sat. lateral permeability | grid, input | LP | LP, Saturated lateral permeability of the soil [m/day] | - |
| Percentage clays | grid, input | PER_C | c, Percentage clays [%] | - |
| Percentage silt | grid, input | PER_Z | z, Percentage silt [%] | - |
| Percentage sand | grid, input | PER_S | s, Percentage sand [%] | - |
| Percentage rock fragments | grid, input | ST | ST, Percentage rock fragments on the soil surface [%] | - |
| Surface roughness | grid, input, optional | RFR | RFR, Surface roughness [cm/m]. In case the surface roughness is not provided as input, v_flow_t is set to 1.0, i.e. natural soil surface roughness is not accounted for. | - |
| Meteorological data | table, input, optional | TAB_METEO | Meteorological data for multiple timestep modelling [model step (day); temperature (Celsius); rainfall (mm), rainfall intensity (mm/h); rainfall duration (day); timespan (days)] | - |
| Output | Interflow | grid, output, optional | IF | IF | - |
| Mean runoff | grid, output | Q | Q, estimation of mean runoff [mm] | - |
| Mean soil loss | grid, output | SL | SL, estimation of mean soil loss [kg] | - |
| Effective Rainfall | grid, output | Rf | Rf | - |
| Total Kinetic Energy | grid, output | KE | KE | - |
| Soil moisture storage capacity | grid, output | Rc | Rc | - |
| Transport Capacity Clay | grid, output | TCc | TCc | - |
| Transport Capacity Silt | grid, output | TCz | TCz | - |
| Transport Capacity Sand | grid, output | TCs | SLs | - |
| Available Clay | grid, output | Gc | Gc | - |
| Available Silt | grid, output | Gz | Gz | - |
| Available Sand | grid, output | Gs | Gs | - |
| Sediment Balance Clay | grid, output | SLc | SLc | - |
| Sediment Balance Silt | grid, output | SLz | SLz | - |
| Sediment Balance Sand | grid, output | SLs | SLs | - |
| Transport Condition Clay | grid, output | TCONDc | Sediment Limited [0], Transport Limited (SL = TC) [1], Transport Limited (SL = G) [2] | - |
| Transport Condition Silt | grid, output | TCONDz | Sediment Limited [0], Transport Limited (SL = TC) [1], Transport Limited (SL = G) [2] | - |
| Transport Condition Sand | grid, output | TCONDs | Sediment Limited [0], Transport Limited (SL = TC) [1], Transport Limited (SL = G) [2] | - |
| Upslope Flow Width | grid, output | W_up | W_up | - |
| Options | Grid System | grid system | PARAMETERS_GRID_SYSTEM | - | - |
| Output file path | file path | OUT_PATH | Full path to the directory for the output grids of each model step | - |
| Simulate Interflow | boolean | INTERFLOW | Simulate interflow | Default: 1 |
| Mean temperature | floating point number | T | T, mean temperature [degree C] | Default: 18.000000 |
| Timespan (days) | integer number | TIMESPAN | The number of days to model. | Minimum: 1
Maximum: 365
Default: 30 |
| Rainfall | floating point number | R | R, height of precipitation in timespan [mm] | Default: 200.000000 |
| Rainfall intensity | floating point number | I | I, rainfall intensity [mm/h] | Default: 20.000000 |
| Rainfall Duration | floating point number | Rn | Rn, number of rain days in timespan [-] | Default: 20.000000 |
| Relationship KE - I | choice | KE_I_METHOD | Relationship between kinetic energy (KE) and rainfall intensity (I) | Available Choices:
[0] North America east of Rocky Mountains (Wischmeier & Smith 1978)
[1] North-western Europe (Marshall & Palmer)
[2] Mediterranean-type climates (Zanchi & Torri 1980)
[3] Western Mediterranean (Coutinho & Tomas 1995)
[4] Tropical climates (Hudson 1965)
[5] Eastern Asia (Onaga et al. 1998)
[6] Southern hemisphere climates (Rosewell 1986)
[7] Bogor, West-Java, Indonesia (McISaac 1990)
Default: 0 |
| Flow Depth (actual flow velocity) | floating point number | FLOWD_VA | The flow depth used to calculate the actual flow velocity [m] (e.g. 0.005 unchannelled flow, 0.01 shallow rills, 0.25 deeper rills. | Minimum: 0.000000
Default: 0.005000 |
| Route Soil along Channel Network | boolean | CHANNELTRANSPORT | Route soil loss along channel network to outlet | Default: 0 |
Command Line
Usage: saga_cmd sim_erosion 0 [-DTM ] [-S ] [-CHANNEL ] [-PI ] [-CC ] [-PH ] [-EtEo ] [-GC ] [-D ] [-NV ] [-MS ] [-BD ] [-EHD ] [-LP ] [-PER_C ] [-PER_Z ] [-PER_S ] [-ST ] [-RFR ] [-TAB_METEO ] [-OUT_PATH ] [-INTERFLOW ] [-T ] [-TIMESPAN ] [-IF ] [-R ] [-I ] [-Rn ] [-KE_I_METHOD ] [-FLOWD_VA ] [-CHANNELTRANSPORT ] [-Q ] [-SL ] [-Rf ] [-KE ] [-Rc ] [-TCc ] [-TCz ] [-TCs ] [-Gc ] [-Gz ] [-Gs ] [-SLc ] [-SLz ] [-SLs ] [-TCONDc ] [-TCONDz ] [-TCONDs ] [-W_up ]
-DTM: Digital Terrain Model
grid, input
-S: Slope
grid, input
-CHANNEL: Channel Network
grid, input, optional
-PI: Permanent Interception
grid, input
-CC: Canopy Cover
grid, input
-PH: Plant Height
grid, input
-EtEo: Ratio Evapotranspiration
grid, input
-GC: Ground cover
grid, input
-D: Diameter plant elements
grid, input
-NV: Number plant elements
grid, input
-MS: Soil moisture (at FC)
grid, input
-BD: Bulk density top layer
grid, input
-EHD: Effective hydrological depth
grid, input
-LP: Sat. lateral permeability
grid, input
-PER_C: Percentage clays
grid, input
-PER_Z: Percentage silt
grid, input
-PER_S: Percentage sand
grid, input
-ST: Percentage rock fragments
grid, input
-RFR: Surface roughness
grid, input, optional
-TAB_METEO: Meteorological data
table, input, optional
-OUT_PATH: Output file path
file path
-INTERFLOW: Simulate Interflow
boolean
Default: 1
-T: Mean temperature
floating point number
Default: 18.000000
-TIMESPAN: Timespan (days)
integer number
Minimum: 1
Maximum: 365
Default: 30
-IF: Interflow
grid, output, optional
-R: Rainfall
floating point number
Default: 200.000000
-I: Rainfall intensity
floating point number
Default: 20.000000
-Rn: Rainfall Duration
floating point number
Default: 20.000000
-KE_I_METHOD:
Relationship KE - I
choice
Available Choices:
[0] North America east of Rocky Mountains (Wischmeier & Smith 1978)
[1] North-western Europe (Marshall & Palmer)
[2] Mediterranean-type climates (Zanchi & Torri 1980)
[3] Western Mediterranean (Coutinho & Tomas 1995)
[4] Tropical climates (Hudson 1965)
[5] Eastern Asia (Onaga et al. 1998)
[6] Southern hemisphere climates (Rosewell 1986)
[7] Bogor, West-Java, Indonesia (McISaac 1990)
Default: 0
-FLOWD_VA:
Flow Depth (actual flow velocity)
floating point number
Minimum: 0.000000
Default: 0.005000
-CHANNELTRANSPORT: Route Soil along Channel Network
boolean
Default: 0
-Q: Mean runoff
grid, output
-SL: Mean soil loss
grid, output
-Rf: Effective Rainfall
grid, output
-KE: Total Kinetic Energy
grid, output
-Rc: Soil moisture storage capacity
grid, output
-TCc: Transport Capacity Clay
grid, output
-TCz: Transport Capacity Silt
grid, output
-TCs: Transport Capacity Sand
grid, output
-Gc: Available Clay
grid, output
-Gz: Available Silt
grid, output
-Gs: Available Sand
grid, output
-SLc: Sediment Balance Clay
grid, output
-SLz: Sediment Balance Silt
grid, output
-SLs: Sediment Balance Sand
grid, output
-TCONDc: Transport Condition Clay
grid, output
-TCONDz: Transport Condition Silt
grid, output
-TCONDs: Transport Condition Sand
grid, output
-W_up: Upslope Flow Width
grid, output