Installation

Get the stable version (2.6.5) from CRAN:

install.packages('soilDB', dependencies = TRUE)

Get the development version (2.6.5) from GitHub:

remotes::install_github("ncss-tech/soilDB", dependencies = FALSE, upgrade = FALSE, build = FALSE)

Website

• http://ncss-tech.github.io/AQP/

Citation

Dylan Beaudette, Jay Skovlin, Stephen Roecker and Andrew Brown (2021). soilDB: Soil Database Interface. R package version 2.6.5. https://CRAN.R-project.org/package=soilDB

soilDB 2.6.5

Notices on Database Interfaces

NASIS

• all NASIS queries now use {DBI}, tested to work with {odbc} or {RSQLite} drivers; replacing {RODBC}
  • Install required R packages with install.packages(c("DBI","odbc","RSQLite"), dependencies = TRUE)
• new methods for connecting to NASIS and querying NASIS data allow for dsn argument to specify a local “static” SQLite file containing NASIS tables.
  • Default argument dsn = NULL uses "nasis_local" ODBC connection to a local NASIS SQL Server instance

Soil Data Access (SDA)

• SDA_query() returns a try-error for queries with invalid syntax or on network error; empty results are an empty data.frame()

• SDA_spatialQuery() can return Soil Survey Area sapolygon data

• fetchSDA_spatial() can return STATSGO gsmmupolygon or Soil Survey Area sapolygon data; and can join to the legend table

• Added several new SDA query methods based on ssurgoOnDemand by Jason Nemecek and Chad Ferguson:

  • get_SDA_property(), get_SDA_interpretation(), get_SDA_muaggatt(), get_SDA_hydric(), get_SDA_pmgroupname(), get_SDA_coecoclass()

MS Access

• AKSite, Montana RangeDB, and PedonPC 6.x get and fetch methods now use {DBI}+{odbc}

SoilWeb API

• ISSR800.wcs() and mukey.wcs() now return a result that inherits from try-error (and a message) if the Web Coverage Service fails

Functions by Data Source

• SDA
  • fetchSDA
  • fetchSDA_spatial
  • SDA_query
  • SDA_spatialQuery
  • “SSURGO On Demand” Queries
    • get_SDA_hydric
    • get_SDA_interpretation
    • get_SDA_muaggatt
    • get_SDA_pmgroupname
    • get_SDA_property
• SSURGO Local Geodatabases
  • fetchGDB
• ROSETTA
  • ROSETTA
• SSURGO/KSSL via SoilWeb
  • fetchKSSL
  • fetchOSD
  • siblings
  • OSDquery
  • seriesExtent
  • taxaExtent
  • mukey.wcs
  • ISSR800.wcs
• NASIS WWW interface
  • parseWebReport
  • fetchNASISWebReport
• SCAN/SNOTEL
  • fetchSCAN
  • SCAN_SNOTEL_metadata
• Henry Mount Soil and Water Database
  • fetchHenry
• NASIS local database
  • fetchNASIS
  • NEW: Argument dsn to specify path to connect to alternate (SQLite) data sources with NASIS schema
  • NEW:dbConnectNASIS (alias NASIS) - create a DBIConnection to local NASIS database
  • NEW:dbQueryNASIS - query NASIS local database (and close connection with close=TRUE)
  • NEW:createStaticNASIS - create list of NASIS tables or write to SQLite
• SoilGrids
  • fetchSoilGrids

Miscellaneous Functions

• estimateSTR
• STRplot
• KSSL_VG_model
• simplfyFragmentData
• simplifyColorData
• uncode
• code
• get_NOAA_GHCND

Tutorials and Demonstrations

• fetchKSSL
• SDA_query
• fetchOSD
• SCAN/SNOTEL Data

Related Packages

• aqp
• sharpshootR

Examples

Load Soil Packages

library(aqp)
library(soilDB)
library(sharpshootR)

Get Official Series Description Data

# get morphology + extended summaries
soils <- c('cecil', 'altavista', 'lloyd', 'wickham', 'woodbridge', 'chewacla', 'congaree')
s <- fetchOSD(soils, extended = TRUE)

Visualize Soil-Landscape Relationships

# viz of hillslope position, from SSURGO component records
res <- vizHillslopePosition(s$hillpos, annotation.cex = 0.9)
print(res$fig)

Make Profile Sketches

par(mar = c(0, 1, 0, 4), xpd = NA)

plotSPC(
  s$SPC,
  plot.order = res$order,
  cex.names = 1,
  axis.line.offset = -0.1,
  width = 0.2
)

Identify Tabular “Siblings”

s <- 'Amador'
amador <- siblings(s, only.major = FALSE, component.data = TRUE)

# limit to named soil series
sib.data <- subset(amador$sib.data, !compkind %in% c('Miscellaneous area', 'Family', 'Taxon above family'))

# get parsed OSD records
sibs <- fetchOSD(c(s, unique(amador$sib$sibling)), extended = TRUE)

Plot Dendrograms with Taxonomic Relationships

# invert colors
par(mar = c(0, 0, 0, 0),
    fg = 'white',
    bg = 'black')

SoilTaxonomyDendrogram(
  sibs$SPC,
  dend.width = 1.5,
  y.offset = 0.4,
  scaling.factor = 0.02,
  width = 0.2,
  cex.taxon.labels = 1,
  cex.names = 1
)

Dependency Graph