compressive_properties5 {soilphysics}

Estimation of compressive properties by O'Sullivan et al. (1999)

Description

It calculates the compressive parameteres N, lambda and kappa using the pedo-transfer function from O'Sullivan et al. (1999)

Usage

 
compressive_properties5(water.content, soil=c("SandyLoam","ClayLoam"))

Arguments

water.content

a numeric vector containing the values of gravimetric water content, \%

soil

the the soil texture group 'SandyLoam' or 'ClayLoam'. See exemples.

Details

See O'Sullivan et al. (1999).

Value

N

the specific volume at p = 1 kPa, N

CI

the compression index, lambda

k

the recompression index, kappa

Author(s)

Renato Paiva de Lima <[e-mail]> Anderson Rodrigo da Silva <[e-mail]>

References

O'sullivan, M. F., Henshall, J. K., Dickson, J. W., 1999. A simplified method for estimating soil compaction. Soil and Tillage Research, 49: 325-335.

See Also

stressTraffic

Examples

# EXAMPLE 1
water.content <- 15
compressive_properties5(water.content=water.content, soil="SandyLoam") 
compressive_properties5(water.content=water.content, soil="ClayLoam") 


# EXAMPLE 2 - SandyLoam

water.content <- seq(from=5,to=20,len=20)
out <- compressive_properties5(water.content=water.content, soil="SandyLoam") 
plot(x=water.content ,y=out$N, 
     ylab="N", xlab="Bulk density") # plot for N
plot(x=water.content ,y=out$lambda, 
    ylab="lambda", xlab="Bulk density") # plot for lambda
plot(x=water.content ,y=out$kappa, 
    ylab="kappa", xlab="Bulk density") # plot for kappa


# EXAMPLE 3 - ClayLoam 

water.content <- seq(from=10,to=25,len=20)
out <- compressive_properties5(water.content=water.content, soil="ClayLoam") 
plot(x=water.content ,y=out$N, 
     ylab="N", xlab="Bulk density") # plot for N
plot(x=water.content ,y=out$lambda, 
     ylab="lambda", xlab="Bulk density") # plot for lambda
plot(x=water.content ,y=out$kappa, 
     ylab="kappa", xlab="Bulk density") # plot for rkappa

# End (not run)