We assume that we can estimate sediment flow at sediment transport capacity as
where R~im, KCP~Kt and LS=Am sin b n (add units here), and m=1.6, n=1.3 for prevailing rill erosion while m=n=1 for prevailing sheet erosion. Then the net erosion/deposition is estimated as
where a [deg] is aspect of the terrain surface.
The model is simple enough to be run as
a sequence of commands provided here for
ArcView3.2,
ArcGIS8.1,
GRASS5.
For routine applications the model can
be run as a script, using Avenue
in ArcView, Visual basic in ArcGIS8, and tcl/tk or UNIX scripting in GRASS.
|