Package net.simplace.sim.components.experimental.erosion

Class ErosionRoseAPSIM

java.lang.Object

net.simplace.sim.model.FWSimComponent

net.simplace.sim.components.experimental.erosion.ErosionRoseAPSIM

All Implemented Interfaces:

net.simplace.sim.util.FWSimFieldContainer

public class ErosionRoseAPSIM
extends net.simplace.sim.model.FWSimComponent

This SimComponent calculates soil erosion using the Rose model which is based on a simplified sediment concentration function described by Rose (1985). The algorithm follows the same formulations used in the APSIM (Agricultural Production Systems sIMulator) platform, which was adapted from the PERFECT model (Littleboy et al., 1992). The general equation is given as:

\[ \begin{eqnarray} A & = & 2700 \cdot S \cdot (100 - COV) \cdot \lambda \cdot \frac{Q}{100} \end{eqnarray} \]

Where A (SoilLoss) is the event soil loss (t/ha); S (cSlopePercent) is the slope angle (%); COV is the soil cover (%) (iCoverFrac * 100); Q (iRunoff) is the event runoff (mm); and \( \lambda \) is a factor approximating efficiency of entrainment (dimensionless). The soil cover fraction (iCoverFrac) (0-1) can be calculated by the SimComponent RunoffSurfaceCover, whereas iRunoff (mm) can be obtained by one of the SimComponents which calculates runoff (e.g., RunoffCurveNumberAPSIM). The efficiency of entrainment \( \lambda \) is calculated as follows:

\[ \begin{eqnarray} \lambda & = & \lambda_{bare} \cdot e^{-\beta \cdot COV} \end{eqnarray} \]

Where \( \lambda_{bare} \) (cEntrainEff) is the efficiency of entrainment of bare surface (dimensionless); and \( \beta \) (cErosRoseB2) is a cover ‘sensitivity’ factor (dimensionless). In special cases where a distinction should be made for the sediment load transport which is maintained in suspension by turbulence in the flowing water, users can prescribe \( \lambda_{bare} \) and \( \beta \) values for such cases respectively (cEntrainEffBed, cEntrainEffSusp, cErosRoseB2Bed, cErosRoseB2Susp). In this case, both cEntrainEff and cErosRoseB2 values should be prescribed as -1. The input soil bulk density (cSoilBD) in g/cm3 of the top layer is used to convert A from Eq. 3 from t/ha to mm (mmSoilLoss):

\[ \begin{eqnarray} mmSoilLoss & = & \frac{[A \cdot (1e^6/1e^8)]}{cSoilBD[0]} \end{eqnarray} \]

References

Rose, C.W., 1985. Developments in soil erosion and deposition models. Adv. soil Sci. 2, 1–63. https://doi.org/10.1007/978-1-4612-5088-3_1

APSIM v7.1 source code. Last time accessed: Dec-14-2022. URL: https://github.com/APSIMInitiative/APSIMClassic/blob/83039e54200481899b4566761f47fd61d53b4fdf/Model/SoilErosion/SoilErosion.cs

Littleboy, M., D.M. Silburn, D.M. Freebairn, D.R. Woodruff, G.L. Hammer, and J.K. Leslie. 1992. Impact of soil erosion on production in cropping systems. I. Development and validation of a simulation model. Aust. J. Soil Res. 30, 757-774.

Holzworth, Dean P., Neil I. Huth, Peter G. deVoil, Eric J. Zurcher, Neville I. Herrmann, Greg McLean, Karine Chenu, et al. “APSIM – Evolution towards a New Generation of Agricultural Systems Simulation.” Environmental Modelling & Software 62 (December 2014): 327–350. https://doi.org/10.1016/j.envsoft.2014.07.009