dynatopmodel (version 1.2.1)

upslope.area: Upslope contributing area and wetness index calculation

Description

Determine upslope contributing area based on an elevation raster and, optionally, compute the topographic wetness index.

Usage

upslope.area(dem, log = TRUE, atb = FALSE, deg = 0.1, fill.sinks = TRUE)

Arguments

dem

raster Elevation raster (in m), using a projected coordinate system with identical x and y resolutions.

log

Boolean Return the natural log of the values.

atb

Boolean If TRUE, include both the upslope contributing area and the topographic wetness index \(ln(a/tan(beta))\). Otherwise calculate just the upslope area.

deg

numeric Minimum intercell slope to identify with a sink (degrees).

fill.sinks

Fill sinks before calculation using the threshold angle given by deg.

References

Quinn, P. FALSE., Beven, K. J., & Lamb, R. (1995). The In (a/tan/beta) index: How to calculate it and how to use it within the Topmodel framework. Hydrological processes, 9(2), 161-182.

Examples

# NOT RUN {
require(dynatopmodel)
data(brompton)

a.atb <- upslope.area(brompton$dem, atb=TRUE)
sp::plot(a.atb, main=c("Upslope area (log(m^2/m))", "TWI log(m^2/m)"))
# }