unitResponse {soilwater}

The water table recharge: the response unit

Description

The water table recharge: the response unit

Usage

unitResponse(t, d = 1, D = 1, H = d, m = 100)

Arguments

t

time coordinate

d

depth of unsaturated zone along the slope-normal direction

D

soil water diffusivity

H

soil depth

m

maximum limit of summary truncation. Default is 100.

Note

This function calcletes the water-table recharge rate in a hillslope assuming:

1. Richards' Equation is linearized and reduced to the form of heat equation;

2. The diffusion water-table rate is connectedwith soil pressure head according with eq. 13 (Cordano and Rigon, 2008);

References

Cordano, E., and R. Rigon (2008), A perturbative view on the subsurface water pressure response at hillslope scale, Water Resour. Res., 44, W05407, doi:10.1029/2006WR005740. http://onlinelibrary.wiley.com/doi/10.1029/2006WR005740/pdf

Examples

library(soilwater)


t <- seq(0,2,by=0.001)
d <- c(1,0.75,0.5,0.25)
val1 <- unitResponse(t, d = d[1], D = 1, H = 1, m = 500)

val2 <- unitResponse(t, d = d[2], D = 1, H = 1, m = 500)

val3 <- unitResponse(t, d = d[3], D = 1, H = 1, m = 500)

val4 <- unitResponse(t, d = d[4], D = 1, H = 1, m = 500)