Week of 09/02/2024 - 09/06/2024

[COBrogan]

DEQ Updates:
Shapefiles of all gage watershed coverages are now available on the online version of the VA Hydro repo as usgsGageWatershedGeofield.csv. Alternatively, these are in media/model/met/ on the server. The easiest way to read these in and use them in plots is shown below.

Plot coverages

library(sf)
library(leaflet)
#Get the watershed coverage from the server
watershedGeo <- read.csv("http://deq1.bse.vt.edu:81/met/usgsGageWatershedGeofield.csv")
#Read in the gageWatershed data as an SF object. Specify the coordinate system
#and the field name in the data frame that contains the well-known text. In this
#case, WKT is the name of the field with the polygon geometries
gageWatershedSF <- st_as_sf(x = watershedGeo,wkt = "WKT",crs = 4326)
#You can plot sf objects like any other R object. You can add "add = TRUE" to
#add them to an existing plot
plot(gageWatershedSF,axes = TRUE)
#Interactive map to help navigate. Harder to share, but may help you moving
#forward
leaflet() %>%  
  addProviderTiles("Esri.WorldStreetMap") %>%
  addPolygons(data = gageWatershedSF,popup = gageWatershedSF$hydrocode)

• On-boarding Nathaniel:
  • Ensure he has access to server, dbase, and can sftp
  • Practice git needed? If so, schedule a call
  • Any issues so far?
• Storm Volume Method Analysis - Storm Volume method now run for all three data sources. Some differences are apparent. Now we need some GIS based analysis:
  • Can we look at a few targeted watersheds where the method performed well for one data source but not others?
  • December has an interesting drop-off in NLDAS. Can we see how this looks at various watersheds compared to other data sources?
  • Do residual plots for NLDAS in those watershed reveal any trends?
  • How does monthly average rainfall compare at these watersheds between methods?
  • Other Qualitative QC storm separation:
    • Are we capturing all the big events?
    • Are we classifying non-storms/partial storms/tiny flow increase as storms? Are multi-day events behaving as intended e.g. do we merge events when they are similar at start/end date?
    • Can we confirm any patterns in the regression results? If we find seasonality, can we visualize or quantify why it's occurring using our basic routines? Or can we establish potential next steps to refine the separation routine?
    • Do we have larger regression errors in summer time? Are we influenced greatly by big events? Look at error and residual plots. Are there mismatches in precipitation volume and storm flow volume e.g. does time of travel exceed 1 day
    • Can we start to find differences in precip timing/magnitude between the datasets and is this connected to the poor or good correlations that result from the linear regressions?
• Simple lm Model Runs
  • run batch of simple method watersheds
  • Can we look at a few targeted watersheds where the method performed well for one data source but not others?
  • Do comparison of storm volume with simple method
• Raster intersection types and performance: Raster Summary: Test st_intersection vs st_clip vs map_algebra model_meteorology#84
  • Establish WDM workflow Workflow wdm model_meteorology#72
    • Add temporal disaggregation Integrate temporal disaggregation step into wdm workflow. model_meteorology#86
    • Add resampling Do resample mode model_meteorology#87
    • Discuss any errors encountered
• Need to make routine to extract coverage rasters, and stash them in a web-accessible directory
• Feedback or questions about Waterjam Poster, with a shareable powerpoint presentation
• Spatial Comparisons
  • Are their trends in our methods over the full geographic area (as well as temporal)
  • Maybe color code watersheds to see what their ratings are compared to one another
• Data model
  • Use om_vahydro_metric_grid() in R to get results of QA and model performance.
  • [ ]

[COBrogan]

@rburghol Could you send Nate an email with his login info to the server?

[rburghol]

Will do @COBrogan !!