High latitudes are predicted to continue warming at higher rates than the global average, with major implications for northern basins where concomitant deglaciation, permafrost thaw and vegetation shifts are expected. The Mackenzie River Basin, a globally significant basin, drains headwaters in the glaciated Canadian Rockies to the Arctic Ocean and is mostly underlain by permafrost. Here, we present scenarios of future change using the MESH distributed hydrological-cryospheric land surface model. MESH was forced with bias-corrected, downscaled RCM forcings and parameterized with a deep subsurface profile, organic soils, and glaciers. The model was validated against discharge, snowpack, and permafrost observations and used to simulate the hydrology and permafrost dynamics over the 21st century under the RCP8.5 climate change scenario with projected land cover change. The results show rapidly increasing rates of permafrost thaw; most of the basin will be permafrost-free by the 2080s. By late century, river discharges shift to earlier and higher peaks in response to projected increases in precipitation, temperature, snowmelt rates, despite increases in evapotranspiration from longer snow-free seasons. Baseflow discharges increase in winter, due to higher precipitation and increased basin connectivity from permafrost thaw resulting in enhanced groundwater flow. Subsurface moisture storage rises slightly but the liquid water fraction increases dramatically, increasing sub-surface runoff and river discharge. Canadian Rockies deglaciation reduces summer and annual discharge in the Athabasca and Peace headwaters. Downstream and northward of the mountain headwaters the direct impacts of climate change on river discharge dominate over those of changing land cover and glaciers.

Despite the proliferation of computer-based research on hydrology and water resources, such research is typically poorly reproducible. Published studies have low reproducibility due to incomplete availability of data and computer code, and a lack of documentation of workflow processes. This leads to a lack of transparency and efficiency because existing code can neither be quality controlled nor re-used. Given the commonalities between existing process-based hydrological models in terms of their required input data and preprocessing steps, open sharing of code can lead to large efficiency gains for the modeling community. Here we present a model configuration workflow that provides full reproducibility of the resulting model instantiations in a way that separates the model-agnostic preprocessing of specific datasets from the model-specific requirements that models impose on their input files. We use this workflow to create large-domain (global, continental) and local configurations of the Structure for Unifying Multiple Modeling Alternatives (SUMMA) hydrologic model connected to the mizuRoute routing model. These examples show how a relatively complex model setup over a large domain can be organized in a reproducible and structured way that has the potential to accelerate advances in hydrologic modeling for the community as a whole. We provide a tentative blueprint of how community modeling initiatives can be built on top of workflows such as this. We term our workflow the “Community Workflows to Advance Reproducibility in Hydrologic Modeling’‘ (CWARHM; pronounced “swarm”).

Cold regions provide water resources for half the global population yet face rapid change. Their hydrology is dominated by snow, ice and frozen soils, and climate warming is having profound effects. Hydrological models have a key role in predicting changing water resources, but are challenged in cold regions. Ground-based data to quantify meteorological forcing and constrain model parameterization are limited, while hydrological processes are complex, often controlled by phase change energetics. River flows are impacted by poorly quantified human activities. This paper reports scientific developments over the past decade of MESH, the Canadian community hydrological land surface scheme. New cold region process representation includes improved blowing snow transport and sublimation, lateral land-surface flow, prairie pothole storage dynamics, frozen ground infiltration and thermodynamics, and improved glacier modelling. New algorithms to represent water management include multi-stage reservoir operation. Parameterization has been supported by field observations and remotely sensed data; new methods for parameter identification have been used to evaluate model uncertainty and support regionalization. Additionally, MESH has been linked to broader decision-support frameworks, including river ice simulation and hydrological forecasting. The paper also reports various applications to the Saskatchewan and Mackenzie River basins in western Canada (0.4 and 1.8 million km2). These basins arise in glaciated mountain headwaters, are partly underlain by permafrost, and include remote and incompletely understood forested, wetland, agricultural and tundra ecoregions. This imposes extraordinary challenges to prediction, including the need to overcoming biases in forcing data sets, which can have disproportionate effects on the simulated hydrology.