The streamflow evaluation tool (ROSET) compares modeled streamflow to in situ USGS monitoring sites, with interactive visualizations supporting an in-depth analysis.
National-scale streamflow modeling remains a modern challenge, as changes in the underlying hydrology from land use and land cover (LULC) change, anthropogentic streamflow modification, and general process components (reach length, hydrogeophysical processes, precipitation, temperature, etc) greatly influence hydrological modeling. In a changing climate, there is a need to anticipate flood intensity, impacts of groundwater depletion on streamflow, western montain low-flow events, eastern rain-on-snow events, storm-induced flooding, and other severe environmental problems that challenge the management of water resources. Given the National Water Model (NWM) bridges the gap between the spatially coarse USGS streamflow observations by providing a near-continuous 2.7 million reach predictions of streamflow, there lies the potential to improve upon the capabilities of the model by characterizing predictive performance across the heterogeneity of processes and land covers present at the national scale. The python-based Streamflow_Evaluator package provides a tool to evaluate national hydrogrphy dataset (nhd) based model outputs with colocated USGS/NWIS streamflow monitorng stations (parameter: 060). The package contains three key methods for evaluation: state-based LULC, HUC level analysis, and USGS station-based analysis. Below is a description of each method and application. While designed to us the NWM version 2.1 retrospective dataset, with minimal modification the tool should be able to take in other model formulations. By using the Streamflow_Evaluator tool, researchers can identify locations where a model may benefit from further training/calibration/parameterization or a need for new model processes/features (e.g., integration of reservoir release operaitons) to ultimately create new post-processing methods and/or hydrological modeling formulations to improve streamflow prediction capabilities with respect to modeling needs (e.g., stormflow, supply, emergency management, flooding, etc).
The Research-Oriented Streamflow Evaluation Tool (ROSET) requires access to USGS/NWIS and colocated model output data for comparision. While all data is publically available through the respective agencies, data download for analysis prevents timely model evaluation. The Alabama Water Institute at the University of Alabama hosts NWM v2.1 for all colocated USGS monitoring stations, and with the package capable of accessing the observed and predicted data, the package supports a fast and repeatable tool for evaluating modeled streamflow performance. We currently provide cloud hosted data via Box, please contact us for access.
Python: Version 3.8 or later
The included .yaml file should set up the correct Streamflow_Evaluator environment. Below is a list of required packages to operate the tool:
| os | hydrotools | pandas |
|---|---|---|
| numpy | matplotlib | sklearn |
| hydroeval | dataretrieval | streamstats |
| geopandas | Ipython | warnings |
| progressbar | datetime | folium |
| mapclassify | time | jenkspy |
| hvplot | holoviews | bokeh |
| branca | vincent | json |
| proplot | pygeohydro | pygeoutils |
| switfter | geopy |
Each streamflow evaluation method requires similar inputs, including a start date, end date, and model. There are currently three different evaluation classes, each providing the user a unique method for evaluating streamflow modeling performance:
- Class Eval_State(): Modeled Streamflow Evaluation by StreamStats
- Class HUC_Eval(): Modeled Streamflow Evaluation by Hydrologic Unit Code (HUC)
- Class Reach_Eval(): - NHD - USGS Streamflow Evaluation
For all examples the prediction are from the NWM v2.1 retrospective. Please see the Examples folder for more information on applying each specific class.
To determine how LULC affects the predictive performance of streamflow models, the Streamflow_Evaluator uses StreamStats to categorize the watershed upstream of each USGS monitoring site by watershed charateristics. Please see the State Land Use - Land Cover Evaluation.md readme to use the tool.
Running the Streamflow_Evaluator LULC_Eval class loads, processes, and visualizes model performance for the state, category, and size of interest
By clicking on a marker a popup of the modeled vs. observed performance at the inputted temporal frequency will appear
The Streamflow_Evaluator supports an interactive engagement with model results
The HUC_Eval class allows the user to evaluate modeled streamflow with observed in situ NWIS monitoring sites for a watershed(s) of interest. The user can input multiple watersheds (e.g., Great Salt lake: ['1601', '1602']) The user must enter a start date, end date, watersheds and model to compare (NWM v2.1 is set up). NWM retrospective data spans from 1980 - 2020, USGS/NWIS data is location dependent.
The HUC_Eval class loads all USGS and colocated modeled NHD reaches for evaluation. Color coding of the markers allows for quick identification of poor and well performing reaches and clicking on the markers will put up a modeled vs. observed graph at the desired temporal resolution.
Similar to the State_Eval class, the HUC_Eval class supports a more in-dpeth graphical analysis of the modeled vs. observed using the holoviews package
The Reach_Eval class allows the user to evaluate modeled streamflow with selected NWIS monitoring sites of interest. The user can input multiple USGS sites (e.g., ['02378780', '02339495', '02342500']) Similar to the other classes, enter a start date, end date, and model to compare (NWM v2.1 is set up). NWM retrospective data spans from 1980 - 2020, USGS/NWIS data is location dependent
The Reach_Eval class quickly loads and maps the selected USGS streamflow monitoring locations, along with the colocated model predictions of NHD reaches. The color code of the marker indicates the model performance at the respective USGS monitoring station site, and by clicking on a marker, the interactive map produces a graph at the desired temporal resolution.
Similar to the State_Eval and HUC_Eval classes, the Reach_Eval class supports a more in-dpeth graphical analysis of the modeled vs. observed using the holoviews package