Present a SoilProfileCollection aligned to a geomorphic summary as cross-section.

plotGeomorphCrossSection(
  x,
  type = c("line", "bar"),
  g = "hillpos",
  clust = TRUE,
  col = c("#377EB8", "#4DAF4A", "#984EA3", "#FF7F00", "#E41A1C"),
  ...
)

Arguments

x

resulting list from soilDB::fetchOSD(..., extended = TRUE)

type

character, 'line' for line plot or 'bar' for barplot of geomorphic proportions

g

character, select a geomorphic summary. Currently 'hillpos' (2D hillslope position) is the only supported choice.

clust

logical, use clustering order of geomorphic proportions (TRUE) or exact hydrologic ordering (FALSE), see hydOrder()

col

character vector of colors

...

additional arguments to iterateHydOrder()

Details

Additional arguments to aqp::plotSPC() can be provided using options(.aqp.plotSPC.args = list(...)). For example, adjustments to maximum depth and profile width can be set via: options(.aqp.plotSPC.args = list(max.depth = 150, width = 0.35). Default arguments can be reset with options(.aqp.plotSPC.args = NULL).

When clust = TRUE, especially for SoilProfileCollections with a wide range in depth, it may be necessary to adjust the scaling.factor argument to aqp::plotSPC() via: options(.aqp.plotSPC.args = list(scaling.factor = 0.01)). Larger values will increase the height of profile sketches.

Author

D.E. Beaudette