Package 'soilDB'

Description

A collection of functions for reading soil data from U.S. Department of Agriculture Natural Resources Conservation Service (USDA-NRCS) and National Cooperative Soil Survey (NCSS) databases

Details

This package provides methods for extracting soils information from local NASIS databases (MS SQL Server), local PedonPC and AKSite databases (MS Access format), Soil Data Access, and other soil-related web services.

Author(s)

J.M. Skovlin, D.E. Beaudette, S.M Roecker, A.G. Brown

See Also

fetchNASIS, SDA_query, loafercreek

Description

The following database types are tested and fully supported:

• SQLite or Geopackage

• DuckDB

• Postgres or PostGIS

In theory any other DBI-compatible data source can be used for output. See conn argument. If you encounter issues using specific DBI connection types, please report in the soilDB issue tracker.

Usage

createSSURGO(
  filename,
  exdir,
  conn = NULL,
  pattern = NULL,
  include_spatial = TRUE,
  overwrite = FALSE,
  header = FALSE,
  quiet = TRUE,
  ...
)

Arguments

Value

Character. Vector of layer/table names in filename.

See Also

downloadSSURGO()

Examples

## Not run: 
 downloadSSURGO("areasymbol IN ('CA067', 'CA077', 'CA632')", destdir = "SSURGO_test")
 createSSURGO("test.gpkg", "SSURGO_test")

## End(Not run)

Description

Create a memory or file-based instance of NASIS database for selected tables.

Usage

createStaticNASIS(
  tables = NULL,
  new_names = NULL,
  SS = TRUE,
  dsn = NULL,
  output_path = NULL,
  verbose = FALSE
)

Arguments

Value

A named list of results from calling dbQueryNASIS for all columns in each NASIS table.

Description

Create a connection to a local NASIS database with DBI

Usage

dbConnectNASIS(dsn = NULL)

NASIS(dsn = NULL)

Arguments

Value

A DBIConnection object, as returned by DBI::dbConnect(). If dsn is a DBIConnection, the attribute isUserDefined of the result is set to TRUE. If the DBIConnection is created by the internal NASIS connection process, isUserDefined is set to FALSE.

Description

Send queries to a NASIS DBIConnection

Usage

dbQueryNASIS(conn, q, close = TRUE, ...)

Arguments

Value

Result of DBI::dbGetQuery

Description

Download ZIP files containing spatial (ESRI shapefile) and tabular (TXT) files with standard SSURGO format; optionally including the corresponding SSURGO Template Database with include_template=TRUE.

Usage

downloadSSURGO(
  WHERE = NULL,
  areasymbols = NULL,
  destdir = tempdir(),
  exdir = destdir,
  include_template = FALSE,
  db = c("SSURGO", "STATSGO"),
  extract = TRUE,
  remove_zip = FALSE,
  overwrite = FALSE,
  quiet = FALSE
)

Arguments

Details

To specify the Soil Survey Areas you would like to obtain data you use a WHERE clause for query of sacatalog table such as areasymbol = 'CA067', "areasymbol IN ('CA628', 'CA067')" or ⁠areasymbol LIKE 'CT%'⁠.

When db="STATSGO" the WHERE argument is not supported. Allowed areasymbols include "US" and two-letter state codes e.g. "WY" for the Wyoming general soils map.

Pipe-delimited TXT files are found in /tabular/ folder extracted from a SSURGO ZIP. The files are named for tables in the SSURGO schema. There is no header / the files do not have column names. See the Soil Data Access Tables and Columns Report: https://sdmdataaccess.nrcs.usda.gov/documents/TablesAndColumnsReport.pdf for details on tables, column names and metadata including the default sequence of columns used in TXT files. The function returns a try-error if the WHERE/areasymbols arguments result in

Several ESRI shapefiles are found in the /spatial/ folder extracted from a SSURGO ZIP. These have prefix soilmu_ (mapunit), soilsa_ (survey area), soilsf_ (special features). There will also be a TXT file with prefix soilsf_ describing any special features. Shapefile names then have an a_ (polygon), l_ (line), p_ (point) followed by the soil survey area symbol.

Value

Character. Paths to downloaded ZIP files (invisibly). May not exist if remove_zip = TRUE.

See Also

createSSURGO()

Description

Estimate color mixtures using weighted average of CIELAB color coordinates

Usage

estimateColorMixture(x, wt = "pct", backTransform = FALSE)

Arguments

Value

A data.frame containing estimated color mixture

Note

See aqp::mixMunsell() for a more realistic (but slower) simulation of subtractive mixing of pigments. An efficient replacement for this function (wt. mean in CIELAB coordinates) is implemented in aqp::mixMunsell(..., mixingMethod = 'estimate').

Author(s)

D.E. Beaudette

Description

Estimate soil temperature regime (STR) based on mean annual soil temperature (MAST), mean summer temperature (MSST), mean winter soil temperature (MWST), presence of O horizons, saturated conditions, and presence of permafrost. Several assumptions are made when O horizon or saturation are undefined.

Usage

estimateSTR(
  mast,
  mean.summer,
  mean.winter,
  O.hz = NA,
  saturated = NA,
  permafrost = FALSE
)

Arguments

Details

Soil Temperature Regime Evaluation Tutorial

Value

Vector of soil temperature regimes.

Author(s)

D.E. Beaudette

References

Soil Survey Staff. 2015. Illustrated guide to soil taxonomy. U.S. Department of Agriculture, Natural Resources Conservation Service, National Soil Survey Center, Lincoln, Nebraska.

See Also

STRplot

Examples

# simple example
estimateSTR(mast=17, mean.summer = 22, mean.winter = 12)

Description

Download soil characterization and morphologic data via BBOX, MLRA, or soil series name query, from the KSSL database.

Usage

fetchKSSL(
  series = NA,
  bbox = NA,
  mlra = NA,
  pedlabsampnum = NA,
  pedon_id = NA,
  pedon_key = NA,
  returnMorphologicData = FALSE,
  returnGeochemicalData = FALSE,
  simplifyColors = FALSE,
  progress = TRUE
)

Arguments

Details

This interface has largely been superseded by the Soil Data Access snapshot of the Laboratory Data Mart, available via fetchLDM().

Series-queries are case insensitive. Series name is based on the "correlated as" field (from KSSL snapshot) when present. The "sampled as" classification was promoted to "correlated as" if the "correlated as" classification was missing.

When returnMorphologicData is TRUE, the resulting object is a list. The standard output from fetchKSSL (SoilProfileCollection object) is stored in the named element "SPC". The additional elements are basic morphologic data: soil color, rock fragment volume, pores, structure, and redoximorphic features. There is a 1:many relationship between the horizon data in "SPC" and the additional dataframes in morph. See examples for ideas on how to "flatten" these tables.

When returnGeochemicalData is TRUE, the resulting object is a list. The standard output from fetchKSSL (SoilProfileCollection object) is stored in the named element "SPC". The additional elements are geochemical and mineralogy analysis tables, specifically: geochemical/elemental analyses "geochem", optical mineralogy "optical", and X-ray diffraction / thermal "xrd_thermal". returnGeochemicalData will include additional dataframes geochem, optical, and xrd_thermal in list result.

Setting simplifyColors=TRUE will automatically flatten the soil color data and join to horizon level attributes.

Function arguments (series, mlra, etc.) are fully vectorized except for bbox.

Value

a SoilProfileCollection object when returnMorphologicData is FALSE, otherwise a list.

Note

SoilWeb maintains a snapshot of these KSSL and NASIS data. The SoilWeb snapshot was developed using methods described here: https://github.com/dylanbeaudette/process-kssl-snapshot. Please use the link below for the live data.

Author(s)

D.E. Beaudette and A.G. Brown

References

http://ncsslabdatamart.sc.egov.usda.gov/

See Also

fetchOSD

Examples

library(aqp)

    # search by series name
    s <- fetchKSSL(series='auburn')

    # search by bounding-box
    # s <- fetchKSSL(bbox=c(-120, 37, -122, 38))

    # how many pedons
    length(s)

    # plot
    plotSPC(s, name='hzn_desgn', max.depth=150)

    ##
    ## morphologic data
    ##

    # get lab and morphologic data
    s <- fetchKSSL(series='auburn', returnMorphologicData = TRUE)

    # extract SPC
    pedons <- s$SPC

   # if (requireNamespace("farver")) {
   #   ## automatically simplify color data (requires farver)
   #   s <- fetchKSSL(series='auburn', returnMorphologicData = TRUE, simplifyColors=TRUE)
   #   # check
   #   par(mar=c(0,0,0,0))
   #   plot(pedons, color='moist_soil_color', print.id=FALSE)
   # }

Description

LDM model diagram: https://jneme910.github.io/Lab_Data_Mart_Documentation/Documents/SDA_KSSL_Data_model.html

Usage

fetchLDM(
  x = NULL,
  what = "pedlabsampnum",
  bycol = "pedon_key",
  tables = c("lab_physical_properties", "lab_chemical_properties",
    "lab_calculations_including_estimates_and_default_values", "lab_rosetta_Key"),
  WHERE = NULL,
  chunk.size = 1000,
  ntries = 3,
  layer_type = c("horizon", "layer", "reporting layer"),
  area_type = c("ssa", "country", "state", "county", "mlra", "nforest", "npark"),
  prep_code = c("S", ""),
  analyzed_size_frac = c("<2 mm", ""),
  dsn = NULL
)

Arguments

Details

If the chunk.size parameter is set too large and the Soil Data Access request fails, the algorithm will re-try the query with a smaller (halved) chunk.size argument. This will be attempted up to 3 times before returning NULL

Currently the lab_area tables are joined only for the "Soil Survey Area" related records.

When requesting data from "lab_major_and_trace_elements_and_oxides", "lab_mineralogy_glass_count", or "lab_xray_and_thermal" multiple preparation codes (prep_code) or size fractions (analyzed_size_frac) are possible. The default behavior of fetchLDM() is to attempt to return a topologically valid (minimal overlaps) SoilProfileCollection. This is achieved by setting prep_code="S" ("sieved") and analyzed_size_frac="<2 mm". You may specify alternate or additional preparation codes or fractions as needed, but note that this may cause "duplication" of some layers where measurements were made with different preparation or on fractionated samples

Value

a SoilProfileCollection for a successful query, a try-error if no site/pedon locations can be found or NULL for an empty lab_layer (within sites/pedons) result

Examples

## Not run: 
  # fetch by ssa_key
  res <- fetchLDM(8297, what = "ssa_key")

  # physical properties correlated as taxonomic subgroup "Typic Argialbolls"
  res <- fetchLDM(x = "Typic Argialbolls", 
                  what = "corr_taxsubgrp", 
                  tables = "lab_physical_properties")

  # fetch by area_code (SSA only)  
  res <- fetchLDM("CA630", what = "area_code")

## End(Not run)

Description

Fetch commonly used site/pedon/horizon or mapunit component data from NASIS, returned as a SoilProfileCollection object.

This function imports data from NASIS into R as a SoilProfileCollection object. It "flattens" NASIS pedon and component tables, including their child tables, into several more manageable data frames. Primarily these functions access the local NASIS database using an ODBC connection. The dsn argument allows you to specify a path or DBIConnection to an SQLite database. The argument from = "pedon_report", data can be read from the NASIS Report 'fetchNASIS', from either text file or URL (specified as url). The primary purpose of fetchNASIS(from = "pedon_report") is importing datasets larger than 8000+ pedons/components.

The value of nullFragsAreZero will have a significant impact on the rock fragment fractions returned by fetchNASIS. Set nullFragsAreZero = FALSE in those cases where there are many data-gaps and NULL rock fragment values should be interpreted as NULL. Set nullFragsAreZero = TRUE in those cases where NULL rock fragment values should be interpreted as 0.

This function attempts to do most of the boilerplate work when extracting site/pedon/horizon or component data from a local NASIS database. Pedon IDs that are missing horizon data, or have errors in their horizonation are printed on the console. Pedons with combination horizons (e.g. B/C) are erroneously marked as errors due to the way in which they are stored in NASIS as two overlapping horizon records.

Tutorials:

• fetchNASIS Columns

• fetchNASIS Pedons Tutorial

• fetchNASIS Components Tutorial

Usage

fetchNASIS(
  from = "pedons",
  url = NULL,
  SS = TRUE,
  rmHzErrors = FALSE,
  nullFragsAreZero = TRUE,
  soilColorState = "moist",
  mixColors = TRUE,
  lab = FALSE,
  fill = FALSE,
  dropAdditional = TRUE,
  dropNonRepresentative = TRUE,
  duplicates = FALSE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_concentrations_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_phfmp_from_NASIS_db(SS = TRUE, stringsAsFactors = NULL, dsn = NULL)

Arguments

Value

a SoilProfileCollection object

Author(s)

D. E. Beaudette, J. M. Skovlin, S.M. Roecker, A.G. Brown

See Also

get_component_data_from_NASIS()

Description

Fetch KSSL laboratory pedon/horizon layer data from a local NASIS database, return as a SoilProfileCollection object.

Usage

fetchNASISLabData(SS = TRUE, dsn = NULL)

Arguments

Details

This function currently works only on Windows, and requires a 'nasis_local' ODBC connection.

Value

a SoilProfileCollection object

Author(s)

J.M. Skovlin and D.E. Beaudette

See Also

get_labpedon_data_from_NASIS_db

Description

Get component tables from NASIS Web Reports

Usage

fetchNASISWebReport(
  projectname,
  rmHzErrors = FALSE,
  fill = FALSE,
  stringsAsFactors = NULL
)

get_component_from_NASISWebReport(projectname, stringsAsFactors = NULL)

get_chorizon_from_NASISWebReport(
  projectname,
  fill = FALSE,
  stringsAsFactors = NULL
)

get_legend_from_NASISWebReport(
  mlraoffice,
  areasymbol,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_lmuaoverlap_from_NASISWebReport(
  areasymbol,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_mapunit_from_NASISWebReport(
  areasymbol,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_projectmapunit_from_NASISWebReport(projectname, stringsAsFactors = NULL)

get_projectmapunit2_from_NASISWebReport(
  mlrassoarea,
  fiscalyear,
  projectname,
  stringsAsFactors = NULL
)

get_project_from_NASISWebReport(mlrassoarea, fiscalyear)

get_progress_from_NASISWebReport(mlrassoarea, fiscalyear, projecttypename)

get_project_correlation_from_NASISWebReport(
  mlrassoarea,
  fiscalyear,
  projectname
)

get_cosoilmoist_from_NASISWebReport(
  projectname,
  impute = TRUE,
  stringsAsFactors = NULL
)

get_sitesoilmoist_from_NASISWebReport(usiteid)

Arguments

Value

A data.frame or list with the results.

Author(s)

Stephen Roecker

Description

This function fetches a variety of data associated with named soil series, extracted from the USDA-NRCS Official Series Description text files and detailed soil survey (SSURGO). These data are updated quarterly and made available via SoilWeb. Set extended = TRUE and see the soilweb.metadata list element for information on when the source data were last updated.

Usage

fetchOSD(soils, colorState = "moist", extended = FALSE)

Arguments

Details

• overview of all soil series query functions

• competing soil series

• siblings

The standard set of "site" and "horizon" data are returned as a SoilProfileCollection object (extended = FALSE). The "extended" suite of summary data can be requested by setting extended = TRUE. The resulting object will be a list with the following elements:

SPC

SoilProfileCollection containing standards "site" and "horizon" data

competing

competing soil series from the SC database snapshot

geog_assoc_soils

geographically associated soils, extracted from named section in the OSD

geomcomp

empirical probabilities for geomorphic component, derived from the current SSURGO snapshot

hillpos

empirical probabilities for hillslope position, derived from the current SSURGO snapshot

mtnpos

empirical probabilities for mountain slope position, derived from the current SSURGO snapshot

terrace

empirical probabilities for river terrace position, derived from the current SSURGO snapshot

flats

empirical probabilities for flat landscapes, derived from the current SSURGO snapshot

shape_across

empirical probabilities for surface shape (across-slope) from the current SSURGO snapshot

shape_down

empirical probabilities for surface shape (down-slope) from the current SSURGO snapshot

pmkind

empirical probabilities for parent material kind, derived from the current SSURGO snapshot

pmorigin

empirical probabilities for parent material origin, derived from the current SSURGO snapshot

mlra

empirical MLRA membership values, derived from the current SSURGO snapshot

ecoclassid

area cross-tabulation of ecoclassid by soil series name, derived from the current SSURGO snapshot, major components only

climate

climate summaries from PRISM stack (CONUS only)

NCCPI

select quantiles of NCCPI and Irrigated NCCPI, derived from the current SSURGO snapshot

metadata

metadata associated with SoilWeb cached summaries

When using extended = TRUE, there are a couple of scenarios in which series morphology contained in SPC do not fully match records in the associated series summary tables (e.g. competing).

1. A query for soil series that exist entirely outside of CONUS (e.g. PALAU).

- Climate summaries are empty data.frames because these summaries are currently generated from PRISM. We are working on a solution that uses DAYMET.

2. A query for data within CONUS, but OSD morphology missing due to parsing error (e.g. formatting, typos).

- Extended summaries are present but morphology missing from SPC. A warning is issued.

These last two cases are problematic for analysis that makes use of morphology and extended data, such as outlined in this tutorial on competing soil series.

Value

a SoilProfileCollection object containing basic soil morphology and taxonomic information.

Note

Requests to the SoilWeb API are split into batches of 100 series names from soils via makeChunks().

Author(s)

D.E. Beaudette, A.G. Brown

References

USDA-NRCS OSD search tools: https://soilseries.sc.egov.usda.gov/

See Also

OSDquery(), siblings()

Examples

library(aqp)
  # soils of interest
  s.list <- c('musick', 'cecil', 'drummer', 'amador', 'pentz',
              'reiff', 'san joaquin', 'montpellier', 'grangeville', 'pollasky', 'ramona')

  # fetch and convert data into an SPC
  s.moist <- fetchOSD(s.list, colorState='moist')
  s.dry <- fetchOSD(s.list, colorState='dry')

  # plot profiles
  # moist soil colors
  par(mar=c(0,0,0,0), mfrow=c(2,1))
  plot(s.moist, name='hzname',
       cex.names=0.85, axis.line.offset=-4)
  plot(s.dry, name='hzname',
       cex.names=0.85, axis.line.offset=-4)

  # extended mode: return a list with SPC + summary tables
  x <- fetchOSD(s.list, extended = TRUE, colorState = 'dry')

  par(mar=c(0,0,1,1))
  plot(x$SPC)
  str(x, 1)

Description

Fetch commonly used site/horizon data from a version 5.x PedonPC database, return as a SoilProfileCollection object.

Usage

fetchPedonPC(dsn)

getHzErrorsPedonPC(dsn, strict = TRUE)

Arguments

Value

a SoilProfileCollection class object

Note

This function attempts to do most of the boilerplate work when extracting site/horizon data from a PedonPC or local NASIS database. Pedons that have errors in their horizonation are excluded from the returned object, however, their IDs are printed on the console. See getHzErrorsPedonPC for a simple approach to identifying pedons with problematic horizonation. Records from the 'taxhistory' table are selected based on 1) most recent record, or 2) record with the least amount of missing data.

Author(s)

D. E. Beaudette and J. M. Skovlin

See Also

get_hz_data_from_pedon_db

Description

NOTICE: The SoilWeb snapshot of the RaCA data has been deprecated. The latest version of the data, including values measured by the Kellogg Soil Survey Laboratory, and supporting documentation, are available here: https://www.nrcs.usda.gov/resources/data-and-reports/rapid-carbon-assessment-raca. Download link on National Agricultural Library Ag Data Commons: https://data.nal.usda.gov/dataset/rapid-carbon-assessment-raca

Get Rapid Carbon Assessment (RaCA) data by state, geographic bounding-box, RaCA site ID, or soil series query from the SoilWeb API. This interface to the data was an experimental delivery service that does not include the latest soil organic carbon (SOC) measurements.

Please use current RaCA distribution if you need lab measured SOC rather than SOC estimated by VNIR.

This interface will be updated sometime calendar year 2022 to include the latest soil morphology, taxonomic classification, and measured SOC values. More detailed coordinates for sample sites should also be available.

Usage

fetchRaCA(
  series = NULL,
  bbox = NULL,
  state = NULL,
  rcasiteid = NULL,
  get.vnir = FALSE
)

Arguments

Details

The VNIR spectra associated with RaCA data are quite large (each gzip-compressed VNIR spectra record is about 6.6kb), so requests for these data are disabled by default. Note that VNIR spectra can only be queried by soil series or geographic BBOX.

Value

pedons:

a SoilProfileCollection object containing site/pedon/horizon data

trees:

a data.frame object containing tree DBH and height

veg:

a data.frame object containing plant species

stock:

a data.frame object containing carbon quantities (stocks) at standardized depths

sample:

a data.frame object containing sample-level bulk density and soil organic carbon values

spectra:

a numeric matrix containing VNIR reflectance spectra from 350–2500 nm

Author(s)

D.E. Beaudette, USDA-NRCS staff

References

https://data.nal.usda.gov/dataset/rapid-carbon-assessment-raca

See Also

fetchOSD

Description

Query soil/climate data from USDA-NRCS SCAN Stations.

Usage

fetchSCAN(
  site.code = NULL,
  year = NULL,
  report = "SCAN",
  timeseries = c("Daily", "Hourly"),
  ...
)

SCAN_sensor_metadata(site.code)

SCAN_site_metadata(site.code = NULL)

Arguments

Details

Possible above and below ground sensor types include: 'SMS' (soil moisture), 'STO' (soil temperature), 'SAL' (salinity), 'TAVG' (daily average air temperature), 'TMIN' (daily minimum air temperature), 'TMAX' (daily maximum air temperature), 'PRCP' (daily precipitation), 'PREC' (daily precipitation), 'SNWD' (snow depth), 'WTEQ' (snow water equivalent),'WDIRV' (wind direction), 'WSPDV' (wind speed), 'LRADT' (solar radiation/langley total).

This function converts below-ground sensor depth from inches to cm. All temperature values are reported as degrees C. Precipitation, snow depth, and snow water content are reported as inches.

SCAN Sensors

All Soil Climate Analysis Network (SCAN) sensor measurements are reported hourly.

SNOTEL Sensors

All Snow Telemetry (SNOTEL) sensor measurements are reported daily.

See the fetchSCAN tutorial for additional usage and visualization examples.

Value

a list of data.frame objects, where each element name is a sensor type, plus a metadata table; different report types change the types of sensor data returned. SCAN_sensor_metadata() and SCAN_site_metadata() return a data.frame. NULL on bad request.

Author(s)

D.E. Beaudette, A.G. Brown

References

See the Soil Climate Analysis Network home page for more information on the SCAN program, and links to other associated programs such as SNOTEL, at the National Weather and Climate Center. You can get information on available web services, as well as interactive maps of snow water equivalent, precipitation and streamflow.

Examples

## Not run: 
    # get data
    x <- try(fetchSCAN(site.code=c(356, 2072), year=c(2015, 2016)))
    str(x)

    # get sensor metadata
    m <- SCAN_sensor_metadata(site.code=c(356, 2072))

    # get site metadata
    m <- SCAN_site_metadata(site.code=c(356, 2072))

    # get hourly data (396315 records)
    # x <- try(fetchSCAN(site.code=c(356, 2072), year=c(2015, 2016), timeseries = "Hourly"))

## End(Not run)

Description

This method facilitates queries to Soil Data Access (SDA) mapunit and survey area geometry. Queries are generated based on map unit key (mukey) and national map unit symbol (nationalmusym) for mupolygon (SSURGO) or gsmmupolygon (STATSGO) geometry OR legend key (lkey) and area symbols (areasymbol) for sapolygon (Soil Survey Area; SSA) geometry).

A Soil Data Access query returns geometry and key identifying information about the map unit or area of interest. Additional columns from the map unit or legend table can be included; see add.fields argument.

Usage

fetchSDA_spatial(
  x,
  by.col = "mukey",
  method = "feature",
  geom.src = "mupolygon",
  db = "SSURGO",
  add.fields = NULL,
  chunk.size = 10,
  verbose = TRUE,
  as_Spatial = getOption("soilDB.return_Spatial", default = FALSE)
)

Arguments

Details

This function automatically "chunks" the input vector (using makeChunks()) of map unit identifiers to minimize the likelihood of exceeding the SDA data request size. The number of chunks varies with the chunk.size setting and the length of your input vector. If you are working with many map units and/or large extents, you may need to decrease this number in order to have more chunks.

Querying regions with complex mapping may require smaller chunk.size. Numerically adjacent IDs in the input vector may share common qualities (say, all from same soil survey area or region) which could cause specific chunks to perform "poorly" (slow or error) no matter what the chunk size is. Shuffling the order of the inputs using sample() may help to eliminate problems related to this, depending on how you obtained your set of MUKEY/nationalmusym to query. One could feasibly use muacres as a heuristic to adjust for total acreage within chunks.

Note that STATSGO data are fetched where CLIPAREASYMBOL = 'US' to avoid duplicating state and national subsets of the geometry.

A prototype interface, geom.src="mlrapolygon", is provided for obtaining Major Land Resource Area (MLRA) polygon boundaries. When using this geometry source x is a vector of MLRARSYM (MLRA Symbols). The geometry source is the MLRA Geographic Database v5.2 (2022) which is not (yet) part of Soil Data Access. Instead of SDA, GDAL utilities are used to read a zipped ESRI Shapefile from a remote URL: https://www.nrcs.usda.gov/sites/default/files/2022-10/MLRA_52_2022.zip. Therefore, most additional fetchSDA_spatial() arguments are not currently supported for the MLRA geometry source. In the future a mlrapolygon table may be added to SDA (analogous to mupolygon and sapolygon), and the function will be updated accordingly at that time.

Value

an sf data.frame corresponding to SDA spatial data for all symbols requested. If as_Spatial=TRUE returns a Spatial*DataFrame from the sp package via sf::as_Spatial() for backward compatibility. Default result contains geometry with attribute table containing unique feature ID, symbol and area symbol plus additional fields in result specified with add.fields.

Author(s)

Andrew G. Brown, Dylan E. Beaudette

Examples

# get spatial data for a single mukey
    single.mukey <- try(fetchSDA_spatial(x = "2924882"))

    # demonstrate fetching full extent (multi-mukey) of national musym
    full.extent.nmusym <- try(fetchSDA_spatial(x = "2x8l5", by = "nmusym"))

    # compare extent of nmusym to single mukey within it
    if (!inherits(single.mukey, 'try-error') && 
        !inherits(full.extent.nmusym, 'try-error')) {
        
        if (requireNamespace("sf")) {
      
         plot(sf::st_geometry(full.extent.nmusym), col = "RED", border = 0)
         plot(sf::st_geometry(single.mukey), add = TRUE, col = "BLUE", border = 0)
       
        }
        
    }

    # demo adding a field (`muname`) to attribute table of result
    head(try(fetchSDA_spatial(x = "2x8l5", by="nmusym", add.fields="muname")))

Description

This function obtains SoilGrids 2.0 properties information (250m raster resolution) given a data.frame containing site IDs, latitudes and longitudes, or a spatial extent (see grid=TRUE argument).

Usage

fetchSoilGrids(
  x,
  loc.names = c("id", "lat", "lon"),
  depth_intervals = c("0-5", "5-15", "15-30", "30-60", "60-100", "100-200"),
  variables = c("bdod", "cec", "cfvo", "clay", "nitrogen", "phh2o", "sand", "silt",
    "soc", "ocd", "wv0010", "wv0033", "wv1500"),
  grid = FALSE,
  filename = NULL,
  overwrite = TRUE,
  target_resolution = c(250, 250),
  summary_type = c("Q0.05", "Q0.5", "Q0.95", "mean"),
  endpoint = ifelse(!grid, "https://rest.isric.org/soilgrids/v2.0/properties/query",
    "https://files.isric.org/soilgrids/latest/data/"),
  ...,
  verbose = FALSE,
  progress = FALSE
)

Arguments

Details

SoilGrids API and maps return values as whole (integer) numbers to minimize the storage space used. These values have conversion factors applied by fetchSoilGrids() to produce conventional units shown in the table below (see Details).

Properties

SoilGrids predictions are made for the six standard depth intervals specified in the GlobalSoilMap IUSS working group and its specifications. The default depth intervals returned are (in centimeters): "0-5", "5-15", "15-30", "30-60", "60-100", "100-200" for the properties "bdod", "cec", "cfvo", "clay", "nitrogen", "phh2o", "sand", "silt", "soc", "ocd", "wv0010", "wv0033", "wv1500"–each with 5th, 50th, 95th, mean and uncertainty values. Soil organic carbon stocks (0-30cm) (variables="ocs") are returned only for depth_intervals="0-30". The uncertainty values are the ratio between the inter-quantile range (90% prediction interval width) and the median : (Q0.95-Q0.05)/Q0.50. All values are converted from "mapped" to "conventional" based on above table conversion factors. Point data requests are made through "properties/query" endpoint of the SoilGrids v2.0 REST API. Please check ISRIC's data policy, disclaimer and citation: https://www.isric.org/about/data-policy.

Find out more information about the SoilGrids and GlobalSoilMap products here:

• https://www.isric.org/explore/soilgrids/faq-soilgrids

• https://www.isric.org/sites/default/files/GlobalSoilMap_specifications_december_2015_2.pdf

Value

A SoilProfileCollection or SpatRaster when grid=TRUE. Returns try-error if all requests fail. Any error messages resulting from parsing will be echoed when verbose=TRUE.

Author(s)

Andrew G. Brown

References

• Common soil chemical and physical properties: Poggio, L., de Sousa, L. M., Batjes, N. H., Heuvelink, G. B. M., Kempen, B., Ribeiro, E., and Rossiter, D.: SoilGrids 2.0: producing soil information for the globe with quantified spatial uncertainty, SOIL, 7, 217–240, 2021. DOI: doi:10.5194/soil-7-217-2021

• Soil water content at different pressure heads: Turek, M.E., Poggio, L., Batjes, N. H., Armindo, R. A., de Jong van Lier, Q., de Sousa, L.M., Heuvelink, G. B. M. : Global mapping of volumetric water retention at 100, 330 and 15000 cm suction using the WoSIS database, International Soil and Water Conservation Research, 11-2, 225-239, 2023. DOI: doi:10.1016/j.iswcr.2022.08.001

Examples

## Not run: 
  library(aqp)
  
  your.points <- data.frame(id  = c("A", "B"), 
                           lat = c(37.9, 38.1), 
                           lon = c(-120.3, -121.5), 
                           stringsAsFactors = FALSE)
  x <- try(fetchSoilGrids(your.points))
 
  if (!inherits(x, 'try-error'))
   aqp::plotSPC(x, name = NA, color = "socQ50")
 
  # organic carbon stocks use 0-30cm interval
  y <- try(fetchSoilGrids(your.points[1, ], 
                          depth_interval = c("0-5", "0-30", "5-15", "15-30"),
                          variables = c("soc", "bdod", "ocd", "ocs")))
                          
  # extract horizons from a SoilProfileCollection where horizon 2 overlaps 1, 3, and 4
  h <- aqp::horizons(y)
  
  # "ocs" (organic carbon stock 0-30cm interval)
  h[2, ]
  
  h$thickness_meters <- ((h$hzdepb - h$hzdept) / 100)

  # estimate "ocs" from modeled organic carbon and bulk density in 0-5, 5-15, 15-30 intervals
  #  (sum the product of soc, bdod, and thickness in meters)
  #  (1 gram per cubic decimeter = 1 kilogram per cubic meter)
  sum(h$socmean * h$bdodmean * h$thickness_meters, na.rm = TRUE)
  
  # estimate "ocs" from modeled organic carbon density in 0-5, 5-15, 15-30 intervals
  #  (sum the product of "ocd" and thickness in meters)
  sum(h$ocdmean * h$thickness_meters, na.rm = TRUE)
 

## End(Not run)

Description

This is a higher level wrapper around the get_SRI and get_SRI_layers functions. This function can fetch multiple File Geodatabases (GDB) and returns all the layers within the GDB.

Usage

fetchSRI(gdb, ...)

Arguments

Value

A list.

Author(s)

Josh Erickson

See Also

get_SRI() get_SRI_layers()

Examples

## Not run: 

# fetch Willamette and Winema SRI

sri <- fetchSRI(gdb = c('will', 'win'), quiet = TRUE)


## End(Not run)

Description

Get vegetation plot data from local NASIS database

Usage

fetchVegdata(SS = TRUE, stringsAsFactors = NULL, dsn = NULL)

get_vegplot_from_NASIS_db(SS = TRUE, stringsAsFactors = NULL, dsn = NULL)

get_vegplot_location_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_trhi_from_NASIS_db(SS = TRUE, stringsAsFactors = NULL, dsn = NULL)

get_vegplot_species_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_transect_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_transpecies_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_transpoints_from_NASIS_db(SS = TRUE, dsn = NULL)

get_vegplot_prodquadrats_from_NASIS_db(SS = TRUE, dsn = NULL)

get_vegplot_tree_si_summary_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_speciesbasalarea_from_NASIS(SS = TRUE, dsn = NULL)

get_vegplot_tree_si_details_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_vegplot_textnote_from_NASIS_db(
  SS = TRUE,
  fixLineEndings = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

Arguments

Value

A named list containing: "vegplot", "vegplotlocation", "vegplotrhi", "vegplotspecies", "vegtransect", "vegtransplantsum", 'vegsiteindexsum', "vegsiteindexdet", and "vegplottext" tables

Description

A function to subset KSSL "geochem" / elemental analysis result table to obtain rows/columns based on: column name, preparation code, major / trace element method.

Usage

filter_geochem(
  geochem,
  columns = NULL,
  prep_code = NULL,
  major_element_method = NULL,
  trace_element_method = NULL
)

Arguments

Value

A data.frame, subset according to the constraints specified in arguments.

Author(s)

Andrew G. Brown.

Description

Concatenate a vector to SQL IN-compatible syntax: letters[1:3] becomes ('a','b','c'). Values in x are first passed through unique().

Usage

format_SQL_in_statement(x)

Arguments

Value

A character vector (unit length) containing concatenated group syntax for use in SQL IN, with unique value found in x.

Note

Only character output is supported.

Examples

format_SQL_in_statement(c(2648889L, 2648890L))

Description

Get, format, mix, and return color data from a NASIS database.

Usage

get_colors_from_NASIS_db(SS = TRUE, mixColors = TRUE, dsn = NULL)

Arguments

Value

A data.frame with the results.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

simplifyColorData, get_hz_data_from_NASIS_db, get_site_data_from_NASIS_db

Description

Get, format, mix, and return color data from a PedonPC database.

Usage

get_colors_from_pedon_db(dsn)

Arguments

Value

A data.frame with the results.

Author(s)

Dylan E. Beaudette and Jay M. Skovlin

See Also

get_hz_data_from_pedon_db, get_site_data_from_pedon_db

Description

Get component month data from a local NASIS Database.

Usage

get_comonth_from_NASIS_db(
  SS = TRUE,
  fill = FALSE,
  stringsAsFactors = NULL,
  dsn = NULL
)

Arguments

Value

A list with the results.

Author(s)

Stephen Roecker

See Also

fetchNASIS

Examples

if(local_NASIS_defined()) {
  # query text note data
  cm <- try(get_comonth_from_NASIS_db())

  # show structure of component month data
  str(cm)
}

Description

Get component data from a local NASIS Database

Usage

get_component_data_from_NASIS_db(
  SS = TRUE,
  nullFragsAreZero = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_component_diaghz_from_NASIS_db(SS = TRUE, dsn = NULL)

get_component_restrictions_from_NASIS_db(SS = TRUE, dsn = NULL)

get_component_correlation_data_from_NASIS_db(
  SS = TRUE,
  dropAdditional = TRUE,
  dropNotRepresentative = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_component_cogeomorph_data_from_NASIS_db(SS = TRUE, dsn = NULL)

get_component_cogeomorph_data_from_NASIS_db2(SS = TRUE, dsn = NULL)

get_component_copm_data_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_component_esd_data_from_NASIS_db(
  SS = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

get_component_otherveg_data_from_NASIS_db(SS = TRUE, dsn = NULL)

get_copedon_from_NASIS_db(SS = TRUE, dsn = NULL)

get_component_horizon_data_from_NASIS_db(
  SS = TRUE,
  fill = FALSE,
  dsn = NULL,
  nullFragsAreZero = TRUE
)

Arguments

Value

a data.frame

Author(s)

Dylan E. Beaudette, Stephen Roecker, and Jay M. Skovlin

See Also

fetchNASIS

Examples

if(local_NASIS_defined()) {
 # query text note data
 fc <- try(get_component_data_from_NASIS_db())

 # show structure of component data returned
 str(fc)
}

Description

Functions to load and flatten commonly used tables and from SSURGO file geodatabases, and create soil profile collection objects (SPC).

Usage

get_component_from_GDB(
  dsn = "gNATSGO_CONUS.gdb",
  WHERE = NULL,
  childs = FALSE,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_legend_from_GDB(
  dsn = "gNATSGO_CONUS.gdb",
  WHERE = NULL,
  droplevels = TRUE,
  stringsAsFactors = NULL,
  stats = FALSE
)

get_mapunit_from_GDB(
  dsn = "gNATSGO_CONUS.gdb",
  WHERE = NULL,
  droplevels = TRUE,
  stringsAsFactors = NULL,
  stats = FALSE
)

fetchGDB(
  dsn = "gNATSGO_CONUS.gdb",
  WHERE = NULL,
  childs = FALSE,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

Arguments

Details

These functions return data from SSURGO file geodatabases with the use of a simple text string that formatted as an SQL WHERE clause (e.g. WHERE = "areasymbol = 'IN001'". Any columns within the target table can be specified (except for fetchGDB() which currently can only target one table (e.g. legend, mapunit or component) at a time with the WHERE clause).

Value

A data.frame or SoilProfileCollection object.

Author(s)

Stephen Roecker

Examples

## replace `dsn` with path to your own geodatabase (SSURGO OR gNATSGO)
##
##
##  download CONUS gNATSGO from here:
##    https://nrcs.app.box.com/v/soils/folder/191790828371
##
# dsn <- "D:/geodata/soils/gNATSGO_CONUS.gdb"

# le <- get_legend_from_GDB(dsn = dsn, WHERE = "areasymbol LIKE '%'")

# mu <- get_mapunit_from_GDB(dsn = dsn, WHERE = "muname LIKE 'Miami%'")

# co <- get_component_from_GDB(dsn, WHERE = "compname = 'Miami'
#                              AND majcompflag = 'Yes'", childs = FALSE)

# f_in_GDB <- fetchGDB(WHERE = "areasymbol LIKE 'IN%'")

Description

Functions to download and flatten commonly used tables and from Soil Data Access, and create soil profile collection objects (SPC).

Usage

get_component_from_SDA(
  WHERE = NULL,
  duplicates = FALSE,
  childs = TRUE,
  droplevels = TRUE,
  nullFragsAreZero = TRUE,
  stringsAsFactors = NULL
)

get_cointerp_from_SDA(
  WHERE = NULL,
  mrulename = NULL,
  duplicates = FALSE,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_legend_from_SDA(WHERE = NULL, droplevels = TRUE, stringsAsFactors = NULL)

get_lmuaoverlap_from_SDA(
  WHERE = NULL,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_mapunit_from_SDA(WHERE = NULL, droplevels = TRUE, stringsAsFactors = NULL)

get_chorizon_from_SDA(
  WHERE = NULL,
  duplicates = FALSE,
  childs = TRUE,
  nullFragsAreZero = TRUE,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

fetchSDA(
  WHERE = NULL,
  duplicates = FALSE,
  childs = TRUE,
  nullFragsAreZero = TRUE,
  rmHzErrors = FALSE,
  droplevels = TRUE,
  stringsAsFactors = NULL
)

get_cosoilmoist_from_SDA(
  WHERE = NULL,
  duplicates = FALSE,
  impute = TRUE,
  stringsAsFactors = NULL
)

Arguments

Details

These functions return data from Soil Data Access with the use of a simple text string that formatted as an SQL WHERE clause (e.g. WHERE = "areasymbol = 'IN001'". All functions are SQL queries that wrap around SDAquery() and format the data for analysis.

Beware SDA includes the data for both SSURGO and STATSGO2. The areasymbol for STATSGO2 is US. For just SSURGO, include WHERE = "areareasymbol != 'US'".

If the duplicates argument is set to TRUE, duplicate components are returned. This is not necessary with data returned from NASIS, which has one unique national map unit. SDA has duplicate map national map units, one for each legend it exists in.

The value of nullFragsAreZero will have a significant impact on the rock fragment fractions returned by fetchSDA. Set nullFragsAreZero = FALSE in those cases where there are many data-gaps and NULL rock fragment values should be interpreted as NULLs. Set nullFragsAreZero = TRUE in those cases where NULL rock fragment values should be interpreted as 0.

Additional examples can be found in the Soil Data Access (SDA) Tutorial

Value

A data.frame or SoilProfileCollection object.

Author(s)

Stephen Roecker

See Also

SDA_query

Description

Read and flatten the component soil moisture month tables from a local NASIS Database.

Usage

get_cosoilmoist_from_NASIS(
  SS = TRUE,
  impute = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

Arguments

Details

The component soil moisture tables within NASIS house monthly data on flooding, ponding, and soil moisture status. The soil moisture status is used to specify the water table depth for components (e.g. status == "Moist").

Value

A data.frame.

Author(s)

S.M. Roecker

See Also

fetchNASIS, get_cosoilmoist_from_NASISWebReport, get_cosoilmoist_from_SDA, get_comonth_from_SDA

Examples

if(local_NASIS_defined()) {
 # load cosoilmoist (e.g. water table data)
 test <- try(get_cosoilmoist_from_NASIS())

 # inspect
 if(!inherits(test, 'try-error')) {
   head(test)
 }
}

Description

Gets the Site Ecological Site History data from local NASIS database. Used by get_extended_data_from_NASIS_db().

Usage

get_ecosite_history_from_NASIS_db(
  best = TRUE,
  SS = TRUE,
  es_classifier = NULL,
  dsn = NULL
)

Arguments

Value

a data.frame, or NULL on error

See Also

get_extended_data_from_NASIS_db()

Description

Data are accessed via Ecological Dynamics Interpretive Tool (EDIT) web services: https://edit.jornada.nmsu.edu/resources/esd. geoUnit refers to MLRA codes, possibly with a leading zero and trailing "X" for two digit MLRA symbols.

Usage

get_EDIT_ecoclass_by_geoUnit(geoUnit, catalog = c("esd", "esg"))

Arguments

Value

A data.frame containing: geoUnit, id, legacyId, name. NULL if no result.

Examples

## Not run: 
   get_EDIT_ecoclass_by_geoUnit(c("018X","022A"))

## End(Not run)

Description

Get accessory tables and summaries from a local NASIS Database

Usage

get_extended_data_from_NASIS_db(
  SS = TRUE,
  nullFragsAreZero = TRUE,
  stringsAsFactors = NULL,
  dsn = NULL
)

Arguments

Value

A list with the results.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

get_hz_data_from_NASIS_db, get_site_data_from_NASIS_db

Examples

if(local_NASIS_defined()) {
 # query extended data
 e <- try(get_extended_data_from_NASIS_db())

 # show contents of extended data
 str(e)
}

Description

Get accessory tables and summaries from a local pedonPC Database.

Usage

get_extended_data_from_pedon_db(dsn)

Arguments

Value

A list with the results.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

get_hz_data_from_pedon_db, get_site_data_from_pedon_db

Description

Get horizon-level data from a local NASIS database.

Usage

get_hz_data_from_NASIS_db(
  SS = TRUE,
  fill = FALSE,
  stringsAsFactors = NULL,
  dsn = NULL
)

Arguments

Value

A data.frame.

Note

NULL total rock fragment values are assumed to represent an absence of rock fragments, and set to 0.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

get_hz_data_from_NASIS_db, get_site_data_from_NASIS_db

Description

Get horizon-level data from a PedonPC database.

Usage

get_hz_data_from_pedon_db(dsn)

Arguments

Value

A data.frame.

Note

NULL total rock fragment values are assumed to represent an absence of rock fragments, and set to 0.

Author(s)

Dylan E. Beaudette and Jay M. Skovlin

See Also

get_colors_from_pedon_db, get_site_data_from_pedon_db

Description

Get lab pedon layer-level (horizon-level) data from a local NASIS database.

Usage

get_lablayer_data_from_NASIS_db(SS = TRUE, dsn = NULL)

Arguments

Value

A data.frame.

Note

This function queries KSSL laboratory site/horizon data from a local NASIS database from the lab layer data table.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

get_labpedon_data_from_NASIS_db

Description

Get lab pedon-level data from a local NASIS database.

Usage

get_labpedon_data_from_NASIS_db(SS = TRUE, dsn = NULL)

Arguments

Details

This function currently works only on Windows, and requires a 'nasis_local' ODBC connection.

Value

A data.frame.

Note

This function queries KSSL laboratory site/horizon data from a local NASIS database from the lab pedon data table.

Author(s)

Jay M. Skovlin and Dylan E. Beaudette

See Also

get_lablayer_data_from_NASIS_db

Description

Get Legend, Mapunit and Legend Mapunit Area Overlap Tables

Usage

get_mapunit_from_NASIS(
  SS = TRUE,
  repdmu = TRUE,
  droplevels = TRUE,
  stringsAsFactors = NULL,
  areatypename = c("Non-MLRA Soil Survey Area", "MLRA Soil Survey Area"),
  dsn = NULL
)

get_legend_from_NASIS(
  SS = TRUE,
  droplevels = TRUE,
  stringsAsFactors = NULL,
  areatypename = c("Non-MLRA Soil Survey Area", "MLRA Soil Survey Area"),
  dsn = NULL
)

get_lmuaoverlap_from_NASIS(
  SS = TRUE,
  droplevels = TRUE,
  stringsAsFactors = NULL,
  areatypename = c("Non-MLRA Soil Survey Area", "MLRA Soil Survey Area"),
  dsn = NULL
)

get_projectmapunit_from_NASIS(SS = TRUE, stringsAsFactors = NULL, dsn = NULL)

Arguments

Description

Retrieve a table containing domain and column names with choice list labels/names/sequences/values from the NASIS 7 metadata tables.

Usage

get_NASIS_metadata(dsn = NULL, include_description = FALSE)

get_NASIS_column_metadata(
  x,
  what = "ColumnPhysicalName",
  include_description = FALSE,
  dsn = NULL
)

Arguments

Details

These data are derived from the MetadataDomainDetail, MetadataDomainMaster, and MetadataTableColumn tables and help with mapping between values stored in the NASIS database and human-readable values. The human-readable values align with the values returned in public facing interfaces such as SSURGO via Soil Data Access and NASIS Web Reports. The data in these tables can also be used to create ordered factors where options for levels of a particular data element follow a logical ChoiceSequence.

If a local NASIS instance is set up, and this is the first time get_NASIS_metadata() has been called, the metadata will be obtained from the NASIS local database. Subsequent runs in the same session will use a copy of the data object NASIS.metadata cached in soilDB.env which can be accessed with get_soilDB_env()$NASIS.metadata.

For users without a local NASIS instance, a cached copy of the NASIS metadata are used (data/metadata.rda).

See ?soilDB::metadata for additional details.

Value

a data.frame containing DomainID, DomainName, DomainRanked, DisplayLabel, ChoiceSequence, ChoiceValue, ChoiceName, ChoiceLabel, ChoiceObsolete, ColumnPhysicalName, ColumnLogicalName and optionally ChoiceDescription when include_description=TRUE.

a data.frame containing selected NASIS metadata sorted first on DomainID and then on ChoiceSequence

Examples

get_NASIS_metadata()
get_NASIS_column_metadata("texcl")

Description

Get a NASIS table key by type and table name

Usage

get_NASIS_table_key_by_name(
  tables,
  keycol = c("all", "fkey", "pkeyref", "pkey")
)

Arguments

Value

The key column name for the specified table name

Examples

## Not run: 
get_NASIS_table_key_by_name(c("site","phorizon_View_1","not_a_table"))

## End(Not run)#'

Description

Retrieve a table containing table and column names with descriptions, help text, units of measure, etc. from NASIS 7 metadata tables.

Usage

get_NASIS_table_metadata(
  table = NULL,
  column = NULL,
  what.table = "TablePhysicalName",
  what.column = "ColumnPhysicalName",
  query_string = FALSE,
  dsn = NULL
)

Arguments

Details

These data are derived from the MetadataTable and MetadataTableColumn tables and describe the expected contents of standard NASIS tables and columns.

For NASIS choice lists based on domain and column names see get_NASIS_metadata() and NASISChoiceList(). This function (get_NASIS_table_metadata()) is intended for higher-level description of the expected contents of a NASIS database instance, rather than the codes/specific values used within columns.

Value

a data.frame

See Also

get_NASIS_metadata() NASISChoiceList() uncode() code()

Examples

if (local_NASIS_defined())
 str(get_NASIS_table_metadata())

Description

Method generalizing concepts of NASIS 7 data model to group tables by "purpose." Most of our more complex queries rely on tables from one or more purposes, so individual higher-level functions might call a function like this to identify the relevant tables from a data source.

Usage

get_NASIS_table_name_by_purpose(
  purpose = c("metadata", "lookup", "nasis", "site", "pedon", "transect", "component",
    "vegetation", "project", "techsoilservice", "area", "soilseries", "legend",
    "mapunit", "datamapunit"),
  SS = FALSE
)

Arguments

Value

character vector of table names

See Also

createStaticNASIS

Examples

## Not run: 
# get the "site" table names
get_NASIS_table_name_by_purpose("site")

# get the pedon table names
get_NASIS_table_name_by_purpose("pedon", SS = TRUE)

# metadata and lookup not affected by SS argument, but site and pedon are
get_NASIS_table_name_by_purpose(c("metadata", "lookup",
                                   "site", "pedon"), SS = TRUE)

## End(Not run)

Description

Obtain daily climatic summary data for a set of station IDs, years, and datatypes.

Note that typically results from the NOAA API are limited to 1000 records. However, by "chunking" up data into individual stationyeardatatypeid combinations, record results generally do not exceed 365 records for daily summaries.

In order to use this function, you must obtain an API token from this website: https://www.ncdc.noaa.gov/cdo-web/token

Usage

get_NOAA_GHCND(stations, years, datatypeids, apitoken)

Arguments

Value

A data.frame containing the GHCND data requested (limit 1000 records)

Examples

#' ## in order to use this function, you must obtain an API token from this website:
##  https://www.ncdc.noaa.gov/cdo-web/token

# get_NOAA_GHCND(c("GHCND:USC00388786", "GHCND:USC00388787"),
#                years = 2017:2020,
#                datatypeids = c("PRCP","SNOW"),
#                apitoken = "yourtokenhere")

Description

Query the NOAA API to get station data (limit 1000 records) near a point. Default extent is plus or minus 0.5 degrees (bounding box) (with bbox = 1) around the specified point [lat, lng].

In order to use this function, you must obtain an API token from this website: https://www.ncdc.noaa.gov/cdo-web/token

Usage

get_NOAA_stations_nearXY(lat, lng, apitoken, bbox = 1, crs = "EPSG:4326")

Arguments

Value

data.frame containing station information for all stations within a bounding box around lat, lng.

Examples

## in order to use this function, you must obtain an API token from this website:
##  https://www.ncdc.noaa.gov/cdo-web/token

# stations <- get_NOAA_stations_nearXY(lat = 37, lng = -120,
#                                      apitoken = "yourtokenhere")

Description

Get Official Series Description Data from JSON, HTML or TXT sources

Usage

get_OSD(
  series,
  base_url = NULL,
  result = c("json", "html", "txt"),
  fix_ocr_errors = FALSE,
  verbose = FALSE
)

get_OSD_JSON(series, base_url = NULL)

Arguments

Details

The default base_url for result="json" is to JSON files stored in a GitHub repository that is regularly updated from the official source of Series Descriptions. Using format: ⁠https://raw.githubusercontent.com/ncss-tech/SoilKnowledgeBase/main/inst/extdata/OSD/{LETTER}/{SERIES}.json⁠ for JSON. And ⁠"https://soilseriesdesc.sc.egov.usda.gov/OSD_Docs/{LETTER}/{SERIES}.html⁠ is for result="html" (official source).

fix_ocr_errors by default is turned off (FALSE). When TRUE, assume that in color data hue/value/chroma lowercase "L" ("l") is a 1, and a capital "O" is interpreted as zero. Also, in horizon designations assume lowercase "L" is a 1, and a string that starts with 0 starts with the capital letter "O".

Value

For JSON result: A data.frame with 1 row per series, and 1 column per "section" in the OSD as defined in National Soil Survey Handbook. For TXT or HTML result a list of character vectors containing OSD text with 1 element per series and one value per line.

Examples

series <- c("Musick", "Hector", "Chewacla")
  get_OSD(series)

Description

Prepare a list of data.frame objects with data from the "phrdxfeatures" and "phredoxfcolor" tables. These tables are related by "phrdxfiid" column, and related to horizon data via "phiid".

Usage

get_RMF_from_NASIS_db(SS = TRUE, dsn = NULL)

Arguments

Value

a list with two data.frame objects:

• RMF: contents of "phrdxfeatures" table, often >1 row per horizon

• RMF_colors: contents of "phredoxfcolor", usually >1 row per record in "phrdxfeatures"

Description

Get mapunit ecological sites from Soil Data Access

Usage

get_SDA_coecoclass(
  method = "None",
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  query_string = FALSE,
  ecoclasstypename = c("NRCS Rangeland Site", "NRCS Forestland Site"),
  ecoclassref = "Ecological Site Description Database",
  not_rated_value = "Not assigned",
  miscellaneous_areas = TRUE,
  include_minors = TRUE,
  threshold = 0,
  dsn = NULL
)

Arguments

Details

When method="Dominant Condition" an additional field ecoclasspct_r is returned in the result with the sum of comppct_r that have the dominant condition ecoclassid. The component with the greatest comppct_r is returned for the component and coecosite level information.

Note that if there are multiple coecoclasskey per ecoclassid there may be more than one record per component.

Description

Get Geomorphic/Surface Morphometry Data from Soil Data Access or a local SSURGO data source and summarize by counts and proportions ("probabilities").

Usage

get_SDA_cosurfmorph(
  table = c("cosurfmorphgc", "cosurfmorphhpp", "cosurfmorphss", "cosurfmorphmr"),
  by = "compname",
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  miscellaneous_areas = FALSE,
  db = c("SSURGO", "STATSGO"),
  dsn = NULL,
  query_string = FALSE
)

Arguments

Details

Default table="cosurfmorphgc" summarizes columns geomposmntn, geomposhill, geomposflats, and geompostrce. table="cosurfmorphhpp" summarizes "hillslopeprof", table="cosurfmorphss" summarizes shapeacross and shapedown, and table="cosurfmorphmr" summarizes geomicrorelief.

Queries are a generalization of now-deprecated functions from sharpshootR package by Dylan Beaudette: geomPosMountainProbability(), geomPosHillProbability(), surfaceShapeProbability(), hillslopeProbability()

Similar summaries of SSURGO component surface morphometry data by series name can be found in fetchOSD(, extended=TRUE) or downloaded from https://github.com/ncss-tech/SoilWeb-data Full component data including surface morphometry summaries at the "site" level can be obtained with fetchSDA().

Value

a data.frame containing the grouping variable (by) and tabular summaries of counts and proportions of geomorphic records.

Author(s)

Dylan E. Beaudette, Andrew G. Brown

See Also

fetchSDA() get_SDA_pmgroupname()

Examples

## Not run: 
 # Summarize by 3D geomorphic components by component name (default `by='compname'`)
 get_SDA_cosurfmorph(WHERE = "areasymbol = 'CA630'")

 # Whole Soil Survey Area summary (using `by = 'areasymbol'`)
 get_SDA_cosurfmorph(by = 'areasymbol', WHERE = "areasymbol = 'CA630'")

 # 2D Hillslope Position summary(using `table = 'cosurfmorphhpp'`)
 get_SDA_cosurfmorph('cosurfmorphhpp', WHERE = "areasymbol = 'CA630'")

 # Surface Shape summary (using `table = 'cosurfmorphss'`)
 get_SDA_cosurfmorph('cosurfmorphss', WHERE = "areasymbol = 'CA630'")

 # Microrelief summary (using `table = 'cosurfmorphmr'`)
 get_SDA_cosurfmorph('cosurfmorphmr', WHERE = "areasymbol = 'CA630'")

## End(Not run)

Description

Assess the hydric soils composition of a map unit.

Usage

get_SDA_hydric(
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  method = "MAPUNIT",
  query_string = FALSE,
  dsn = NULL
)

Arguments

Details

The default classes for method="MAPUNIT" are as follows:

• 'Nonhydric' - no hydric components

• 'Hydric' - all hydric components

• 'Predominantly Hydric' - hydric component percentage is 50% or more

• 'Partially Hydric' - one or more of the major components is hydric

• 'Predominantly Nonhydric' - hydric component percentage is less than 50%

The default result will also include the following summaries of component percentages: total_comppct, hydric_majors and hydric_inclusions.

Default method "Mapunit" produces aggregate summaries of all components in the mapunit. Use "Dominant Component" and "Dominant Condition" to get the dominant component (highest percentage) or dominant hydric condition (similar conditions aggregated across components), respectively. Use "None" for no aggregation (one record per component).

Value

a data.frame

Author(s)

Jason Nemecek, Chad Ferguson, Andrew Brown

Description

Get map unit interpretations from Soil Data Access by rule name

Usage

get_SDA_interpretation(
  rulename,
  method = c("Dominant Component", "Dominant Condition", "Weighted Average", "None"),
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  query_string = FALSE,
  not_rated_value = NA_real_,
  wide_reason = FALSE,
  dsn = NULL
)

Arguments

Details

Rule Names in cointerp table

• AGR - Avocado Root Rot Hazard (CA)

• AGR - California Revised Storie Index (CA)

• AGR - Hops Site Suitability (WA)

• AGR - Map Unit Cropland Productivity (MN)

• AGR - Nitrate Leaching Potential, Nonirrigated (WA)

• AGR - No Till (TX)

• AGR - Pesticide Loss Potential-Soil Surface Runoff (NE)

• AGR - Ridge Till (TX)

• AGR - Selenium Leaching Potential (CO)

• AGR - Water Erosion Potential (NE)

• AGR - Wind Erosion Potential (TX)

• AGR - Winter Wheat Yield (MT)

• AGR-Pesticide and Nutrient Runoff Potential (ND)

• AGR-Rooting Depth (ND)

• American Wine Grape Varieties Site Desirability (Long)

• American Wine Grape Varieties Site Desirability (Medium)

• American Wine Grape Varieties Site Desirability (Very Long)

• AWM - Animal Mortality Disposal (Catastrophic) (MO)

• AWM - Irrigation Disposal of Wastewater (OH)

• AWM - Irrigation Disposal of Wastewater (VT)s

• AWM - Land Application of Municipal Biosolids, summer (OR)

• AWM - Manure and Food Processing Waste (MD)

• AWM - Manure and Food Processing Waste (OH)

• AWM - Overland Flow Process Treatment of Wastewater (VT)

• AWM - Rapid Infil Disposal of Wastewater (DE)

• AWM - Sensitive Soil Features (MN)

• AWM - Sensitive Soil Features (WI)

• BLM - Fencing

• BLM - Fire Damage Susceptibility

• BLM - Mechanical Treatment, Rolling Drum

• BLM - Rangeland Drill

• BLM - Rangeland Seeding, Colorado Plateau Ecoregion

• BLM - Rangeland Seeding, Great Basin Ecoregion

• BLM-Reclamation Suitability (MT)

• CLASS RULE - Depth to lithic bedrock (5 classes) (NPS)

• CLASS RULE - Soil Inorganic Carbon kg/m2 to 2m (NPS)

• CLASS RULE - Soil Organic Carbon kg/m2 to 2m (NPS)

• CLR-pastureland limitation (IN)

• Commodity Crop Productivity Index (Soybeans) (TN)

• CPI - Alfalfa Hay, NIRR - Palouse, Northern Rocky Mtns. (WA)

• CPI - Barley, IRR - Eastern Idaho Plateaus (ID)

• CPI - Grass Hay, IRR - Klamath Valleys and Basins (OR)

• CPI - Small Grains, IRR - Snake River Plains (ID)

• CPI - Wheat, IRR - Eastern Idaho Plateaus (ID)

• CZSS - Salinization due to Coastal Saltwater Inundation (CT)

• DHS - Catastrophic Event, Large Animal Mortality, Burial

• DHS - Catastrophic Mortality, Large Animal Disposal, Pit

• DHS - Catastrophic Mortality, Large Animal Disposal, Trench

• DHS - Potential for Radioactive Bioaccumulation

• DHS - Potential for Radioactive Sequestration

• DHS - Suitability for Composting Medium and Final Cover

• ENG - Construction Materials; Gravel Source

• ENG - Construction Materials; Gravel Source (AK)

• ENG - Construction Materials; Gravel Source (ID)

• ENG - Construction Materials; Gravel Source (OH)

• ENG - Construction Materials; Gravel Source (VT)

• ENG - Construction Materials; Gravel Source (WA)

• ENG - Construction Materials; Roadfill (OH)

• ENG - Construction Materials; Sand Source (OR)

• ENG - Construction Materials; Sand Source (WA)

• ENG - Construction Materials; Topsoil (GA)

• ENG - Construction Materials; Topsoil (MD)

• ENG - Daily Cover for Landfill

• ENG - Daily Cover for Landfill (AK)

• ENG - Disposal Field Suitability Class (NJ)

• ENG - Dwellings W/O Basements (OH)

• ENG - Dwellings with Basements (AK)

• ENG - Large Animal Disposal, Pit (CT)

• ENG - Lawn, landscape, golf fairway (CT)

• ENG - Lined Retention Systems

• ENG - Local Roads and Streets (OH)

• ENG - On-Site Waste Water Absorption Fields (MO)

• ENG - Septic Tank Absorption Fields

• ENG - Septic Tank Absorption Fields (MD)

• ENG - Septic Tank Absorption Fields (TX)

• ENG - Septic Tank, Gravity Disposal (TX)

• ENG - Sewage Lagoons

• ENG - Small Commercial Buildings (OH)

• ENG - Soil Potential Ratings of SSDS (CT)

• FOR (USFS) - Road Construction/Maintenance (Natural Surface)

• FOR - Compaction Potential (WA)

• FOR - Conservation Tree/Shrub Groups (MT)

• FOR - Damage by Fire (OH)

• FOR - General Harvest Season (VT)

• FOR - Hand Planting Suitability

• FOR - Hand Planting Suitability, MO13 (DE)

• FOR - Hand Planting Suitability, MO13 (MD)

• FOR - Log Landing Suitability

• FOR - Log Landing Suitability (ME)

• FOR - Log Landing Suitability (VT)

• FOR - Log Landing Suitability (WA)

• FOR - Mechanical Planting Suitability (CT)

• FOR - Mechanical Planting Suitability, MO13 (MD)

• FOR - Mechanical Site Preparation (Deep)

• FOR - Mechanical Site Preparation (Deep) (DE)

• FOR - Mechanical Site Preparation (Surface) (DE)

• FOR - Mechanical Site Preparation (Surface) (MI)

• FOR - Mechanical Site Preparation; Surface (ME)

• FOR - Potential Erosion Hazard, Road/Trail, Spring Thaw (AK)

• FOR - Potential Seedling Mortality (PIA)

• FOR - Potential Seedling Mortality(ME)

• FOR - Puddling Hazard

• FOR - Road Suitability (Natural Surface) (ME)

• FOR - Road Suitability (Natural Surface) (WA)

• FOR - Soil Rutting Hazard (OH)

• FOR - Soil Sustainability Forest Biomass Harvesting (CT)

• FOR - White Oak Suitability (MO)

• FOR-Biomass Harvest (WI)

• FOTG - Indiana Corn Yield Calculation (IN)

• GRL - Excavations to 24 inches for Plastic Pipelines (TX)

• GRL - Fencing, 24 inch Post Depth (MT)

• GRL - NV range seeding (Wind C = 100) (NV)

• GRL - NV range seeding (Wind C = 40) (NV)

• GRL - NV range seeding (Wind C = 60) (NV)

• GRL - NV range seeding (Wind C = 80) (NV)

• GRL - NV range seeding (Wind C >= 160) (NV)

• GRL - Rangeland Planting by Mechanical Seeding (TX)

• GRL - Rangeland Root Plowing (TX)

• Hybrid Wine Grape Varieties Site Desirability (Long)

• Low Pressure Pipe Septic System (DE)

• MIL - Bivouac Areas (DOD)

• MIL - Excavations Crew-Served Weapon Fighting Position (DOD)

• MIL - Excavations for Individual Fighting Position (DOD)

• MIL - Trafficability Veh. Type 1 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 2 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 4 1-pass wet season (DOD)

• MIL - Trafficability Veh. Type 4 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 6 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 7 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 7 dry season (DOD)

• NCCPI - Irrigated National Commodity Crop Productivity Index

• Nitrogen Loss Potential (ND)

• Potential Windthrow Hazard (TN)

• REC - Foot and ATV Trails (AK)

• REC - Playgrounds (AK)

• Reclamation Suitability (ND)

• RSK-risk assessment for manure application (OH)

• SAS - CMECS Substrate Origin

• SAS - CMECS Substrate Subclass/Group/Subgroup

• SAS - Mooring Anchor - Deadweight

• Septic System A/B Soil System (Alternate) (PA)

• Septic System CO-OP RFS III w/Spray Irrigation (PA)

• Septic System Dual Field Trench (conventional) (WV)

• Septic System Elevated Field (alternative) (WV)

• Septic System In Ground Trench (conventional) (PA)

• Septic System In Ground Trench (conventional) (WV)

• AGR - Filter Strips (TX)

• AGR - Hops Site Suitability (ID)

• AGR - Mulch Till (TX)

• AGR - Nitrate Leaching Potential, Nonirrigated (MT)

• AGR - Nitrate Leaching Potential, Nonirrigated (WV)

• AGR - No Till (VT)

• AGR - Oats Yield (MT)

• AGR - Pesticide Loss Potential-Leaching

• AGR - Pesticide Loss Potential-Leaching (NE)

• AGR - Rutting Hazard =< 10,000 Pounds per Wheel (TX)

• AGR - S. Highbush Blueberry Suitability MLRA 153 (SC)

• AGR - Wind Erosion Potential (NE)

• AGR-Available Water Capacity (ND)

• AGR-Physical Limitations (ND)

• AGR-Sodicity (ND)

• AGR-Surface Crusting (ND)

• AGR-Wind Erosion (ND)

• AWM - Irrigation Disposal of Wastewater (DE)

• AWM - Land App of Municipal Sewage Sludge (DE)

• AWM - Land App of Municipal Sewage Sludge (MD)

• AWM - Land Application of Milk (CT)

• AWM - Land Application of Municipal Biosolids, spring (OR)

• AWM - Land Application of Municipal Sewage Sludge

• AWM - Land Application of Municipal Sewage Sludge (OH)

• AWM - Land Application of Municipal Sewage Sludge (VT)

• AWM - Large Animal Disposal, Pit (MN)

• AWM - Manure and Food Processing Waste

• AWM - Manure and Food Processing Waste (VT)

• AWM - Rapid Infil Disposal of Wastewater (MD)

• AWM - Rapid Infiltration Disposal of Wastewater (VT)

• AWM - Slow Rate Process Treatment of Wastewater (VT)

• BLM - Chaining Suitability

• BLM - Fugitive Dust Resistance

• BLM - Soil Restoration Potential

• BLM - Yellow Star-thistle Invasion Susceptibility

• CLASS RULE - Depth to non-lithic bedrock (5 classes) (NPS)

• CLR-cropland limitation for corn and soybeans (IN)

• Commodity Crop Productivity Index (Corn) (WI)

• CPI - Grass Hay, NIRR - Klamath Valleys and Basins (OR)

• CPI - Potatoes Productivity Index (AK)

• CPI - Potatoes, IRR - Eastern Idaho Plateaus (ID)

• CPI - Small Grains, NIRR - Palouse Prairies (ID)

• DHS - Emergency Animal Mortality Disposal by Shallow Burial

• DHS - Rubble and Debris Disposal, Large-Scale Event

• ENG - Aquifer Assessment - 7081 (MN)

• ENG - Construction Materials - Gravel Source (MN)

• ENG - Construction Materials; Gravel Source (MI)

• ENG - Construction Materials; Gravel Source (OR)

• ENG - Construction Materials; Reclamation

• ENG - Construction Materials; Reclamation (OH)

• ENG - Construction Materials; Sand Source

• ENG - Construction Materials; Sand Source (AK)

• ENG - Construction Materials; Sand Source (ID)

• ENG - Construction Materials; Sand Source (IN)

• ENG - Construction Materials; Sand Source (OH)

• ENG - Construction Materials; Topsoil

• ENG - Construction Materials; Topsoil (WA)

• ENG - Ground-based Solar Arrays, Soil-based Anchor Systems

• ENG - Local Roads and Streets

• ENG - New Ohio Septic Rating (OH)

• ENG - Sanitary Landfill (Trench) (OH)

• ENG - Septic Tank Absorption Fields (AK)

• ENG - Septic Tank Absorption Fields (DE)

• ENG - Septic Tank Absorption Fields (NY)

• ENG - Sewage Lagoons (OH)

• ENG - Shallow Excavations (AK)

• ENG - Shallow Excavations (MI)

• ENG - Unpaved Local Roads and Streets

• FOR - Black Walnut Suitability Index (MO)

• FOR - Conservation Tree and Shrub Groups (TX)

• FOR - Construction Limitations - Haul Roads/Log Landing (OH)

• FOR - Construction Limitations For Haul Roads (MI)

• FOR - Hand Planting Suitability (ME)

• FOR - Harvest Equipment Operability (MD)

• FOR - Harvest Equipment Operability (OH)

• FOR - Harvest Equipment Operability (VT)

• FOR - Mechanical Planting Suitability

• FOR - Mechanical Planting Suitability (ME)

• FOR - Mechanical Planting Suitability, MO13 (DE)

• FOR - Potential Erosion Hazard (Off-Road/Off-Trail)

• FOR - Potential Erosion Hazard (Road/Trail) (PIA)

• FOR - Potential Seedling Mortality (VT)

• FOR - Potential Windthrow Hazard (NY)

• FOR - Potential Windthrow Hazard (VT)

• FOR - Puddling Potential (WA)

• FOR - Road Suitability (Natural Surface)

• FOR - Road Suitability (Natural Surface) (OH)

• FOR - Road Suitability (Natural Surface) (OR)

• FOR - Rutting Hazard by Season

• FOR - Shortleaf pine littleleaf disease susceptibility

• FOR - Soil Compactibility Risk

• FOR - Soil Rutting Hazard (ME)

• FOR - Windthrow Hazard

• FOR-Construction Limitations for Haul Roads/Log Landings(ME)

• FOTG - Indiana Slippage Potential (IN)

• Gravity Full Depth Septic System (DE)

• GRL - Fencing, Post Depth =<36 inches

• GRL - NV range seeding (Wind C = 50) (NV)

• GRL - Ranch Access Roads (TX)

• GRL - Rangeland Roller Chopping (TX)

• Ground Penetrating Radar Penetration

• Ground-based Solar Arrays_bedrock(ME)

• Ground-based Solar Arrays_bedrock_slope_ballast(ME)

• Hybrid Wine Grape Varieties Site Desirability (Short)

• ISDH Septic Tank Interpretation (IN)

• Land Application of Municipal Sewage Sludge (PA)

• MIL - Helicopter Landing Zones (DOD)

• MIL - Trafficability Veh. Type 2 1-pass wet season (DOD)

• MIL - Trafficability Veh. Type 5 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 5 dry season (DOD)

• MIL - Trafficability Veh. Type 7 1-pass wet season (DOD)

• NCCPI - National Commodity Crop Productivity Index (Ver 3.0)

• REC - Camp and Picnic Areas (AK)

• REC - Picnic Areas (CT)

• REC - Playgrounds (CT)

• SAS - CMECS Substrate Subclass

• Septic System Drip Irrigation (Alternate) (PA)

• Septic System Free Access Sand Filter w/Drip Irrigation (PA)

• Septic System In Ground Bed (conventional) (PA)

• Septic System Peat Based Option1 (UV & At-Grade Bed)Alt (PA)

• Septic System Peat Sys Opt3 w/Subsurface Sand Filter (PA)

• Septic System Sand Mound Bed or Trench (PA)

• Septic System Shallow Placement Pressure Dosed (Alt.) (PA)

• SOH - Aggregate Stability (ND)

• SOH - Agricultural Organic Soil Subsidence

• SOH - Dynamic Soil Properties Response to Biochar

• SOH - Organic Matter Depletion

• SOIL HEALTH ASSESSMENT (NJ)

• URB - Commercial Brick Bldg; w/Reinforced Concrete Slab (TX)

• URB - Reinforced Concrete Slab (TX)

• URB/REC - Camp Areas

• URB/REC - Camp Areas (OH)

• URB/REC - Off-Road Motorcycle Trails (OH)

• URB/REC - Paths and Trails (OH)

• URB/REC - Picnic Areas

• URB/REC - Playgrounds

• URB/REC - Playgrounds (GA)

• Vinifera Wine Grape Site Desirability (Short to Medium)

• WLF - Irr. Domestic Grasses & Legumes for Food & Cover (TX)

• WLF - Upland Coniferous Trees (TX)

• WLF - Upland Deciduous Trees (TX)

• WLF - Upland Desertic Shrubs & Trees (TX)

• WLF - Upland Native Herbaceous Plants (TX)

• WLF - Upland Shrubs & Vines (TX)

• WLF-Soil Suitability - Karner Blue Butterfly (WI)

• WMS - Drainage (IL)

• WMS - Drainage - (MI)

• WMS - Embankments, Dikes, and Levees

• WMS - Embankments, Dikes, and Levees (OH)

• WMS - Grassed Waterways - (MI)

• AGR - Air Quality; PM10 (TX)

• AGR - Air Quality; PM2_5 (TX)

• AGR - Aronia Berry Suitability (SD)

• AGR - Farmland of Statewide Importance (TX)

• AGR - Index for alfalfa hay, irrigated (NV)

• AGR - Nitrate Leaching Potential, Nonirrigated (MA)

• AGR - Rangeland Grass/Herbaceous Productivity Index (TX)

• AGR - Rutting Hazard > 10,000 Pounds per Wheel (TX)

• AGR - Water Erosion Potential (TX)

• AGR - Wine Grape Site Suitability (WA)

• AGR-Natural Fertility (ND)

• AGR-Subsurface Salinity (ND)

• AWM - Filter Group (OH)

• AWM - Irrigation Disposal of Wastewater

• AWM - Land Application of Dry and Slurry Manure (TX)

• AWM - Land Application of Municipal Biosolids, winter (OR)

• AWM - Overland Flow Process Treatment of Wastewater

• AWM - Rapid Infiltration Disposal of Wastewater

• AWM - Vegetated Treatment Area (PIA)

• AWM - Waste Field Storage Area (VT)

• BLM - Mechanical Treatment, Shredder

• BLM - Medusahead Invasion Susceptibility

• BLM - Soil Compaction Resistance

• Capping Fill Gravity Septic System (DE)

• CLASS RULE - Depth to any bedrock kind (5 classes) (NPS)

• CPI - Alfalfa Hay, IRR - Eastern Idaho Plateaus (ID)

• CPI - Alfalfa Hay, IRR - Klamath Valley and Basins (OR)

• CPI - Alfalfa Hay, IRR - Snake River Plains (ID)

• CPI - Alfalfa Hay, NIRR- Eastern Idaho Plateaus (ID)

• CPI - Grass Hay, NIRR - Palouse, Northern Rocky Mtns. (WA)

• CPI - Small Grains Productivity Index (AK)

• DHS - Catastrophic Event, Large Animal Mortality, Incinerate

• DHS - Emergency Land Disposal of Milk

• DHS - Site for Composting Facility - Subsurface

• DHS - Suitability for Clay Liner Material

• ENG - Cohesive Soil Liner (MN)

• ENG - Construction Materials - Sand Source (MN)

• ENG - Construction Materials; Gravel Source (CT)

• ENG - Construction Materials; Gravel Source (NY)

• ENG - Construction Materials; Reclamation (DE)

• ENG - Construction Materials; Roadfill

• ENG - Construction Materials; Roadfill (AK)

• ENG - Construction Materials; Sand Source (NY)

• ENG - Construction Materials; Sand Source (VT)

• ENG - Construction Materials; Topsoil (AK)

• ENG - Construction Materials; Topsoil (DE)

• ENG - Construction Materials; Topsoil (MI)

• ENG - Construction Materials; Topsoil (OR)

• ENG - Conventional On-Site Septic Systems (TN)

• ENG - Deep Infiltration Systems

• ENG - Disposal Field Gravity (DE)

• ENG - Dwellings With Basements (OH)

• ENG - Ground-based Solar Arrays, Ballast Anchor Systems

• ENG - Large Animal Disposal, Trench (CT)

• ENG - Lawn, Landscape, Golf Fairway (MI)

• ENG - Lawn, Landscape, Golf Fairway (VT)

• ENG - Sanitary Landfill (Area) (OH)

• ENG - Sanitary Landfill (Trench)

• ENG - Sanitary Landfill (Trench) (AK)

• ENG - Septage Application - Surface (MN)

• ENG - Septic Tank Absorption Fields - At-Grade (MN)

• ENG - Septic Tank Absorption Fields - Mound (MN)

• ENG - Septic Tank Leaching Chamber (TX)

• ENG - Septic Tank, Subsurface Drip Irrigation (TX)

• ENG - Shallow Excavations

• ENG - Shallow Infiltration Systems

• ENG - Small Commercial Buildings

• ENG - Soil Potential of Road Salt Applications (CT)

• ENG - Source of Caliche (TX)

• ENG - Stormwater Management / Ponds (NY)

• ENG - Unlined Retention Systems

• Farm and Garden Composting Facility - Surface

• FOR - Biomass Harvest (MA)

• FOR - Black Walnut Suitability Index (KS)

• FOR - Displacement Potential (WA)

• FOR - Drought Vulnerable Soils

• FOR - General Harvest Season (ME)

• FOR - Harvest Equipment Operability

• FOR - Mechanical Site Preparation (Deep) (MD)

• FOR - Mechanical Site Preparation (Surface)

• FOR - Mechanical Site Preparation; Deep (CT)

• FOR - Potential Erosion Hazard (Road/Trail)

• FOR - Potential Fire Damage Hazard

• FOR - Potential Seedling Mortality

• FOR - Potential Seedling Mortality (MI)

• FOR - Potential Windthrow Hazard (ME)

• FOR - Potential Windthrow Hazard (MI)

• FOR - Road Suitability (Natural Surface) (ID)

• FOR - Rutting Hazard by Month

• FOR - Windthrow Hazard (WA)

• FOTG - NLI Interp Calculation - (IN)

• Fragile Soil Index

• GRL - Juniper Encroachment Potential (NM)

• GRL - NV range seeding (Wind C = 20) (NV)

• GRL - Pasture and Hayland SG (OH)

• GRL - Rangeland Prescribed Burning (TX)

• GRL - Rangeland Soil Seed Bank Suitability (NM)

• GRL-FSG-NP-W (MT)

• GRL-SHSI Soil Health Sustainability Index (MT)

• Ground-based Solar Arrays_saturationt(ME)

• Ground-based Solar Arrays_slope(ME)

• Inland Wetlands (CT)

• IRR-restrictive features for irrigation (OH)

• MIL - Excavations for Vehicle Fighting Position (DOD)

• MIL - Trafficability Veh. Type 1 1-pass wet season (DOD)

• MIL - Trafficability Veh. Type 2 dry season (DOD)

• MIL - Trafficability Veh. Type 3 50-passes wet season (DOD)

• MIL - Trafficability Veh. Type 6 1-pass wet season (DOD)

• MIL - Trafficability Veh. Type 6 dry season (DOD)

• Muscadine Wine Grape Site Desirability (Very Long)

• NCCPI - NCCPI Cotton Submodel (II)

• Permafrost Sensitivity (AK)

• Pressure Dose Capping Fill Septic System (DE)

• REC - Camp Areas (CT)

• REC - Off-Road Motorcycle Trails (CT)

• SAS - CMECS Substrate Class

• SAS - CMECS Substrate Subclass/Group

• SAS - Eelgrass Restoration Suitability

• SAS - Land Utilization of Dredged Materials

• SAS - Northern Quahog (Hard Clam) Habitat Suitability

• Septic System At Grade Shallow Field (alternative) (WV)

• Septic System At-Grade Bed (Alternate) (PA)

• Septic System CO-OP RFS III w/Drip Irrigation (PA)

• Septic System Drip Irrigation (alternative) (WV)

• Septic System Free Access Sand Filterw/Spray Irrigation (PA)

• Septic System Peat Based Option1 w/At-Grade Bed (Alt.) (PA)

• Septic System Spray Irrigation (PA)

• Septic System Steep Slope Sand Mound (Alternate) (PA)

• Shallow Infiltration Systems

• SOH - Organic Matter Depletion Potential, Irrigated (CA)

• SOH - Soil Surface Sealing

• TROP - Plantains Productivity

• URB/REC - Camp Areas (GA)

• URB/REC - Camp Areas (MI)

• URB/REC - Golf Fairways (OH)

• URB/REC - Off-Road Motorcycle Trails

• URB/REC - Paths and Trails (MI)

• URB/REC - Playgrounds (OH)

• Vinifera Wine Grape Site Desirability (Long to Medium)

• WLF - Chufa for Turkey Forage (LA)

• WLF - Food Plots for Upland Wildlife < 2 Acres (TX)

• WLF - Freshwater Wetland Plants (TX)

• WLF - Irrigated Saline Water Wetland Plants (TX)

• WLF - Riparian Herbaceous Plants (TX)

• WLF - Riparian Shrubs, Vines, & Trees (TX)

• WLF - Saline Water Wetland Plants (TX)

• WLF - Upland Mixed Deciduous & Coniferous Trees (TX)

• WMS - Constructing Grassed Waterways (TX)

• WMS - Constructing Terraces and Diversions (OH)

• WMS - Embankments, Dikes, and Levees (VT)

• WMS - Irrigation, Sprinkler (close spaced outlet drops)

• WMS - Irrigation, Sprinkler (general)

• WMS - Pond Reservoir Area (GA)

• WMS-Subsurface Water Management, Installation (ND)

• WMS-Subsurface Water Management, Outflow Quality (ND)

• AGR - Barley Yield (MT)

• AGR - Conventional Tillage (TX)

• AGR - Grape non-irrigated (MO)

• AGR - Industrial Hemp for Fiber and Seed Production

• AGR - Nitrate Leaching Potential, Irrigated (WA)

• AGR - Pasture hayland (MO)

• AGR - Pesticide Loss Potential-Soil Surface Runoff

• AGR - Prime Farmland (TX)

• AGR - Spring Wheat Yield (MT)

• AGR-Agronomic Concerns (ND)

• AGR-Pesticide and Nutrient Leaching Potential, NIRR (ND)

• AGR-Surface Salinity (ND)

• AGR-Water Erosion Potential (ND)

• Alaska Exempt Wetland Potential (AK)

• American Wine Grape Varieties Site Desirability (Short)

• AWM - Irrigation Disposal of Wastewater (MD)

• AWM - Manure and Food Processing Waste (DE)

• AWM - Manure Stacking - Site Evaluation (TX)

• AWM - Phosphorus Management (TX)

• AWM - Slow Rate Process Treatment of Wastewater

• BLM - Pygmy Rabbit Habitat Potential

• BLM - Rangeland Tillage

• BLM - Site Degradation Susceptibility

• CA Prime Farmland (CA)

• CLASS RULE - Depth to root limiting layer (5 classes) (NPS)

• Commodity Crop Productivity Index (Corn) (TN)

• CPI - Alfalfa Hay, NIRR - Palouse, Northern Rocky Mtns. (ID)

• CPI - Barley, NIRR - Eastern Idaho Plateaus (ID)

• CPI - Grass Hay, IRR - Eastern Idaho Plateaus (ID)

• CPI - Grass Hay, NIRR - Palouse, Northern Rocky Mtns. (ID)

• CPI - Potatoes, IRR - Snake River Plains (ID)

• CPI - Small Grains, NIRR - Palouse Prairies (OR)

• CPI - Small Grains, NIRR - Palouse Prairies (WA)

• CPI - Small Grains, NIRR - Snake River Plains (ID)

• CPI - Wheat, NIRR - Eastern Idaho Plateaus (ID)

• CPI - Wild Hay, NIRR - Eastern Idaho Plateaus (ID)

• CPI - Wild Hay, NIRR - Palouse, Northern Rocky Mtns. (ID)

• CPI - Wild Hay, NIRR - Palouse, Northern Rocky Mtns. (WA)

• Deep Infiltration Systems

• DHS - Site for Composting Facility - Surface

• Elevated Sand Mound Septic System (DE)

• ENG - Animal Disposal by Composting (Catastrophic) (WV)

• ENG - Application of Municipal Sludge (TX)

• ENG - Closed-Loop Horizontal Geothermal Heat Pump (CT)

• ENG - Construction Materials; Gravel Source (IN)

• ENG - Construction Materials; Gravel Source (NE)

• ENG - Construction Materials; Reclamation (MD)

• ENG - Construction Materials; Reclamation (MI)

• ENG - Construction Materials; Roadfill (GA)

• ENG - Construction Materials; Sand Source (CT)

• ENG - Construction Materials; Sand Source (GA)

• ENG - Construction Materials; Topsoil (ID)

• ENG - Construction Materials; Topsoil (OH)

• ENG - Daily Cover for Landfill (OH)

• ENG - Disposal Field (NJ)

• ENG - Disposal Field Type Inst (NJ)

• ENG - Dwellings W/O Basements

• ENG - Dwellings With Basements

• ENG - Dwellings without Basements (AK)

• ENG - Lawn and Landscape (OH)

• ENG - Lawn, Landscape, Golf Fairway

• ENG - Local Roads and Streets (AK)

• ENG - Local Roads and Streets (GA)

• ENG - On-Site Waste Water Lagoons (MO)

• ENG - Pier Beam Building Foundations (TX)

• ENG - Sanitary Landfill (Area)

• ENG - Sanitary Landfill (Area) (AK)

• ENG - Septage Application - Incorporation or Injection (MN)

• ENG - Septic System; Disinfection, Surface Application (TX)

• ENG - Septic Tank Absorption Fields (FL)

• ENG - Septic Tank Absorption Fields (OH)

• ENG - Septic Tank Absorption Fields - Trench (MN)

• ENG - Sewage Lagoons (AK)

• ENG - Shallow Excavations (OH)

• ENG - Soil Suitability for SLAMM Marsh Migration (CT)

• ENG - Stormwater Management / Infiltration (NY)

• ENG - Stormwater Management / Wetlands (NY)

• FOR - Black Walnut Suitability (WI)

• FOR - Black Walnut Suitability (WV)

• FOR - Construction Limitations for Haul Roads/Log Landings

• FOR - Displacement Hazard

• FOR - Harvest Equipment Operability (DE)

• FOR - Harvest Equipment Operability (ME)

• FOR - Harvest Equipment Operability (MI)

• FOR - Log Landing Suitability (ID)

• FOR - Log Landing Suitability (MI)

• FOR - Log Landing Suitability (OR)

• FOR - Mechanical Planting Suitability (OH)

• FOR - Mechanical Site Preparation (Surface) (MD)

• FOR - Mechanical Site Preparation (Surface) (OH)

• FOR - Mechanical Site Preparation; Surface (CT)

• FOR - Potential Erosion Hazard (Off-Road/Off-Trail) (MI)

• FOR - Potential Erosion Hazard (Off-Road/Off-Trail) (OH)

• FOR - Potential Seedling Mortality (FL)

• FOR - Potential Seedling Mortality (OH)

• FOR - Road Suitability (Natural Surface) (VT)

• FOR - Soil Rutting Hazard

• FOTG - Indiana Soy Bean Yield Calculation (IN)

• FOTG - Indiana Wheat Yield Calculation (IN)

• FOTG - NLI report Calculation - (IN)

• GRL - Fencing, Post Depth =<24 inches

• GRL - Fencing, Post Depth Less Than 24 inches (TX)

• GRL - Fencing, Post Depth Less Than 36 inches (TX)

• GRL - NV range seeding (Wind C = 10) (NV)

• GRL - NV range seeding (Wind C = 30) (NV)

• GRL - Rangeland Chaining (TX)

• GRL - Rangeland Disking (TX)

• GRL - Rangeland Dozing/Grubbing (TX)

• GRL - Utah Juniper Encroachment Potential

• GRL - Western Juniper Encroachment Potential (OR)

• Ground-based Solar Arrays_bedrock_slope_anchor(ME)

• Ground-based Solar Arrays_saturation_flooding_Frost(ME)

• Hybrid Wine Grape Varieties Site Desirability (Medium)

• Lined Retention Systems

• MIL - Trafficability Veh. Type 1 dry season (DOD)

• MIL - Trafficability Veh. Type 3 1-pass wet season (DOD)

• MIL - Trafficability Veh. Type 3 dry season (DOD)

• MIL - Trafficability Veh. Type 4 dry season (DOD)

• MIL - Trafficability Veh. Type 5 1-pass wet season (DOD)

• NCCPI - NCCPI Corn Submodel (I)

• NCCPI - NCCPI Small Grains Submodel (II)

• NCCPI - NCCPI Soybeans Submodel (I)

• Peony Flowers Site Suitability (AK)

• Pressure Dose Full Depth Septic System (DE)

• REC - Camp Areas; Primitive (AK)

• REC - Paths and Trails (CT)

• Salinity Risk Index (ND)

• SAS - Eastern Oyster Habitat Restoration Suitability

• SAS - Mooring Anchor - Mushroom

• Septic System CO-OP RFS III w/At-Grade Bed (PA)

• Septic System Free Access Sand Filter w/At-Grade Bed (PA)

• Septic System Modified Subsurface Sand Filter (Alt.) (PA)

• Septic System Shallow In Ground Trench (conventional) (WV)

• Septic System Subsurface Sand Filter Bed (conventional) (PA)

• Septic System Subsurface Sand Filter Trench (standard) (PA)

• SOH - Limitations for Aerobic Soil Organisms

• URB - Concrete Driveways and Sidewalks (TX)

• URB - Dwellings on Concrete Slab (TX)

• URB - Lawns and Ornamental Plantings (TX)

• URB/REC - Paths and Trails

• URB/REC - Paths and Trails (GA)

• URB/REC - Playgrounds (MI)

• Vinifera Wine Grape Site Desirability (Long)

• WLF - Crawfish Aquaculture (TX)

• WLF - Desertic Herbaceous Plants (TX)

• WLF - Gopher Tortoise Burrowing Suitability

• WLF - Grain & Seed Crops for Food and Cover (TX)

• WMS - Constructing Grassed Waterways (OH)

• WMS - Irrigation, Surface (graded)

• WMS - Subsurface Drains - Installation (VT)

• WMS - Subsurface Water Management, System Performance

• WMS - Surface Drains (TX)

• WMS - Surface Irrigation Intake Family (TX)

• Septic System Low Pressure Pipe (alternative) (WV)

• Septic System Mound (alternative) (WV)

• Septic System Peat Based Option2 w/Spray Irrigation (PA)

• Septic System Steep Slope Mound (alternative) (WV)

• SOH - Concentration of Salts- Soil Surface

• SOH - Soil Susceptibility to Compaction

• Soil Habitat for Saprophyte Stage of Coccidioides

• Unlined Retention Systems

• URB - Commercial Metal Bldg; w/Reinforced Concrete Slab (TX)

• URB/REC - Picnic Areas (GA)

• URB/REC - Picnic Areas (MI)

• URB/REC - Picnic Areas (OH)

• Vinifera Wine Grape Site Desirability (Short)

• WLF - Burrowing Mammals & Reptiles (TX)

• WLF - Desert Tortoise (CA)

• WLF - Domestic Grasses & Legumes for Food and Cover (TX)

• WLF - Irrigated Grain & Seed Crops for Food & Cover (TX)

• WMS - Excavated Ponds (Aquifer-fed)

• WMS - Excavated Ponds (Aquifer-fed) (VT)

• WMS - Irrigation, General

• WMS - Irrigation, Micro (above ground)

• WMS - Irrigation, Micro (above ground) (VT)

• WMS - Irrigation, Micro (subsurface drip)

• WMS - Irrigation, Sprinkler (general) (VT)

• WMS - Pond Reservoir Area

• WMS - Pond Reservoir Area (OH)

• WMS - Subsurface Water Management, System Installation

• WMS - Constructing Terraces & Diversions (TX)

• WMS - Drainage (OH)

• WMS - Excavated Ponds (Aquifer-fed) (OH)

• WMS - Grape Production with Drip Irrigation (TX)

• WMS - Irrigation, Micro (subsurface drip) (VT)

• WMS - Irrigation, Surface (level)

• WMS - Pond Reservoir Area (MI)

• WMS - Pond Reservoir Area (VT)

• WMS - Sprinkler Irrigation (MT)

• WMS - Sprinkler Irrigation RDC (IL)

• WMS - Subsurface Drains - Performance (VT)

• WMS - Subsurface Water Management, Outflow Quality

• WMS - Surface Water Management, System

• WMS-Subsurface Water Management, Performance (ND)

Value

a data.frame

Author(s)

Jason Nemecek, Chad Ferguson, Andrew Brown

Examples

# get two forestry interpretations for CA630
  get_SDA_interpretation(c("FOR - Potential Seedling Mortality",
                           "FOR - Road Suitability (Natural Surface)"),
                         method = "Dominant Condition",
                         areasymbols = "CA630")

Description

Obtain pre-calculated tabular reports of usage, activities, areas of interest (AOI), exports, ecological sites, ratings and reports for specific areas, times and intervals.

Usage

get_SDA_metrics(query_name, query_frequency, query_year, state = NULL)

Arguments

Value

A data.frame containing query results

Author(s)

Jason Nemecek

Examples

## Not run: 
get_SDA_metrics('SDA_Usage', 'CY', 2019:2021)

## End(Not run)

Description

Get map unit aggregate attribute information from Soil Data Access

Usage

get_SDA_muaggatt(
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  query_string = FALSE,
  dsn = NULL
)

Arguments

Value

a data.frame

Author(s)

Jason Nemecek, Chad Ferguson, Andrew Brown

Description

Get map unit parent material group information from Soil Data Access

Usage

get_SDA_pmgroupname(
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  method = "DOMINANT COMPONENT",
  simplify = TRUE,
  miscellaneous_areas = FALSE,
  query_string = FALSE,
  dsn = NULL
)

Arguments

Details

Default method is "Dominant Component" to get the dominant component (highest percentage). Use "Dominant Condition" or dominant parent material condition (similar conditions aggregated across components). Use "None" for no aggregation (one record per component).

Value

a data.frame

Author(s)

Jason Nemecek, Chad Ferguson, Andrew Brown

Description

Get map unit properties from Soil Data Access

Usage

get_SDA_property(
  property,
  method = c("Dominant Component (Category)", "Weighted Average", "Min/Max",
    "Dominant Component (Numeric)", "Dominant Condition", "None"),
  areasymbols = NULL,
  mukeys = NULL,
  WHERE = NULL,
  top_depth = 0,
  bottom_depth = 200,
  FUN = NULL,
  include_minors = FALSE,
  miscellaneous_areas = FALSE,
  query_string = FALSE,
  dsn = NULL
)

Arguments

Details

The property argument refers to one of the property names or columns specified in the tables below. Note that property can be specified as either a character vector of labeled properties, such as "Bulk Density 0.33 bar H2O - Rep Value", OR physical column names such as "dbthirdbar_r". To get "low" and "high" values for a particular property, replace the ⁠_r⁠ with ⁠_l⁠ or ⁠_h⁠ in the physical column name; for example property = c("dbthirdbar_l","dbthirdbar_r","dbthirdbar_h"). You can view exhaustive lists of component and component horizon level properties in the Soil Data Access "Tables and Columns Report".

Selected Component-level Properties

Selected Horizon-level Properties

Value

a data.frame with result

Author(s)

Jason Nemecek, Chad Ferguson, Andrew Brown

Examples

# get 1/3 bar bulk density [0,25] centimeter depth weighted average from dominant component
 get_SDA_property(property = c("dbthirdbar_l","dbthirdbar_r","dbthirdbar_h"),
                  method = "Dominant Component (Numeric)",
                  areasymbols = "CA630",
                  top_depth = 0,
                  bottom_depth = 25)