Development of an Automatic Calibration Framework for Hydrologic Modelling Using Approximate Bayesian Computation

Authors: A. Chowdhury, P. Egodawatta, J. M. McGree, A. Goonetilleke

Abstract:

Hydrologic models are increasingly used as tools to predict stormwater quantity and quality from urban catchments. However, due to a range of practical issues, most models produce gross errors in simulating complex hydraulic and hydrologic systems. Difficulty in finding a robust approach for model calibration is one of the main issues. Though automatic calibration techniques are available, they are rarely used in common commercial hydraulic and hydrologic modelling software e.g. MIKE URBAN. This is partly due to the need for a large number of parameters and large datasets in the calibration process. To overcome this practical issue, a framework for automatic calibration of a hydrologic model was developed in R platform and presented in this paper. The model was developed based on the time-area conceptualization. Four calibration parameters, including initial loss, reduction factor, time of concentration and time-lag were considered as the primary set of parameters. Using these parameters, automatic calibration was performed using Approximate Bayesian Computation (ABC). ABC is a simulation-based technique for performing Bayesian inference when the likelihood is intractable or computationally expensive to compute. To test the performance and usefulness, the technique was used to simulate three small catchments in Gold Coast. For comparison, simulation outcomes from the same three catchments using commercial modelling software, MIKE URBAN were used. The graphical comparison shows strong agreement of MIKE URBAN result within the upper and lower 95% credible intervals of posterior predictions as obtained via ABC. Statistical validation for posterior predictions of runoff result using coefficient of determination (CD), root mean square error (RMSE) and maximum error (ME) was found reasonable for three study catchments. The main benefit of using ABC over MIKE URBAN is that ABC provides a posterior distribution for runoff flow prediction, and therefore associated uncertainty in predictions can be obtained. In contrast, MIKE URBAN just provides a point estimate. Based on the results of the analysis, it appears as though ABC the developed framework performs well for automatic calibration.

Keywords: Automatic calibration framework, approximate Bayesian computation, hydrologic and hydraulic modelling, MIKE URBAN software, R platform.

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1338892

References:

[1] C. Zoppou, “Review of urban storm water models”, Environmental Modelling & Software, vol. 16, no. 3, pp. 195-231, 2001.
[2] H. Medsen, G. Wilson, H. C. Ammentorp, “Comparison of different automated strategies for calibration of rainfall-runoff models”, Journal of Hydrology, Vol. 261, pp. 48-59, 2002.
[3] D. Pilgrim, “Model evaluation, testing and parameter estimation in hydrology”, National Symposium on Hydrology, Australian Academy of Science, Canberra, 1975.
[4] M. URBAN, “Mike Urban Model User Manual”, Danish Hydraulic, 2012.
[5] R. Rasmussen, and G. Hamilton, “An approximate Bayesian computation approach for estimating parameters of complex environmental processes in a cellular automata”, Environmental Modelling & Software, vol. 29, no. 1, pp. 1-10, 2012.
[6] K. Csilléry, M. G. Blum, O. E. Gaggiotti, and O. François, “Approximate Bayesian computation (ABC) in practice”, Trends in ecology & evolution, vol. 25, no. 7, pp. 410-418, 2010.
[7] P. Egodawatta, “Translation of small-plot scale pollutant build-up and wash-off measurements to urban catchment scale”, Queensland University of Technology, Brisbane, PhD Thesis, 2007.
[8] F. Jabot, T. Faure and N. Dumoulin, “EasyABC: performing efficient approximate Bayesian computation sampling schemes using R”, Methods in Ecology and Evolution, vol. 4, no. 7, pp. 684-687, 2013.
[9] P. Torfs, and C. Brauer, “A (very) short introduction to R”, Hydrology and Quantitative Water Management Group, Wageningen University, The Netherlands, 2014.
[10] G. O’Loughlin and B. Stack, “DRAINS user manual”, Watercom Pty Ltd, 2014.
[11] S. C. Service, “Time of Concentration”, In National Engineering Handbook: Natural Resources Conservation Service, 2010.
[12] K. Csilléry, M. G. Blum, O. E. Gaggiotti, and O. François, “Approximate Bayesian computation (ABC) in practice”, Trends in ecology & evolution, vol. 25, no. 7, pp. 410-418, 2010.
[13] A. E. Raftery, and J. Gill, “Bayesian statistics: one-day course for the American Sociological Association”, Bayesian statistics, pp. 1-61, 2002.
[14] M. A. Beaumont, “Approximate Bayesian computation in evolution and ecology” Annual review of ecology, evolution, and systematics, vol. 41, pp. 379-406, 2010.
[15] M. A. Beaumont, W. Zhang, and D. J. Balding, “Approximate Bayesian computation in population genetics”, Genetics, vol. 162, no. 4, pp. 2025- 2035, 2002
[16] P. Marjoram, J. Molitor, V. Plagnol, and S. Tavaré, “Markov chain Monte Carlo without likelihoods”, Proceedings of the National Academy of Sciences, Vol. 100, No. 26, pp. 15324-15328, 2003.
[17] J. K. Pritchard, M. T. Seielstad, A. Perez-Lezaun, and M. W. Feldman, “Population growth of human Y chromosomes: a study of Y chromosome microsatellites”, Molecular Biology and Evolution, vol. 16, no. 12, pp. 1791-1798, 1999.
[18] K. Loague and R. E. Green, “Statistical and graphical methods for evaluating solute transport models: overview and application”, Journal of contaminant hydrology, vol. 7, no. 1, pp. 51-73, 1991.
[19] I. Green and D. Stephenson, “Criteria for comparison of single eventmodels”, Hydrological Sciences Journal, vol. 31, no. 3, pp. 395- 411, 1986.
[20] H. J. Henriksen, L. Troldborg, P. Nyegaard, T. O. Sonnenborg, J. C. Refsgaard and B. Madsen, “Methodology for construction, calibration and validation of a national hydrological model for Denmark”, Journal of Hydrology, vol. 280, pp. 52-71, 2003.
[21] G. Maps, Download link: https://www.google.com.au/maps/place/ 28%C2%B000%2754.4%22S+153%C2%B020%2710.1%22E/@- 28.0146414,153.3353743,17z/data=!3m1!1e3