Calculate the "Profile Darkness Index" by the method of Thompson & Bell (1996) "Color index for identifying hydric conditions for seasonally saturated mollisols in Minnesota" DOI: 10.2136/sssaj1996.03615995006000060051x. The Thompson-Bell Index has been shown to reflect catenary relationships in some Mollisols of Minnesota (generally: wetter landscape positions = thicker, darker surfaces).
thompson.bell.darkness(
p,
name = NULL,
pattern = "^A",
value = "m_value",
chroma = "m_chroma"
)
A single-profile SoilProfileCollection (e.g. via profileApply())
Column name containing horizon designations used to find A horizons (default: first column name containing 'name')
Regular expression to match A horizons (default: "^A" which means horizon designation _starts with_ A)
Column name containing horizon color values (default: "m_value")
Column name containing horizon color chromas (default: "m_chroma")
A numeric vector reflecting horizon redness (higher values = redder).