GRASS logo


r.divergence - Computes divergence of a vector field defined by magnitude and direction


raster, divergence


r.divergence --help
r.divergence magnitude=name direction=name output=name [--overwrite] [--help] [--verbose] [--quiet] [--ui]


Allow output files to overwrite existing files
Print usage summary
Verbose module output
Quiet module output
Force launching GUI dialog


magnitude=name [required]
Name of input raster map representing magnitude
direction=name [required]
Name of input raster map representing direction
Direction is in degrees ccw from the east
output=name [required]
Name of output divergence raster

Table of contents


Module r.divergence computes the divergence of a vector field given by magnitude and direction raster maps. Direction is in degrees counterclockwise from the east and can be computed using r.slope.aspect. This module can be used for estimating erosion and deposition rates for a steady state overland flow using USPED (Unit Stream Power based Erosion Deposition) model. Net erosion/deposition is estimated as a change in sediment flow rate expressed by a divergence in sediment flow.


In North Carolina sample dataset, we compute net erosion/deposition.
g.region raster=elev_lid792_1m -p
r.slope.aspect elevation=elev_lid792_1m slope=slope aspect=aspect
r.flow elevation=elev_lid792_1m flowaccumulation=flowacc

# exponents m=1.3 and n=1.2
# multiply flowaccumulation by cell area/resolution to get contributing area per unit width
r.mapcalc "sflowtopo = pow(flowacc * 1.,1.3) * pow(sin(slope),1.2)"

# Compute sediment flow by combining the rainfall, soil and land cover factors
# with the topographic sediment transport factor.
# We use a constant value of 270. for rainfall intensity factor.
r.mapcalc "sedflow = 270. * soils_Kfactor * cfactorbare_1m * sflowtopo"
r.divergence magnitude=sedflow direction=aspect output=erosion_deposition

# set suitable color table
r.colors map=erosion_deposition rules=- << EOF
0% 100 0 100   #dark magenta
-100 magenta
-10 red
-1 orange
-0.1 yellow
0 200 255 200     #light green
0.1 cyan
1 aqua
10 blue
100 0 0 100       #dark blue
100% black




Anna Petrasova, NCSU OSGeoREL
Helena Mitasova, NCSU OSGeoREL


Available at: r.divergence source code (history)

Latest change: Thursday Feb 03 09:32:35 2022 in commit: f17c792f5de56c64ecfbe63ec315307872cf9d5c

Main index | Raster index | Topics index | Keywords index | Graphical index | Full index

© 2003-2024 GRASS Development Team, GRASS GIS 8.3.3dev Reference Manual