In USGS-R/powstreams: powstreams

Powell Streams data report

library(powstreams)
library(dplyr)
library(choroplethr)
library(dataRetrieval)
library(choroplethrMaps)
data(state.regions)

sites <- list_sites(with_timeseries = 'doobs')
df <- state.regions[,c(1,3)]
names(df)[2] <- 'value'
df[,2] <- 0
nwis_sites <- mda.streams:::parse_site_name(sites)
site_data <- readNWISsite(nwis_sites)
for (i in 1:length(sites)){

  region <- state.regions[as.numeric(site_data$state_cd[i])==state.regions[,3], 1]
  u.i = df[,1] == region
  if (sum(u.i) == 1){
    df[u.i, 2]<- df[u.i, 2]+1
  }
}

sites_w_either <- unique(c(list_sites(with_timeseries = c('doobs','wtr','disch')), list_sites(with_timeseries = c('doobs','wtr','stage'))))

state_choropleth(df)

sites with DO data: r length(list_sites(with_timeseries = 'doobs'))

sites with DO, temp, and discharge data: r length(list_sites(with_timeseries = c('doobs','wtr','disch')))

sites with DO, temp, discharge AND/OR stage: r length(sites_w_either)

sites with DO and nitrate: r length(list_sites(with_timeseries = c('doobs','nitrate')))

Range of DO for sites with DO and temperature: