Skip to contents

Argic / argillic horizon via aqp::getArgillicBounds()

Source: R/aqp-interop-v0962.R
argic_aqp.Rd

Wraps aqp::getArgillicBounds() (Beaudette et al.) in soilKey's DiagnosticResult contract. The aqp implementation is the canonical NRCS R port and uses the tiered USDA-NRCS clay-increase thresholds:

  • Eluvial clay < 15\

  • Eluvial clay 15-40\

  • Eluvial clay \>= 40\

(vs. soilKey's hand-coded argic which uses the WRB 6/1.4/20 thresholds). For BDsolos / FEBR / KSSL profiles the aqp rule is closer to KST 13ed and BDsolos field practice.

Usage

argic_aqp(pedon, require_t = FALSE, ...)

Arguments

pedon

A PedonRecord.

require_t

Whether to require an explicit "t" suffix in the horizon designation (default FALSE for BDsolos / FEBR; TRUE matches the strict KST 13ed text).

...

Reserved for future arguments.

Value

A DiagnosticResult with name = "argic_aqp". $layers are the row indices of horizons in the argillic / argic depth interval. $evidence carries the raw aqp c(ubound, lbound) bounds for traceability.

Details

By default aqp requires a "t" suffix in the horizon designation (require_t = TRUE); we expose this so callers can be permissive on datasets where designation is missing or non-conforming (BDsolos exports often drop the "t").

See also

argic (soilKey hand-coded; WRB 6/1.4/20), aqp::getArgillicBounds.