
Returns the top depth of the argillic horizon as a numeric vector.
argillic.clay.increase.depth(p, clay.attr = "clay")
A single-profile SoilProfileCollection
object.
OPTIONAL: horizon attribute name referring to clay content.
default: clay
A numeric vector containing top depth of argillic horizon, if present, or NA.
Uses crit.clay.argillic
to determine threshold clay increase, and
get.increase.matrix
to determine where increase is met within a
vertical distance of 30 cm.
getArgillicBounds
, get.increase.matrix
,
crit.clay.argillic
# NOT RUN {
data(sp1, package = 'aqp')
depths(sp1) <- id ~ top + bottom
site(sp1) <- ~ group
p <- sp1[1]
attr <- 'prop' # clay contents
foo <- argillic.clay.increase.depth(p, clay.attr = attr)
foo
# }
Run the code above in your browser using DataLab