Skip to content

The sensitivity of hydrologic parameters and processes across North America

License

Notifications You must be signed in to change notification settings

julemai/xSSA-North-America

Folders and files

NameName
Last commit message
Last commit date

Latest commit

 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 
 

Repository files navigation

The Sensitivity of Simulated Streamflow to Individual Hydrologic Processes Across North America

by Juliane Mai 1, James R. Craig 1, Bryan A. Tolson 1, and Richard Arsenault 2

1 Dept. of Civil and Environmental Engineering, University of Waterloo, Waterloo, ON, Canada.
2 Dept. of Construction Engineering, École de technologie supérieure, Montreal, QC, Canada.

Abstract

Streamflow sensitivity to different hydrologic processes varies in both space and time. This sensitivity is traditionally evaluated for the parameters specific to a given hydrologic model simulating streamflow. In this study, we apply a novel analysis over more than 3000 basins across North America considering a blended hydrologic model structure, which includes not only parametric, but also structural uncertainties. This enables seamless quantification of model process sensitivities and parameter sensitivities across a continuous set of models. It also leads to high-level conclusions about the importance of water cycle components on streamflow predictions, such as quickflow being the most sensitive process for streamflow simulations across the North American continent. The results of the 3000 basins are used to derive an approximation of sensitivities based on physiographic and climatologic data without the need to perform expensive sensitivity analyses. Detailed spatio-temporal inputs and results are shared through an interactive website.

Usage

The provided scripts enable the users to run their own extended Sobol' sensitivity analysis (xSSA). The forcing data for the basin need to be provided under data_in/data_obs/<basin-id>. See the provided examples for details on the format (Raven standard forcing format). The physical characteristics such as forest cover and basin-average elevation need to be provided under data_in/basin_metadata/basin_physical_characteristics.txt which currently contains the information for the 3826 HYSETS basins used in this study. The user also needs to provide an executable of the Raven modeling framework annd place it under data_in/data_model/Raven.exe. To run an analysis for one basin (e.g., USGS gauge station 03054500) use the following commands. The final results will be stored under data_out/<basin-id>/.

cd scripts

bb='03054500'           # an example basin
nsets='1000'            # number of iterations for analysis
                        # nsets='1000' takes about 22h
                        # number of model runs = (43+2) x nsets + (27+2) x nsets + (11+2) x nsets
tmpdir='/tmp/test-xssa' # a temporary directory; use ${SLURM_TMPDIR} for jobs on super-computers such as Graham

python raven_sa-usgs-canopex.py -i ${bb} -n ${nsets} -t ${tmpdir} -o nc

Website

The results of this study are available on an interactive webpage. The website includes the download of figures and model setups as well as the visualization of sensitivities of parameters, process options and processes.

Creating Plots

This GitHub contains all scripts and data to reproduce the plots in the paper and Supplementary Material. Please see the main plotting script here and select the figures you want to plot. The final figures will then be found in the folder figures/. All the data used to produce those figures can be found in the folder scripts/data/.

Citation

Journal Publication

Mai, J., Craig, J. R., Tolson, B. A., and R. Arsenault (2021).
The Sensitivity of Simulated Streamflow to Individual Hydrologic Processes Across North America.
Nature Communications, 13(1), 455.
https://doi.org/10.1038/s41467-022-28010-7

See the following publication for details on the introduction of the xSSA method and the Blended Model:
Mai, J., Craig, J. R., and Tolson, B. A. (2020).
Simultaneously determining global sensitivities of model parameters and model structure.
Hydrol. Earth Syst. Sci., 24, 5835–5858.
https://doi.org/10.5194/hess-24-5835-2020

See the following publication for details on calibration of the Blended Model:
Chlumsky, R., Mai, J., Craig, J. R., & Tolson, B. A. (2021).
Simultaneous Calibration of Hydrologic Model Structure and Parameters Using a Blended Model.
Water Resources Research, 57(5), e2020WR029229.
http://doi.org/10.1029/2020WR029229

See the following publication for details about the HYSETS database:
Arsenault, R., Brissette, F., Martel, J.-L., Troin, M., Lévesque, G., Davidson-Chaput, J., et al. (2020).
A comprehensive, multisource database for hydrometeorological modeling of 14,425 North American watersheds.
Scientific Data, 7, 1–12.
http://doi.org/10.1038/s41597-020-00583-2

Code and Data Publication

J. Mai, J. R. Craig, B. A. Tolson, and R. Arsenault (2021).
The Sensitivity of Simulated Streamflow to Individual Hydrologic Processes Across North America.
Zenodo
DOI

About

The sensitivity of hydrologic parameters and processes across North America

Resources

License

Stars

Watchers

Forks

Packages

No packages published