compressive_properties4: Estimation of compressive properties by de Lima et al. (2020)

Description Usage Arguments Details Value Author(s) References See Also Examples

View source: R/compressive_properties4.R

Description

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

Usage

1
2
 
compressive_properties4(matric.suction, soil=c("PloughLayer","PloughPan"))

Arguments

matric.suction

a numeric vector containing the values of matric suction, hPa.

soil

the soil compaction state 'PloughLayer' or 'PloughPan'. See the examples.

Details

Pedo-transfer function developed for a sandy loam soil texture. 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 <renato_agro_@hotmail.com> Anderson Rodrigo da Silva <anderson.agro@hotmail.com>

References

de Lima, R. P., Rolim, M. M., da C. Dantas, D., da Silva, A. R., Mendonca, E. A., 2020. Compressive properties and least limiting water range of plough layer and plough pan in sugarcane fields. Soil Use and Management, 00: 1-12.

See Also

stressTraffic

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
# EXAMPLE 1

compressive_properties4(matric.suction=100, soil="PloughLayer")
compressive_properties4(matric.suction=100, soil="PloughPan")



# EXAMPLE 2 for "PloughLayer"

matric.suction <- seq(from=10,to=10000,len=100)
out <- compressive_properties4(matric.suction=matric.suction, soil="PloughLayer")
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 "PloughPan"

matric.suction <- seq(from=10,to=10000,len=100)
out <- compressive_properties4(matric.suction=matric.suction, 
                               soil="PloughPan")
# plot for N
plot(x=matric.suction,y=out$N, 
     ylab="N", xlab="Matric suction (hPa)", log="x") 
# 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)

soilphysics documentation built on Jan. 13, 2021, 7:06 p.m.