Exploring the operational impacts of climate change and glacier loss in the upper Columbia River Basin, Canada

The mountainous watersheds of western Canada are generally thought to be in a state of transition from snow-dominated to hybrid regimes. In stream networks that are regulated, the effects of this transition on streamflow can have compelling operational consequences. Seasonal magnitude changes may impact spill-risk management, while changes in the composition of summer runoff may increase its variability and reduce the forecasting capabilities of state variables like peak snow water equivalent. Though glacier loss can have a considerable impact on summer runoff, few studies explicitly model the ongoing glacier recession in conjunction with other primary hydrological processes. In this study, we incorporate glacier dynamics from a previous run of the Regional Glaciation Model into the University of British Columbia Watershed Model via the Raven modelling framework. We use this modelling system to explore potential changes under Representative Concentration Pathways 4.5 and 8.5 to the hydrology of the similar to 20000 km(2) Mica Basin, a regulated watershed containing the headwaters of the Columbia River. Our results project statistically significant increases in spring flow in future eras, which may force lower reservoir drafting in late winter, creating potential for energy shortfalls in early spring. We project the coefficient of variation of summer runoff generally goes unchanged in future eras as does the summer runoff forecasting capability of April 1st SWE. Hence, despite modelled glacier loss and reduced snowmelt contribution, our study does not reject the null hypothesis that the predictability of the Mica Basin's summer runoff is unchanged in future eras. We explore these results in detail because they superficially appear to contrast the conventional conceptualization that reduced snowmelt negatively affects the predictive powers of snowpack and glacier loss increases the variability of runoff. We argue that our results' apparent discordance from convention displays the complexities inherent in isolating the effects of changes to a single water balance component when other components are also non-stationary and highlights the benefits of using modelling to more explicitly explore such implications.