We pull this characteristic, but I think it is the temperature of a water grab sample. Note really cold temperatures throughout the year:
sample_data <- readRDS('1_wqp_pull/out/wqp_data.rds') %>%
filter(CharacteristicName %in% 'Water, sample') %>%
mutate(doy = lubridate::yday(ActivityStartDate)
ggplot(sample_data, aes(x = doy, y = ResultMeasureValue)) +
geom_point()

For now will remove these after the WQP pull (in munge_wqp_files.R) but we should remove this from the WQP pull params.