soilphysics (version 2.1)

fitsoilwater4: Self-starting Nls Power Models for Soil Water Retention

Description

Function to self start the following nonlinear power models for soil water retention: $$\theta = \exp(a + b*Bd) \psi^c$$ (Silva et al., 1994) $$\theta = a \psi^c$$ (Ross et al., 1991) where $\theta$ is the soil water content. fitsoilwater() creates initial estimates (by log-linearization) of the parameters a, b and c and uses them to provide its least-squares estimates through nls.

Usage

fitsoilwater4(theta, psi, Bd, model = c("Silva", "Ross"))

Arguments

theta

a numeric vector containing values of soil water content.

psi

a numeric vector containing values of water potential (Psi).

Bd

a numeric vector containing values of dry bulk density.

model

a character; the model to be used for calculating the soil water content. It must be one of the two: "Silva" (default) or "Ross".

Value

• A "nls" object containing the fitted model.

References

Ross et al. (1991). Equation for extending water-retention curves to dryness. Soil Science Society of America Journal, 55:923-927. Silva et al. (1994). Characterization of the least limiting water range of soils. Soil Science Society of America Journal, 58:1775-1781.

See Also

fitsoilwater4, soilwater, soilwater2, soilwater3

Examples

# Example 1
ex1 <- with(skp1994,
    fitsoilwater4(theta = Water, psi = Potential, model = "Ross"))
ex1
summary(ex1)

# Example 2
ex2 <- with(skp1994,
    fitsoilwater4(theta = Water, psi = Potential, Bd = BD, model = "Silva"))
ex2
summary(ex2)

# Not run