Annals of Glaciology (55)66 2014 doi: 10.3189/2014AoG66A102

25

Initial field observations on Qaanaaq ice cap, northwestern
Greenland
Shin SUGIYAMA,1* Daiki SAKAKIBARA,1,2 Satoshi MATSUNO,1,2
Satoru YAMAGUCHI,3 Sumito MATOBA,1 Teruo AOKI4
1

Institute of Low Temperature Science, Hokkaido University, Sapporo, Japan
E-mail: sugishin@lowtem.hokudai.ac.jp
2
Graduate School of Environmental Science, Hokkaido University, Sapporo, Japan
3
Snow and Ice Research Center, National Research Institute for Earth Science and Disaster Prevention, Suyoshi-machi,
Nagaoka, Japan
4
Climate Research Department, Meteorological Research Institute, Nagamine, Tsukuba, Japan
ABSTRACT. To study the glaciological processes controlling the mass budget of Greenland’s peripheral
glaciers and ice caps, field measurements were carried out on Qaanaaq ice cap, a 20 km long ice cap innorthwestern Greenland. In the summer of 2012, we measured surface melt rate, ice flow velocity andice thickness along a survey route spanning the ice margin (200 m a.s.l.) to the ice-cap summit(1110 m a.s.l.). Melt rates in the ablation area were clearly influenced by dark materials covering theice surface, where degree-day factors varied from 5.44 mm w.e. K–1 d–1 on a clean surface to8.26 mm w.e. K–1 d–1 in the dark regions. Ice velocity showed diurnal variations, indicating the presenceof surface-meltwater induced basal sliding. Mean ice thickness along the survey route was 120 m, with amaximum thickness of 165 m. Ice velocity and temperature fields were computed using athermomechanically coupled numerical glacier model. Modelled ice temperature, obtained by imposing
estimated annual mean air temperature as the surface boundary condition, was substantially lower than
implied by the observed ice velocity. This result suggests that the ice dynamics and thermodynamics of
the ice cap are significantly influenced by heat transfer from meltwater and changing ice geometry.
KEYWORDS: Arctic glaciology, glacier flow, glacier modelling, ice cap, surface melt

1. INTRODUCTION
Coastal regions of Greenland contain a number of glaciersand ice caps (GICs), which are physically separated fromthe Greenland ice sheet. The total area of these ice bodies(1.3  105 km2) accounts for 7% of the entire ice cover
over Greenland (Rastner and others, 2012), so their volume
change has substantial influence on sea-level rise. Because
the GICs are located at relatively lower elevations, they aresusceptible to the recent strong warming trend in Greenland (e.g. Hanna and others, 2012). Recent estimates of
total ice mass loss of Greenland’s GICs are 28–41 Gt a–1,
depending on the definition of the GICs in terms of thedegree of connection with the ice sheet (Bolch andothers, 2013).
Despite the importance and the large mass loss revealedby satellite remote-sensing techniques, only a few fieldobservations have been reported from Greenland’s GICs.
One of the most detailed field datasets is available for
Mittivakkat Gletscher, southeast Greenland. The massbalance record of this glacier shows increasing mass loss
since 1995 (Mernild and others, 2011), and the mass loss in
2009/10 was more than twice the mean over the previous
15 years. Importantly, recent changes in the GICs are notuniformly distributed over Greenland. In northeast Greenland, the surface elevation change over the Flade Isblink
ice cap from 2002 to 2009 was near zero (Rinne andothers, 2011), suggesting spatial variability in the massbalance trend. Thus, observational data from different
regions are needed to better understand the drivers of thechanging mass balance of the GICs as well as the massbalance variations over Greenland.

A likely reason for the recent mass loss of the GICs is rising
air temperature, but other processes (e.g. surface albedo andice velocity) may be playing a crucial role in the changing
mass balance. Reduction in the ice and snow surface albedo
is one such potentially important process. Over the ablation
zone of the western Greenland ice sheet, a large area of theice surface is progressively darkening, and the area of thedark region is expanding (e.g. Wientjes and Oerlemans,
2010). Similar changes in the ice surface appearance are also
expected on peripheral GICs. It is believed that melt
increases in the dark regions because of low albedo, resulting
in more negative mass balance of the ice sheet and the GICs.
The effect of a darkening ice surface on the mass balance has
been studied and confirmed in other regions (e.g. Oerlemans
and others, 2009), but few in situ data exist in Greenland.
Another important process is the influence of meltwater on
the ice dynamics. The magnitude of surface melt and its
temporal variability have the potential to control the flow
velocity of the predominantly cold ice sheet (Zwally andothers, 2002; Schoof, 2010; Colgan and others, 2012), whichindirectly affects the mass budget by changing the rate of icetransport to lower elevations. Seasonal flow velocity changes
are reported in the Flade Isblink ice cap (Palmer and others,2010), but field data on short-term velocity variations arelacking in Greenland’s GICs.
In this paper, we report melt rates, ice velocity and ice
thickness of Qaanaaq ice cap (QIC) as the results of theinitial glaciological investigation on this ice cap in northwestern Greenland. We also present preliminary results ofnumerical experiments performed with a glacier flow model
based on the ice thickness data.

26

Sugiyama and others: Field observations on Qaanaaq ice cap

2. FIELD MEASUREMENTS
2.1. Qaanaaq ice cap
Qaanaaq ice cap (77828’ N, 69814’ W) is located in northwestern Greenland, covering an area of 289 km2 and anelevation range of 30–1110 m a.s.l. over the central part of apeninsula (Fig. 1). The peninsula is populated by about 600
people in the village of Qaanaaq located on the southerncoast. Climatic data since 1996 are available from regional
weather stations, but no glaciological observations havebeen reported on the ice cap. A number of outlet glaciers
protrude from the ice cap for distances of 1–10 km. Our
study site in the 2012 summer field campaign was Qaanaaq
Gletscher, which flows southwest from the southern part ofthe ice cap (Fig. 1c). We performed field observations along
a survey route, which follows a glacier flowline from one of
the highest peaks of the ice cap to the terminus of
Qaanaaq Gletscher.

2.2. Melt rate
We installed seven 2 m long aluminum poles along thesurvey route to measure ice/snow melt rates (Fig. 1c). The
installation sites ranged from 243 to 1079 m a.s.l. (Table 1).
When the poles were installed on 18 and 20 July 2012, ice
was exposed at the lower five sites (Q1201–Q1205),
whereas the upper two sites (Q1206 and Q1207) werecovered by a snow layer. The snow cover was wet andclassified as granular snow (Fierz and others, 2009). The
poles were drilled into the surface using a mechanical ice
drill (Kovacs) equipped with an electric drill driver (Makita
HR262DRDX). The distance from the top of the pole to theice/snow surface was measured with an accuracy of 5 mm.
We repeated the measurements with intervals of 1–9 days
from the dates of the installation until 29 July. Water
equivalent melt rates were computed assuming ice densitiesof  = 910 kg m–3 at Q1201–Q1205, and firn densities of = 500 kg m–3 at Q1206 and Q1207. Because firn density
was not measured on site, it was estimated from within arange suggested for granular snow (Paterson, 1994), and its
uncertainty is about 20%.

2.3. Ice velocity

Fig. 1. (a) Satellite image (Landsat, 5 February 1999) showing
northwestern Greenland including Qaanaaq. The inset shows thelocation of the region in Greenland, and the box indicates the area
shown in (b). (b) Satellite image of QIC (ALOS PRISM, 25 August2009). The locations of Qaanaaq Airport and the temperature
measurement site on the ice cap (SIGMA-B) are indicated byasterisks. The box indicates the area covered by (c). (c) Satellite
image (ALOS PRISM, 25 August 2009) of the study site, showing thelocations of the measurement sites for surface melt (+), ice velocity(*), ice radar (*) and the GPS reference station (*). The arrows arehorizontal surface flow vectors from 18 to 29 July 2012.

The poles were also used for ice velocity measurements.
The three-dimensional (3-D) coordinates of the poles weresurveyed twice, on the date of installation and on 29 July.
We mounted the antenna of a dual-frequency GPS receiver(GNSS Technologies Inc., GEM-1) on the top of the poles torecord GPS satellite signals every 1 s for at least 30 min.
These GPS data were post-processed in static mode withthose collected by a reference GPS station establishedwithin the village of Qaanaaq (Fig. 1c). The same device
and recording intervals were used for the reference station.
The coordinates were computed with the GPS processing
software RTKLIB (http://www.rtklib.com/rtklib.htm). According to the estimate of the software, the standard
deviations of the GPS positioning are generally <5 mm in
the horizontal direction and <20 mm in the verticaldirection. Vertical velocity was not obtained at Q1207
because the GPS solution was not accurate enough due topoor satellite configuration.
The survey pole at Q1203 was equipped with apermanent GPS station from 18 to 26 July. The satellite
signals were recorded continuously with the same type ofGPS receiver powered by a battery (13 Ah) and a solar panel

Sugiyama and others: Field observations on Qaanaaq ice cap

27

Table 1. Surface melt rates, horizontal and vertical ice velocities measured at sites Q1201–Q1207. Elevations, measurement periods andsurface conditions based on in situ visual observations are listed for the sites

Elevation (m a.s.l.)Dates (July 2012)
Surface condition
Melt rate (mm w.e. d–1)
Horizontal velocity (mm d–1)
Vertical velocity (mm d–1)

Q1201

Q1202

Q1203

Q1204

Q1205

Q1206

Q1207

243
18–29
ice/clean
46
4.0
0.9

427
18–29
ice/clean
38
34.0
–0.7

584
18–29
ice/dark
39
67.5
–1.7

739
18–29
ice/dark
38
69.0
–4.2

839
18–29
ice/dark
30
44.8
–3.6

968
18–29
snow
14
18.6
–4.4

1079
20–29
snow
7
1.1
–

(20 W). The GPS data were split into 30 min segments tocompute 3-D displacements. Each of the 30 min data
segments was processed in a static mode with data recordedat the reference station as described above. Ice velocity wascomputed after filtering the half-hourly positioning datausing a Gaussian smoothing routine with a time window of
3 hours.

2.4. Ice thickness and bed elevation
Ice thickness was measured approximately every 200 m
along the survey route (Fig. 1c). We used an ice-penetratingradar system manufactured by The Ohio State University,
USA, which consists of a transmitter, receiver and 20 m long
dipole antennas. The central wave frequency was 5 MHz.
Vertical resolution expected for ice thickness measurementswith this frequency is 8.4 m. The transmitter and receiver
were placed 20 m apart along the survey route, and theantenna cables were stretched perpendicular to the route. To
compute ice thickness from two-way travel time of bedreflected radar signals, wave propagation speed through theglacier was assumed to be 168 m ms–1. We took this valuethat has been proposed for the speed in solid ice (Glen andParen, 1975), because a surface firn layer was observed onlyin the uppermost regions (>970 m a.s.l.). Glacier thickness
may have been slightly underestimated in the firn-covered
region because of this assumption.
Ice-cap surface elevation was surveyed by a kinematic
GPS positioning technique at the sites of the radarsounding. We used a dual-frequency GPS receiver (Leica
System 1200) as a rover station and post-processed the datawith those from the aforementioned reference station. The
ice thickness measured by the radar sounding was thensubtracted from the surface elevation to obtain ice-cap
bed elevation.

3. RESULTS OF FIELD MEASUREMENTS
3.1. Melt rate
Melt rates measured at Q1201–Q1207 ranged from 7 mm
w.e. d–1 at the ice-cap summit Q1207 to 46 mm w.e. d–1 atthe lowermost site Q1201 (Table 1). The mean rate of theseven sites was 30 mm w.e. d–1. The rates were greater
down-glacier in general, but the relationship to elevation is
not simple (Fig. 2). The melt rate is relatively less sensitive toelevation from about 400 to 700 m a.s.l. (e.g. same rates at
Q1202 and Q1204), and more sensitive to it from 700 to1000 m a.s.l. (e.g. –10.5 mm w.e. d–1 (100 m)–1 between
Q1203 and Q1205) as compared to the mean gradient for
the entire elevation range (–4.7 mm w.e. d–1 (100 m)–1 between Q1201 and Q1207).

3.2. Ice velocity
Horizontal and vertical components of the ice motion aresummarized in Table 1, and flow vectors are shown in
Figure 1c. Ice flows faster in the middle part of the survey
route, where ice is expected to be thick. The greatest
horizontal velocities were observed at Q1203 (67.5 mm d–1)
and Q1204 (69.0 mm d–1), and from that region the velocity
progressively decreases down- and up-glacier. Flow velocity
was only 4.0 mm d–1 at the lowermost site Q1201 500 m
from the terminus. Vertical ice motion was downward except
for Q1201, where ice emerged at a rate of 0.9 mm d–1.

2.5. Air temperature data
Air temperature recorded by a weather station at Qaanaaq
Airport (778280 N, 698140 W; z = 16 m a.s.l.) located at
3.5 km northwest of Qaanaaq (Fig. 1b) was used in thisstudy. Daily mean, minimum and maximum temperatures atthis station since 1996 are distributed by the US National
Oceanic and Atmospheric Administration. Daily mean
temperature from 2010 to 2013 was used to study temperature conditions and estimate surface melt.
We also measured air temperature on the ice cap at
1 km to the east of the survey route (778310 N, 698040 W;
z = 944 m a.s.l.) (SIGMA-B in Fig. 1b). Hourly temperature
measurement has been in operation since July 2012 using asensor (Vaisala, HMP155) installed at 3.0 m above thesurface. This record was used to compute the temperature
lapse rate over the ice cap.

Fig. 2. Mean melt rate at Q1201–Q1207 measured over the periods
indicated in Table 1. The dashed line is the extrapolation of the line
connecting the data at Q1201 and Q1202.

28

Sugiyama and others: Field observations on Qaanaaq ice cap

Fig. 3. Horizontal ice velocity variations at Q1203 and daily mean temperature with its variation range measured at Qaanaaq Airport.

The continuous GPS measurement at Q1203 has shown
significant ice velocity variations within the range 36–
94 mm d–1, equivalent to 50–140% of the mean velocity(Fig. 3). The variations were diurnal, with minimum andmaximum velocities in early morning and late afternoon,
respectively. More detailed discussion on the timings of thevelocity peaks is difficult as the GPS positioning is notaccurate enough. The diurnal variations were clearlyrecognized until 23 July, but became more irregular for the
rest of the period. The daily mean air temperature and its
range of variation at Qaanaaq Airport were relatively smaller
during the latter half of the observation period.

3.3. Ice thickness and bed elevation
Ice thickness, bed and surface elevations were obtained by
the ice radar and kinematic GPS measurements along thesurvey route (Fig. 4). Mean ice thickness along the route was
120 m. Ice was thickest (165 m) at 3.2 and 4.4 km from theinterior of the ice cap (Fig. 4a), where the bed deepens
locally (Fig. 4b). Relatively large bedrock bumps anddepressions (amplitude of several tens of metres) are
observed in the middle reaches of the ice cap. These

irregular bed geometries affect the surface elevation, as
represented by a small surface bulge above a bedrock bumpat 5.9 km from the summit.

4. DEGREE-DAY ANALYSISThe temperature record at Qaanaaq Airport (Fig. 5a) was
utilized to estimate total summer melt at our measurementsites by the degree-day method. First, daily mean temperature at the measurement sites over summer 2012 was
calculated by assuming a temperature lapse rate of–7.80  10–3 K m–1 (results for Q1201 and Q1207 are shownin Fig. 5b). This lapse rate was calculated from temperaturesat Qaanaaq Airport and SIGMA-B during the period of ourmelt measurement (20–29 July 2012). The temperature series
were then used to compute the positive degree-days (PDD)
over the melt season (Fig. 5b). Daily melt expected for one
degree of positive daily mean temperature (1 K d), so-called
degree-day factor k, was obtained at each elevation of themeasurement sites askðzÞ ¼

MðzÞ
MðzÞ
¼P
,
PDDðzÞ
Tm 0 Tm ðzÞ tð1Þ

where M(z) and Tm(z) are observed melt during the fieldstudy period and estimated mean temperature over a time
period t = 1 day, respectively, at elevation z. Total meltover the 2012 summer melt season was estimated by
Mt ðzÞ ¼ kðzÞ  PDDt ðzÞ,ð2Þ

where suffix t denotes the total values for the whole summer
melt season. The computations were performed for theelevations of the measurement sites, and the results are
summarized in Table 2 and Figure 6. The degree-day factor
at the lower two sites is 5.4–5.5 mm w.e. K–1 d–1, which are
lower bounds of the previously reported values for an ice
surface (Hock, 2003). The factors are greater at Q1203–
Q1205, with a peak value of 8.3 mm w.e. K–1 d–1 at Q1204(739 m a.s.l.). The degree-day factor is smaller at the upper
two sites (Q1206 and Q1207), where the ice-cap surface
was covered with a snow layer.

5. NUMERICAL MODELLING
5.1. Model
Fig. 4. (a) Ice thickness and (b) surface and bed elevations along thesurvey route shown in Figure 1c. Ice thicknesses and bed elevations
after the ice radar measurements are indicated by (*).

To investigate ice flow and temperature conditions of QIC, athermomechanically coupled finite-element glacier model
was applied to the longitudinal cross section of the ice cap
along the survey route. The aim of the modelling was toSugiyama and others: Field observations on Qaanaaq ice cap

29Fig. 5. (a) Daily mean temperature measured at Qaanaaq Airport for 2010 (blue), 2011 (green) and 2012 (red). (b) Daily mean temperatures(solid lines) and PDD (dashed line) computed for Q1201 (red) andcompare the measured ice velocities with those diagnosed
from the present geometry and climate conditions. Also ofinterest was basal temperature, which is relevant to thediurnal ice velocity variations observed at Q1203.
The model is based on the previously developed twodimensional ice flow model (Sugiyama and others, 2003),
which is coupled with a thermodynamic model in this study.
Taking x- and z-coordinates in the along-flow and verticaldirections, momentum and mass conservation equations in
the x–z plane are@xx @xz
þ
¼0
@x
@z

ð3Þ

@zx @zz
þ
¼ g
@x
@z

ð4Þ

and
@ux @uz
þ
¼ 0,
ð5Þ
@x
@z
where ij (i, j = x, z) are components of the Cauchy stresstensor,  = 910 kg m–3 is the density of ice, g = 9.81 m s–2 isthe vertical component of the gravitational acceleration
vector and ux and uz are the horizontal and verticalcomponents of the velocity vector. We link the stresses inFig. 6. (a) PDD (&), degree-day factors (*) and (b) total melt amountat Q1201–Q1207 computed for the entire summer melt season in
2012. Blue line in (b) is the brightness intensity of the ALOS PRISM
image (25 August 2009; see Fig. 1) along the survey route.

Table 2. PDD over the study period and for the year 2012 (PDDt), positive degree-day factors (k) and total melt estimated for the summer
melt season in 2012 (Mt) at sites Q1201–Q1207

PDD (K d)
PDDt (K d)k (mm w.e. K–1 d–1)
Mt (mm w.e.)

Q1201

Q1202

Q1203

Q1204

Q1205

Q1206

Q1207

92.7
425.0
5.45
2315

76.9
301.3
5.44
1640

63.4
214.0
6.74
1443

50.1
143.5
8.26
1185

41.6
107.9
7.88
851

30.5
70.6
4.97
351

16.3
45.7
3.69
169

30

Sugiyama and others: Field observations on Qaanaaq ice cap

Table 3. Parameter values used in the numerical model. Values arefrom Greve and Blatter (2009)
Parameter

SymbolRate factor constant

A0

Activation energyQ

Gas constant
Tm at 1 atm
p dependence of Tm
Heat capacity
Heat conductivity

R
T0
c


Value

Unit

1.916  103 (T0  263.15 K) s–1 Pa–3
3.985  10–13 (T0 < 263.15 K)
139 (T0  263.15 K)
kJ mol–1
60 (T0 < 263.15 K)
8.314
J mol–1 K–1
273.15K
9.8  10–8
K Pa–1146.3 + 7.253T
J kg–1 K–1
9.828 exp (–0.0057T)
W m–1 K–1

Eqns (3) and (4) to the strain rate via the constitutive
equation given by Glen’s flow law (Glen, 1955)"_ ij ¼ Aen 1 ij ,ð6Þwhere "_ ij and ij are the components of the strain-rate anddeviatoric-stress tensors, respectively, and  e is the effective
stress. We assume the flow-law exponent n = 3 and the rate
factor A is a function of T 0 , ice temperature relative to thepressure-melting point Tm (e.g. Greve and Blatter, 2009).

Q
A ¼ A0 exp
ð7Þ
RT 0
T0 ¼ TTm ¼ T

ðT0

pÞ

ð8Þ

The values used for the constant A0, activation energy Q, gasconstant R, ice melting point at ambient pressure T0 and its
pressure (p) dependence are listed in Table 3. The steadystate energy-balance equation states




@T
@T
@
@T
c ux
¼
þ uz

@x
@z
@x
@x


ð9Þ
@
@T
þ

þ 2Ae nþ1@z@zwhere c is the specific heat capacity and  is the heat
conductivity of ice with the values listed in Table 3.
The modelled domain was discretized into 800 triangle
elements (Fig. 7a) to solve Eqns (1–4) for ux, uz and Eqn (9) for
T with the finite-element method. The horizontal resolutionof the velocity field was 100 m, and the vertical resolutions
were less than several metres near the glacier bed and 10–
15 m near the surface. The resolutions of the temperaturefield were twice as coarse as those for the velocity. This isbecause linear triangle elements were used for the temperature field, whereas quadratic elements were applied for thevelocity field. To solve the velocity field, a stress-freeboundary condition was applied to the upper surface andzero displacement was assumed at the bed. We also assumed
horizontal velocity is zero at x = 0. For the temperature fieldcomputation, surface temperature Ts was prescribed by
Ts ðzÞ ¼ T ref þ ðz

zref Þ,ð10Þ

where T ref = –8.08C is annual mean temperature in 2012
measured at Qaanaaq Airport (z ref = 16 m). We took
¼ 4:42  10–3 K m–1 for the temperature lapse rate,
which was calculated from air temperature records from
July 2012 to June 2013 at Qaanaaq Airport and SIGMA-B.
Geothermal heat flux into the basal boundary was assumed
to be 50 mW m–2, which is a value suggested for this region(Greve, 2005).

Fig. 7. (a) Finite-element mesh used for modelling the longitudinalcross section of QIC. (b) Ice flow vectors and (c) temperature
distribution computed by the numerical model.

We first solved the ice dynamics equations for a velocity
field. The velocities were then used to compute a temperature field, and the velocity computation was subsequently
updated by the temperature using Eqn (7). When ice
temperature reached the pressure-melting point, temperature
was kept constant at this value and melting/freezing
processes were not taken into account. This iteration was
repeated until the velocity and temperature fields converged.

5.2. Model results
Modelling results show ice velocity maxima (up to 6 m a–1)
in the lower middle reaches of the glacier over the relatively
steep surface (Fig. 7b). Ice motion is slow in the upper
reaches, where surface velocity is <2.5 m a–1 within 3000 mof the ice-cap summit. Ice temperature is within the range–13 to –6.58C (Fig. 7c). The location of the greatest ice
motion is influenced by the temperature field, because
relatively warm ice in the lower reaches favors deformation.
The magnitude of the computed velocity is significantly
smaller than the measurements (control experiment in
Fig. 8a). For example, the modelled surface horizontal
velocity is 5.9 m a–1 at Q1203, whereas field measurements
showed 24.6 m a–1 at the same location. The measured
summer velocity is possibly greater than the annual mean,

Sugiyama and others: Field observations on Qaanaaq ice cap

31

6. DISCUSSION
6.1. Summer temperature and melt conditions in
2012

Fig. 8. (a) Horizontal components of the surface ice motion. Black
curve is computed with the numerical model with the standard
conditions described in the text (control experiment). Green, blue
and red curves are computed with the assumptions geothermal heat
flux = 60 mW m–2 (experiment 1), +10 m ice thickness (experiment2) and T ref = –2.08C, (experiment 3), respectively. Observational
data are indicated by the asterisks. (b) Surface (dashed) and bed(solid) temperature computed with the numerical model. Experimental conditions for the black and coloured curves are as in (a).

but the discrepancy is very large. Sensitivity tests were
performed to investigate the influence of uncertainties in thegeothermal heat flux and ice thickness. The surface velocity
increases by taking a greater geothermal heat flux (60 mWm–2) (experiment 1: u ¼ þ11% on average), and byassuming the glacier is 10 m thicker than estimated in
Figure 4 (experiment 2: u ¼ þ45% on average) (Fig. 8a).
However, the effect is too small to explain the inconsistency
with the field data. The changes in glacier bed temperature
are also small, i.e. T = +0.68C and +0.38C for experiments
1 and 2, respectively (Fig. 8b).
The most likely reason for the inconsistency with the fielddata is that the modelled temperature is too low. The
computed glacier bed temperature is below –6.58C (control
experiment in Fig. 8b), which also contradicts the diurnal ice
flow variations observed at Q1203. The diurnal velocity
variations strongly suggest the changes in the basal velocity,i.e. part of the glacier bed is melting and lubricated by
diurnal input of surface meltwater.
To investigate the influence of temperature on the ice
dynamics, an additional experiment was performed byassuming an annual mean temperature T ref = –2.08C at thereference elevation. The model slightly under/overestimatesthe observations in the upper/lower reaches, but relatively
good agreement was achieved (experiment 3 in Fig. 8a).
Glacier bed temperature reaches the pressure-melting pointin the lower reaches at x = 5800–7500 m, which is consistent
with the observed diurnal flow variations. Overall, the model
better reproduces the observations under an unrealisticallyhigher surface temperature than the actual air temperatures.

The melt rate obtained near the terminus of Qaanaaq
Gletscher (46 mm d–1 at Q1201) is comparable to the values
reported for mid-latitude glaciers. For example, mean melt
rate within the lower 1 km of Rhonegletscher, Switzerland(2200–2350 m a.s.l.), from July to September was reported to
be 56 mm w.e. d–1 (Sugiyama and others, 2011). At 4.5 km
from the terminus of Glaciar Perito Moreno, Southern
Patagonia Icefield (360 m a.s.l.), mean melt rate from November to March was 55 mm w.e. d–1 (Stuefer and others,2007). Climate conditions on QIC during such an intensive
melt period were analysed by air temperature recorded at
Qaanaaq Airport. Daily mean temperature variations showthat our study period corresponds to the warmest period in
the year (Fig. 5a). Mean temperature from 18 to 29 July 2012
was 10.28C, and one of the peak temperatures of the year wasobserved on 27 July. In July 2012, an extraordinary hightemperature event hit Greenland and melt occurred overnearly the entire surface of the ice sheet (Nghiem and others,2012). However, temperature at Qaanaaq Airport shows noclear indication of this melt event, and the summer
temperature in 2012 differs little from those in 2010 and2011. Summer (June/July/August) mean temperature in 2012(6.18C) was higher than that in 2010 by 0.38C, but lower thanthat in 2011 by 0.68C. Therefore, the observed melt rates
represent midsummer melt conditions on QIC over the last
several years.

6.2. Enhanced melt on dark surface
The degree-day factors computed at the measurement sites
show an interesting dependence on elevation (Fig 6a). The
greater degree-day factors in the middle reach of the ice cap
can be attributed to ice surface conditions. The ice surface
over the region including Q1203–Q1205 was covered bydark materials, and thus absorbed more solar radiation as aresult of low albedo. The spatial variation in surface
reflectance was clearly observed by in situ visual observations (Fig. 9). Moreover, the brightness intensity of the
ALOS (Advanced Land Observing Satellite) PRISM (Panchromatic Remote-sensing Instrument for Stereo Mapping) satellite image acquired on 25 August 2009 (Fig. 1c) gives aquantitative measure of the surface reflectance along thesurvey route (Fig. 6b). The intensity decreases up-glacier of
the terminus and shows consistently smaller values between
Q1203 and Q1205. The smaller degree-day factors at theupper two sites (Q1206 and Q1207) can be explained byhigher albedo due to the snow cover.
Darkening of the glacier surface has been reported inmany glaciers in the world (e.g. Takeuchi and others, 2001;
Oerlemans and others, 2009), including marginal regions of
the Greenland ice sheet (e.g. Van de Wal and Oerlemans
1994; Wientjes and Oerlemans 2010; Wientjes and others,2011). Such changes in ice surface conditions result in melt
increase due to the reduction in albedo, and thus acceleratemass loss of glaciers and ice sheets. A field study has shownthat the dark appearance over QIC is caused by materials
consisting of mineral dust, snow algae and other organic
matter (Takeuchi and others, 2013). The degree-day factors
obtained in this study confirm the link between thisdarkening surface and melt rate increase. For example, thefactor at Q1204 was 50% greater than that at Q1202,
suggesting that the ice in the dark regions melts 50% more32

Sugiyama and others: Field observations on Qaanaaq ice cap

others, 2008), as well as for McCall Glacier, Alaska, USA(798170 N, 1438500 W) (Pattyn and others, 2005).
Several mechanisms have been proposed for the inadequate modelling results in these glaciers, which are locatedin similar latitude and elevation ranges. First, the ice
temperature near the surface is not simply controlled by
air temperature (Blatter and Hutter, 1991; Bingham andothers, 2008). A large amount of latent heat is generated
when meltwater percolates and refreezes in a cold firn layer.
The insulation effect of a firn layer also plays a role in
keeping the near-surface ice warmer than the annual mean
air temperature. Second, surface meltwater transports sensible heat into a glacier if it penetrates into englacial andsubglacial environments (Bingham and others, 2008; Wohlleben and others, 2009; Phillips and others, 2010). The
diurnal velocity variations confirm penetration of surface
water to the bed of QIC, and we also observed water-filledcrevasses in the accumulation area. Finally, recent changes
in ice thickness may have influenced the ice temperaturefield. The ice temperature may not be in equilibrium withthe present geometry, so warm ice remains as a remnant ofearlier thicker ice conditions. These complex processes must
be investigated further in the field and accurately taken into
account in the model to better simulate the actual flow
velocity and temperature fields of QIC.
Fig. 9. Photographs showing (a) clean and (b) dark ice surface
conditions on QIC. Photographs were taken (a) at Q1201 viewing
up-glacier on 18 July 2012 and (b) at the middle reach of Q1203
and Q1204 viewing down-glacier on 22 July 2012.

than the magnitude expected if the surface was clean, as it isin the lower reaches. Because of the recent negative massbalance in the Greenland coastal regions (Bolch and others,2013), ablation areas are expanding and ice surface is
exposed for longer periods in the summer. This is a favorable
condition for biological activities because they are more
active on a bare ice surface rather than snow (Takeuchi,2013). Thus, further darkening and melt rate increases arelikely to occur over the GICs and ice-sheet margins
in Greenland.

6.3. Ice flow and temperature conditions
Diurnal flow velocity variations at Q1203 (Fig. 3) imply that
surface meltwater drained to the bed through cold ice andenhanced basal ice motion. Thus, at least a part of theglacier bed is melting in this region. This observation
contradicts the wholly frozen ice temperature computedby the numerical model (Fig. 7c). Moreover, the measured
surface velocities are substantially greater than thosecomputed by the model (Fig. 8a). These observations suggestthat the modelled ice temperature is significantly lower than
the actual temperature, and some important thermodynamical processes are missing from the model.
Discrepancies between the modelling results and observational data have previously been reported for other Arctic
glaciers. Blatter and Hutter (1991) applied a thermomechanical glacier model to Laika Glacier, Coburg Island, Canadian
Arctic (758500 N, 798100 W). In contrast to the boreholemeasurements showing the existence of a temperate bed in
the ablation area (Blatter and Kappenberger, 1988), themodel predicted that the glacier is entirely cold. Similar
results were obtained for John Evans Glacier, Ellesmere
Island, Canadian Arctic (798400 N, 748000 W) (Bingham and7. CONCLUSION
To better understand the GICs in Greenland, which arephysically disconnected from the ice sheet, we carried out
initial field observations on QIC. Surface melt rate, ice flow
velocity, and ice thickness were measured along an 8 km
survey route from the ice-cap summit to the terminus of anoutlet glacier.
Melt rates were clearly influenced by dark materials
covering the middle elevation range of the ice cap. The
degree-day factor was greater on the dark surface(8.26 mm w.e. K–1 d–1 at 739 m a.s.l.) than on a clean surface(5.44 mm w.e. K–1 d–1 at 427 m a.s.l.). This result strongly
suggests the important role of darkening ice surface in themelt increase on Greenland’s GICs. Ice velocity showeddiurnal variations, indicating surface meltwater is affecting
the basal ice motion. Meltwater should have been routed
through ice to the bed, although surface annual mean air
temperature is far below the melting point (below –88C at
the terminus). Ice-cap surface and bed geometries measured
in the field were used to compute ice velocity andtemperature using a numerical glacier model. The model
assumed that ice surface temperature is equal to the annual
mean air temperature, and the ice cap is in equilibrium withpresent climate and geometrical conditions. Computed ice
temperature was too low to reproduce the observed ice
velocity, implying the important role of latent and sensible
heat transport by meltwater as well as recent changes in theice thickness. These processes will have to be accurately
taken into account in the model to understand the current
status and future evolution of the GICs in Greenland.

ACKNOWLEDGEMENTS
We thank the members of the 2012 summer field campaign
in Qaanaaq. Special thanks are due to S. Daorana andT. Yamasaki for providing logistical support in Qaanaaq. We
also thank N. Takeuchi and T. Tanikawa for providing

Sugiyama and others: Field observations on Qaanaaq ice cap

helpful comments on the paper, and M. Niwano forprocessing the temperature data at SIGMA-B. The qualityof the paper was improved by critical comments from two
anonymous reviewers and the scientific editor, S. Adhikari.
This research was funded by MEXT Japan (Japanese Ministryof Education, Culture, Sports, Science and Technology)
through the Green Network of Excellence (GRENE) Arctic
Climate Change Research Project and Grant-in-Aid
23221004 (2011–2015). We thank H. Enomoto for his
contribution as leader of the research project.

REFERENCES
Bingham RG, Hubbard AL, Nienow PW and Sharp MJ (2008) An
investigation into the mechanisms controlling seasonal speed-up
events at a High Arctic glacier. J. Geophys. Res., 113(F2),F02006 (doi: 10.1029/2007JF000832)
Blatter H and Hutter K (1991) Polythermal conditions in Arctic
glaciers. J. Glaciol., 37(126), 261–269
Blatter H and Kappenberger G (1988) Mass balance and thermal
regime of Laika ice cap, Coburg Island, N.W.T., Canada.
J. Glaciol., 34(116), 102–110
Bolch T and 6 others (2013) Mass loss of Greenland’s glaciers andice caps 2003–2008 revealed from ICESat laser altimetry data.
Geophys. Res. Lett., 40(5), 875–881 (doi: 10.1002/grl.50270)
Colgan W and 6 others (2012) The annual glaciohydrology cycle in
the ablation zone of the Greenland ice sheet: Part 2. Observed
and modeled ice flow. J. Glaciol., 58(207), 51–64 (doi: 10.3189/
2012JoG11J081)
Fierz C and 8 others. (2009) The international classification forseasonal snow on the ground. (IHP Technical Documents in
Hydrology 83) UNESCO–International Hydrological Programme, Paris
Glen JW (1955) The creep of polycrystalline ice. Proc. R. Soc.
London, Ser. A, 228(1175), 519–538 (doi: 10.1098/rspa.1955.
0066)
Glen JW and Paren JG (1975) The electrical properties of snow andice. J. Glaciol., 15(73), 15–38
Greve R (2005) Relation of measured basal temperatures and thespatial distribution of the geothermal heat flux for the Greenland
ice sheet. Ann. Glaciol., 42(1), 424–432 (doi: 10.3189/
172756405781812510)
Greve R and Blatter H (2009) Dynamics of ice sheets and glaciers.
Springer, Dordrecht
Hanna E, Mernild SH, Cappelen J and Steffen K (2012) Recent
warming in Greenland in a long-term instrumental (1881–2012)
climatic context: I. Evaluation of surface air temperature records.
Environ. Res. Lett., 7(4), 045404 (doi: 10.1088/1748-9326/7/4/
045404)
Hock R (2003) Temperature index melt modelling in mountainareas. J. Hydrol., 282(1–4), 104–115 (doi: 10.1016/S00221694(03)00257-9)
Mernild SH and 6 others (2011) Increasing mass loss fromGreenland’s Mittivakkat Gletscher. Cryosphere, 5(2), 341–348(doi: 10.5194/tc-5-341-2011)
Nghiem SV and 8 others (2012) The extreme melt across theGreenland ice sheet in 2012. Geophys. Res. Lett., 39(20),
L20502 (doi: 10.1029/2012GL053611)
Oerlemans J, Giesen RH and Van den Broeke MR (2009) Retreating
alpine glaciers: increased melt rates due to accumulation of dust
(Vadret da Morteratsch, Switzerland). J. Glaciol., 55(192),
729–736 (doi: 10.3189/002214309789470969)
Palmer SJ, Shepherd A, Sundal A, Rinne E and Nienow P (2010)
InSAR observations of ice elevation and velocity fluctuations at

33

the Flade Isblink ice cap, eastern North Greenland. J. Geophys.
Res., 15(F4), F04037 (doi: 10.1029/2010JF001686)
Paterson WSB (1994) The physics of glaciers, 3rd edn. Elsevier,OxfordPattyn F, Nolan M, Rabus B and Takahashi S (2005) Localizedbasal motion of a polythermal Arctic glacier: McCall Glacier,
Alaska, USA. Ann. Glaciol., 40, 47–51 (doi: 10.3189/
172756405781813537)
Phillips T, Rajaram H and Steffen K (2010) Cryo-hydrologic
warming: a potential mechanism for rapid thermal response ofice sheets. Geophys. Res. Lett., 37(20), L20503 (doi: 10.1029/
2010GL044397)
Rastner P, Bolch T, Mölg N, Machguth H, Le Bris R and Paul F(2012) The first complete inventory of the local glaciers and ice
caps on Greenland. Cryosphere, 6(6), 1483–1495 (doi: 10.5194/
tc-6-1483-2012)
Rinne EJ and 6 others (2011) On the recent elevation changes at the
Flade Isblink Ice Cap, northern Greenland. J. Geophys. Res.,
116(F3), F03024 (doi: 10.1029/2011JF001972)
Schoof C (2010) Ice-sheet acceleration driven by melt supply
variability. Nature, 468(7325), 803–806 (doi: 10.1038/
nature09618)
Stuefer M, Rott H and Skvarca P (2007) Glaciar Perito Moreno,
Patagonia: climate sensitivities and glacier characteristics preceding the 2003/04 and 2005/06 damming events. J. Glaciol.,53(180), 3–16 (doi: 10.3189/172756507781833848)
Sugiyama S, Gudmundsson GH and Helbing J (2003) Numericalinvestigation of the effects of temporal variations in basal
lubrication on englacial strain-rate distribution. Ann. Glaciol.,37, 49–54 (doi: 10.3189/172756403781815618)
Sugiyama M, Yoshizawa T, Huss M, Tsutaki S and Nishimura D(2011) Spatial distribution of surface ablation in the terminus ofRhonegletscher, Switzerland. Ann. Glaciol., 52(58), 1–8 (doi:10.3189/172756411797252185)
Takeuchi N (2013) Seasonal and altitudinal variations in snow algal
communities on an Alaskan glacier (Gulkana Glacier in the
Alaska range). Environ. Res. Lett., 8(3), 035002 (doi: 10.1088/
1748-9326/8/3/035002)
Takeuchi N, Kohshima S and Seko K (2001) Structure, formation,and darkening process of albedo-reducing material (cryoconite)
on a Himalayan glacier: a granular algal mat growing on theglacier. Arct. Antarct. Alp. Res., 33(2), 115–122
Takeuchi N, Nagatsuka N, Shimada R and Uetake J (2013) Biogenic
impurities darkening the Greenland ice sheet. In Detecting theChange in the Arctic System and Searching the Global Influence.
3rd International Symposium on the Arctic Research (ISAR-3),
14–17 January 2013, Tokyo, Japan. http://www.jcar.org/isar-3/
abstracts/ISAR3_abstract.pdf
Van de Wal RSW and Oerlemans J (1994) An energy balance modelfor the Greenland ice sheet. Global Planet. Change, 9(1–2),
115–131
Wientjes IGM and Oerlemans J (2010) An explanation for the dark
region in the western melt zone of the Greenland ice sheet.
Cryosphere, 4(3), 261–268 (doi: 10.5194/tc-4-261-2010)
Wientjes IGM, Van de Wal RSW, Reichart GJ, Sluijs A andOerlemans J (2011) Dust from the dark region in the western
ablation zone of the Greenland ice sheet. Cryosphere, 5(3),589–601 (doi: 10.5194/tc-5-589-2011)
Wohlleben T, Sharp M and Bush A (2009) Factors influencing thebasal temperatures of a High Arctic polythermal glacier. Ann.
Glaciol., 50(52), 9–16 (doi: 10.3189/172756409789624210)
Zwally HJ, Abdalati W, Herring T, Larson K, Saba J and Steffen K(2002) Surface melt-induced acceleration of Greenland icesheet flow. Science, 297(5579), 218–222 (doi: 10.1126/science.
1072708)