Multiple model combination methods for annual maximum water level prediction during river ice breakup

The Athabasca River is the largest unregulated river in Alberta, Canada, with ice jams frequently occurring in the vicinity of Fort McMurray. Modelling tools are desired to forecast ice‐related flood events. Multiple model combination methods can often obtain better predictive performances than any member models due to possible variance reduction of forecast errors or correction of biases. However, few applications of this method to river ice forecasting are reported. Thus, a framework of multiple model combination methods for maximum breakup water level (MBWL) Prediction during river ice breakup is proposed. Within the framework, the member models describe the relations between the MBWL (predicted variable) and their corresponding indicators (predictor variables); the combining models link the relations between the predicted MBWL by each member model and the observed MBWL. Especially, adaptive neuro‐fuzzy inference systems, artificial neural networks, and multiple linear regression are not only employed as member models but also as combining models. Simple average methods (SAM) are selected as the basic combining model due to simple calculations. In the SAM, an equal weight (1/n) is assigned to n member models. The historical breakup data of the Athabasca River at Fort McMurray for the past 36 years (1980 to 2015) are collected to facilitate the comparison of models. These models are examined using the leave‐one‐out cross validation and the holdout validation methods. A SAM, which is the average output from three optimal member models, is selected as the best model as it has the optimal validation performance (lowest average squared errors). In terms of lowest average squared errors, the SAM improves upon the optimal artificial neural networks, adaptive neuro‐fuzzy inference systems, and multiple linear regression member models by 21.95%, 30.97%, and 24.03%, respectively. This result sheds light on the effectiveness of combining different forecasting models when a scarce river ice data set is investigated. The indicators included in the SAM may indicate that the MBWL is affected by water flow conditions just after freeze‐up, overall freezing conditions during winter, and snowpack conditions before breakup.