HBV Rainfall-Runoff Model (C++)
By Matteo Giuliani, Josh Kollat, Jon Herman, and others.
HBV Rainfall-runoff model, based on the work by (Bergstrom 1995). Runs on a daily timestep and saves all states and fluxes from each day for further analysis.
Both simulation and optimization (calibration) are available. Simulation mode is currently configured to read multiple parameter sets from
stdin and evaluate them in order. Calibration is currently configured to work with MOEAFramework, but may be easily modified for use with another application.
example_data/: Example forcing data files showing the input format
hbv_model.h: Defines the
HBVclass to store all states and fluxes at each timestep over the course of the evaluation.
hbv_model.cpp: Defines the functions for the processes in the model: degree-day snow, PDM soil moisture, Hamon PE, and the water balance between reservoirs.
main_HBV.cpp: Defines the initialization function (called once), the calculation function (called for each model evaluation), and the main function
CalHBV.java: Example Java class for calibration with MOEAFramework (optional).
moeaframework.c/h: Required libraries for communication with stdin/out
utils.cpp/h: Utilities for vector operations
To compile and run:
maketo compile. Modify the makefile first to use a different compiler or flags.
./SimHBV my_forcing_data.txt my_output_file.txt < my_parameter_samples.txtto perform simulation
- For calibration using MOEAFramework, follow the instructions for connecting an external optimization problem here. More detailed instructions are available from the MOEAFramework Setup Guide.
- Note that the second argument (the output filename) is only available in simulation mode.
my_forcing_data.txt: see the
example_data/directory for the format being used.
my_output_file.txt: name of file to output performance metric(s) (simulation mode only)
my_parameter_samples.txt: parameter sets to be evaluated in the model, with one parameter per column (e.g., hbv_param.txt). Currently there are 12 parameters being read into the model, which would correspond to 12 columns per row of this file. The parameters are read from
stdin, hence the
<operator to pipe the contents of the file to the executable. The order of parameters to be read in can be modified at
In its current form, the model will output (or optimize) three objectives: the relative variability (alpha), absolute value of the relative bias (beta) and the correlation (r) between the simulated and observed flows over the simulated time period. These objectives represent three components of the Nash Sutcliffe Efficiency (see Gupta et al. (2009)). However, the output can easily be modified to include any combination of states/fluxes or error metrics from any time during the simulation.
Based on work from the following paper: Herman, J.D., P.M. Reed, and T. Wagener (2013), Time-varying sensitivity analysis clarifies the effects of watershed model formulation on model behavior, Water Resour. Res., 49, doi:10.1002/wrcr.20124. (Link to Paper)
Copyright (C) 2010-2017 Matteo Giuliani, Josh Kollat, Jon Herman, and others.
HBV is free software: you can redistribute it and/or modify it under the terms of the GNU Lesser General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
HBV is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU Lesser General Public License for more details.
You should have received a copy of the GNU Lesser General Public License along with HBV. If not, see http://www.gnu.org/licenses/.