Search

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.

Glacier mass balance models are needed at sites with scarce long-term observations to reconstruct past glacier mass balance and assess its sensitivity to future climate change. In this study, North American Regional Reanalysis (NARR) data were used to force a physically based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess its sensitivity to climate change. A 2-year record (2014–2016) from an on-glacier automatic weather station (AWS) and historical precipitation records from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the NARR to the station sites. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multitemporal digital elevation model dataset. The model showed a good performance against recent (2012–2016) direct glaciological mass balance observations as well as with cumulative geodetic mass balance estimates. The simulated mass balance was not very sensitive to the NARR spatial interpolation method, as long as station data were used for bias correction. The simulated mass balance was however sensitive to the biases in NARR precipitation and air temperature, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these two parameters with ancillary data. The glacier-wide simulated energy balance regime showed a large contribution (57 %) of turbulent (sensible and latent) heat fluxes to melting in summer, higher than typical mid-latitude glaciers in continental climates, which reflects the local humid “icefield weather” of the Columbia Icefield. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 and 7 ∘C and precipitation between −20 % and +20 %, which comprise the spread of ensemble Representative Concentration Pathway (RCP) climate scenarios for the mid (2041–2070) and late (2071–2100) 21st century. The climate sensitivity experiments showed that future changes in precipitation would have a small impact on glacier mass balance, while the temperature sensitivity increases with warming, from −0.65 to −0.93 m w.e. a−1 ∘C−1. The mass balance response to warming was driven by a positive albedo feedback (44 %), followed by direct atmospheric warming impacts (24 %), a positive air humidity feedback (22 %) and a positive precipitation phase feedback (10 %). Our study underlines the key role of albedo and air humidity in modulating the response of winter-accumulation type mountain glaciers and upland icefield-outlet glacier settings to climate.

Glacier mass balance models are needed at sites with scarce long-term observations to reconstruct past glacier mass balance and assess its sensitivity to future climate change. In this study North American Regional Reanalysis (NARR) data are used to force a physically-based, distributed glacier mass balance model of Saskatchewan Glacier for the historical period 1979–2016 and assess it sensitivity to climate change. A two-year record (2014–2016) from an on-glacier automatic weather station (AWS) and a homogenized historical precipitation record from nearby permanent weather stations were used to downscale air temperature, relative humidity, wind speed, incoming solar radiation and precipitation from the nearest NARR gridpoint to the glacier AWS site. The model was run with fixed (1979, 2010) and time-varying (dynamic) geometry using a multi-temporal digital elevation model (DEM) dataset. The model showed a good performance against recent (2012–2016) direct glaciological mass balance observations as well as with cumulative geodetic mass balance estimates. The simulated mass balance showed a large sensitivity to the biases in NARR precipitation and solar radiation, as well as to the prescribed precipitation lapse rate and ice aerodynamic roughness lengths, showing the importance of constraining these parameters with ancillary data. The difference between the static (1979) and dynamic simulations showed small differences (mean = 0.06 m w.e. a−1 or 1.5 m w.e. over 37 yrs), indicating minor effects of elevation changes on the glacier specific mass balance. The static mass balance sensitivity to climate was assessed for prescribed changes in regional mean air temperature between 0 to 7 °C and precipitation between −20 to +20 %, which comprise the spread of ensemble IPCC representative concentration pathways climate scenarios for the mid (2041–2070) and late (2071–2100) 21st century. The climate sensitivity experiments showed that future changes in precipitation would have a small impact on glacier mass-balance, while the temperature sensitivity increases with warming, from −0.65 to −0.93 m w.e. °C−1. Increased melting accounted for 90 % of the temperature sensitivity while precipitation phase feedbacks accounted for only 10 %. Roughly half of the melt response to warming was driven by a positive albedo feedback, in which glacier albedo decreases as the snow cover on the glacier thins and recedes earlier in response to warming, increasing net solar radiation fluxes. About one quarter of the melt response to warming was driven by latent heat energy gains (positive humidity feedback). Our study underlines the key role of albedo and air humidity in modulating the response of winter-accumulation type mountain glaciers and upland icefield-outlet glacier settings to climate.

Glacial mass balance estimated through the geodetic method requires glacial surface coordinate observations from historical and contemporary sources. Contemporary observations and historical topographic maps are typically referenced to separate horizontal and vertical datums and observed with different sampling intervals. This research demonstrates the sensitivity of glacial change detection to the datum considerations and sampling schemes through case studies of Andrei, Bridge and Peyto glaciers in Western Canada. To simulate the procedure of observing the glacial surfaces, profile lines were sampled from Digital Elevation Model (DEMs) on contour intervals for historical data and horizontal intervals for contemporary data. Profile lines from the following scenarios were compared: (1) different horizontal and vertical sampling schemes; (2) the horizontal datum was correctly reconciled but the vertical datum was not; (3) the vertical datum was correctly reconciled but the horizontal datum was not; (4) both the horizontal and vertical datums were correctly reconciled; and (5) both the horizontal and vertical datums were incorrectly reconciled. Vertical errors of up to 6.9 m, 6.0 m and 5.0 m were observed due to sampling effects and vertical errors of 22.2 m, 9.9 m and 55.0 m were observed due to datum inconsistencies on Bridge, Andrei and Peyto glacier respectively. Horizontal datum inconsistencies manifested as erratic levels of growth or downwasting along the glacial surface profile and vertical datum errors manifested as a consistent vertical offset. Datum inconsistencies were identified to contribute errors of up to 257.2 × 106 m3 (or 87%) and 54.6 × 106 m3 (or 580%) of estimated volume change below and above the equilibrium line respectively on Peyto Glacier. The results of this study provide an estimate of typical errors due to sampling constraints or datum inconsistencies as well as guidance for identifying where these error sources have contaminated mass balance results.

Ground-penetrating radar (GPR) surveys at a center frequency of 500 MHz were used to determine winter (2007/08) and net annual (2005–07) snow water equivalent (SWE) patterns across the upper parts of Belcher Glacier, Devon Ice Cap, Nunavut, Canada. The GPR measurements were validated against snow depths determined from avalanche probe measurements, and converted to SWE values using densities measured with a down-borehole neutron density probe and in shallow snow pits. Distinct internal reflection horizons (IRHs) in the GPR record were formed during warm summers in 2007 and 2005, and a large rain event in summer 2006 which caused ice to accumulate above the 2005 melt surface. Elevation provides the dominant control on winter SWE distribution across the basin, with surface topography (e.g. gullies) also being locally important. Based on the location where IRHs intersected the ice-cap surface, the basin-wide firn line occurred at an altitude of 1260–1300 m over the period 2005–08. Net mass balance across the accumulation area of Belcher Glacier averaged 0.24 m w.e. a−1 over the period 2005–07, mainly dependent on altitude. This is a little higher than most previous estimates for the period since the 1960s, although the differences lie within error limits.

Recognizing that climate influences both annual tree-ring growth and glacier mass balance, changes in the mass balance of Place Glacier, British Columbia, were documented from increment core records. Annually resolved ring-width (RW), maximum (MXD), and mean density (MD) chronologies were developed from Engelmann spruce and Douglas-fir trees sampled at sites within the surrounding region. A snowpack record dating to AD 1730 was reconstructed using a multivariate regression of spruce MD and fir RW chronologies. Spruce MXD and RW chronologies were used to reconstruct winter mass balance (Bw) for Place Glacier to AD 1585. Summer mass balance (Bs) was reconstructed using the RW chronology from spruce, and net balance was calculated from Bw and Bs. The reconstructions provide insight into the changes that snowpack and mass balance have undergone in the last 400 years, as well as identifying relationships to air temperature and circulation indices in southern British Columbia. These changes are consistent with other regional mass-balance reconstructions and indicate that the persistent weather systems characterizing large scale climate-forcing mechanisms play a significant glaciological role in this region. A comparison to dated moraine surfaces in the surrounding region substantiates that the mass-balance shifts recorded in the proxy data are evident in the response of glaciers throughout the region.

The influence of digital elevation model (DEM) resolution to modelled glacier melt during peak melt production was evaluated by performing a clear sky GIS radiation simulation over the Peyto Glacier in the Canadian Rockies. DEMs were generated at eight resolutions ranging from 1 m to 1000 m grid spacing from airborne lidar data. When applied to the planar area (PA) of the terrain, it was found that total melt increased with DEM resolution (r2 = 0·63) by 4% over 3 orders of magnitude. This systematic scaling-effect was mitigated at the basin scale, however, when the DEM slope variant area (SVA) was used to account for the increased divergence from PA as resolution increases. However, even after the inclusion of SVA in glacier surface melt simulations, localized melt variations with scale were still evident in the ablation and accumulation zone observations. In the ablation zone, there was a systematic increase in simulated melt (∼4%) as resolution decreased from 1 m to 1000 m (r2 = 0·89), with the opposite effect in the accumulation zone (r2 = 0·81). DEM resolution also affected the diurnal melt cycle, such that for the entire glacier there was a tendency for a morning over-estimation and afternoon underestimation of melt rate with decreasing resolution. For the accumulation zone, there was an increased melt rate at low resolutions occurring in the afternoon, while in the ablation zone there was a tendency for increasing melt rates with decreasing resolution throughout the day. These localized spatio-temporal variations in simulated melt are largely due to the lowering of ridges and raising of valley floors that occur as resolution decreases. This scale dependence in the representation of terrain morphology directly controls the pattern and relative proportion of direct beam shadowing over actively melting surfaces and thereby has a systematic influence on the grid cell-level hydrological balance. It is recommended that GIS-based glacier melt modelling routines take into account the slope area of grid cells, while noting that the choice of DEM scale can have a discernible and systematic influence on modelled runoff magnitude. It is important to note that while higher grid resolutions mitigate the effect of terrain smoothing on spatio-temporal melt patterns, lower resolutions actually mitigate the systematic error associated with assuming all surface areas are planar. Copyright © 2010 John Wiley & Sons, Ltd and Her Majesty the Queen in right of Canada.

Due to access difficulties in active alpine moraine environments, it can be challenging to accurately map and quantify debris cover and ice-core extent. To aid in identifying the presence and extent of ice-cored moraine, a non-invasive method of mapping spatial and temporal moraine temperature patterns using a light detection and ranging (lidar) digital elevation model (DEM) and sequences of oblique thermal imagery was evaluated. A procedure of lidar DEM-based orthorectification of thermal images collected through time from different locations enabled maps of temperature change to be generated and thermal signatures plotted. Although no exposed ice was visible on the moraine slope studied, the presence of shallow ice core beneath the debris-covered surface was inferred in areas of cooler temperatures during daylight solar heating and rapid thermal decay after sunset. It is presumed that this apparent increased heat loss in some areas of the moraine is being used to drive internal melt processes. It is bel...

The hydrological model WATFLOOD and a separate volume–area scaling relationship are applied to estimate glacier wastage and seasonal Melt contribution to the North and South Saskatchewan Rivers originating in the Canadian Rocky Mountains (1975–1998). Wastage is the ice melt volume that exceeds the volume of snow accumulation into the glacier system in a hydrological year, causing an annual net loss of glacier volume. Melt is the ice melt volume that is equal to, or less than, the volume of snow that accumulates into the glacier system in a hydrological year. By our definition then, glacier Melt is a storage term and does not contribute to increased total annual streamflow. Water is stored as snow on accumulation into the glacier system, and the water equivalent runoff is delayed until ice melts in the late summer months of the otherwise low streamflow. Wastage varied between basins with similar glacierized areas reflecting the individual response of glaciers to climate, contributing over 10% to July-to-September streamflow in some headwater basins, but under 3% annually to the regulated flow at Edmonton and Calgary. Melt was positively correlated with basin glacierized area and contributed over 27% to July-to-September flow from basins with greater than 1% glacierized area, and over double the wastage volume at Edmonton and Calgary. Future glacier decline is therefore expected to result mainly in an advancement of peak flow towards a non-glacierized snowmelt regime hydrograph, resulting in significantly reduced late summer flows further reduced by decreasing wastage contributions. Copyright © 2009 John Wiley & Sons, Ltd and Her Majesty the Queen in right of Canada.

(2009). Invited Commentary: A Framework for Integrated Research and Monitoring (FIRM) Canadian Water Resources Journal / Revue canadienne des ressources hydriques: Vol. 34, No. 1, pp. 1-6.

A large rock and ice avalanche occurred on the north face of Mount Steele, southwest Yukon Territory, Canada, on July 24, 2007. In the days and weeks preceding the landslide, several smaller avalanches initiated from the same slope. The ice and rock debris traveled a maximum horizontal distance 5.76 km with a maximum vertical descent of 2,160 m, leaving a deposit 3.66 km2 in area on Steele Glacier. The seismic magnitude estimated from long-period surface waves (M s) is 5.2. Modeling of the waveforms suggests an estimated duration of approximately 100 s and an average velocity of between 35 and 65 m/s. This landslide is one of 18 large rock avalanches known to have occurred since 1899 on slopes adjacent to glaciers in western Canada. We describe the setting, reconstruct the event chronology and present a preliminary characterization of the Mount Steele ice and rock avalanches based on field reconnaissance, analysis of seismic records and an airborne LiDAR survey. We also present the results of a successful dynamic simulation for the July 24 event.

Interpretation of ice mass elevation changes observed by satellite altimetry demands quantification of the proportion of elevation change which is attributable to variations in firn densification. Detailed stratigraphic logging of snowpack structure and density was carried out at ~1km intervals along a 47 km transect on Devon Ice Cap, Canada, in spring (pre-melt) and autumn (during/ after melt) 2004 and 2006 to characterize seasonal snowpack variability across the full range of snow facies. Simultaneous meteorological measurements were gathered. Spring (pre-melt) snowpacks show low variability over large spatial scales, with low-magnitude changes in density. The end-of-summer/ autumn density profiles show high variability in both 2004 and 2006, with vastly different melt regimes generating dissimilar patterns of ice-layer formation over the two melt seasons. Dye-tracing experiments from spring to autumn 2006 reveal that vertical and horizontal distribution of meltwater flow within and below the annual snowpack is strongly affected by the pre-existing, often subtle stratigraphic interfaces in the snowpack, rather than its bulk properties. Strong interannual variability suggests that using a simple relationship between air temperature, elevation and snowpack densification to derive mass change from measurements of elevation change across High Arctic ice caps may be misguided. Melt timing and duration are important extrinsic factors governing snowpack densification and ice-layer formation in summer, rather than averaged air temperatures.

Three ice cores were recovered on or near Mount Logan, Yukon, Canada, at 3017, 4135 and 5340 ma.s.l. in 2002. Prior to ice-core drilling, we collected snow-pit and shallow core samples from Mount Logan in 2001 to study seasonal and spatial variations of snow chemistry. We dug snow pits at six sites between 2420 and 5340 m a.s.l. before the beginning of the melt season, with the exception of a pit at 3180 m a.s.l., where the melt season had just started but had affected only the near-surface stratigraphy. Three of the pits were extended deeper with a shallow core. The snow-pit and core samples were analyzed for ion chemistry and δ18O. A series of depth profiles of ions and δ18O shows spatial variations, though characteristic peaks can usually be traced across all the profiles. Concentrations and deposition fluxes of Na+ and Cl−, which are mainly of sea-salt origin, decrease with altitude. On the other hand, deposition fluxes of NO3−, SO42–, Ca2+ and NH4+ show a weak positive relationship with elevation below the summit plateau. Stable isotopes (δ18O) decrease with altitude, with a distinctive jump between 3200 and 4500 m a.s.l., as was reported previously. Stable isotopes (δ18O), Cl−, CH3SO3− (MSA), Na+ and Ca2+ show clear seasonal variations, which would enable us to date the cores by annual-layer counting.

Knowledge of the changing dimensions of alpine glacier surfaces is critical from both a water resources and climate change indication perspective. With the development of airborne light detection and ranging (lidar) technologies with the capability to rapidly map large areas of topography at high resolutions, there is a need to assess the utility of this technology for glacier surface change detection and water resources assessment. The study presented here compares two lidar digital elevation models (DEMs) collected 23 months apart in 2000 and 2002 over Peyto Glacier, Canadian Rocky Mountains, for the purposes of intensity image feature recognition and surface downwasting assessment. The 2002 DEM was subtracted from the 2000 DEM to quantify the total and spatial variability in surface downwasting (or growth) within the glacial and periglacial environments. It was found that there was a reduction in volume totaling 33 × 106 m3 over the Peyto Glacier surface and surrounding ice-cored moraines. This downwas...

A new glacier inventory was developed for the Ragged Range, an area only recently added to Nahanni National Park Reserve, Northwest Territories. Federal mapping from aerial photography in 1982 was compared with Landsat satellite imagery from 2008. Glacier cover decreased in area by 30 % over that period; the greatest percentage loss was in glaciers with areas between 1 and 10 km2, while the largest number of disappearances was of glaciers 0.1–1 km2 in area. The smallest glaciers less than 0.1 km2 in area showed little loss to modest growth in topographic niches protected from solar radiation. Recommendations for further work in this environmentally significant frontier region are provided.

Although a great deal of research has focused on the hydrologic effects of climate variability and change, relatively little research has examined the effects on streamflow of interactions between climate variability and change and resulting glacier response. Place Glacier, in the southern Coast Mountains of British Columbia, Canada, has been monitored for mass balance since 1965, and a stream gauge was operated just below the glacier terminus from 1969 to 1989. This paper presents analyses of the mass balance history and streamflow variations in relation to recorded climatic variability. Place Glacier's winter and net balances are correlated with the Pacific Decadal Oscillation (PDO). Summer balance is positively correlated with summer temperature and negatively with the preceding winter balance, which enhances the effects of changes in winter balance on net balance. The well-documented post-1976 shift from the PDO cold phase to the present warm phase initiated a significant and persistent period of more negative net balance and terminal retreat. A reconstruction of net balance extending back to the 1890s, based on a regression with winter precipitation and summer temperature, displays decadal-scale fluctuations consistent with the PDO. Summer streamflow responded to interannual variations in winter snow accumulation and summer temperatures, which control the rate of rise of the glacier snowline and melt rates. After accounting for these influences via regression analysis, August streamflow displayed a negative trend in total runoff. Examination of air photographs and the reconstructed mass balance history suggest that significant firn depletion had occurred prior to 1965, such that the dominant effect of glacier changes was a reduction in ice area, resulting in decreased meltwater production. Copyright © 2001 John Wiley & Sons, Ltd.

Maps of glacier area in western Canada have recently been generated for 1985 and 2005 (Bolch et al., 2010), providing the first complete inventory of glacier cover in Alberta and British Columbia. Western Canada lost about 11% of its glacier area over this period, with area loss exceeding 20% on the eastern slopes of the Canadian Rockies. Glacier area is difficult to relate to glacier volume, which is the attribute of relevance to water resources and global sea level rise. We apply several possible volume-area scaling relations and glacier slope-thickness relations to estimate the volume of glacier ice in the headwater regions of rivers that spring from the eastern slopes of the Canadian Rocky Mountains, arriving at an estimate of 55 ± 15 km3. We cannot preclude higher values, because the available data indicate that large valley glaciers in the Rocky Mountains may be anomalously thick relative to what is typical in the global database that forms the basis for empirical volume-area scaling relations. Inco...

Discusses logistics and applications of the first complete airborne scanning LiDAR survey of a glacierised mountain basin. The survey was conducted over a 25/spl times/6 km area with ground height variation from 1800 m asl to 3400 m asl. The resolution of survey points on the ground varied from 1 per sq metre to around 1 per 4 sq metres, depending on elevation. Correspondence from swath to swath is very good despite the extreme nature of the topography and individual ground laser heights are accurate to within 20 cm. Glaciological applications for the high-resolution digital elevation data are outlined. Special attention is paid to radiation loading model-scaling issues.