# Soil erosion processes

The soil erosion model is based on the Modified MMF model (Morgan and Duzant, 2008), which determines total soil erosion from detachment by raindrop impact and detachment by runoff and considers cell deposition. The next paragraphs provide a detailed description of all these processes.

#### Estimation of rainfall energy <a href="#toc534724182" id="toc534724182"></a>

The total kinetic energy of the effective precipitation $$KE (J m^{-2})$$ is used to determine the detachment of soil particles by raindrop impact and is defined as:

_Equation 53_

$$
KE=KE_{LD}+KE_{DT}
$$

Where $$KE_{LD}$$ is the kinetic energy of the leaf drainage $$(Jm^{-2})$$ and $$KE_{DT}$$ is the kinetic energy of the direct throughfall $$(Jm^{-2})$$.

The kinetic energy of the leaf drainage $$KE_{LD}$$ is based on Brandt (1990):

_Equation 54_

$$
KE_{LD}=\begin{cases}
   0 &\text{for } PH < 0.15 \\
   LD(15.8\sqrt{PH}-5.87) &\text{for } PH \ge 0.15
\end{cases}
$$

Where _LD_ is the leaf drainage (mm) and _PH_ is the plant height (m), specified for each landuse class.

The kinetic energy of the direct throughfall is based on a relationship described by Marshall & Palmer (1948), which is representative of a wide range of environments (Morgan, 2005):

_Equation 55_

$$
KE_{DT}=DT(8.95+8.44 log_{10}I)
$$

Where _DT_ is the direct throughfall (mm) and _I_ is the intensity of the erosive precipitation $$(mm \space h^{-1})$$. The intensity of the erosive precipitation is a model parameter and varies according to geographical location. Morgan & Duzant (2008) proposes $$10 \space mm \space h^{-1}$$ for temperate climates, $$25 \space mm \space h^{-1}$$ for tropical climates and $$30 \space mm \space h^{-1}$$ for strongly seasonal climates (e.g. Mediterranean, tropical monsoon).

The leaf drainage _LD_, i.e. precipitation that reaches the soil surface as flow or drips from the leaves and stems of the vegetation, and direct throughfall _DT_ i.e. precipitation that reaches the soil surface directly through gaps in the vegetation cover, from Equation 54 and 55, are obtained from the effective precipitation $$Pe_t (mm)$$, also known as throughfall. The effective precipitation from the hydrological model is first corrected for the slope angle, following Choi et al. (2017):

_Equation 56_

$$
Pe_t=Pe_t\space cos \space S
$$

Where $$Pe_t$$ is the effective precipitation (mm) and _S_ is the slope (°).

Leaf drainage is determined as:

_Equation 57_

$$
LD=Pe_t\space CC
$$

Where _CC_ is the canopy cover (proportion between zero and unity). The canopy cover is either introduced by a landuse-class specific tabular value or determined by the vegetation module. When the vegetation module is used, the canopy cover is obtained from the _LAI_ (Equation 6), maximized by 1.

Direct throughfall becomes:

_Equation 58_

$$
DT=Pe_t-LD
$$

#### Detachment of soil particles <a href="#toc534724183" id="toc534724183"></a>

Detachment of soil particles is determined separately for raindrop impact and accumulated runoff and is subsequently summed. The detachment of soil particles by raindrop impact $$(F, kg \space m^{-2})$$ and the detachment of soil particles by runoff   $$(H, kg \space m^{-2})$$ are determined for each of the soil texture classes separately and subsequently summed. The detachment of soil particles by raindrop impact is calculated as:

_Equation 59_

$$
F_i=K_i \space \frac{\% i}{100}(1-GC)KE\cdot 10^{-3}
$$

With _K_ the detachability of the soil by raindrop impact $$(g J^{-1})$$, _i_ the textural class, with _c_ for clay, _z_ for silt and _s_ for sand, and _GC_ the ground cover (-). The detachability of the soil for each texture class is included as a model parameter, for which Quansah (1982) proposed $$K_c = 0.1$$, $$K_z = 0.5$$, and $$K_s = 0.3 \space g J^{-1}$$. The ground cover _GC_, expressed as a proportion between zero and unity, protects the soil from detachment and is determined by the proportion of vegetation and rocks covering the surface. The ground cover is set to 1 in case the surface is covered with snow, which is determined from the snow module.

The detachment of soil particles by runoff   $$(H, kg \space m^{-2})$$ is calculated as:

_Equation 60_

$$
H_i=DR_i \frac{\%i}{100}Q_{rout,t}^{1.5}(1-GC)\space sin^{0.3}S \cdot10^{-3}
$$

Where _DR_ the detachability of the soil by runoff$$(g \space mm^{-1})$$ and $$Q_{rout,t}$$ is the routed streamflow on day $$t$$ (mm). The detachability of the soil _DR_ for each texture class is included as a model parameter for which Quansah (1982) proposed $$DR_c =1.0, DR_z=1.6$$, and $$DR_s = 1.5 g \space mm^{-1}$$.

#### Immediate deposition of detached particles <a href="#toc534724184" id="toc534724184"></a>

A proportion of the detached soil is deposited in the cell of its origin as a function of the abundance of vegetation and the surface roughness. The percentage of the detached sediment that is deposited (_DEP_) is estimated from the relationship obtained by Tollner et al. (1976) and calculated separately for each texture class:

_Equation 61_

$$
DEP=44.1N_{fi}^{0.29}
$$

Where $$N_f$$ is the particle fall number (-), defined as:

_Equation 62_

$$
N_{fi}=\frac{l \upsilon_{si}}{\upsilon d}
$$

Where $$l$$ is the length of a grid cell (m),  $$\upsilon_s$$ the particle fall velocity $$(m s^{-1})$$, $$\upsilon$$ the flow velocity$$(m s^{-1})$$, and _d_ the depth of flow (m). Particle fall velocities $$\upsilon_s$$ are estimated from:

_Equation 63_

$$
\upsilon_s=\frac{\frac{1}{18}\delta^2(\rho_s+\rho)g}{\eta}
$$

Where $$\delta$$ is the diameter of the particle (m), $$\rho_s$$ the sediment density$$(=2650\space kg \space m^{-3})$$, $$\rho$$ the flow density (typically $$1100 \space kg \space m^{-3}$$ for runoff on hillslopes; Abrahams et al., 2001), _g_ gravitational acceleration (taken as  $$9.81 \space m \space s^{-2}$$) and $$\eta$$ the fluid viscosity (nominally $$0.001 \space kg \space m^{-1} \space s^{-1}$$, but taken as 0.0015 to allow for the effects of the sediment in the flow; Morgan & Duzant, 2008). When Equation 63 is applied to the three texture sizes of 2 μm for clay, 60 μm for silt and 200 μm for sand, this gives respective values of $$2 \cdot 10^{-6}\space m \space s^{-1}$$ for clay, $$2 \cdot 10^{-3}\space m \space s^{-1}$$ for silt and $$0.02 \space m \space s^{-1}$$for sand.

The flow velocity $$\upsilon$$ from Equation 62 is obtained by the Manning formula:

_Equation 64_

$$
\upsilon = \frac{1}{n'}d^{\frac{2}{3}}S^{\frac{1}{2}}
$$

Where _n'_ is the modified Manning's roughness coefficient ($$s \space m^{-1/3}$$), which is a combination of the Manning's roughness coefficient for the soil surface and vegetation, defined as (Petryk and Bosmajian, 1975):

_Equation 65_

$$
n'=\sqrt{n^2_{soil}+n^2_{vegetation}}
$$

The Manning's roughness coefficient for bare soil $$n_{soil}$$ is set to $$0.015 \space s \space m^{-1/3}$$, as suggested by Morgan and Duzant (2008) (Figure 4a). For tilled conditions (Figure 4b) the following equation is applied:

_Equation 66_

$$
n_{soil}=exp(-2.1132+0.0349RFR)
$$

Where _RFR_ is the surface roughness parameter $$(cm \space m^{-1})$$. Values for _RFR_ are tillage implementation specific and can be obtained from Morgan & Duzant (2008) Table IV.

The Manning's roughness coefficient for regular spaced vegetation (Figure 4c) $$n_{vegetation}$$ is obtained from the following equation (Jin et al., 2000):

_Equation 67_

$$
n_{vegetation}=\frac{d^{\frac{2}{3}}}{\sqrt{\frac{2g}{D NV}}}
$$

Where _D_ is the stem diameter (m) and _NV_ the stem density $$(stems \space m^{-2})$$. Stem diameter and stem density may be difficult to obtain for certain landuse classes with irregular spaced vegetation (e.g. forest, shrubland), therefore, users may opt to use tabular values for $$n_{vegetation}$$, e.g. from Chow (1959) (Figure 4d). Preliminary model runs showed that Equation 67 results in unrealistically high flow velocity values for landuse classes where the stem density is very low, such as in orchards. Therefore, in these conditions where the influence of vegetation on flow velocity is negligible, $$n_{vegetation}$$ can be set to 0.

Equations 62, 64 and 67 require a flow depth _d_, a model parameter that can be used in the model calibration. The value _d_ for should be taken such that it corresponds to a water depth from runoff generated within the cell margins, i.e. without accumulation of flow from upstream located cells.

![Figure 4. Surface and vegetation roughness options: (a) bare soil, (b) tilled soil (Equation 66), (c) regular vegetation (Eq. 67), and (d) irregular vegetation.](https://4210278315-files.gitbook.io/~/files/v0/b/gitbook-x-prod.appspot.com/o/spaces%2F68y0QNUHmAMyqeH7yW67%2Fuploads%2FWRqE9qPDxaVHIRyrseFc%2F5?alt=media)

#### Sediment taken into transport <a href="#toc534724185" id="toc534724185"></a>

The amount of sediment that is routed to downstream cells is determined from the sum of the detached sediment from raindrop impact (Equation 59) and accumulated runoff (Equation 60), subtracting the proportion of the sediment that is deposited within the cell of its origin (Equation 61):

_Equation 68_

$$
G_i=(F_i+H_i)\Big(1-\frac{DEP_i}{100}\Big)
$$

The amount of sediment that is routed to downstream cells is the summation of the individual amounts for clay, silt and sand.
