compressive_properties {soilphysics}

Estimation of compressive properties by Defossez et al. (2003)

Description

It calculates the compressive parameters N and lambda using the pedo-transfer function from Defossez et al. (2003)

Usage

 
compressive_properties(water.content, soil=c("Loess", "Calcareous"))

Arguments

water.content

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

soil

the soil group 'Loess' or 'Calcareous'. See examples

Details

In Defossez et al. (2003), the recompression index, kappa, was assumed as 0.0058 for both soil group.

Value

N

the specific volume at p = 1 kPa, N

CI

the compression index, lambda

Author(s)

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

References

Defossez, P., Richard, G., Boizard, H., & O'Sullivan, M. F., 2003. Modeling change in soil compaction due to agricultural traffic as function of soil water content. Geoderma, 116: 89–105.

See Also

stressTraffic

Examples

# EXAMPLE 1 - For Loess and Calcareous soil

water.content <- 25
compressive_properties(water.content=water.content, soil="Loess") 
compressive_properties(water.content=water.content, soil="Calcareous") 

# EXAMPLE 2 - For Loess soil

water.content <- seq(from=5,to=30,len=20)
out <- compressive_properties(water.content=water.content, soil="Loess") 
plot(x=water.content ,y=out$N, ylab="N", xlab="Bulk density") # plot for N
plot(x=water.content ,y=out$CI, ylab="Lambda", xlab="Bulk density") # plot for compression index


# EXAMPLE 3 - For Calcareous soil

water.content <- seq(from=5,to=30,len=20)
out <- compressive_properties(water.content=water.content, soil="Calcareous") 
plot(x=water.content ,y=out$N, ylab="N", xlab="Bulk density") # plot for N
plot(x=water.content ,y=out$CI, ylab="Lambda", xlab="Bulk density") # plot for compression index

# End (not run)