Skip to contents

aqp's getArgillicBounds() requires an NRCS texture class column (e.g. "SCL", "C", "CL", "FS"). soilKey horizons only carry the percent fractions; this helper derives the class from the standard USDA texture triangle.

Usage

texture_class_from_pct(clay, silt, sand)

Arguments

clay

Numeric vector of clay percent (0-100).

silt

Numeric vector of silt percent.

sand

Numeric vector of sand percent. (clay + silt + sand should sum to ~100; mild deviations are tolerated.)

Value

Character vector of NRCS texture class abbreviations.

Details

Returns the standard NRCS abbreviation:

COSCoarse sand
SSand
FSFine sand
VFSVery fine sand
LSLoamy sand
LFSLoamy fine sand
SLSandy loam
FSLFine sandy loam
LLoam
SILSilt loam
SISilt
SCLSandy clay loam
CLClay loam
SICLSilty clay loam
SCSandy clay
SICSilty clay
CClay

Implementation follows the canonical USDA texture triangle; vector- ised over the input. NA in / NA out.