Abstract Reliable, long-term records of glacier mass change are invaluable to the glaciological and climate-change communities and used to assess the importance of glacier wastage on streamflow. Here we evaluate the in-situ observations of glacier mass change for Place (1982–2020) and Peyto glaciers (1983–2020) in western Canada. We use geodetic mass balance to calibrate a physically-based mass-balance model coupled with an ice dynamics routine. We find large discrepancies between the glaciological and geodetic records for the periods 1987–1993 (Place) and 2001–2006 (Peyto). Over the period of observations, the exclusion of ice dynamics in the model increased simulated cumulative mass change by ~10.6 (24%) and 7.1 (21%) m w.e. for Place and Peyto glacier, respectively. Cumulative mass loss using geodetic, modelled and glaciological approaches are respectively − 30.5 ± 4.5, − 32.0 ± 3.6, − 29.7 ± 3.6 m w.e. for Peyto Glacier (1982–2017) and − 45.9 ± 5.2, − 43.1 ± 3.1, − 38.4 ± 5.1 m w.e. for Place Glacier (1981–2019). Based on discrepancies noted in the mass-balance records for certain decades (e.g. 1990s), we caution the community if these data are to be used for hydrological model development.

Abstract. This paper presents hydrometeorological, glaciological and geospatial data from the Peyto Glacier Research Basin (PGRB) in the Canadian Rockies. Peyto Glacier has been of interest to glaciological and hydrological researchers since the 1960s, when it was chosen as one of five glacier basins in Canada for the study of mass and water balance during the International Hydrological Decade (IHD, 1965–1974). Intensive studies of the glacier and observations of the glacier mass balance continued after the IHD, when the initial seasonal meteorological stations were discontinued, then restarted as continuous stations in the late 1980s. The corresponding hydrometric observations were discontinued in 1977 and restarted in 2013. Datasets presented in this paper include high-resolution, co-registered digital elevation models (DEMs) derived from original air photos and lidar surveys; hourly off-glacier meteorological data recorded from 1987 to the present; precipitation data from the nearby Bow Summit weather station; and long-term hydrological and glaciological model forcing datasets derived from bias-corrected reanalysis products. These data are crucial for studying climate change and variability in the basin and understanding the hydrological responses of the basin to both glacier and climate change. The comprehensive dataset for the PGRB is a valuable and exceptionally long-standing testament to the impacts of climate change on the cryosphere in the high-mountain environment. The dataset is publicly available from Federated Research Data Repository at https://doi.org/10.20383/101.0259 (Pradhananga et al., 2020).

Snowmelt contributions to streamflow in mid-latitude mountain basins typically dominate other runoff sources on annual and seasonal timescales. Future increases in temperature and changes in precipitation will affect both snow accumulation and seasonal runoff timing and magnitude, but the underlying and fundamental roles of mountain basin geometry and hypsometry on snowmelt sensitivity have received little attention. To investigate the role of basin geometry in snowmelt sensitivity, a linear snow accumulation model and the Cold Regions Hydrological Modeling (CRHM) platform driven are used to estimate how hypsometry affects basin-wide snow volumes and snowmelt runoff. Area-elevation distributions for fifty basins in western Canada were extracted, normalized according to their elevation statistics, and classified into three clusters that represent top-heavy, middle, and bottom-heavy basins. Prescribed changes in air temperature alter both the snow accumulation gradient and the total snowmelt energy, leading to snowpack volume reductions (10–40%), earlier melt onsets (1–4 weeks) and end of melt season (3 weeks), increases in early spring melt rates and reductions in seasonal areal melt rates (up to 50%). Basin hypsometry controls the magnitude of the basin response. The most sensitive basins are bottom-heavy, and have a greater proportion of their area at low elevations. The least sensitive basins are top-heavy, and have a greater proportion of their area at high elevations. Basins with similar proportional areas at high and low elevations fall in between the others in terms of sensitivity and other metrics. This work provides context for anticipating the impacts of ongoing hydrological change due to climate change, and provides guidance for both monitoring networks and distributed modeling efforts.

Abstract. The interaction of mountain terrain with meteorological processes causes substantial temporal and spatial variability in snow accumulation and ablation. Processes impacted by complex terrain include large-scale orographic enhancement of snowfall, small-scale processes such as gravitational and wind-induced transport of snow, and variability in the radiative balance such as through terrain shadowing. In this study, a multi-scale modelling approach is proposed to simulate the temporal and spatial evolution of high-mountain snowpacks. The multi-scale approach combines atmospheric data from a numerical weather prediction system at the kilometre scale with process-based downscaling techniques to drive the Canadian Hydrological Model (CHM) at spatial resolutions allowing for explicit snow redistribution modelling. CHM permits a variable spatial resolution by using the efficient terrain representation by unstructured triangular meshes. The model simulates processes such as radiation shadowing and irradiance to slopes, blowing-snow transport (saltation and suspension) and sublimation, avalanching, forest canopy interception and sublimation, and snowpack melt. Short-term, kilometre-scale atmospheric forecasts from Environment and Climate Change Canada's Global Environmental Multiscale Model through its High Resolution Deterministic Prediction System (HRDPS) drive CHM and are downscaled to the unstructured mesh scale. In particular, a new wind-downscaling strategy uses pre-computed wind fields from a mass-conserving wind model at 50 m resolution to perturb the mesoscale HRDPS wind and to account for the influence of topographic features on wind direction and speed. HRDPS-CHM was applied to simulate snow conditions down to 50 m resolution during winter 2017/2018 in a domain around the Kananaskis Valley (∼1000 km2) in the Canadian Rockies. Simulations were evaluated using high-resolution airborne light detection and ranging (lidar) snow depth data and snow persistence indexes derived from remotely sensed imagery. Results included model falsifications and showed that both wind-induced and gravitational snow redistribution need to be simulated to capture the snowpack variability and the evolution of snow depth and persistence with elevation across the region. Accumulation of windblown snow on leeward slopes and associated snow cover persistence were underestimated in a CHM simulation driven by wind fields that did not capture lee-side flow recirculation and associated wind speed decreases. A terrain-based metric helped to identify these lee-side areas and improved the wind field and the associated snow redistribution. An overestimation of snow redistribution from windward to leeward slopes and subsequent avalanching was still found. The results of this study highlight the need for further improvements of snowdrift-permitting models for large-scale applications, in particular the representation of subgrid topographic effects on snow transport.

Abstract. This paper presents hydrometeorological, glaciological and geospatial data of the Peyto Glacier Research Basin (PGRB) in the Canadian Rockies. Peyto Glacier has been of interest to glaciological and hydrological researchers since the 1960s, when it was chosen as one of five glacier basins in Canada for the study of mass and water balance during the International Hydrological Decade (IHD, 1965–1974). Intensive studies of the glacier and observations of the glacier mass balance continued after the IHD, when the initial seasonal meteorological stations were discontinued, then restarted as continuous stations in the late 1980s. The corresponding hydrometric observations were discontinued in 1977 and restarted in 2013. Data sets presented in this paper include: high resolution, co-registered DEMs derived from original air photos and LiDAR surveys; hourly off-glacier meteorological data recorded from 1987 to present; precipitation data from nearby Bow Summit; and long-term hydrological and glaciological model forcing datasets derived from bias-corrected reanalysis products. These data are crucial for studying climate change and variability in the basin, and to understanding the hydrological responses of the basin to both glacier and climate change. The comprehensive data set for the PGRB is a valuable and exceptionally long-standing testament to the impacts of climate change on the cryosphere in the high mountain environment. The dataset is publicly available from Federated Research Data Repository at https://doi.org/10.20383/101.0259 (Pradhananga et al., 2020).

Abstract. The interaction of mountain terrain with meteorological processes causes substantial temporal and spatial variability in snow accumulation and ablation. Processes impacted by complex terrain include large-scale orographic enhancement of snowfall, small-scale processes such as gravitational and wind-induced transport of snow, and variability in the radiative balance such as through terrain shadowing. In this study, a multi-scale modeling approach is proposed to simulate the temporal and spatial evolution of high mountain snowpacks using the Canadian Hydrological Model (CHM), a multi-scale, spatially distributed modelling framework. CHM permits a variable spatial resolution by using the efficient terrain representation by unstructured triangular meshes. The model simulates processes such as radiation shadowing and irradiance to slopes, blowing snow redistribution and sublimation, avalanching, forest canopy interception and sublimation and snowpack melt. Short-term, km-scale atmospheric forecasts from Environment and Climate Change Canada's Global Environmental Multiscale Model through its High Resolution Deterministic Prediction System (HRDPS) drive CHM, and were downscaled to the unstructured mesh scale using process-based procedures. In particular, a new wind downscaling strategy combines meso-scale HRDPS outputs and micro-scale pre-computed wind fields to allow for blowing snow calculations. HRDPS-CHM was applied to simulate snow conditions down to 50-m resolution during winter 2017/2018 in a domain around the Kananaskis Valley (~1000 km2) in the Canadian Rockies. Simulations were evaluated using high-resolution airborne Light Detection and Ranging (LiDAR) snow depth data and snow persistence indexes derived from remotely sensed imagery. Results included model falsifications and showed that both blowing snow and gravitational snow redistribution need to be simulated to capture the snowpack variability and the evolution of snow depth and persistence with elevation across the region. Accumulation of wind-blown snow on leeward slopes and associated snow-cover persistence were underestimated in a CHM simulation driven by wind fields that did not capture leeside flow recirculation and associated wind speed decreases. A terrain-based metric helped to identify these lee-side areas and improved the wind field and the associated snow redistribution. An overestimation of snow redistribution from windward to leeward slopes and subsequent avalanching was still found. The results of this study highlight the need for further improvements of snowdrift-permitting models for large-scale applications, in particular the representation of subgrid topographic effects on snow transport.

Western North American (WNA) glaciers outside of Alaska cover 14,384 km2 of mountainous terrain. No comprehensive analysis of recent mass change exists for this region. We generated over 15,000 multisensor digital elevation models from spaceborne optical imagery to provide an assessment of mass change for WNA over the period 2000–2018. These glaciers lost 117 ± 42 gigatons (Gt) of mass, which accounts for up to 0.32 ± 0.11 mm of sea level rise over the full period of study. We observe a fourfold increase in mass loss rates between 2000–2009 [−2.9 ± 3.1 Gt yr−1] and 2009–2018 [−12.3 ± 4.6 Gt yr−1], and we attribute this change to a shift in regional meteorological conditions driven by the location and strength of upper level zonal wind. Our results document decadal‐scale climate variability over WNA that will likely modulate glacier mass change in the future.

Abstract Horizontal and altitudinal redistribution of snow by wind transport and avalanches can be important controls on small- and large-scale snow accumulation patterns that control meltwater supply in alpine environments. Redistribution processes control the spatial variability of snow accumulation, which not only controls meltwater supply, but also regulates snowmelt timing, duration, and rates, as well as snow-covered area depletion and the variable contributing area for meltwater runoff generation. However, most hydrological models and land surface schemes do not consider snow redistribution processes, and those that do are difficult to verify without spatially distributed snow depth measurements. These are rarely available in both high resolution and covering large scales. As an increased number of hydrological models include snow redistribution processes there is a need for additional snowcover metrics to verify snow redistribution schemes over large areas using readily available data. This study develops novel high-resolution (20 m), snowcover indices from remotely sensed imagery (Landsat-8 and Sentinel-2) to evaluate snow redistribution models over alpine areas without in-situ or airborne snow observations. A snowcover absence (SA) index, calculated from snow-free areas during the winter, identifies areas of wind erosion or avalanche source areas. A snowcover persistence (SP) index, calculated from snow-covered areas during the summer, identifies snow deposition in drifts and avalanche deposits. The snowcover indices captured the relative differences in surface observations of snow presence and absence between exposed and sheltered sites on an intensely instrumented ridge in the Canadian Rockies Hydrological Observatory. Within the Tuolumne River Basin in central California (1100 km2), the SP index captured roughly half of the spatial variability (R2 = 0.49 to 0.56) in peak SWE as estimated from airborne LiDAR-derived snow depths. At the individual mountain ridge scale (~800 m), variability in both ablation and snow redistribution controlled the SP patterns over 7979 ridges. Differences in shortwave irradiance explained 76% of the SP variance across ridges, but could not explain smaller-scale (~100 m) SP peaks that are associated with snowdrifts and avalanche deposits. The snowcover indices can be used to evaluate snow redistribution models of the finer scale impacts of snow redistribution by wind and gravity as long as the larger scale influences of spatially variable solar irradiance effects are also simulated.