-- -
NAME
-r.landscape.evol2.4.1.py - Simulates the cumulative effect of erosion and deposition on a landscape over time. -KEYWORDS
-raster, hydrology, erosion modeling, landscape evolution -SYNOPSIS
--
Flags:
--
-
- -m -
- -m Apply smoothing (useful to mitigate possible unstable conditions in streams) - -
- -p -
- -p Run a sampling procedure to generate a vector points map with scaled flow accumulation values suitable for determining transport equation thresholds. Overrides all other output. - -
- -k -
- -k Keep ALL temporary maps (overides flags -drst). This will make A LOT of maps! - -
- -d -
- -d Don't output yearly soil depth maps - -
- -r -
- -r Don't output yearly maps of the erosion/deposition rates ("ED_rate" map, in vertical meters) - -
- -s -
- -s Keep all slope maps - -
- -t -
- -t Keep yearly maps of the Transport Capacity at each cell ("Qs" maps) - -
- -e -
- -e Keep yearly maps of the Excess Transport Capacity (divergence) at each cell ("DeltaQs" maps) - -
- --overwrite -
- Allow output files to overwrite existing files -
- --help -
- Print usage summary -
- --verbose -
- Verbose module output -
- --quiet -
- Quiet module output -
- --ui -
- Force launching GUI dialog -
Parameters:
--
-
- elev=name [required] -
- Input elevation map (DEM of surface) - -
- initbdrk=name [required] -
- Bedrock elevations map (DEM of bedrock) -
- Default: - -
- prefx=basename [required] -
- Name for output basename raster map(s) -
- Default: levol_ - -
- outdem=basename [required] -
- Name stem for output elevation map(s) (preceded by prefix and followed by numerical suffix if more than one iteration) -
- Default: elevation - -
- outsoil=basename [required] -
- Name stem for the output soil depth map(s) (preceded by prefix and followed by numerical suffix if more than one iteration) -
- Default: soildepth - -
- number=integer [required] -
- Number of iterations (cycles) to run -
- Default: 1 - -
- k=name -
- Soil erodability index (K factor) map or constant (values <= 0.09 [t.ha.h /ha.MJ.mm]) -
- Default: 0.05 - -
- c=name -
- Landcover index (C factor) map or constant (values <=1.0 [unitless]) -
- Default: 0.005 - -
- p=name -
- Landuse practices factor (P factor) map or constant (values <=1.0 [unitless]) -
- Default: 1.0 - -
- sdensity=name -
- Soil density map or constant for conversion from mass to volume (values typically >=1000 [kg/m3]) -
- Default: 1218.4 - -
- transp_eq=string -
- The sediment transport equation to use (USPED: Tc=R*K*C*P*A^m*B^n, Stream power: Tc=Kt*gw*1/N*h^m*B^n, or Shear stress: Tc=Kt*tau^m ). -
- Options: StreamPower, ShearStress, USPED -
- Default: StreamPower +
- exp_m=string -
- Exponent m relates to the influence of upslope area (and thus flow depth, discharge) on transport capacity. Values generally thought to scale inversely with increasing depth of flow between the two cutoff thresholds specified: "thresh1,m1,thresh2,m2" -
- Default: 10,2,100,1 - -
- exp_n=string -
- Exponent n relates to the influence of local topographic slope on transport capacity. Default values set to scale inversely with increasing local slope between the two slope cutoff thresholds specified: "thresh1,n1,thresh2,n2" -
- Default: 10,2,45,0.5 +
- r=name [required] -
- Rainfall (R factor) map or constant (Employed only in the USPED equation) (values typically between 500 and 10000 [MJ.mm/ha.h.yr]) -
- Default: 720 +
- rain=name [required] -
- Precip total for the average erosion-causing storm map (Employed in stream power and shear stress equations) (values typically >=30.0 [mm]) -
- Default: 30 +
- storms=name [required] -
- Number of erosion-causing storms per year map or constant (Employed in stream power and shear stress equations) (values >=0 [integer]) -
- Default: 2 +
- stormlength=name [required] -
- Average length of the storm map or constant (Employed in stream power and shear stress equations) (values >=0.0 [h]) -
- Default: 24.0 - -
- stormi=name [required] -
- Proportion of the length of the storm where the storm is at peak intensity map or constant (Employed in stream power and shear stress equations) (values typically ~0.05 [unitless proportion]) -
- Default: 0.05 - -
- climfile=string -
- Path to climate file of comma separated values of "rain,R,storms,stormlength,stormi", with a new line for each year of the simulation. This option will override values or maps entered above. - -
- manningn=name -
- Map or constant of the value of Manning's "N" value for channelized flow. (Employed in stream power and shear stress equations) (0.03 = clean/straight stream channel, 0.035 = major river, 0.04 = sluggish stream with pools, 0.06 = very clogged streams [unitless]) -
- Default: 0.03 - -
- flowcontrib=name -
- Map or constant indicating how much each cell contributes to downstream flow (this typically relates to vegetation or conservation practices). If no map or value entered, routine will assume 100% downstream contribution (values between 0 and 100 [unitless percentage]) -
- Default: 100 - -
- convergence=integer -
- Value for the flow convergence variable in r.watershed. Small values make water spread out, high values make it converge in narrower channels. -
- Options: 1, 2, 3, 4, 5, 6, 7, 8, 9, 10 -
- Default: 5 - -
- statsout=name -
- Name for the statsout text file (optional, if none provided, a default name will be used) - -
- Tc = Transport Capacity [kg/meters.second] -
- K*C*P ~ Kt = mitigating effects of soil type, vegetation cover, and land-use practices. [unitless] -
- gw = Hydrostatic pressure of water 9810 [kg/m2.second] -
- N = Manning's coefficient (0.3-0.6 for different types of stream channels) [unitless] -
- i = rainfall intensity [m/rainfall event] -
- A = upslope accumulated area per contour (cell) width [m2/m] = [m] -
- 0.595 = constant for time-lagged peak flow (assumes symmetrical unit-hydrograph) -
- t = length of rainfall event [seconds] -
- S = topographic slope [degrees] -
- m = transport coefficient for upslope area [unitless] -
- n = transport coefficient for slope [unitless] -
- Tc = Transport Capacity [kg/meters.second]
- K*C*P ~ +Kt = mitigating effects of soil type, vegetation cover, and land-use +practices. [unitless]
- gw = Hydrostatic pressure of water 9810 +[kg/m2.second]
- N = Manning's coefficient (0.3-0.6 for different +types of stream channels) [unitless]
- i = rainfall intensity +[m/rainfall event]
- A = upslope accumulated area per contour (cell) +width [m2/m] = [m]
- 0.595 = constant for time-lagged peak flow +(assumes symmetrical unit-hydrograph)
- t = length of rainfall event +[seconds]
- S = topographic slope [degrees]
- m = transport +coefficient for upslope area [unitless]
- n = transport coefficient +for slope [unitless]
- This is likely the best of the three equations for simulating erosion at the scale of small watersheds, including overland flow on hillslopes and channelized flow in gullies and streams. -
- It is likely not appropriate for simulating erosion and deposition processes in larger rivers, especially meandering flood plains. -
K*C*Pshould equal an appropriate value ofKt: 0.001 for a soft substrate, 0.0001 for a normal substrate, 0.00001 for a hard substrate, 0.000001 for a very hard substrate. See note below about methods for scaling these values.
-Nshould likely scale with channel vegetation so that 0.03 = clean/straight stream channel, 0.035 = major free-flowing river, 0.04 = sluggish stream with pools, 0.06 = very clogged streams. See below for methods to scale these values.
-- This is likely the best of the three equations for simulating erosion +at the scale of small watersheds, including overland flow on hillslopes +and channelized flow in gullies and streams.
- It is likely not +appropriate for simulating erosion and deposition processes in larger rivers, +especially meandering flood plains.
K*C*Pshould equal +an appropriate value ofKt: 0.001 for a soft substrate, 0.0001 +for a normal substrate, 0.00001 for a hard substrate, 0.000001 for a very +hard substrate. See note below about methods for scaling these values.
+Nshould likely scale with channel vegetation so that 0.03 += clean/straight stream channel, 0.035 = major free-flowing river, 0.04 = +sluggish stream with pools, 0.06 = very clogged streams. See below for +methods to scale these values.- Tc = Transport Capacity [kg/meters.second] -
- K*C*P ~ Kt = mitigating effects of soil type, vegetation cover, and land-use practices. [unitless] -
- gw = Hydrostatic pressure of water 9810 [kg/m2.second] -
- N = Manning's coefficient ~0.3-0.6 for different types of stream channels [unitless] -
- i = rainfall intensity [m/rainfall event] -
- A = upslope accumulated area per contour (cell) width [m2/m] = [m] -
- 0.595 = constant for time-lagged peak flow (assumes symmetrical unit-hydrograph) -
- t = length of rainfall event [seconds] -
- B = topographic slope [degrees] -
- n = transport coefficient (here assumed to be scaled to slope) [unitless] -
- Tc = Transport Capacity [kg/meters.second]
- K*C*P ~ +Kt = mitigating effects of soil type, vegetation cover, and land-use +practices. [unitless]
- gw = Hydrostatic pressure of water 9810 +[kg/m2.second]
- N = Manning's coefficient ~0.3-0.6 for different +types of stream channels [unitless]
- i = rainfall intensity [m/rainfall +event]
- A = upslope accumulated area per contour (cell) width [m2/m] = +[m]
- 0.595 = constant for time-lagged peak flow (assumes symmetrical +unit-hydrograph)
- t = length of rainfall event [seconds]
- B = +topographic slope [degrees]
- n = transport coefficient (here assumed +to be scaled to slope) [unitless]
- This implementation of the Shear Stress equation assumes the critical shear stress is 0. -
- This means the equation is likely to over predict erosion in situations where shear stress is less than the actual critical shear stress, such as on vegetated hillslopes. -
K*C*Pshould equal an appropriate value ofKt: 0.001 for a soft substrate, 0.0001 for a normal substrate, 0.00001 for a hard substrate, 0.000001 for a very hard substrate. See note below about methods for scaling these values.
-Nshould likely scale with channel vegetation so that 0.03 = clean/straight stream channel, 0.035 = major free-flowing river, 0.04 = sluggish stream with pools, 0.06 = very clogged streams. See below for methods to scale these values.
-- This implementation of the Shear Stress equation assumes the +critical shear stress is 0.
- This means the equation is likely +to over predict erosion in situations where shear stress is less than +the actual critical shear stress, such as on vegetated hillslopes. +
K*C*Pshould equal an appropriate value ofKt: +0.001 for a soft substrate, 0.0001 for a normal substrate, 0.00001 for a +hard substrate, 0.000001 for a very hard substrate. See note below about +methods for scaling these values.Nshould likely +scale with channel vegetation so that 0.03 = clean/straight stream channel, +0.035 = major free-flowing river, 0.04 = sluggish stream with pools, 0.06 = +very clogged streams. See below for methods to scale these values.- Tc = Transport Capacity [kg/meters.second] -
- R = Rainfall intensity factor [MJ.mm/ha.h.yr] -
- K*C*P ~ Kt = mitigating effects of soil type, vegetation cover, and land-use practices. [unitless] -
- A = upslope accumulated area per contour (cell) width [m2/m] = [m] -
- S = topographic slope [degrees] -
- m = transport coefficient for upslope area [unitless] -
- n = transport coefficient for slope [unitless] -
- Tc = Transport Capacity [kg/meters.second]
- R = Rainfall +intensity factor [MJ.mm/ha.h.yr]
- K*C*P ~ Kt = mitigating effects of +soil type, vegetation cover, and land-use practices. [unitless]
- A = +upslope accumulated area per contour (cell) width [m2/m] = [m]
- S = +topographic slope [degrees]
- m = transport coefficient for upslope area +[unitless]
- n = transport coefficient for slope [unitless]
- The USPED equation is best suited for modeling erosion and deposition on hillslopes and small gullies. -
- It will vastly over predict erosion/deposition in channels and streams. -
- The USPED equation is best suited for modeling erosion and deposition +on hillslopes and small gullies.
- It will vastly over predict +erosion/deposition in channels and streams.
DESCRIPTION
-r.landscape.evol takes as input a raster digital elevation model +(DEM) of surface topography and an input raster DEM of bedrock elevations, +as well as several environmental variables, and computes the net change in +elevation due to erosion and deposition Stream Power equation, the Shear +Stress equation, or the USPED equation.
-NOTES
-Transport capacity equations.
-Users may select to use the Stream Power equation, the Shear Stress
+equation, or the USPED equations with variable transp_eq. All
+three equations estimate transport capacity as [kg/m.s], and
+thus eventually erosion/deposition rate as [kg/m2.s], which
+is transformed to [vertical meters/cell] using the variable
+sdensity (see below for details of these conversions).
Transport capacity equations.
- -Users may select to use the Stream Power equation, the Shear Stress equation, or the USPED equations with variable transp_eq. All three equations estimate transport capacity as [kg/m.s], and thus eventually erosion/deposition rate as [kg/m2.s], which is transformed to [vertical meters/cell] using the variable sdensity (see below for details of these conversions).
It is important to note that in this new version of r.landscape.evol, only one transport equation will be used to model sediment flux across the entire landscape. Chane in process will be simulated through scalar m and n exponents (see below for details).
It is important to note that in this new version of
+r.landscape.evol, only one transport equation will be
+used to model sediment flux across the entire landscape. Chane in process
+will be simulated through scalar m and n exponents
+(see below for details).
1) Stream power equation:
Tc=Kt*gw*1/N*h^m*B^n
- where: h = depth of flow = (i*A)/(0.595*t)
- and: B = slope (rise over run)
+ where: h = depth of flow = (i*A)/(0.595*t) and: B = slope (rise over run)
a) GIS Implementation:
-Tc=K*C*P*gw*(1/N)*((i*A)/(0.595*t))^m*(tan(S)^n)
-Tc=K*C*P*gw*(1/N)*((i*A)/(0.595*t))^m*(tan(S)^n) b) Variables:
--
-
c) Converted to Map Algebra:
- ${K}*${C}*${P} * exp(${manningn}, -1) * 9810. * exp((((${rain}/1000.)*\
+ ${K}*${C}*${P} * exp(${manningn}, -1) * 9810. *
+exp((((${rain}/1000.)*\
${flowacc})/(0.595*${stormtimet})), graph(${flowacc}, ${exp_m1a},\
${exp_m1b}, ${exp_m2a},${exp_m2b}) ) * exp(tan(${slope}), graph(${slope},\
${exp_n1a},${exp_n1b}, ${exp_n2a},${exp_n2b}))
@@ -200,57 +60,65 @@ 1) Stream power equation:
d) NOTES:
-
-
+
2) Shear stress equation:
Tc=Kt*tau^m
- where: tau = shear stress = gw*h*B
- and: B = slope (rise over run)
- and: h = depth of flow = (i*A)/(0.595*t)
+ where: tau = shear stress = gw*h*B and: B = slope (rise over run) and:
+ h = depth of flow = (i*A)/(0.595*t)
a) GIS Implmentation:
-Tc=K*C*P*(gw*((i*A)/(0.595*t)*(tan(S))))^n
-
+Tc=K*C*P*(gw*((i*A)/(0.595*t)*(tan(S))))^n
b) Variables:
-
-
+
c) Converted to Map Algebra:
- ${K}*${C}*${P} * exp(9810.*(((${rain}/1000)*${flowacc})/(0.595*\
+ ${K}*${C}*${P} *
+exp(9810.*(((${rain}/1000)*${flowacc})/(0.595*\
${stormtimet}))*tan(${slope}), graph(${flowacc}, ${exp_n1a},${exp_n1b},\
${exp_n2a},${exp_n2b}))
d) NOTES:
-
-
+
3) USPED equation:
@@ -261,214 +129,580 @@ 3) USPED equation:
a) GIS Implementation:
-Tc=R*K*C*P*A^m*tan(S)^n
-
+Tc=R*K*C*P*A^m*tan(S)^n
b) Variables:
-
-
+
c) Converted to Map Algebra:
- (${R}*${K}*${C}*${P}*exp((${flowacc}*${res}),graph(${flowacc}, ${exp_m1a},\
- ${exp_m1b}, ${exp_m2a},${exp_m2b}))*exp(sin(${slope}), graph(${slope}, \
- ${exp_n1a},${exp_n1b}, ${exp_n2a},${exp_n2b})))
+ (${R}*${K}*${C}*${P}*exp((${flowacc}*${res}),graph(${flowacc},
+${exp_m1a},\
+ ${exp_m1b}, ${exp_m2a},${exp_m2b}))*exp(sin(${slope}), graph(${slope},
+ \ ${exp_n1a},${exp_n1b}, ${exp_n2a},${exp_n2b})))
d) NOTES:
-
-
-
-Scalar m and n exponents to simulate changing process across landscapes.
-
-Exponents m and n are used to influence the behavior of the transport equations by differentially weighting the influence of upslope accumulated area (and thus depth of flow) (m) or the influence of local slope (n). Depending on how these are each weighted, transport estimates can be made for overland flow processes, rilling and gullying, or channelized flow (see references below, but in particular Peckham 2003, Mathier et al 1989, and Kwang and Parker 2017). Following a suggestion in Peckham 2003, this new version of r.landscape.evol simulates change in process across the landscape by scaling m and n to changes in topography and flow accumulation. As this is largely an experimental process, the specifics of this scaling are exposed to the user via the m and n variables. The user can define the scalar relationship of m to surface flow accumulation, and n to local slope. Sensible default values are included to help the user know where to start.
-
-Exponent m relates to the influence of upslope area (and thus flow depth, discharge) on transport capacity in the Stream Power and USPED, but is not used in the Shear Stress equation. Values of m are generally thought to be between 2 and 1, and experimentation suggests that they should scale inversely with increasing depth of flow. Exponent m will scale linearly with the value of flow accumulation between the two cutoff thresholds specified: "thresh1,m1,thresh2,m2". So, for example, if you would like the value of exponent m to scale from 1.2 to 1 between a flow accumulation value of 5 and 50, enter the following into the variable m: "5,1.2,50,1". The exponent m will remain 1.2 for all cells where flow accumulation is below 5, and will remain 1 for all cells with flow accumulation above 50. It will scale linearly between 1.2 and 1 for all cells with values of flow accumulation between 5 and 50.
-
-A literature search indicates that maximum values of m should be less than or equal to 2, and that scaling between 1.2 and 1 is probably a good range to start with.
-
-Exponent n relates to the influence of local topographic slope on transport capacity, and is used in the Stream Power, Shear Stress, and USPED equations. Values of n are generally thought to be between 2 and 1, and experimentation suggests that they should scale inversely with increasing local slope. Exponent n will scale linearly with slope between the two slope cutoff thresholds specified: "thresh1,n1,thresh2,n2". So, for example, if you would like the value of exponent n to scale from 1.3 to 1 between a slope value of 10 and 30, enter the following into the variable n: "10,1.3,30,1". The exponent n will remain 1.3 for all cells where slope is below 5, and will remain 1 for all cells with slope above 30. It will scale linearly between 1.3 and 1 for all cells with values of slope between 5 and 30.
-
-A literature search indicates that maximum values values of n should be less than or equal to 2, and that scaling between 1.3 and 1 is probably a good range to start with.
-
-Scaling other input values.
-
-To ensure proper behavior for landscape evolution simulation over long periods, it is important that most of the important variables be allowed to vary spatially as they would on a real landscape. The three most important sets of variables are a) Soil, vegetation cover, and land use factors k, c, p, which together approximate erodibility factor Kt, b) Manning's N manningn which is used to estimate stream power/shear stress of flowing water in different types of channels and surface conditions, and c) flowcontrib, the rainfall excess rate (percentage of direct precipitation that will flow off of a cell), which is used to estimate the flow depth (see below).
-
-Because upslope accumulated area A is a major influencing factor in each of the three equations, transport capacity (and thus erosion/deposition rate) will be inordinately governed by A as values of flow accumulation approach very large numbers (e.g., >> 10,000). This will be partially mitigated with scalar m and n (see above), but will need additional dampening by scaling Kt, N, and rainfall excess.
-
-Kt is composed of the K, C, and P factors. If empirical patterns of K, C, and P are known (e.g., digitized or classified from remotely senses data products), these should be entered as maps in input variables k, c, and p.
-
-If empirically determined maps of these variables are not available, it is possible to use constants in their place, but it will be much better to create maps using some theoretical concepts. The simplest way is to scale C to a wetness index using the principle that the more water accumulation, the denser the vegetation. From a DEM, it is possible to calculate the TCI topographic wetness index using r.watershed with output parameter tci. Here is an example set of r.recode rules to create a C map from TCI to enter in input variable c:
-
-0:3:0.1:0.01
-3:7:0.01:0.005
-7:10:0.005:0.004
-10:*:0.004
+
+
+Scalar m and n exponents to simulate changing process across
+landscapes.
+
+Exponents m and n are used to influence the
+behavior of the transport equations by differentially weighting the influence
+of upslope accumulated area (and thus depth of flow) (m) or the
+influence of local slope (n). Depending on how these are each
+weighted, transport estimates can be made for overland flow processes, rilling
+and gullying, or channelized flow (see references below, but in particular
+Peckham 2003, Mathier et al 1989, and Kwang and Parker 2017). Following a
+suggestion in Peckham 2003, this new version of r.landscape.evol
+simulates change in process across the landscape by scaling m
+and n to changes in topography and flow accumulation. As
+this is largely an experimental process, the specifics of this scaling
+are exposed to the user via the m and n
+variables. The user can define the scalar relationship of m
+to surface flow accumulation, and n to local slope. Sensible
+default values are included to help the user know where to start.
+
+Exponent m relates to the influence of upslope area
+(and thus flow depth, discharge) on transport capacity in the Stream
+Power and USPED, but is not used in the Shear Stress equation. Values
+of m are generally thought to be between 2 and 1, and
+experimentation suggests that they should scale inversely with
+increasing depth of flow. Exponent m will scale linearly with
+the value of flow accumulation between the two cutoff thresholds specified:
+"thresh1,m1,thresh2,m2". So, for example, if you
+would like the value of exponent m to scale from 1.2 to 1
+between a flow accumulation value of 5 and 50, enter the following into
+the variable m: "5,1.2,50,1". The
+exponent m will remain 1.2 for all cells where flow accumulation
+is below 5, and will remain 1 for all cells with flow accumulation above
+50. It will scale linearly between 1.2 and 1 for all cells with values of
+flow accumulation between 5 and 50.
+
+A literature search indicates that maximum values of m
+should be less than or equal to 2, and that scaling between 1.2 and 1 is
+probably a good range to start with.
+
+Exponent n relates to the influence of local topographic
+slope on transport capacity, and is used in the Stream Power, Shear Stress,
+and USPED equations. Values of n are generally thought to
+be between 2 and 1, and experimentation suggests that they should scale
+inversely with increasing local slope. Exponent n
+will scale linearly with slope between the two slope cutoff thresholds
+specified: "thresh1,n1,thresh2,n2". So, for example,
+if you would like the value of exponent n to scale from 1.3 to
+1 between a slope value of 10 and 30, enter the following into the variable
+n: "10,1.3,30,1". The exponent
+n will remain 1.3 for all cells where slope is below 5, and
+will remain 1 for all cells with slope above 30. It will scale linearly
+between 1.3 and 1 for all cells with values of slope between 5 and 30.
+
+A literature search indicates that maximum values values of n
+should be less than or equal to 2, and that scaling between 1.3 and 1 is
+probably a good range to start with.
+
+Scaling other input values.
+
+To ensure proper behavior for landscape evolution simulation over long
+periods, it is important that most of the important variables be allowed
+to vary spatially as they would on a real landscape. The three most
+important sets of variables are a) Soil, vegetation cover, and land use
+factors k, c, p, which
+together approximate erodibility factor Kt, b) Manning's N
+manningn which is used to estimate stream power/shear stress
+of flowing water in different types of channels and surface conditions,
+and c) flowcontrib, the rainfall excess rate (percentage
+of direct precipitation that will flow off of a cell), which is used to
+estimate the flow depth (see below).
+
+Because upslope accumulated area A is a major influencing
+factor in each of the three equations, transport capacity (and thus
+erosion/deposition rate) will be inordinately governed by A
+as values of flow accumulation approach very large numbers (e.g., >>
+10,000). This will be partially mitigated with scalar m and
+n (see above), but will need additional dampening by scaling
+Kt, N, and rainfall excess.
+
+Kt is composed of the K, C,
+and P factors. If empirical patterns of K,
+C, and P are known (e.g., digitized or classified
+from remotely senses data products), these should be entered as maps in input
+variables k, c, and p.
+
+If empirically determined maps of these variables are not available,
+it is possible to use constants in their place, but it will be much better
+to create maps using some theoretical concepts. The simplest way is to
+scale C to a wetness index using the principle that the more
+water accumulation, the denser the vegetation. From a DEM, it is possible
+to calculate the TCI topographic wetness index using r.watershed
+with output parameter tci. Here is an example set of
+r.recode rules to create a C map from TCI to enter in
+input variable c:
+
+0:3:0.1:0.01 3:7:0.01:0.005 7:10:0.005:0.004 10:*:0.004
-Here, low values of TCI will be coded as shrubs or open woodlands. Moderate values of TCI will become wooded, and high values of TCI will coded as dense riparian vegetation. It's important to note that this should be done with a TCI map created with r.watershed on the same DEM that will be used as the initial DEM for the simulation.
+Here, low values of TCI will be coded as shrubs or open woodlands. Moderate
+values of TCI will become wooded, and high values of TCI will coded as dense
+riparian vegetation. It's important to note that this should be done
+with a TCI map created with r.watershed on the same DEM that will
+be used as the initial DEM for the simulation.
-From here, it is possible to map rainfall excess to values of C. The following recode rules will achieve a reasonable mapping:
+From here, it is possible to map rainfall excess to values of
+C. The following recode rules will achieve a reasonable
+mapping:
-0.1:0.05:85:80
-0.05:0.01:80:60
-0.01:0.005:60:45
-0.005:0.001:45:35
+0.1:0.05:85:80 0.05:0.01:80:60 0.01:0.005:60:45 0.005:0.001:45:35
-Here, as vegetation becomes more protective of detachment, it is also scaled to become more conducive to water infiltration, and thus more prohibitive to excess water escaping from the cell. The resulting map should be entered into input variable flowcontrib.
+Here, as vegetation becomes more protective of detachment, it is also scaled
+to become more conducive to water infiltration, and thus more prohibitive
+to excess water escaping from the cell. The resulting map should be entered
+into input variable flowcontrib.
-Finally, Manning's N can be scaled to flow accumulation (i.e., computed with r.watershed) using the following recode rules to create an input map for variable manningn:
+Finally, Manning's N can be scaled to flow accumulation
+(i.e., computed with r.watershed) using the following recode rules
+to create an input map for variable manningn:
-0:10:0.03:0.04
-10:100:0.04:0.05
-100:10000:0.05:0.06
-10000:*:0.06
+0:10:0.03:0.04 10:100:0.04:0.05 100:10000:0.05:0.06 10000:*:0.06
-Here, the assumption is that as flow accumulation increases, the channel will become more complex. These particular rules assume that the scale of analysis is at the level of small watershed feeding into a small trunk stream, not a large free-flowing river. If some empirical data about channel conditions are known, then the values used in the recode statement should be adjusted to reflect this. Again, it's important to note that this should be done with a flow accumulation map created with r.watershed on the same DEM that will be used as the initial DEM for the simulation. Further, the -a flag in r.watershed should be checked so that the output flow accumulation will contain only positive numbers
+Here, the assumption is that as flow accumulation increases, the channel
+will become more complex. These particular rules assume that the scale
+of analysis is at the level of small watershed feeding into a small trunk
+stream, not a large free-flowing river. If some empirical data about channel
+conditions are known, then the values used in the recode statement should be
+adjusted to reflect this. Again, it's important to note that this should
+be done with a flow accumulation map created with r.watershed on the
+same DEM that will be used as the initial DEM for the simulation. Further,
+the -a flag in r.watershed should be checked so that the output
+flow accumulation will contain only positive numbers
Creating a hydrologically appropriate base DEM.
-It is vitally important the the input starting DEM be hydrologically valid and at an appropriate raster resolution. Resolution should be scaled to the size of the region being modeled, with the caveat that the assumptions of the way the transport equations are implemented will start to break down at larger cell resolutions. As a general rule of thumb, cell resolution should be <= 10m. This can be achieved through resampling/interpolation from coarser data sets (e.g., a 30m SRTM DEM). If interpolation is used, it is best to use an interpolation procedure that will result in relatively smooth interpolated DEM with minimal depressions. Generally, v.surf.bspline achieves good results when the spline step is double to triple the cell resolution of the coarser input map, and the smoothing parameter is set to provide some additional smoothing (e.g., ~0.1). This results in an interpolated DEM with a smooth surface and minimal localized depressions caused by over-fitting to localized surface trends. Although v.surf.rst can also be used, it often produces rectilinear artifacts from it's segmentation procedure that can adversely affect simulation of water flow on the interpolated DEM.
-
-The DEM should be clipped to a contiguous watershed boundary (e.g., extracted with r.watershed or r.water.outlet). Rectilinear input maps will produce erroneous results outside of internally contiguous watersheds leading to faulty statistics, so it is more useful to clip to the watershed of interest (e.g., using r.mapcalc). Finally, in order to assure that water will flow naturally across the DEM, it is important to ensure that the DEM is depressionless. This could be achieved with r.fill.dir to fill any interior basins to an elevation level with their spill point, but doing so creates many flat areas where otherwise channelized flow will diverge (and thus deposit). This can be partially addressed by adjusting convergence to a low value, which forces the flow accumulation routine in r.watershed to send a higher proportion of the flow to the most downstream cell. Perhaps a better approach is to create a depressionless DEM by carving the main streams through any blockages using r.carve. To do so, extract streams using r.watershed or r.streams.extract and an appropriate interior basin threshold parameter to isolate main trunk streams from smaller branches. This streams vector is used as input into r.carve, along with a small value for additional carving (e.g. > 0.5m) to provide a “pre-carved” stream channel for the flow to converge in. Used in this way r.carve will cut through any locally high blockages between depressions in the channel, and will optionally carve an additional depth at all cells crossed by the input stream vector. This can be especially useful if flow will tend to diverge in low lying areas. Further, you can specify the width of the stream, which will carve adjacent cells to create a wider stream path. Optionally, the addon module r.stream.order can be used to separate streams by their network order (e.g., Strahler order), which can be use to iteratively carve narrower and shallower portions of the stream network, starting with the highest order streams, and proceeding to the lowest order, creating interim carved DEM's along the way.
-
-Estimating soil depth.
-
-Soil depth is important in the routine, as it provides a depth-based limitation on the amount of erosion that can occur at any particular cell (see below). The depth of soil available to erode is the difference between the current surface elevations (DEM) and the bedrock elevation map initbdrk. The simplest way to estimate the bedrock elevation map is to subtract a constant from the starting DEM map used for elev using r.mapcalc. A more complex bedrock topography can be estimated using the addon module r.soildepth. In either case, it is important to use the same DEM to derive the bedrock elevations as you will use for the initial starting topography in the simulation.
-
-Climate data file.
-
-Users can use constants for climate data, or can use an input climate file with columns of comma separated values arranged in order of: "R,rain,storms,stormlength,stormsi" A new line should be used for each year of the simulation. The file can have a one-line header or no header. Do not included a column containing dates, but ensure that the number of rows matches the value you input for number.
-
-Note that only the USPED equation needs a value for R factor, and USPED does not need the remaining climate variables. In the case of using USPED, only the first column needs to contain data (for R factor), but you still need to include all columns (the remaining columns can be with zeros or NaN's).
-
-In the case of using the Stream Power or Shear Stress equations, you still must create a CSV file with 5 columns, but the first column (for R Factor) can be filled with zeros or NaN's.
-
-When using a climate file, you enter the path to the text file as variable climfile. This will override values or maps entered into variables r, rain, storms, stormlength, or stormsi. A fatal error message will be raised if the number of rows in the input climate file does not match the value entered for the variable number.
+It is vitally important the the input starting DEM be hydrologically valid
+and at an appropriate raster resolution. Resolution should be scaled to the
+size of the region being modeled, with the caveat that the assumptions of the
+way the transport equations are implemented will start to break down at larger
+cell resolutions. As a general rule of thumb, cell resolution should be <=
+10m. This can be achieved through resampling/interpolation from coarser data
+sets (e.g., a 30m SRTM DEM). If interpolation is used, it is best to use an
+interpolation procedure that will result in relatively smooth interpolated
+DEM with minimal depressions. Generally, v.surf.bspline achieves
+good results when the spline step is double to triple the cell resolution
+of the coarser input map, and the smoothing parameter is set to provide some
+additional smoothing (e.g., ~0.1). This results in an interpolated DEM with
+a smooth surface and minimal localized depressions caused by over-fitting
+to localized surface trends. Although v.surf.rst can also be used,
+it often produces rectilinear artifacts from it's segmentation procedure
+that can adversely affect simulation of water flow on the interpolated DEM.
+
+The DEM should be clipped to a contiguous watershed boundary (e.g.,
+extracted with r.watershed or r.water.outlet). Rectilinear
+input maps will produce erroneous results outside of internally contiguous
+watersheds leading to faulty statistics, so it is more useful to clip
+to the watershed of interest (e.g., using r.mapcalc). Finally,
+in order to assure that water will flow naturally across the DEM, it is
+important to ensure that the DEM is depressionless. This could
+be achieved with r.fill.dir to fill any interior basins to an
+elevation level with their spill point, but doing so creates many flat areas
+where otherwise channelized flow will diverge (and thus deposit). This can
+be partially addressed by adjusting convergence to a low
+value, which forces the flow accumulation routine in r.watershed
+to send a higher proportion of the flow to the most downstream cell. Perhaps
+a better approach is to create a depressionless DEM by carving
+the main streams through any blockages using r.carve. To do so,
+extract streams using r.watershed or r.streams.extract
+and an appropriate interior basin threshold parameter to isolate main trunk
+streams from smaller branches. This streams vector is used as input into
+r.carve, along with a small value for additional carving (e.g. >
+0.5m) to provide a “pre-carved” stream channel for the flow to
+converge in. Used in this way r.carve will cut through any locally
+high blockages between depressions in the channel, and will optionally carve an
+additional depth at all cells crossed by the input stream vector. This can be
+especially useful if flow will tend to diverge in low lying areas. Further, you
+can specify the width of the stream, which will carve adjacent cells to create
+a wider stream path. Optionally, the addon module r.stream.order
+can be used to separate streams by their network order (e.g., Strahler order),
+which can be use to iteratively carve narrower and shallower portions of the
+stream network, starting with the highest order streams, and proceeding to
+the lowest order, creating interim carved DEM's along the way.
+
+Estimating soil depth.
+
+Soil depth is important in the routine, as it provides a depth-based
+limitation on the amount of erosion that can occur at any particular
+cell (see below). The depth of soil available to erode is the difference
+between the current surface elevations (DEM) and the bedrock elevation
+map initbdrk. The simplest way to estimate the bedrock
+elevation map is to subtract a constant from the starting DEM map used
+for elev using r.mapcalc. A more complex bedrock
+topography can be estimated using the addon module r.soildepth. In
+either case, it is important to use the same DEM to derive the bedrock
+elevations as you will use for the initial starting topography in the
+simulation.
+
+Climate data file.
+
+Users can use constants for climate data, or can use an input
+climate file with columns of comma separated values arranged in order of:
+"R,rain,storms,stormlength,stormsi" A new line should
+be used for each year of the simulation. The file can have a one-line header
+or no header. Do not included a column containing dates, but ensure that the
+number of rows matches the value you input for number.
+
+Note that only the USPED equation needs a value for R factor, and USPED
+does not need the remaining climate variables. In the case of using USPED,
+only the first column needs to contain data (for R factor), but you still
+need to include all columns (the remaining columns can be with zeros or
+NaN's).
+
+In the case of using the Stream Power or Shear Stress equations, you still
+must create a CSV file with 5 columns, but the first column (for R Factor)
+can be filled with zeros or NaN's.
+
+When using a climate file, you enter the path to the text file as variable
+climfile. This will override values or maps entered into
+variables r, rain, storms,
+stormlength, or stormsi. A fatal error
+message will be raised if the number of rows in the input climate file does
+not match the value entered for the variable number.
Rainfall excess and flow accumulation.
-This module will take rainfall totals into account when calculating the value of flow accumulation. It does so using r.watershed and the value of flowcontrib to calculate flow accumulation scaled by the percentage of rain that will flow off the cell (i.e., rainfall - infiltration). See above for a method to scale flowcontrib to C factor.
-
-Temporal Interval
-
-The USPED equation relies on the value of R from the RUSLE equation to define the temporal interval for landscape evolution. Typically, R is estimated at a yearly temporal interval, so it is important to understand the time step of your R input data before simulation with the USPED equation.
-
-The Stream Power and Shear Stress equations, on the other hand, accept storm-level data. This can be aggregated at any time step (per-storm, daily, weekly, monthly, yearly, decadal, etc.). The time step does not need to be an even interval; this means you can model on a per-storm basis where the interval between storms is not the same. To do so, you would use the option to enter a climate file where each line would detail the timing and intensity of each storm. You would then run the simulation with variable number equal to the total number of storms in your study interval.
-
-Approximation of depth of flow for Stream Power and Shear Stress equations.
-
-Flow depth is an important component for estimating stream power or shear stress. Here, it is estimated using upslope accumulated area (as modified by rainfall excess), rain fall in a typical erosion causing event (e.g., greater than ~30mm), and the length of the typical erosion causing event. Depth at peak flow is then estimated by assuming a symmetrical unit-hydrograph where total flow is the area below the hydrograph curve, and the total length equal to duration of the storm. The constant 0.595 is used to estimate the depth at peak flow under a symmetrical hydrograph where the area under the graph equals A (upslope accumulated area), and the horizontal width of the base of the hydrograph is equal to the length of the storm in seconds (stormlength).
-
-One of the benefits of this approach is that it is not tied to any specific time scale; any amount of time equal to or greater than 1 second can be modeled. For example, hourly rainfall totals can be entered as rain in sequence, with stormtime set to 3600 seconds, storms set to 1, and stormi set to 1. Hourly data could be aggregated to the level of the individual storm with the total for each storm entered as rain, stormtime equal to the total number of seconds each storm lasted, storms set to 1, and stormi set to some proportion of the storm where flow was at or near peak depths (e.g. 0.05). Daily rainfall totals can be entered as rain in sequence, with stormtime set to 84600 seconds, storms set to 1, and stormi set to some proportion where flow is at peak (e.g., 0.05). Monthly totals can be broken up into proportions per rain day, entered as rain with stormtime set to 84600 seconds, storms set to the number of storms that occurred that month, and stormi set to some proportion where flow is at or near peak depth (e.g., 0.05). Weekly, yearly, decadal, etc., totals can be entered in the same manner.
-
-This approach is more flexible than using R factor to encapsulate rainfall intensivity, as with USPED, as often R factor can only be estimated from rainfall totals at the timescale of the year or decade.
-
-Conversion of output of divergence to calculated erosion and deposition in vertical meters of elevation change.
-
-In order to convert the changes in transport capacity into the amount of elevation gained or lost by deposition or erosion, first the divergence in transport capacity is calculated in the EW and NS directions. These are then added back together to calculate the divergence in transport capacity (flux) in the direction of flow across the cell. Once this is done, the units are in kg/m2.s of sediment gained or lost. This is converted to meters of elevation gained or lostby dividing by soil density [kg/m3]. For USPED, which is tied to the temporal interval of R factor, this typically provides [m/year] as the output units. For the shear stress and stream power equations, however, this first comes out in units of [m/s]. It is then necessary to multiply by the number of seconds at peak flow depth (stormi * stormtime) and then by the number of erosive storms (storms) per year to get [m/year] elevation change.
-
-Computing elevation changes from one year to next.
-
-To compute the new surface elevation after erosion and deposition have occurred, it is necessary to add this year's ED map to last year's DEM, checking first if the amount of erodible soil in a given cell is less than the amount of erosion calculated. The cell will be prevented from eroding past this amount. If there is some soil depth remaining in the cell, then if the amount of erosion is more than the amount of soil, the routine will remove all the remaining soil and stop. Otherwise it will remove the amount of calculated erosion. If there is deposition, then it will be added on top of current depth of sediment (even if no sediment is currently in the cell).
-
-Finally, this routine is sensitive to edge effects carried forward from calculation of slope or other neighborhood routines used earlier in the module. To prevent null cells at the edges of maps, (the edge cells have no upstream cell, so get turned null), the initial DEM is patched underneath. Thus, the perimeter cells will never change in elevation throughout the simulation. Users are therefore strongly suggested to use a watershed boundary for their input maps (e.g., extracted from r.watershed, and then clipped with the map calculator), as cells at the watershed boundary should not change in elevation much in real world scenarios over the time spans of landscape evolution intended to be modeled with this module (100's to 1000's of years).
-
-References
-
-Aiello, A., Adamo, M., Canora, F., 2015. Remote sensing and GIS to assess soil erosion with RUSLE3D and USPED at river basin scale in southern Italy. CATENA 131, 174–185. https://doi.org/10.1016/j.catena.2015.04.003
-
-Aksoy, H., Kavvas, M.L., 2005. A review of hillslope and watershed scale erosion and sediment transport models. CATENA 64, 247–271. https://doi.org/10.1016/j.catena.2005.08.008
-
-Ayala, G., French, C., 2005. Erosion modeling of past land-use practices in the Fiume di Sotto di Troina river valley, north-central Sicily. Geoarchaeology 20, 149–167.
-
-Benavidez, R., Jackson, B., Maxwell, D., Norton, K., 2018. A review of the (Revised) Universal Soil Loss Equation (R/USLE): with a view to increasing its global applicability and improving soil loss estimates. Hydrology and Earth System Sciences Discussions 1–34. https://doi.org/10.5194/hess-2018-68
-
-Bosco, C., de Rigo, D., Dewitte, O., Poesen, J., Panagos, P., 2015. Modelling soil erosion at European scale: towards harmonization and reproducibility. Natural Hazards and Earth System Science 15, 225–245. https://doi.org/10.5194/nhess-15-225-2015
-
-Davy, P., Crave, A., 2000. Upscaling local-scale transport processes in large-scale relief dynamics. Physics and Chemistry of the Earth, Part A: Solid Earth and Geodesy 25, 533–541. https://doi.org/10.1016/S1464-1895(00)00082-X
-
-Dietrich, W.E., Bellugi, D.G., Sklar, L.S., Stock, J.D., Heimsath, A.M., Roering, J.J., 2003. Geomorphic Transport Laws for Predicting Landscape form and Dynamics, in: Wilcock, P.R., Iverson, R.M. (Eds.), Prediction in Geomorphology, Geophysical Monograph. American Geophysical Union, pp. 103–132.
-
-Diodato, N., 2006. Predicting RUSLE (Revised Universal Soil Loss Equation) Monthly Erosivity Index from Readily Available Rainfall Data in Mediterranean Area. The Environmentalist 26, 63–70. https://doi.org/10.1007/s10669-006-5359-x
-
-Hammad, A.A., Lundekvam, H., Børresen, T., 2004. Adaptation of RUSLE in the Eastern Part of the Mediterranean Region. Environmental Management 34, 829–841.
-
-Hancock, G.R., 2004. Modelling soil erosion on the catchment and landscape scale using landscape evolution models – a probabilistic approach using digital elevation model error, in: Super Soil 2004:3rd Australian New Zealand Soils Conference. University of Sydney, Australia.
-
-Kelley, A.D., Malin, M.C., Nielson, G.M., 1988. Terrain simulation using a model of stream erosion. ACM SIGGRAPH Computer Graphics 22, 263–268.
-
-Koko, Š., 2011. Simulation of gully erosion using the SIMWE model and GIS. Landform Analysis 17, 81–86.
-
-Kwang, J.S., Parker, G., 2017. Landscape evolution models using the stream power incision model show unrealistic behavior when <i>m</i> ∕ <i>n</i> equals 0.5. Earth Surface Dynamics 5, 807–820. https://doi.org/10.5194/esurf-5-807-2017
-
-Martínez-Casasnovas, J.A., Sánchez-Bosch, I., 2000. Impact assessment of changes in land use/conservation practices on soil erosion in the Penedès-Anoia vineyard region (NE Spain). Soil and Tillage Research 57, 101–106.
-
-Mathier, L., Roy, A.G., Paré, J.P., 1989. The effect of slope gradient and length on the parameters of a sediment transport equation for sheetwash. CATENA 16, 545–558. https://doi.org/10.1016/0341-8162(89)90041-6
-
-Mitasova, H., Barton, C.M., Ullah, I.I., Hofierka, J., Harmon, R.S., 2013. GIS-based soil erosion modeling, in: Shroder, J., Bishop, M.P. (Eds.), Remote Sensing and GIScience in Geomorphology, Treatise in Geomorphology. Academic Press, San Diego, pp. 228–258.
-
-Mitasova, H., Brown, W.M., Johnston, D., 2002. Terrain Modeling and Soil Erosion Simulation Final Report. Geographic Modeling Systems Lab, University of Illinois at Urbana-Champaign.
-
-Mitasova, H., Hofierka, J., Zlocha, M., Iverson, L.R., 1996a. Modelling topographic potential for erosion and deposition using GIS. International journal of geographical information systems 10, 629–641. https://doi.org/10.1080/02693799608902101
-
-Mitasova, H., Mitas, L., Brown, W.M., 2001. Multiscale Simulation of Land Use Impact on Soil Erosion and Deposition Patterns, in: Stott, D.E., Mohtar, R.H., Steinhardt, G.C. (Eds.), Sustaining the Global Farm: 10th International Soil Conservation Organization Meeting Held May 24-29, 1999. Purdue University and the USDA-ARS National Soil Erosion Research Laboratory, pp. 1163–1169.
-
-Mitasova, H., Mitas, L., Brown, W.M., Johnston, D., 1996b. Multidimensional Soil Erosion/Deposition Modeling Part III: Process based erosion simulation. Geographic Modeling and Systems Laboratory, University of Illinois at Urban-Champaign.
-
-Mitasova, H., Mitas, L., Brown, W.M., Johnston, D.M., 1999. Terrain modeling and Soil Erosion Simulations for Fort Hood and Fort Polk test areas. Geographic Modeling and Systems Laboratory, University of Illinois at Urbana-Champaign.
-
-Onori, F., De Bonis, P., Grauso, S., 2006. Soil erosion prediction at the basin scale using the revised universal soil loss equation (RUSLE) in a catchment of Sicily (southern Italy). Environmental Geology 50, 1129–1140.
-
-Panagos, P., Ballabio, C., Borrelli, P., Meusburger, K., Klik, A., Rousseva, S., Tadić, M.P., Michaelides, S., Hrabalíková, M., Olsen, P., Aalto, J., Lakatos, M., Rymszewicz, A., Dumitrescu, A., Beguería, S., Alewell, C., 2015a. Rainfall erosivity in Europe. Science of The Total Environment 511, 801–814. https://doi.org/10.1016/j.scitotenv.2015.01.008
-
-Panagos, P., Borrelli, P., Meusburger, K., Alewell, C., Lugato, E., Montanarella, L., 2015b. Estimating the soil erosion cover-management factor at the European scale. Land Use Policy 48, 38–50. https://doi.org/10.1016/j.landusepol.2015.05.021
-
-Panagos, P., Borrelli, P., Meusburger, K., van der Zanden, E.H., Poesen, J., Alewell, C., 2015c. Modelling the effect of support practices (P-factor) on the reduction of soil erosion by water at European scale. Environmental Science & Policy 51, 23–34. https://doi.org/10.1016/j.envsci.2015.03.012
-
-Panagos, P., Meusburger, K., Ballabio, C., Borrelli, P., Alewell, C., 2014. Soil erodibility in Europe: A high-resolution dataset based on LUCAS. Science of The Total Environment 479–480, 189–200. https://doi.org/10.1016/j.scitotenv.2014.02.010
-
-Peckham, S.D., 2003. Fluvial landscape models and catchment-scale sediment transport. Global and Planetary Change 39, 31–51. https://doi.org/10.1016/S0921-8181(03)00014-6
-
-Peeters, I., Rommens, T., Verstraeten, G., Govers, G., Van Rompaey, A., Poesen, J., Van Oost, K., 2006. Reconstructing ancient topography through erosion modelling. Geomorphology 78, 250–264.
-Pistocchi, A., Cassani, G., Zani, O., n.d. Use of the USPED model for mapping soil erosion and managing best land conservation practices 7.
-
-Renard, K.G., Foster, G.R., Weesies, G.A., McCool, D.K., Yoder, D.C., 1997. Predicting soil erosion by water: a guide to conservation planning with the Revised Universal Soil Loss Equation (RUSLE), in: Agriculture Handbook. US Department of Agriculture, Washington, DC, pp. 1–251.
-
-Renard, K.G., Foster, G.R., Weesies, G.A., Porter, J.P., 1991. RUSLE: Revised Universal Soil Loss Equation. Journal of Soil and Water Conservation 46, 30–33.
-
-Renard, K.G., Freimund, J.R., 1994. Using monthly precipitation data to estimate the R-factor in the revised USLE. Journal of Hydrology 157, 287–306.
-
-Sklar, L.S., Riebe, C.S., Marshall, J.A., Genetti, J., Leclere, S., Lukens, C.L., Merces, V., 2017. The problem of predicting the size distribution of sediment supplied by hillslopes to rivers. Geomorphology 277, 31–49. https://doi.org/10.1016/j.geomorph.2016.05.005
-
-Terranova, O., Antronico, L., Coscarelli, R., Iaquinta, P., 2009. Soil erosion risk scenarios in the Mediterranean environment using RUSLE and GIS: An application model for Calabria (southern Italy). Geomorphology 112, 228–245. https://doi.org/10.1016/j.geomorph.2009.06.009
-
-Tucker, G.E., Whipple, K.X., 2002. Topographic outcomes predicted by stream erosion models: Sensitivity analysis and intermodel comparison. J. Geophys. Res 107, 1–1.
-
-Warren, S.D., Mitasova, H., Hohmann, M.G., Landsberger, S., Skander, F.Y., Ruzycki, T.S., Senseman, G.M., 2005. Validation of a 3-D enhancement of the Universal Soil Loss Equation for preediction of soil erosion and sediment deposition. Catena 64, 281–296.
-
-Whipple, K.X., Tucker, G.E., 2002. Implications of sediment-flux-dependent river incision models for landscape evolution. Journal of Geophysical Research 107.
-
-Whipple, K.X., Tucker, G.E., 1999. Dynamics of the stream-power river incision model; implications for height limits of mountain ranges, landscape response timescales, and research needs. Journal of Geophysical Research 104, 17,661-17,674.
-
-Willgoose, G., 2005. Mathematical Modeling of Whole Landscape Evolution. Annual Review of Earth and Planetary Sciences 33, 443–459.
-
-
-
This module will take rainfall totals into account when calculating the +value of flow accumulation. It does so using r.watershed and the +value of flowcontrib to calculate flow accumulation scaled +by the percentage of rain that will flow off the cell (i.e., rainfall - +infiltration). See above for a method to scale flowcontrib +to C factor.
+ +Temporal Interval
+ +The USPED equation relies on the value of R from the RUSLE equation +to define the temporal interval for landscape evolution. Typically, R is +estimated at a yearly temporal interval, so it is important to understand the +time step of your R input data before simulation with the USPED equation.
+ +The Stream Power and Shear Stress equations, on the other hand, accept +storm-level data. This can be aggregated at any time step (per-storm, daily, +weekly, monthly, yearly, decadal, etc.). The time step does not need to be an +even interval; this means you can model on a per-storm basis where the interval +between storms is not the same. To do so, you would use the option to enter +a climate file where each line would detail the timing and intensity of each +storm. You would then run the simulation with variable number +equal to the total number of storms in your study interval.
+ +Approximation of depth of flow for Stream Power and Shear Stress +equations.
+ +Flow depth is an important component for estimating stream power or shear +stress. Here, it is estimated using upslope accumulated area (as modified by +rainfall excess), rain fall in a typical erosion causing event (e.g., greater +than ~30mm), and the length of the typical erosion causing event. Depth at +peak flow is then estimated by assuming a symmetrical unit-hydrograph where +total flow is the area below the hydrograph curve, and the total length +equal to duration of the storm. The constant 0.595 is used to estimate +the depth at peak flow under a symmetrical hydrograph where the area under +the graph equals A (upslope accumulated area), and the horizontal width of +the base of the hydrograph is equal to the length of the storm in seconds +(stormlength).
+ +One of the benefits of this approach is that it is not tied to any +specific time scale; any amount of time equal to or greater than 1 second +can be modeled. For example, hourly rainfall totals can be entered as +rain in sequence, with stormtime set to +3600 seconds, storms set to 1, and stormi +set to 1. Hourly data could be aggregated to the level of the individual +storm with the total for each storm entered as rain, +stormtime equal to the total number of seconds each storm +lasted, storms set to 1, and stormi +set to some proportion of the storm where flow was at or near peak depths +(e.g. 0.05). Daily rainfall totals can be entered as rain +in sequence, with stormtime set to 84600 seconds, +storms set to 1, and stormi set to some +proportion where flow is at peak (e.g., 0.05). Monthly totals can be broken +up into proportions per rain day, entered as rain with +stormtime set to 84600 seconds, storms set to +the number of storms that occurred that month, and stormi set +to some proportion where flow is at or near peak depth (e.g., 0.05). Weekly, +yearly, decadal, etc., totals can be entered in the same manner.
+ +This approach is more flexible than using R factor to encapsulate rainfall +intensivity, as with USPED, as often R factor can only be estimated from +rainfall totals at the timescale of the year or decade.
+ +Conversion of output of divergence to calculated erosion and deposition +in vertical meters of elevation change.
+ +In order to convert the changes in transport capacity into the amount of +elevation gained or lost by deposition or erosion, first the divergence in +transport capacity is calculated in the EW and NS directions. These are then +added back together to calculate the divergence in transport capacity (flux) +in the direction of flow across the cell. Once this is done, the units are in +kg/m2.s of sediment gained or lost. This is converted to meters of elevation +gained or lostby dividing by soil density [kg/m3]. For USPED, which is tied +to the temporal interval of R factor, this typically provides [m/year] as +the output units. For the shear stress and stream power equations, however, +this first comes out in units of [m/s]. It is then necessary to multiply +by the number of seconds at peak flow depth (stormi +* stormtime) and then by the number of erosive storms +(storms) per year to get [m/year] elevation change.
+ +Computing elevation changes from one year to next.
+ +To compute the new surface elevation after erosion and deposition have +occurred, it is necessary to add this year's ED map to last year's DEM, +checking first if the amount of erodible soil in a given cell is less than +the amount of erosion calculated. The cell will be prevented from eroding +past this amount. If there is some soil depth remaining in the cell, then +if the amount of erosion is more than the amount of soil, the routine will +remove all the remaining soil and stop. Otherwise it will remove the amount +of calculated erosion. If there is deposition, then it will be added on top of +current depth of sediment (even if no sediment is currently in the cell).
+ +Finally, this routine is sensitive to edge effects carried forward +from calculation of slope or other neighborhood routines used earlier in +the module. To prevent null cells at the edges of maps, (the edge cells +have no upstream cell, so get turned null), the initial DEM is patched +underneath. Thus, the perimeter cells will never change in elevation throughout +the simulation. Users are therefore strongly suggested to use a watershed +boundary for their input maps (e.g., extracted from r.watershed, +and then clipped with the map calculator), as cells at the watershed +boundary should not change in elevation much in real world scenarios over +the time spans of landscape evolution intended to be modeled with this module +(100's to 1000's of years).
+ +KNOWN ISSUES
+ +This module is sensitive to the geometry of the input DEM. False flat +areas and very steep slope transitions that are in the path of the flowlines +will result in erroneous values, and perhaps even lead to instability in the +landscape evolution algorithms that will exhibit as large “spikes” +and “pits” in the output DEM's after several iterations, +and may lead to numerical instability and NULL values in the various output +maps. Preconditioning the input DEM to reduce these issues, which can be +introduced during initial interpolation, or by the process of filling basins +with r.fill.basins or carving streams with r.carve.
+ +The module is also sensitive to input climate parameters and the exponents +of flow and how they are scaled. It is important to test these out extensively +before use.
+ +At this time, this module should be considered to be at a robust alpha +stage. It appears stable enough, but needs to be tested more extensively +before it can be considered stable and ready for production use.
+ +SEE ALSO
+ +The MEDLAND project at Arizona +State University
+ +r.watershed, +r.terraflow, r.mapcalc
+ +Mitasova, H., C. M. Barton, I. I. Ullah, J. Hofierka, and R. S. Harmon +2013 GIS-based soil erosion modeling. In Remote Sensing and GIScience in +Geomorphology, edited by J. Shroder and M. P. Bishop. 3:228-258. San Diego: +Academic Press.
+ +REFERENCES
+ +Aiello, A., Adamo, M., Canora, F., 2015. Remote sensing and +GIS to assess soil erosion with RUSLE3D and USPED at river +basin scale in southern Italy. CATENA 131, 174–185. https://doi.org/10.1016/j.catena.2015.04.003
+ +Aksoy, H., Kavvas, M.L., 2005. A review of hillslope and watershed +scale erosion and sediment transport models. CATENA 64, 247–271. https://doi.org/10.1016/j.catena.2005.08.008
+ +Ayala, G., French, C., 2005. Erosion modeling of past land-use practices in +the Fiume di Sotto di Troina river valley, north-central Sicily. Geoarchaeology +20, 149–167.
+ +Benavidez, R., Jackson, B., Maxwell, D., Norton, K., 2018. A +review of the (Revised) Universal Soil Loss Equation (R/USLE): with +a view to increasing its global applicability and improving soil loss +estimates. Hydrology and Earth System Sciences Discussions 1–34. https://doi.org/10.5194/hess-2018-68
+ +Bosco, C., de Rigo, D., Dewitte, O., Poesen, J., Panagos, P., +2015. Modelling soil erosion at European scale: towards harmonization and +reproducibility. Natural Hazards and Earth System Science 15, 225–245. https://doi.org/10.5194/nhess-15-225-2015
+ +Davy, P., Crave, A., 2000. Upscaling local-scale transport +processes in large-scale relief dynamics. Physics and Chemistry +of the Earth, Part A: Solid Earth and Geodesy 25, 533–541. https://doi.org/10.1016/S1464-1895(00)00082-X
+ +Dietrich, W.E., Bellugi, D.G., Sklar, L.S., Stock, J.D., Heimsath, A.M., +Roering, J.J., 2003. Geomorphic Transport Laws for Predicting Landscape +form and Dynamics, in: Wilcock, P.R., Iverson, R.M. (Eds.), Prediction +in Geomorphology, Geophysical Monograph. American Geophysical Union, +pp. 103–132.
+ +Diodato, N., 2006. Predicting RUSLE (Revised Universal Soil Loss +Equation) Monthly Erosivity Index from Readily Available Rainfall +Data in Mediterranean Area. The Environmentalist 26, 63–70. https://doi.org/10.1007/s10669-006-5359-x
+ +Hammad, A.A., Lundekvam, H., Børresen, T., 2004. Adaptation of RUSLE in +the Eastern Part of the Mediterranean Region. Environmental Management 34, +829–841.
+ +Hancock, G.R., 2004. Modelling soil erosion on the catchment and landscape +scale using landscape evolution models – a probabilistic approach using +digital elevation model error, in: Super Soil 2004:3rd Australian New Zealand +Soils Conference. University of Sydney, Australia.
+ +Kelley, A.D., Malin, M.C., Nielson, G.M., 1988. Terrain simulation using +a model of stream erosion. ACM SIGGRAPH Computer Graphics 22, 263–268.
+ +Koko, Š., 2011. Simulation of gully erosion using the SIMWE model and +GIS. Landform Analysis 17, 81–86.
+ +Kwang, J.S., Parker, G., 2017. Landscape evolution models using the stream +power incision model show unrealistic behavior when <i>m</i> +∕ <i>n</i> equals 0.5. Earth Surface Dynamics 5, 807–820. https://doi.org/10.5194/esurf-5-807-2017
+ +Martínez-Casasnovas, J.A., Sánchez-Bosch, I., 2000. Impact assessment +of changes in land use/conservation practices on soil erosion in the +Penedès-Anoia vineyard region (NE Spain). Soil and Tillage Research 57, +101–106.
+ +Mathier, L., Roy, A.G., Paré, J.P., 1989. The effect of +slope gradient and length on the parameters of a sediment +transport equation for sheetwash. CATENA 16, 545–558. https://doi.org/10.1016/0341-8162(89)90041-6
+ +Mitasova, H., Barton, C.M., Ullah, I.I., Hofierka, J., Harmon, +R.S., 2013. GIS-based soil erosion modeling, in: Shroder, J., Bishop, +M.P. (Eds.), Remote Sensing and GIScience in Geomorphology, Treatise in +Geomorphology. Academic Press, San Diego, pp. 228–258.
+ +Mitasova, H., Brown, W.M., Johnston, D., 2002. Terrain Modeling and Soil +Erosion Simulation Final Report. Geographic Modeling Systems Lab, University +of Illinois at Urbana-Champaign.
+ +Mitasova, H., Hofierka, J., Zlocha, M., Iverson, L.R., 1996a. Modelling +topographic potential for erosion and deposition using GIS. International +journal of geographical information systems 10, 629–641. https://doi.org/10.1080/02693799608902101
+ +Mitasova, H., Mitas, L., Brown, W.M., 2001. Multiscale Simulation +of Land Use Impact on Soil Erosion and Deposition Patterns, in: Stott, +D.E., Mohtar, R.H., Steinhardt, G.C. (Eds.), Sustaining the Global Farm: +10th International Soil Conservation Organization Meeting Held May 24-29, +1999. Purdue University and the USDA-ARS National Soil Erosion Research +Laboratory, pp. 1163–1169.
+ +Mitasova, H., Mitas, L., Brown, W.M., Johnston, D., 1996b. Multidimensional +Soil Erosion/Deposition Modeling Part III: Process based erosion +simulation. Geographic Modeling and Systems Laboratory, University of Illinois +at Urban-Champaign.
+ +Mitasova, H., Mitas, L., Brown, W.M., Johnston, D.M., 1999. Terrain +modeling and Soil Erosion Simulations for Fort Hood and Fort Polk test +areas. Geographic Modeling and Systems Laboratory, University of Illinois +at Urbana-Champaign.
+ +Onori, F., De Bonis, P., Grauso, S., 2006. Soil erosion prediction at +the basin scale using the revised universal soil loss equation (RUSLE) +in a catchment of Sicily (southern Italy). Environmental Geology 50, +1129–1140.
+ +Panagos, P., Ballabio, C., Borrelli, P., Meusburger, K., Klik, +A., Rousseva, S., Tadić, M.P., Michaelides, S., Hrabalíková, +M., Olsen, P., Aalto, J., Lakatos, M., Rymszewicz, A., Dumitrescu, +A., Beguería, S., Alewell, C., 2015a. Rainfall erosivity in +Europe. Science of The Total Environment 511, 801–814. https://doi.org/10.1016/j.scitotenv.2015.01.008
+ +Panagos, P., Borrelli, P., Meusburger, K., Alewell, C., Lugato, E., +Montanarella, L., 2015b. Estimating the soil erosion cover-management +factor at the European scale. Land Use Policy 48, 38–50. https://doi.org/10.1016/j.landusepol.2015.05.021
+ +Panagos, P., Borrelli, P., Meusburger, K., van der Zanden, E.H., +Poesen, J., Alewell, C., 2015c. Modelling the effect of support +practices (P-factor) on the reduction of soil erosion by water at +European scale. Environmental Science & Policy 51, 23–34. https://doi.org/10.1016/j.envsci.2015.03.012
+ +Panagos, P., Meusburger, K., Ballabio, C., Borrelli, P., Alewell, +C., 2014. Soil erodibility in Europe: A high-resolution dataset based +on LUCAS. Science of The Total Environment 479–480, 189–200. https://doi.org/10.1016/j.scitotenv.2014.02.010
+ +Peckham, S.D., 2003. Fluvial landscape models and catchment-scale +sediment transport. Global and Planetary Change 39, 31–51. https://doi.org/10.1016/S0921-8181(03)00014-6
+ +Peeters, I., Rommens, T., Verstraeten, G., Govers, G., Van Rompaey, A., +Poesen, J., Van Oost, K., 2006. Reconstructing ancient topography through +erosion modelling. Geomorphology 78, 250–264. Pistocchi, A., Cassani, G., +Zani, O., n.d. Use of the USPED model for mapping soil erosion and managing +best land conservation practices 7.
+ +Renard, K.G., Foster, G.R., Weesies, G.A., McCool, D.K., Yoder, D.C., +1997. Predicting soil erosion by water: a guide to conservation planning +with the Revised Universal Soil Loss Equation (RUSLE), in: Agriculture +Handbook. US Department of Agriculture, Washington, DC, pp. 1–251.
+ +Renard, K.G., Foster, G.R., Weesies, G.A., Porter, J.P., 1991. RUSLE: +Revised Universal Soil Loss Equation. Journal of Soil and Water Conservation +46, 30–33.
+ +Renard, K.G., Freimund, J.R., 1994. Using monthly precipitation data +to estimate the R-factor in the revised USLE. Journal of Hydrology 157, +287–306.
+ +Sklar, L.S., Riebe, C.S., Marshall, J.A., Genetti, J., Leclere, S., Lukens, +C.L., Merces, V., 2017. The problem of predicting the size distribution of +sediment supplied by hillslopes to rivers. Geomorphology 277, 31–49. https://doi.org/10.1016/j.geomorph.2016.05.005
+ +Terranova, O., Antronico, L., Coscarelli, R., Iaquinta, +P., 2009. Soil erosion risk scenarios in the Mediterranean +environment using RUSLE and GIS: An application model for +Calabria (southern Italy). Geomorphology 112, 228–245. https://doi.org/10.1016/j.geomorph.2009.06.009
+ +Tucker, G.E., Whipple, K.X., 2002. Topographic outcomes predicted by stream +erosion models: Sensitivity analysis and intermodel comparison. J. Geophys. Res +107, 1–1.
+ +Warren, S.D., Mitasova, H., Hohmann, M.G., Landsberger, S., Skander, F.Y., +Ruzycki, T.S., Senseman, G.M., 2005. Validation of a 3-D enhancement of the +Universal Soil Loss Equation for preediction of soil erosion and sediment +deposition. Catena 64, 281–296.
+ +Whipple, K.X., Tucker, G.E., 2002. Implications of sediment-flux-dependent +river incision models for landscape evolution. Journal of Geophysical +Research 107.
+ +Whipple, K.X., Tucker, G.E., 1999. Dynamics of the stream-power river +incision model; implications for height limits of mountain ranges, landscape +response timescales, and research needs. Journal of Geophysical Research +104, 17,661-17,674.
+ +Willgoose, G., 2005. Mathematical Modeling of Whole Landscape +Evolution. Annual Review of Earth and Planetary Sciences 33, 443–459.
+ +AUTHORS
+ +Isaac I. Ullah, C. Michael Barton, and Helena Mitasova