SoilWaterStorage: Soil Moisture Mapping

Description Usage Arguments Examples

Description

Vertical or Slope-Normal Integration/Aggregation of soil water content

Usage

1
2
3
SoilWaterStorage(theta, psi, layers = 1, fun = sum, na.rm = TRUE,
  psi_thres = 0, comparison = c("<", ">", "<=", ">=", "==", "!="),
  indices = 1, bedrock.depth = NULL, ...)

Arguments

theta

soil water content RasterBrick-class object

psi

soil pressure head content RasterBrick-class object

layers

layer thickness. It is a scalar or a vector whose length is equal to the number of theta layers.

fun

function aggregation

na.rm

a logical value for fun indicating whether NA values should be stripped before the computation proceeds.

psi_thres

threshold value. It is 0 for separating satureted and unsatureted zones, or NA if the two zones are not distinguished.

comparison

comparison operator for psi .See Comparison. Default is the first element of c("<",">","<=",">=","==","!=")

indices

see stackApply

bedrock.depth

map of deprock depth. If it is not NULL, theta cells under bedrock are not considered.

...

further aguments for stackApply

Examples

1
2
3
4
5
6
data(PsiTheta)

 t <- 5
 dz <- 50 ## 50 millimiters layer depth
 unsatWaterVolume <- SoilWaterStorage(theta[[t]],psi[[t]],layer=dz,comparison="<",psi_thres=0)
 satWaterVolume <- SoilWaterStorage(theta[[t]],psi[[t]],layer=dz,comparison=">=",psi_thres=0)

ecor/geotopsim documentation built on May 15, 2019, 10:05 p.m.