compressive_properties3 {soilphysics}

Estimation of compressive properties by de Lima et al. (2018)

Description

It calculates the compressive parameters N, lambda and kappa using the pedo-transfer function from de Lima et al. (2018)

Usage

 
compressive_properties3(bulk.density, matric.suction, soil=c("SandyLoam","SandyClayLoam"))

Arguments

bulk.density

a numeric vector containing the values of bulk density, Mg m^{-3}

matric.suction

a numeric vector containing the values of matric suction, hPa

soil

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

Details

Pedo-transfer function deveploed under no-till condition. See de Lima et al. (2018)

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

de Lima, R. P., da Silva, A. P., Giarola, N. F., da Silva, A. R., Rolim, M. M., Keller, T., 2018. Impact of initial bulk density and matric suction on compressive properties of two Oxisols under no-till. Soil and Tillage Research, 175: 168-177.

See Also

stressTraffic

Examples

# EXAMPLE 1

compressive_properties3(bulk.density=1.5, matric.suction=100, soil="SandyLoam")
compressive_properties3(bulk.density=1.5, matric.suction=100, soil="SandyClayLoam")

# EXAMPLE 2 for SandyLoam soil

matric.suction <- seq(from=30,to=1000,len=100)
out <- compressive_properties3(bulk.density=1.5, 
	matric.suction=matric.suction, soil="SandyLoam")
plot(x=matric.suction,y=out$N, ylab="N", 
	xlab="Matric suction (hPa)", log="x") # plot for N
# plot for lambda
plot(x=matric.suction,y=out$lambda, ylab="lambda", 
	xlab="Matric suction (hPa)", log="x") 
# plot for kappa
plot(x=matric.suction,y=out$k, ylab="kappa", 
	xlab="Matric suction (hPa)", log="x") 


# EXAMPLE 3 for SandyClayLoam soil

matric.suction <- seq(from=30,to=1000,len=100)
out <- compressive_properties3(bulk.density=1.5, 
				matric.suction=matric.suction, 
				soil="SandyClayLoam")
plot(x=matric.suction,y=out$N, ylab="N", 
     xlab="Matric suction (hPa)", log="x") # plot for N
# plot for lambda
plot(x=matric.suction,y=out$lambda, 
     ylab="lambda", xlab="Matric suction (hPa)", log="x")
# plot for kappa
plot(x=matric.suction,y=out$k, 
     ylab="kappa", xlab="Matric suction (hPa)", log="x")

# End (not run)