Estimation of compressive properties by O'Sullivan et al. (1999)
It calculates the compressive parameteres N, lambda and kappa using the pedo-transfer function from O'Sullivan et al. (1999)
compressive_properties5(water.content, soil=c("SandyLoam","ClayLoam"))water.content |
a numeric vector containing the values of gravimetric water content, \% |
soil |
the the soil texture group 'SandyLoam' or 'ClayLoam'. See exemples. |
See O'Sullivan et al. (1999).
N |
the specific volume at p = 1 kPa, N |
CI |
the compression index, lambda |
k |
the recompression index, kappa |
Renato Paiva de Lima <renato_agro_@hotmail.com> Anderson Rodrigo da Silva <anderson.agro@hotmail.com>
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.
# 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)Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.