1. Introduction
Glaciers in the southwest Yukon and Alaska represent ∼14% of Earth’s glaciated area outside Greenland and Antarctica (ACIA, 2004), and like most mountain glaciers, they are generally thinning and retreating. Reference Arendt, Echelmeyer, Harrison, Lingle and ValentineArendt and others (2002) measured the average rate of thinning of glaciers in this area between the mid-1950s and mid-1990s to be ∼0.52 m a−1. Between the mid-1990s and 2001, the average rate of thinning increased to ∼1.8 m a−1. This change amounts to 96 ± 35 km3 a−1 of water loss, which raises global sea level by 0.27 ± 0.10 mm a−1, or ∼9% of the total observed rise over the past 50 years (Reference Arendt, Echelmeyer, Harrison, Lingle and ValentineArendt and others, 2002). Further altimeter-based measurements showed that losses in this region between the 1950s and 2002 may be greater by a factor of two (Reference ArendtArendt and others, 2006; Reference Larsen, Motyka, Arendt, Echelmeyer and GeisslerLarsen and others, 2007). These values are closely corroborated by Gravity Recovery and Climate Experiment (GRACE) data which indicate that glaciers in the Gulf of Alaska region lost mass at an average rate of 101 ± 22 km3 a−1 between 2002 and 2005 (Reference Chen, Tapley and WilsonChen and others, 2006). Using mascon (mass concentration) corrections for GRACE, Reference Luthcke, Arendt, Rowlands, McCarthy and LarsenLuthcke and others (2008) determined that glaciers in the Gulf of Alaska lost volume at a rate of 84 ± 5 km3 a−1 between 2003 and 2006. Reference Arendt, Luthcke, Larsen, Abdalati, Krabill and BeedleArendt and others (2008) narrowed GRACE mascon calculations to the St Elias Mountains and by combining these data with laser altimetry measurements, losses were calculated to be between 20.6 ± 3.0 and 21.1 ± 3.8 km3 a−1, a rate that had been approximately steady for several decades.
Recently, temperatures in the southwest Yukon have been increasing, with a warming of 2.0°C between 1948 and 2008 (Environment Canada, http://www.ec.gc.ca/adsccmda/default.asp?lang=en&n=4CC724DA-1). Warming is predicted to continue, with mean annual temperature in the southwest Yukon expected to increase by ∼3–3.5°C by 2099 (Reference SolomonSolomon and others, 2007). As a consequence of this temperature increase and a projected intensification of the Aleutian low, average annual precipitation in this region is expected to increase by up to 20% by 2099, with up to a 30% increase in winter (Reference SolomonSolomon and others, 2007). There is therefore a need for improved understanding of how climate change is affecting glaciers in the southwest Yukon, particularly as most existing studies focus on the Alaskan side of the St Elias Mountains. This study aims to address this need through calculations of changes in the volume of Kaskawulsh Glacier.
Kaskawulsh Glacier is a large, temperate valley glacier located within Kluane National Park in the St Elias Mountains (60°43′ N, 138°51′ W; Fig. 1). It flows generally northeast, is ∼70 km long and is not known to surge (although some tributaries do). The extensive Kluane Icefields, located to the west of the glacier, supply ice to the north and central arms of Kaskawulsh Glacier. The South Arm of the glacier is supplied by a catchment basin located to the south and southeast. The terminus is currently at ∼830 m a.s.l., and since 1836 it has been generally retreating, with increased wastage since 1980 (Reference Borns and GoldthwaitBorns and Goldthwait, 1966; Reference Wiles, Jacoby, Davi and McAllisterWiles and others, 2002; Reference Reyes, Luckman, Smith, Clague and van DorpReyes and others, 2006).
Kaskawulsh Glacier (base image: Landsat 2, 3 September 1977).
2. Methods
Changes of Kaskawulsh Glacier are examined via changes in surface height and extent, with this information combined to quantify changes in water equivalent volume. Elevation data were obtained for four periods between 1977 and 2007 from stereo aerial photography and airborne laser altimetry (Table 1). Glacier extent was obtained for nine periods between 1956 and 2007 from aerial photographs and satellite imagery (Table 1). All data were projected in Universal Transverse Mercator (UTM) zone 7N (NAD83 datum, GRS80 ellipsoid) in ESRI ArcGIS 9.2 for analysis.
Elevation and area data used in this study
2.1. Elevation data
The 1977 surface consisted of a Canadian Digital Elevation Data (CDED) digital elevation model (DEM) derived from stereo matching of aerial photographs collected in 1976 and 1977. Cumulative vertical errors in the CDED DEM are assumed to be ±15 m below the equilibrium-line altitude (ELA), and ±30–45 m above the ELA due to difficulties with matching surface features in snow-covered terrain (Reference EchelmeyerEchelmeyer and others, 1996; Reference Larsen, Motyka, Arendt, Echelmeyer and GeisslerLarsen and others, 2007). One area of anomalously steep surface slope was observed in this dataset in the accumulation area, and when this DEM is compared to all other elevation data a prominent depression is observed, the shape and size of which suggested that a low-level cloud may have been interpreted locally as the glacier surface. This region was removed from all analyses. Additionally, a tributary glacier on the central arm of Kaskawulsh Glacier surged between 1977 and 1994 (Fig. 1), which resulted in a rapid advance of its terminus by ∼2 km and a local gain in surface height. Portions of surface profiles affected by this surge event were removed from analysis as their local changes in surface height are controlled by surge dynamics rather than glacier mass balance.
Airborne laser altimetry (light detection and ranging; lidar) is currently one of the best methods to monitor the mass balance of glaciers in mountainous areas (Reference EchelmeyerEchelmeyer and others, 1996), as current technology allows small changes to be quantified over large featureless areas where derivation of topographic information from stereo imagery can be unreliable (Reference Hopkinson, Demuth, Sitar and ChasmerHopkinson and others, 2001; Reference Hopkinson and DemuthHopkinson and Demuth, 2006). Surface elevations for 1995 and 2000 were derived from flight-lines undertaken by the University of Alaska Fairbanks with a non-scanning lidar mounted beneath a Piper PA12 aircraft flying at ∼100–500 m above the glacier surface. This system only samples the points directly beneath the aircraft, with a surveyed area (laser footprint) of 0.18 m at a flight altitude of 100 m and measurements collected every ∼1.2 m along track (Reference EchelmeyerEchelmeyer and others, 1996). Positioning was provided with a differential GPS (dGPS). The range accuracy of this laser altimeter is ∼0.3 m (Reference EchelmeyerEchelmeyer and others, 1996; Reference Arendt, Echelmeyer, Harrison, Lingle and ValentineArendt and others, 2002), although errors increase with increasing slope.
Surface elevations in 2007 were determined using an airborne scanning laser altimeter operated by C-CLEAR (Canadian Consortium for Lidar Environmental Applications Research). The sensor used was an Optech Inc. airborne laser terrain mapper (ALTM) model 3100 flown on a Twin Otter at 400–1500 m above the glacier surface. The data were collected at a pulse repetition frequency of 33 kHz across a swath of ±25° from nadir, which resulted in cross-track and along-track point spacing of 2–10 m at the glacier surface. All ground points were referenced to the aircraft trajectory, which was determined from an inertial measurement unit and dGPS mounted inside the Twin Otter. The dGPS data were differentially corrected to a base station located at Kluane Lake Research Station (<30 km from the glacier terminus). The ALTM 3100 was calibrated in Calgary prior to and following the Kluane mission. Elevation root-meansquare error (RMSE) across the full 25° swath was <15 cm when flown at 1500 m a.g.l. For Kaskawulsh Glacier survey configuration, system-related errors are anticipated to be within 0.5 m vertically and 2 m horizontally (Reference Goulden and HopkinsonGoulden and Hopkinson, 2010). Coverage of the main glacier was achieved by surveying two adjacent flight-lines in an up-and down-glacier direction. To best match the coarsely spaced 1995/2000 lidar data, a DEM was created using a nearest-neighbour routine, and the data were resampled to 70 m grid spacing. One of the 2007 flight-lines deviated moderately from the 1995/2000 lidar lines over the glacier terminus (below 1353 m), and these non-overlapping regions could not be analyzed for change.
2.2. Areal data
Partial coverage of Kaskawulsh Glacier for 1956 was provided by mosaicking 11 aerial photographs, each encompassing an area of ∼15.2 × 15.2 km2. These photos covered the terminus and central arm of the glacier almost up to the equilibrium line, together with the entire South Arm. The images were acquired during the summer at flight altitudes of 3400–4900 m, and were scanned at 150 dpi, which yielded a resolution of 10 m.
In the recent period, Landsat images were acquired for 1977, 1986, 1990, 1994, 2001 and 2007 (Table 1). Images from 2007 have bands with missing data due to failure of the Landsat 7 scan-line corrector. Complete 2007 coverage was therefore provided by supplementing a 9 June 2007 image with a 12 August 2007 image in the accumulation area and a 28 August 2007 image over the terminus. Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER) scenes covering the entirety of Kaskawulsh Glacier were also used from 2001 and 2006, with the 2001 scene used to digitize the glacier terminus (except where there was cloud over the eastern terminus, where Landsat was used instead). Finally, a SPOT-5 image was acquired on 3 September 2007, via a collaborative agreement with SPIRIT (Spot 5 stereoscopic survey of Polar Ice: Reference Images and Topographies; Reference Korona, Berthier, Bernard, Rémy and ThouvenotKorona and others, 2009).
2.3. Image correction and data validation
Fieldwork was conducted in July 2008 to collect ground control points (GCPs) to georectify and verify the air photographs, satellite imagery and DEMs. These GCPs were collected on distinctive rock nunataks and medial moraines with a Trimble R7 dGPS. The data were processed using the online Precise Point Positioning (PPP) service provided by Natural Resources Canada (http://www.geod.nrcan.gc.ca/products-produits/ppp_e.php), which returned positions accurate to a few decimeters both horizontally and vertically. Alignment between the GCPs and 2007 SPOT-5 image was ±10 m (2 pixels) horizontally, and therefore the SPOT-5 scene was used as a base image to georectify the other scenes in locations away from the GCPs.
For georectification, 20–70 tie points were selected across each satellite image in ESRI ArcGIS 9.2. The central moraine GCPs were used for ∼6 tie points on average, while tie points on the SPOT-5 base image were selected based on recognizable features that were stationary in time and therefore not located on the ice or snow surface. An ‘adjust’ transformation was then applied, which optimizes local accuracy and the global least-square fitting by using polynomial transformation and triangular irregular network interpolation techniques depending on how well the image initially aligns with the GCPs (ESRI ArcGIS 9.2 Desktop Help). The rectified images were saved when visual investigation revealed pixel alignment of all major features within the area of interest and the RMSE was <10 m. For the air photographs, 35–80 points were selected for each photograph and a RMSE of <10 m was deemed acceptable. The registered aerial photographs were mosaicked together to provide maximum coverage across the glacier.
Comparison of the 1977 CDED DEM with the digitized extent of the georectified 1977 Landsat image reveals an estimated horizontal error of ±1 pixel (60 m). This DEM was therefore not adjusted after its initial creation. The laser altimetry-derived DEMs did not require any warping or adjustment due to the use of dGPS for positioning during data acquisition. The accuracy of the ALTM data (collected in August 2007) was determined by comparing the surface elevation of 50 points with ground measurements made with the Trimble R7 dGPS (in July 2008) along a relatively flat portion of the medial moraine (Fig. 1). The mean difference between these measurements was 2.0 m, with part of this likely due to unquantified interannual ice melt.
2.4. Determination of changes in glacier extent
The terminus position of Kaskawulsh Glacier was digitized for all area datasets, with measurements taken along 27 planes in the direction of ice flow to quantify its change in length over time. The difference in terminus area between images was also calculated. Most of the rest of the ablation area is bordered by steep valley side-walls that would show little change in glacier area even if there were substantial changes in ice thickness. Consequently, changes in glacier width in this region were deemed negligible and were excluded from area change analysis.
In the accumulation area, changes in extent were examined through changes in the area of five exposed rock nunataks. An effort was made to select rock complexes on each of the glacier arms to maximize representativeness of the findings across the accumulation zone. Furthermore, only images from summer and early fall were selected for analysis to reduce the impact of seasonal snowfall. A ratio between the area of the five rock complexes (7.93 km2) and the total exposed rock in the accumulation zone (244.20 km2) was calculated from the 1977 Landsat satellite imagery. This ratio was then used to produce an estimate of total change in rock area over time based on measurements at just the five rock complexes in later satellite images (Table 1). This method likely introduces errors in the area change estimates for the upper glacier, but still provides the best available areal data for this region in the absence of any ground control. Errors in the area change calculations can also arise from factors such as summer snowfall shortly before image acquisition, and from the image acquisition not coinciding with the minimum annual snow extent. Therefore, area changes for the accumulation area should be regarded as minimum estimates.
2.5. Determination of changes in glacier height and volume
The spatial distribution of absolute vertical surface-height changes along Kaskawulsh Glacier was derived by differencing the DEMs (Table 2). By combining these measurements with the area changes, the change in water equivalent ice volume for the entire glacier (Δv all) over a given time interval was calculated:
where a is the area of remaining glacier ice, Δh is the change in surface height, c is the change in area of the glacier and d is near-surface density. For Δh and c, losses are indicated with negative values, and gains with positive values. The subscript ‘acc’ represents values from the accumulation zone, while ‘abl’ represents values from the ablation zone. In certain cases, the change in height of the accumulation area is positive while the change in area is negative (or vice versa). In these cases the volume change accounted for by the change in area (c accΔh acc) is negligible and a value of 0 is used.
Area, height and volume changes of Kaskawulsh Glacier, 1977–2007
The ELA divides the accumulation and ablation zones and was determined from visual classification of the end-ofsummer snowline (ESS) from the 3 September 2007 SPOT image. This ESS occurred at 1958 m, the highest over the 30 year period, although there was little variation observed in the earlier Landsat imagery (Table 1). Therefore an ELA of 1958 m is used throughout this study.
Densities in the accumulation zone were determined from snow samples taken at 10 cm increments in an 8.1 m deep snow pit dug at an elevation of 2606 m near the summit of Kaskawulsh Glacier in July 2007 (60°40′ 44″ N, 139°47′ 44″ W; Fig. 1). The mean density in this pit was 0.51 g cm −3 (d acc), with little variability with depth (std dev. 0.11 g cm−3). This depth is representative of the order of thickness changes observed, and was therefore used as a mean value for the accumulation area. Glacier ice dominates the ablation zone, so a density of 0.90 g cm−3 (dabl) was used for this area.
The change in volume of the glacier was calculated for the time intervals over which elevation data are available (1977–95, 1995–2000, 2000–07 and 1977–2007). These intervals are not exactly the same intervals over which data on changes in area are available (Table 1), so a linear rate of area change was assumed between known dates to adjust the area calculations.
3. Results
3.1. Changes in surface height by individual period for 1977–2007
Changes in surface height along the centre line of the glacier over the period 1977–95 showed maximum thinning near the terminus (Table 2; Figs 2 and 3). Thinning is dominant throughout the ablation zone, although some areas of thickening (up to 1 m a−1) do occur. For the entire length of the glacier the mean thinning was 1.0 m over the 18 year period, an average rate of 0.06 m a−1. Within the ablation zone the average change was a thinning of 2.0 m, while the mean change in the accumulation zone was a thickening of 1.5 m.
Annual thinning/thickening along the central profile (shown in Fig. 1) for the periods 1977–95, 1995–2000, 2000–07 and 1977–2007.
Elevation profile along the centre line of Kaskawulsh Glacier (purple line in Fig. 1) for 1977, 1995, 2000 and 2007 for (a) terminus, (b) lower ablation zone, (c) upper ablation zone and (d) accumulation zone. Error bars shown for 1977 (±15 m in ablation zone, ±30 m in accumulation zone); error bars too small to be shown for 1995 (±0.3 m), 2000 (±0.3 m) or 2007 (±0.5 m).
Between 1995 and 2000, clear thinning of the glacier continued, with maximum losses of 31.1 m (6.2 m a−1) occurring near the terminus (Figs 2 and 3). There were a few small, isolated areas that appear to have thickened, although the localized nature of these changes suggests that they could be due to measurement error. For the ablation zone the mean change was a thinning of 6.3 m, whereas in the accumulation zone there was a mean thinning of 2.3 m. Overall, for the entire glacier length between 1995 and 2000 the mean thinning was 6.0 m, or 1.2 m a−1.
Between 2000 and 2007, thinning continued at a reduced rate (Figs 2 and 3). The terminus underwent an average thinning of 9.2 m (1.3 m a−1), with maximum thinning in excess of 29 m (4.1 m a−1).In the accumulation zone, average thinning of 2.7 m was observed. The average thinning for the glacier centre line as a whole was 4.8 m, a rate of 0.7 m a−1.
3.2. Overall changes in surface height for 1977–2007
Due to the fact that the 1977 DEM covers the entire Kaskawulsh Glacier and the 2007 DEM was created using swath mapping, a greater area of the glacier can be compared for thickness change by comparing these datasets than by plotting centre-line changes for the individual periods. Using this increased coverage, the ablation zone thinned by an average of 10.8 m between 1977 and 2007 (0.36 m a−1; Table 2; Figs 2 and 3). Thinning was particularly extensive across the glacier terminus, with maximum losses (up to 88 m) near the glacier centre. The accumulation zone underwent modest overall thinning of 1.1 m, although some areas did experience localized thickening. Overall the mean thinning of Kaskawulsh Glacier between 1977 and 2007 is 6.1 m, which equates to 0.20 ma−1.
When the thickness changes from the individual periods 1977–95, 1995–2000 and 2000–07 are added together, the accumulation area is shown to have thinned by an average of 3.44 m, and the ablation area by an average of 17.53 m (Table 2; Figs 2 and 3). This produces an overall thinning of the glacier of 13.80 m (0.46 m a−1). These values differ somewhat from the 0.20 m a−1 thinning measured directly from 1977 to 2007, because narrow centre-line profiles neglect differential thinning across the glacier width. Because the changes in thickness at the glacier centre are assumed to be constant across the glacier width for the periods 1977–95, 1995–2000 and 2000–07, between 1 and 3 m a−1 of vertical error at the glacier margins may have been introduced, with errors likely greatest at low elevations (Reference EchelmeyerEchelmeyer and others, 1996; Reference Arendt, Luthcke, Larsen, Abdalati, Krabill and BeedleArendt and others, 2008). Reference Barrand, James and MurrayBarrand and others (2010) discuss this issue in detail, and conclude that significant differences in mass-balance estimates may occur between centre-line and full-glacier surveys due to the high complexity of surface height changes within elevation bands. The error introduced by ignoring this non-uniformity cannot be quantified in this study due to data only being available from centre-line surveys in 1995 and 2000, although Figure 4 illustrates the variability in ice loss that occurs across the glacier width when comparing the 1977 and 2007 surveys directly. To address this limitation we calculated ice losses using both methods, and use these to bracket the likely range of mass losses over the study period.
(a) Location of sampled cross sections; (b) cross section at 1340 m (a–a1); (c) cross section at 1850 m (b–b1); (d) cross section at 2330 m (c–c1).
3.3. Changes in area
Figure 5 shows the changes in terminus position of Kaskawulsh Glacier between 1956 and 2007, with the details provided in Table 3. Over this period the glacier retreated by an average of 655 m, with a decrease in terminus area of 8.20 km2. Over the entire 51 year period, there was only one brief interval of readvance, between 1986 and 1990, when the terminus moved forward by 62 m, resulting in an increase in terminus area of 0.54 km2. The greatest rate of retreat occurred in 2006/07, when the glacier receded 79 m a−1, resulting in a terminus area loss of 0.87 km2. This was accommodated largely by the expansion of terminal lakes on the eastern and western lobes. Over the period 1977–2007, glacier retreat amounted to 471 m, accounting for an area loss of 5.46 km2.
Temporal variations in terminus position of Kaskawulsh Glacier displayed over a 1956 air photograph.
Change in terminus position of Kaskawulsh Glacier, 1956–2007
The rate of area change in the accumulation zone has varied over time (Fig. 6). By extrapolating the changes in the five representative nunatak complexes across the entire accumulation zone, an estimated increase in exposed rock of 40.05 km2 occurred between 1977 and 2001, while the area of exposed rock decreased by 28.47 km2 between 2001 and 2007. Thus, between 1977 and 2007 there was an overall ice loss of 11.58 km2, which equates to an ice area decrease of 0.39 km2 a−1 (0.05% a−1). Variability in snow accumulation and melt patterns across the accumulation zone of the glacier appears to be relatively low as fluctuations are relatively uniform across all nunatak complexes.
Area changes over time of nunataks across the accumulation area of Kaskawulsh Glacier. NA represents exposed rock complexes on the northern arm of the glacier, CA the central arm and SA the southern arm.
In 1977, the area of the glacier as a whole was 1111.81 km2. In 2001, the area had decreased to 1069.24 km2. By 2007, the glacier had increased in area by 25.55 km2, to 1094.79 km2, largely due to increases in the accumulation area. Therefore, over the 30 year period, the glacier decreased in area by 17.02 km2, or 1.53%.
3.4. Changes in volume
The total reduction in water equivalent volume of Kaskawulsh Glacier, calculated by summing the centre-line data from the separate measurement intervals and using the appropriate densities from the accumulation and ablation zones, is a decrease of 5.94 km3 w.e. (−0.20 km3 w.e. a−1) between 1977 and 2007 (Table 2). Losses in volume occurred prominently in the ablation zone (−4.53 km3 w.e.), and more moderately in the accumulation zone (−1.42 km3 w.e.).
Volume change was also calculated directly from the overall area changes between 1977 and 2007, and the overall surface height changes between the 1977 DEM and the full 2.3 km wide swath DEM from 2007. These data show that between 1977 and 2007 the glacier lost a total of 3.27 km3 w.e., at an overall rate of −0.11 km3 w.e. a−1.There were small volume losses in the accumulation zone (−0.46 km3 w.e.), and moderate volume losses in the ablation zone (−2.81 km3 w.e.). Any inaccuracy in estimation of the ELA will affect calculated volume losses as different density values are associated with the ablation and accumulation zones, although the effect of this error will be minimal given the total volume change. Changes in volume cannot be directly interpreted as changes in mass balance, however, as at least some of the volume change could be attributed to internal accumulation, particularly in the accumulation zone. Although difficult to quantify, internal accumulation can account for 5–100% of annual net accumulation (Reference Kaser, Cogley, Dyurgerov, Meier and OhmuraKaser and others, 2006).
4. Discussion and Conclusion
Our results clearly indicate that Kaskawulsh Glacier has lost mass over the period 1977–2007, with average overall thinning of 0.20–0.46 m a−1 and a decrease in overall volume of 3.27–5.94 km3 w.e. These changes have largely occurred as a function of changes in surface height, rather than area. These changes are also variable with altitude; negative changes have been prominent in the ablation zone (average thinning 0.4–0.5 m w.e.a−1), but minimal in the accumulation zone (average thinning 0.04–0.11 m a−1) (Fig. 7). At 2300 m a.s.l., net change in surface height is zero, with net gains above this altitude over the period 1977–2007 (although the deviation about this trend is large). This is consistent with surface height changes on other glaciers in this region, which have demonstrated slight to no thinning at higher elevations, and increasingly prominent thinning towards the glacier terminus (Reference ArendtArendt and others, 2006; Reference MolniaMolnia 2007). For example, Reference ArendtArendt and others (2006) found that at least 23 glaciers in the western Chugach Mountains displayed this pattern of surface height change between the 1950s and the early 2000s.
Relationship between changes in surface height and altitude for 1977–2007 along the swath profile shown in inset map.
The glacier reductions reported in this paper and by Reference Barrand and SharpBarrand and Sharp (2010) are in step with observed decadal warming trends in interior Alaska and northwestern Canada (Reference Stafford, Wendler and CurtisStafford and others, 2000; Environment Canada, http://www.ec.gc.ca/adsc-cmda/default.asp?lang=en&n=4CC724DA-1). The impact of such massive glacier changes is particularly significant for watersheds such as that of the Yukon River, where perennial ice and snow cover, while only accounting for ∼1% of the total catchment area, play a major role in regulating discharge year-round. Already, glacier losses have altered the hydrological regime of rivers from Alberta to Alaska (Reference DemuthDemuth and others, 2008; Reference JanowiczJanowicz, 2008; Reference MooreMoore and others, 2009). What remains uncertain, however, is whether the recent warming is accompanied by increases in accumulation on the glacerized regions of Alaska and the Yukon.
There exist no long-term (>10 year) instrumental records of precipitation from the St Elias or Wrangell Mountains. However, reconstructions of historical snowfall have been developed from winter balance (accumulation) measurements on southeastern Alaskan glaciers, and from ice cores drilled on or near Mount Logan (60°35′ N, 140°30′ W; 5300 m a.s.l.) in the central St Elias Mountains (Fig. 8). The winter balance of the maritime Wolverine and Gulkana Glaciers, southeast Alaska, is influenced by the Pacific Decadal Oscillation (PDO; Hartmann and Wedler 2005), and Wolverine Glacier experienced an increase in winter accumulation after the 1976 PDO modal shift (Reference Hodge, Trabant, Krimmel, Heinrichs, March and JosbergerHodge and others, 1998). However, since 1989 both glaciers have had sustained negative net balance trends, with summer melt losses surpassing winter nourishment (Reference Josberger, Bidlake, March and KennedyJosberger and others, 2007). In the southern Yukon, Reference Moore, Holdsworth and AlversonMoore and others (2002) reported a positive trend in snow accumulation of 1.2 cm w.e. a−1 over the period 1976–91 using a composite record of cores from Northwest Col (NWC) on Mount Logan. However, a separate core, drilled in 2001/02 from nearby Prospector-Russell Col (PRC; 5340 m a.s.l.; Reference FisherFisher and others, 2008), also on Mount Logan, shows no such trend over the same period (Fig. 8). Meanwhile, cores drilled in 1996 and 2002 from the Eclipse Icefield (60°51′ N, 139°47′ W; 3017 m a.s.l.), 45 km northeast of Kaskawulsh Glacier’s summit area (Reference Wake, Yalcin and GundestrupWake and others, 2002), show post-1976 trends of 2.5 and 1.1 cm a−1, respectively, but neither trend is statistically significant (p = 0.05). The trend in the stacked 1996+2002 Eclipse record is only 0.04 cm a−1. Elsewhere, Yasunari and others (2007) reported on a 10 year ice-core record from Mount Wrangell, Alaska, (62° N, 144° W; 4100 m a.s.l.) which appears to lack any clear accumulation trend. Finally, a >500 year long record of accumulation was developed from a core drilled on Mount Bona–Churchill, Alaska (61°24′ N, 42°00′ W; 4420 m a.s.l.) (Reference UrmannUrmann, 2009), but its significance is yet to be tested.
Historical variations in snow accumulation, 1950–2000, reconstructed from ice cores drilled in the Mount Logan area of the St Elias Mountains, and from winter balance measurements on two southeastern Alaska glaciers (Wolverine and Gulkana). Mount Logan ice cores: NWC = Northwest Col; PRC = Prospector-Russell Col. Thin curves denote annual variations, bold curves are 5 year running averages. The data were normalized to ease comparison. The PDO index for 1950–2000 is shown for comparison (data from Joint Institute for the Study of the Atmosphere and Ocean, University of Washington). The vertical line highlights the 1976 PDO modal shift. Of all the accumulation series shown, only two (Logan NWC and Wolverine Glacier) have post-1976 trends that are statistically significant (p = 0.05; trends shown in parentheses).
In summary, evidence for increasing snowfall in the central St Elias Mountains over the period of interest (post-1976) remains ambiguous. In the absence of such evidence, it must be anticipated that the mass balance of glaciers in the region will continue to become increasingly negative in the foreseeable future, with melt-induced losses far exceeding snowfall nourishment. Regional impact assessments on hydrological consequences should take this into consideration.
Acknowledgements
We thank the Canadian Foundation for Innovation (CFI), Ontario Research Fund (ORF), Natural Sciences and Engineering Research Council of Canada (NSERC), Northern Scientific Training Program (NSTP), University of Ottawa and Parks Canada for financial and logistical support. We also thank the Arctic Institute of North America, C. Wong and S. Pope for contributions to fieldwork. SPOT-5 data access was provided by the SPIRIT project. Special thanks to K. Echelmeyer, A. Arendt and the University of Alaska Fairbanks for supplying data. Research students from the Applied Geomatics Research Group are acknowledged for assisting with lidar data collection activities. Lidar data collection was facilitated by the Canadian Consortium for Lidar Environmental Applications Research (C-CLEAR) and a funding grant to M.N.D. and C.H. from the Canadian Space Agency. This is Natural Resources Canada Earth Sciences Sector (ESS) contribution 20100139.
