The SNOTEL network is composed of over 800 automated data collection sites located in remote, high-elevation mountain watersheds in the western U.S. They are used to monitor snowpack, precipitation, temperature, and other climatic conditions. The data collected at SNOTEL sites are transmitted to a central database. This package queries this centralized database to provide easy access to these data and additional seasonal metrics of snow accumulation (snow phenology).
The SNOTEL network consists of a vast number of observation sites,
all of them listed together with their meta-data on the SNOTEL website.
snotel_info() function allows you to query this table
and import it as a neat table into
R. Some of the
meta-data, in particular the site id (
site_id), you will
need of you want to download the data for a site. You can save this
table to disk using the
path variable to specify a location
on your computer where to store the data as a csv. If this parameter is
missing the data is returned as an
# download and list site information <- snotel_info() site_meta_data head(site_meta_data)
If you downloaded the meta-data for all sites you can make a
selection using either geographic coordinates, or
columns. For the sake of brevity I’ll only query data for one site using
site_id below. By default the data, reported in
imperial values, are converted to metric measurements.
# downloading data for a random site <- snotel_download( snow_data site_id = 670, internal = TRUE ) # show the data head(snow_data)
# A plot of snow accummulation through the years plot(as.Date(snow_data$date), $snow_water_equivalent, snow_datatype = "l", xlab = "Date", ylab = "SWE (mm)" )
Although the main function of the package is to provide easy access
to the SNOTEL data a function
provided to calculate seasonal metrics of snow deposition.
# calculate snow phenology <- snotel_phenology(snow_data)phenology
# subset data to the first decade of the century <- subset(snow_data, as.Date(date) > as.Date("2000-01-01") & snow_data_subset as.Date(date) < as.Date("2010-01-01")) # plot the snow water equivalent time series plot(as.Date(snow_data_subset$date), $snow_water_equivalent, snow_data_subsettype = "l", xlab = "Date", ylab = "SWE (mm)" ) # plot the dates of first snow accumulation as a red dot points(phenology$first_snow_acc, rep(1,nrow(phenology)), col = "red", pch = 19, cex = 0.5 )
A list of all provided snow phenology statistics is provided below.
|year||The year in which the an event happened|
|first_snow_melt||day of first full snow melt (in DOY)|
|cont_snow_acc||start of continuous snow accumulation / retention (in DOY)|
|last_snow_melt||day on which all snow melts for the remaining year (in DOY)|
|first_snow_acc||day on which the first snow accumulates (in DOY)|
|max_swe||maximum snow water equivalent value during a given year (in mm)|
|max_swe_doy||day on which the maximum snow water equivalent value is reached (in DOY)|
Please use the proper Zenodo DOI when using this software for research purposes.