Devon Ice Cap Research Articles

Variability in ice motion and dynamic discharge from Devon Ice Cap, Nunavut, Canada

ABSTRACTFeature tracking of approximately annually separated Landsat-7 ETM+ imagery acquired from 1999 to 2010 and speckle tracking of 24-day separated RADARSAT-2 imagery acquired from 2009 to 2015 reveal that motion of the major tidewater glaciers of Devon Ice Cap is more variable than previously described. The flow of almost half (six of 14) of the outlet glaciers slowed over the observation period, while that of the terminus regions of three of 14 of the glaciers sped up in the most recent years of observation. The North Croker Bay Glacier of southern Devon Ice Cap showed the greatest variability in motion, oscillating between multi-year (three or more) periods of slower and faster flow and exhibited a pattern of velocity variability that is different from that of the rest of the ice cap's outlet glaciers. Comparisons between areas of dynamic variability and glacier bed topography indicate that velocity variability is largely restricted to regions where the glacier bed is grounded below sea level. Derived velocities are combined with measurements of ice thickness at the fronts of tidewater glacier to determine a mean annual (2009; 2011–15) dynamic ice discharge of 0.41 ± 0.11 Gt a−1for Devon Ice Cap. The Belcher Glacier is becoming a larger source of mass loss via ice discharge.

Characterizing near‐surface firn using the scattered signal component of the glacier surface return from airborne radio‐echo sounding

AbstractWe derive the scattered component (hereafter referred to as the incoherent component) of glacier surface echoes from airborne radio‐echo sounding measurements over Devon Ice Cap, Arctic Canada, and compare the scattering distribution to firn stratigraphy observations from ground‐based radar data. Low scattering correlates to laterally homogeneous firn above 1800 m elevation containing thin, flat, and continuous ice layers and below 1200 m elevation where firn predominantly consists of ice. Increased scattering between elevations of 1200–1800 m corresponds to firn with inhomogeneous, undulating ice layers. No correlation was found to surface roughness and its theoretical incoherent backscattering values. This indicates that the scattering component is mainly influenced by the near‐surface firn stratigraphy, whereas surface roughness effects are minor. Our results suggest that analyzing the scattered signal component of glacier surface echoes is a promising approach to characterize the spatial heterogeneity of firn that is affected by melting and refreezing processes.

A tidewater glacier landform assemblage in Belcher Inlet, Canadian Arctic

Belcher Glacier, a 35 km long tidewater outlet glacier of the 12 000 km2 ice cap on Devon Island (Dowdeswell et al. 2004), is one of the fastest-flowing glaciers in the Canadian Arctic (Van Wychen et al. 2014) (Fig. 1). Belcher Glacier and neighbouring Fitzroy Glacier to the southeast account for about 55% of the iceberg calving loss from the Devon Ice Cap (Van Wychen et al. 2014). The terminus of Belcher Glacier remained relatively stable between the 1960s (light blue dashed line in Fig. 1a) and 2000 (Landsat 7 satellite image in Fig. 1a). In contrast, the unnamed glacier immediately to the north retreated 2 km during this period (Fig. 1a). Belcher Glacier and the unnamed glacier retreated around 500 and 250 m, respectively, between 2000 and 2014 (dark blue dashed line in Fig. 1a). The bed topography of Belcher Glacier, which is around 250 m below sea level at the present-day glacier margin (Fig. 1c) and remains below sea level in the lower 11 km of the glacier, suggests that its terminus region may become unstable in the event of future retreat. Seafloor mapping of Belcher Inlet beyond the termini of Belcher Glacier and the unnamed glacier (Fig. 1a), together with sub-bottom profiling, provide information about the dynamic behaviour of tidewater glaciers. Fig. 1. ( a ) Swath bathymetry of Belcher Inlet, Arctic Canada (acquired in 2006), superimposed on 2000 Landsat 7 satellite imagery. Multibeam acquisition systems Kongsberg EM300 and Kongsberg EM3002S. Frequency 30 and 300 kHz, respectively. Grid-cell size 10 m. Light blue and dark blue dashed lines are terminus positions derived from 1960s and 2014 aerial photographs, respectively. ( b ) Multibeam-bathymetric image of the seafloor immediately in front of Belcher Glacier. ( c ) 2014 airborne radar profile of Belcher Glacier, showing the depth of the …

CryoSat-2 delivers monthly and inter-annual surface elevation change for Arctic ice caps

Abstract. We show that the CryoSat-2 radar altimeter can provide useful estimates of surface elevation change on a variety of Arctic ice caps, on both monthly and yearly timescales. Changing conditions, however, can lead to a varying bias between the elevation estimated from the radar altimeter and the physical surface due to changes in the ratio of subsurface to surface backscatter. Under melting conditions the radar returns are predominantly from the surface so that if surface melt is extensive across the ice cap estimates of summer elevation loss can be made with the frequent coverage provided by CryoSat-2. For example, the average summer elevation decreases on the Barnes Ice Cap, Baffin Island, Canada were 2.05 ± 0.36 m (2011), 2.55 ± 0.32 m (2012), 1.38 ± 0.40 m (2013) and 1.44 ± 0.37 m (2014), losses which were not balanced by the winter snow accumulation. As winter-to-winter conditions were similar, the net elevation losses were 1.0 ± 0.20 m (winter 2010/11 to winter 2011/12), 1.39 ± 0.20 m (2011/12 to 2012/13) and 0.36 ± 0.20 m (2012/13 to 2013/14); for a total surface elevation loss of 2.75 ± 0.20 m over this 3-year period. In contrast, the uncertainty in height change from Devon Ice Cap, Canada, and Austfonna, Svalbard, can be up to twice as large because of the presence of firn and the possibility of a varying bias between the true surface and the detected elevation due to changing year-to-year conditions. Nevertheless, the surface elevation change estimates from CryoSat for both ice caps are consistent with field and meteorological measurements.

Linking surface hydrology to flow regimes and patterns of velocity variability on Devon Ice Cap, Nunavut

AbstractSupraglacial meltwater reaching a glacier bed can increase ice surface velocities via hydraulic jacking and enhanced basal sliding. However, the relationships between the structure of supraglacial drainage systems, sink-point distributions, glacier flow processes and the magnitude of interannual velocity variability are poorly understood. To explore the hypothesis that spatial variations in the rate and mechanisms of glacier flow are linked to variations in supraglacial drainage system structure and sink-point distribution across an ice cap, we mapped supraglacial drainage systems on Devon Ice Cap from Landsat-7 ETM+ imagery. Spatial patterns of surface velocity and interannual velocity variability were determined using gradient correlation applied to Landsat-7 ETM+ images. Velocity variability is greater in areas close to sink-point locations, presumably because hydrologically forced basal sliding and/or bed deformation are enhanced in such areas. The distribution and characteristics of supraglacial drainage systems may play an important role in determining the distribution of regions of basal sliding, highlighting the need for knowledge of the supraglacial drainage system structure and sink-point distribution to inform efforts to model the dynamic response of Arctic ice caps to future climate warming.

Field‐calibrated model of melt, refreezing, and runoff for polar ice caps: Application to Devon Ice Cap

AbstractUnderstanding the controls on the amount of surface meltwater that refreezes, rather than becoming runoff, over polar ice masses is necessary for modeling their surface mass balance and ultimately for predicting their future contributions to global sea level change. We present a modified version of a physically based model that includes an energy balance routine and explicit calculation of near‐surface meltwater refreezing capacity, to simulate the evolution of near‐surface density and temperature profiles across Devon Ice Cap in Arctic Canada. Uniquely, our model is initiated and calibrated using high spatial resolution measurements of snow and firn densities across almost the entire elevation range of the ice cap for the summer of 2004 and subsequently validated with the same type of measurements obtained during the very different meteorological conditions of summer 2006. The model captures the spatial variability across the transect in bulk snowpack properties although it slightly underestimates the flow of meltwater into the firn of previous years. The percentage of meltwater that becomes runoff is similar in both years; however, the spatial pattern of this melt‐runoff relationship is different in the 2 years. The model is found to be insensitive to variation in the depth of impermeable layers within the firn but is very sensitive to variation in air temperature, since the refreezing capacity of firn decreases with increasing temperature. We highlight that the sensitivity of the ice cap's surface mass balance to air temperature is itself dependent on air temperature.

How well is firn densification represented by a physically based multilayer model? Model evaluation for Devon Ice Cap, Nunavut, Canada

AbstractEvaluation of how accurately snowpack models can capture hydrological processes in firn is critical to determining how well they will simulate future glacier mass-balance changes. Here we compare simulations using the Crocus snowpack model with the evolving firn stratigraphy recorded in 14 cores drilled at four elevations in the accumulation zone of Devon Ice Cap, Canada, during the 2004–12 period of rapid summer warming. Simulations were forced with a combination of surface observations and reanalysis data. Simulations resulted in positive model bias in near-surface density, and negative bias in density at depth compared to observations. Results point to the importance of incorporating heterogeneous percolation in firn in order to improve the representation of meltwater flow, better reproduce observed firn density and temperature profile evolution, and improve simulations of glacier mass balance during periods of climate warming.

Interferometric swath processing of Cryosat data for glacial ice topography

Abstract. We have derived digital elevation models (DEMs) over the western part of the Devon Ice Cap in Nunavut, Canada, using "swath processing" of interferometric data collected by Cryosat between February 2011 and January 2012. With the standard ESA (European Space Agency) SARIn (synthetic aperture radar interferometry) level 2 (L2) data product, the interferometric mode is used to map the cross-track position and elevation of the "point-of-closest-approach" (POCA) in sloping glacial terrain. However, in this work we explore the extent to which the phase of the returns in the intermediate L1b product can also be used to map the heights of time-delayed footprints beyond the POCA. We show that there is a range of average cross-track slopes (~ 0.5 to ~ 2°) for which the returns will be dominated by those beneath the satellite in the main beam of the antenna so that the resulting interferometric phase allows mapping of heights in the delayed range window beyond the POCA. In this way a swath of elevation data is mapped, allowing the creation of DEMs from a sequence of L1b SARIn Cryosat data takes. Comparison of the Devon results with airborne scanning laser data showed a mean difference of order 1 m with a standard deviation of about 1 m. The limitations of swath processing, which generates almost 2 orders of magnitude more data than traditional radar altimetry, are explored through simulation, and the strengths and weaknesses of the technique are discussed.

Changes in accumulation-area firn stratigraphy and meltwater flow during a period of climate warming: Devon Ice Cap, Nunavut, Canada

[1] To document climate-driven changes in firn stratigraphy and their implications for meltwater flow patterns within firn on the Devon Ice Cap, Nunavut, during the 21st century summer warming, 500 MHz ground-penetrating radar (GPR) surveys were conducted along a 40 km transect in each spring from 2007 to 2012. These linear GPR surveys were supplemented by four 190 m by 100 m GPR grid surveys and 36 firn cores. Increased meltwater percolation and infiltration ice formation associated with high surface melt rates since 2005 modified the firn stratigraphy substantially over a horizontal distance of nearly 30 km. The most dramatic change involved the growth of a thick ice layer within the firn body. This layer grew primarily by upward accretion over an initial widespread ice layer formed during summer 2005. It thickened by between 0.5 and 4.5 m over the study period and filled much of the pore volume in the upper part of the firn, reducing vertical percolation of meltwater into deeper sections of the firn and thus the water storage potential of much of the firn reservoir. Heterogeneous percolation of surface meltwater promoted by rolling topography played an important role in meltwater infiltration and drainage, encouraging lateral flow at the tops of small hills and ponding and refreezing of meltwater beneath surface depressions.

Spatial patterns of snow accumulation across Belcher Glacier, Devon Ice Cap, Nunavut, Canada

AbstractGround-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.

Changes in melt season characteristics on Devon Ice Cap, Canada, and their association with the Arctic atmospheric circulation

AbstractUsing June–August surface meteorological measurements for 2004–10 we found that the duration of the summer melt season at elevations of 1800, 1300 and 1000 m a.s.l. on Devon Ice Cap, Canada, increased at a rate of 3.4 (p < 0.1; weak significance), 6.1 (p< 0.01) and 8.8 (p < 0.05) days a−1, respectively. The calculated surface melt rate at 1300m a.s.l. increased from 74 to 133 cm w.e. a−1 between 2007 and 2010. These changes are linked to two types of change in the Arctic atmospheric circulation. Strengthening of the 500 hPa ridge over the Arctic in June–July resulted in increases in both the advection of warm air into the region and the occurrence of cloud-free conditions over the ice cap, causing the available melt energy to increase by 4–24% relative to June–July 2007–10 mean conditions. More frequent southwesterly low-pressure systems in August after 2004 accounted for a 12–38% increase in available melt energy relative to the August 2007–10 daily mean due to advection of warm air into the Arctic, reduced incoming shortwave radiation and increased net longwave radiation. This delayed the timing of freeze-up by an average of 5.5 days a–1 (p < 0.05) at the three sites over 2004–10.

Firn profile changes in response to extreme 21st-century melting at Devon Ice Cap, Nunavut, Canada

AbstractEstimating regional-scale rates of glacier mass change from repeat airborne and satellite altimetry requires assumptions about the density of the material being added to, or removed from, the glacier surface. This may vary in both space and time, especially in a period of warming climate such as has occurred in Arctic Canada since 2005. Here we compare firn properties measured on 20 shallow cores recovered from elevations of 1400–1900 m a.s.l. on Devon Ice Cap in 2012 with equivalent measurements from 34 cores recovered from the same locations between 2004 and 2011. The average density of the top 2.5 m of the firn column has increased by 13–80%. As a result, measured rates of density increase in the top 1.25 mw.e. of the firn layer have likely thinned this layer at rates of 0.021–0.168 m a−1 over the period 2004–12. These thinning rates are significant relative to local rates of surface height change measured by repeat airborne laser altimetry between 2004 and 2012 (−0.077 to −0.252 m a−1) and indicate that, during periods of climate warming or cooling, firn densification rate changes must be taken into account when estimating mass change rates from measured height changes.

Atmospheric deposition of current use pesticides in the Arctic: Snow core records from the Devon Island Ice Cap, Nunavut, Canada

Current use pesticides (CUPs) have been detected in the Arctic, even though there are no direct sources and their long range atmospheric transport potential is generally lower than that of legacy pesticides. Data on the deposition of CUPs in the Arctic are required to assess the impact of their global usage and emission. In this study, selected CUPs were measured in the layers of a snow pit sampled on the Devon Ice Cap, Nunavut, Canada. The oldest sampled layers correspond to deposition from the early 1990s. Dacthal and endosulfan sulfate were most frequently detected, with peak deposition fluxes of 1.0 and 0.4 pg cm(-2) per year. While endosulfan sulfate was more abundant than its parent compounds in most years, endosulfan (sum of α and β isomers) was predominant in 2003 and 2006, which together with air mass backward trajectories suggests a possible origin from ongoing use in Eurasia. The interannual variation in CUP deposition fluxes could not be explained with annual variations in the extent of air mass origin over agricultural lands, suggesting that other factors, such as the interannual variation in pesticide use, play a role in affecting the long range transport of CUPs to the Arctic. The very high variability in the concentrations of CUPs in the horizontal layers of Arctic ice caps is most plausibly explained by the highly episodic nature of long range atmospheric transport and deposition. While this strong influence of rare events limits the suitability of ice caps as reliable records of historical trends in Arctic contaminant deposition with annual resolution, the presence of concentration peaks in the ice record is proof of the possibility of such transport and deposition.

Spatial and temporal variation of ice motion and ice flux from Devon Ice Cap, Nunavut, Canada

AbstractSpeckle tracking of repeat RADARSAT-2 fine-beam imagery acquired over 24 day periods in March 2009 allowed the creation of updated surface motion maps for the entire Devon Ice Cap, Canada. Error analysis indicates that speckle tracking can determine ice motion to an accuracy of ~5 ma-1. Comparisons with earlier velocity maps from the mid-1990s and 2000 reveal velocity patterns that largely agree with flow regimes described previously. However, motion determined along East5 Glacier indicates an increase in surface velocities between the studies. Additionally, Southeast2 Glacier has significantly accelerated over the past decade, with velocities greater in 2009 than in the early 1990s along almost the entire length of the glacier. This is likely indicative of a surge. Present-day total mass loss from Devon Ice Cap due to iceberg calving is calculated as 0.40 ± 0.09 Gta-1, similar to that reported by Burgess and others (2005), with Belcher Glacier accounting for ~42% of the entire loss.

Modelling the delivery of supraglacial meltwater to the ice/bed interface: application to southwest Devon Ice Cap, Nunavut, Canada

AbstractThe transfer of surface-generated meltwater to the subglacial drainage system through full ice thickness crevassing may lead to accelerated glacier velocities, with implications for ice motion under future climatic scenarios. Accurate predictions of where surface meltwater accesses the ice/bed interface are therefore needed in fully coupled hydrodynamic ice-sheet models. We present a spatially distributed modelling routine for predicting the location and timing of delivery of surface-derived meltwater to the ice/bed interface through moulins and supraglacial lake drainage. The model is explained as it is applied to the Croker Bay glacial catchment of Devon Ice Cap, Canada. The formation of moulins, drainage of lakes, and the transfer of meltwater through the full ice thickness are modelled for the 2004 and 2006 ablation seasons. Through this case study we assess the model’s sensitivity to degree-day factors, fracture toughness, tensile strength and crevasse width, and confirm that parameters influencing the rate at which water fills a crevasse are the most significant controls on the ability of a crevasse to reach the bed. Increased surface melt production, therefore, has the potential to significantly influence the spatial and temporal transfer of meltwater through surface-to-bed connections in a warmer climate

Relationships between iceberg plumes and sea-ice conditions on northeast Devon Ice Cap, Nunavut, Canada

AbstractThis study investigates the impact of sea-ice conditions on the production of iceberg plumes from two tidewater glaciers on Devon Ice Cap, Nunavut, Canada. These effects are quantified using a 12 year RADARSAT-1 satellite record from 1997–2008 that contains imagery from approximately every 1–2 weeks in the winter and every 1–4 days in the summer. Iceberg plumes identified in this record are verified against terrestrial time-lapse photography of Belcher Glacier from 2007–08. Results suggest a strong relationship between iceberg plumes and the retreat of sea ice from the glacier termini, with the plumes caused by both the release of previously calved icebergs (ice melange) and new glacier calving. Iceberg plumes are also sometimes observed at other times in the summer and in midwinter (occasionally on both glaciers simultaneously), with these events likely due to new glacier calving alone. Analysis of tides and air temperatures suggests that they provide a minor influence on the timing of iceberg plumes. Instead, it appears that changes in the presence of sea ice are dominant on seasonal timescales, although internal glacier dynamics likely play a significant role for winter plume events that occur when substantial thicknesses of landfast sea ice are present.

Deposition of Brominated Flame Retardants to the Devon Ice Cap, Nunavut, Canada

Brominated flame retardants (BFRs) can be transported to Arctic regions via atmospheric long-range transport, however, relatively little is known about their deposition to terrestrial environments. Snow cores from the Devon Ice Cap in Nunavut, Canada served to determine the recent depositional trends of BFRs. Snow pits were dug in 2005, 2006, and 2008. Dating using annual snow accumulation data, ion chemistry, and density measurements established that the pits covered the period from approximately 1993 to spring 2008. Samples were extracted under clean room conditions, and analyzed using GC-negative ion MS for 26 tri- to decabromodiphenyl ethers (BDEs), as well as other BFRs, nonbrominated flame retardants, and industrial chemicals. Decabromodiphenyl ether (BDE-209) was the major congener present in all samples followed by nona-BDEs (BDE-207, BDE-206, and BDE-208), both accounting for 89% and 7% of total BDE, respectively. BDE-209 concentrations were in most cases significantly correlated (P < 0.05) to tri- to nona-BDE homologues, and the strength of the correlations increased with increasing degree of bromination. Prior to or after deposition BDE-209 may be subject to debromination to lighter congeners. Deposition fluxes of BDE-209 show no clear temporal trend and range between 90 and 2000 pg·cm(-2)·year(-1). Back trajectory origin in densely populated areas of northeastern North America is significantly correlated (P < 0.005) with the BDE-209 deposition flux. Several other high production volume and/or alternative BFRs such as hexabromocyclododecane (HBCD), 1,2-bis(2,4,6-dibromophenoxy)ethane (BTBPE), pentabromo ethyl benzene (PBEBz), and pentabromobenzene (PBBz), as well as the industrial chemical 1,3,5-tribromobenzene (135-TBBz) were found consistently in the snow pits.

A Comparison of Recent Elevation Change Estimates of the Devon Ice Cap as Measured by the ICESat and EnviSAT Satellite Altimeters

We have used surface elevation measurements acquired by the Ice, Cloud,and land Elevation Satellite Geoscience Laser Altimeter System (GLAS) and EnviSAT Radar Altimeter 2 (RA-2) satellite altimeters to assess the elevation change of the 13 700-km2 Devon Ice Cap (DIC) in Arctic Canada between 2002 and 2008. We present algorithms for the retrieval of elevation change rates over ice caps using data acquired from these satellites. A comparison of GLAS elevation data to those acquired by the RA-2 shows reasonable agreement between the two instruments; the root mean square elevation change difference was 56 cm, and the correlation coefficient between the two data sets was 0.68. Using only RA-2 elevation measurements, which are spatially and temporally more continuous, we determined the elevation change rate of the areas of the DIC where the surface geometry allows the RA-2 retracker to maintain lock. This includes most of the DIC, excluding large parts of the eastern half of the ice cap. The elevation change rate was found to be insignificant given a statistical estimate of the measurement error (-0.09 ± 0.29 m/a). We also present an assessment of the regional variations of the DIC elevation change, including a significant -0.71 ± 0.49 m/a elevation change rate of the 1980-km2 western arm. Furthermore, we present evidence of a localized 2-m drop in the surface elevation of the South Croker Bay Glacier during summer 2007. This drop is apparent within both satellite data sets, and we interpret this signal to reflect a sudden speedup of the glacier.

Estimating the equilibrium line of Devon Ice Cap, Nunavut, from RADARSAT-1 ScanSAR wide imagery

C-band synthetic aperture radar (SAR) imagery from RADARSAT-1 is used to estimate the equilibrium line of Devon Ice Cap, Nunavut, over 3 years (2004–2006). The equilibrium line is determined from the interpretation of radar glacier zones and the transient snow line (TSL). Winter SAR data cannot be used to accurately estimate the equilibrium line because the firn line obscures the end-of-summer TSL, and accumulated superimposed ice cannot be identified. A methodology to estimate the equilibrium line from frequent summer SAR acquisitions is proposed. The TSL is tracked through the ablation season, and the equilibrium line is estimated as being equal to its maximum elevation. The summer SAR data are compared with near-surface air temperatures, QuikSCAT backscatter, optical imagery, and in situ mass balance measurements. The summer SAR-derived equilibrium line altitude (ELA) is always estimated as being above the in situ measured ELA because accumulated superimposed ice cannot be identified. Uncertainty in the equilibrium line estimation is reduced where slopes face the radar look direction. The summer SAR data provide a more accurate estimation of the equilibrium line than winter data, and the estimated equilibrium line can be used to infer relative changes in net balance between years.

Sensitivity of net mass-balance estimates to near-surface temperature lapse rates when employing the degree-day method to estimate glacier melt

AbstractGlacier mass-balance models that employ the degree-day method of melt modeling are most commonly driven by surface air temperatures that have been downscaled over the area of interest, using digital elevation models and assuming a constant free air lapse rate that is often taken to be the moist adiabatic lapse rate (MALR: –6.5°Ckm–1). Air-temperature lapse rates measured over melting glacier surface are, however, consistently less steep than free air values and have been shown to vary systematically with lower-tropospheric temperatures. In this study, the implications of including a variable near-surface lapse rate in a 26 year (1980–2006) degree-day model simulation of the surface mass balance of Devon Ice Cap, Nunavut, Canada, are examined and compared with estimates derived from surface air temperatures downscaled using a constant near-surface lapse rate equal to the measured summer mean (–4.9°Ckm–1) and the MALR. Our results show that degree-day models are highly sensitive to the choice of lapse rate. When compared with 23 years of surface mass-balance measurements from the northwest sector of the ice cap, model estimates are significantly better when surface air temperatures are downscaled using a modeled daily lapse rate rather than a constant lapse equal to either the summer mean or the MALR.