Learn R Programming

HydroMe (version 2.1.1)

Tani: A Tani water retention model

Description

This is a three- parameter function to determine water retention hydraulic parameters contained in the Tani water retention model

Usage

Tani(x, thr, ths, alp)

Value

A list of objects returned by standard R regression functions such as nls and lm

Arguments

x

Suction potential/head as contained in the x-column of the xy water retention table or data

thr

This is the residual moisture content. It's the moisture content when suction potential is very high (almost at the drying point)

ths

This is the saturated moisture content. It's the moisture content when suction potential is very low (almost at the saturation point)

alp

It's the inverse of air-entry potential or bubbling pressure

Author

Christian Thine Omuto

Details

Tani model is a type of nonlinear curve fitting model for fitting water retention characteristics using experimental data. It requires initial parameter estimates for the model to work. This can be done using Dstart function (e.g. thr=Dstart(data)[1]). If warnings are given during the estimation process, try to increase the iteration limit using control function. Sometimes the warnings can be suppressed using (warn=-1)/ignored since they do not terminate the process

References

Tani M. 1982. The properties of water-table rise produced by a one-dimensional, vertical, unsaturated flow. Journal of Japan Forestry Society 64: 409-418

See Also

Brook, Ruso, Campbel

Examples

Run this code
data(isric)
pf <- subset(isric, Sample=="Italy4")
library(minpack.lm)
tani.ns <- nlsLM(y ~ Tani(x,thr,ths,alp), data = pf,
                 control = nls.lm.control(maxiter=200),
                 start = c(thr=Dstart(pf)[1],ths=Dstart(pf)[2], alp=Dstart(pf)[3]))
coef(tani.ns)

Run the code above in your browser using DataLab