1. Introduction
Mass loss from the Greenland ice sheet has contributed to global sea level rise at a mean rate of 0.47 mm a−1 during the period from 1991 to 2015 (van den Broeke and others, Reference van den Broeke2016). Since ~40% of the mass loss was due to changes in ice discharge from marine-terminating outlet glaciers (van den Broeke and others, Reference van den Broeke2016), a future projection of ice discharge from the Greenland ice sheet is crucial for understanding the magnitude and timing of global sea level rise. Ice discharge from the ice sheet is dependent on the ice speed of marine-terminating outlet glaciers. Glaciers in Greenland had accelerated substantially in the early 2000s, and ice speed has been kept at the accelerated level since 2005 (e.g. Howat and others, Reference Howat, Joughin, Fahnestock, Smith and Scambos2008; Moon and others, Reference Moon, Joughin, Smith and Howat2012; King and others, Reference King2018; Mankoff and others, Reference Mankoff2019). Observed speed changes are inhomogeneous in space and time; thus further investigations and insights are needed to project future changes in marine-terminating outlet glaciers in view of our changing climate.
Recent progress in satellite remote sensing techniques contributes to the studies of seasonal speed patterns of marine-terminating glaciers in Greenland (Howat and others, Reference Howat, Box, Ahn, Herrington and McFadden2010; Moon and others, Reference Moon2014, Reference Moon, Joughin and Smith2015; Vijay and others, Reference Vijay2019). Reported data indicate significant summer acceleration, suggesting the importance of high-frequency ice-speed measurements to the quantification of annual ice discharge. Sole and others (Reference Sole2011) reported the ~1% contribution of summer speedup to annual ice motion of a marine-terminating outlet glacier in southwest Greenland. However, as the study site was far inland from the glacier terminus at >1200 m a.s.l., the importance of summer speedup to glacier dynamics in the vicinity of the terminus remains poorly understood.
Seasonal speed variations are essential also for understanding the response of ice dynamics to external forcings (Howat and others, Reference Howat, Box, Ahn, Herrington and McFadden2010; Carr and others, Reference Carr, Stokes and Vieli2013; Moon and others, Reference Moon2014; Lemos and others, Reference Lemos, Shepherd, McMillan and Hogg2018; Vijay and others, Reference Vijay2019). Moon and others (Reference Moon2014) reported the seasonal speed patterns of 55 marine-terminating glaciers along the coast of northwest, west, southeast and southwest Greenland. This study classifies seasonal speed patterns into three types, depending on their sensitivity to the terminus position, meltwater production and the subglacial hydrological system. Seasonal speed patterns in northwestern Greenland were correlated with meltwater production, suggesting a higher sensitivity of the glacier dynamics in this region to climate variability (Moon and others, Reference Moon2014). More recently, higher frequency ice-speed measurements were performed on 45 glaciers over broader areas in Greenland (Vijay and others, Reference Vijay2019). This study shows the concurrence of summer speedup and the onset of glacier surface melting, which essentially agrees with the conclusion by Moon and others (Reference Moon2014). These studies indicate the importance of expanding the spatial and temporal coverage of these data to include the investigation of future responses of Greenlandic outlet glaciers to changing climate.
Seasonal speed variations in marine-terminating glaciers are influenced by several factors such as subglacial hydraulic conditions, surface melting and drainage into the glacier bed, terminus variations and sea-ice/ice-mélange conditions (e.g. Howat and others, Reference Howat, Box, Ahn, Herrington and McFadden2010; Carr and others, Reference Carr, Stokes and Vieli2013; Moon and others, Reference Moon2014; Pimentel and others, Reference Pimentel2017). Numerical modeling and observational studies indicate the importance of the glacier terminus position in ice dynamics (Howat and others, Reference Howat, Joughin, Fahnestock, Smith and Scambos2008; Joughin and others, Reference Joughin2008a, Reference Joughin2008b; Nick and others, Reference Nick, Vieli, Howat and Joughin2009). Observations in southeast Greenland have shown that resistive stress acting on the terminus is diminished when a marine-terminating glacier retreats, a process which is significant enough to cause flow acceleration (Howat and others, Reference Howat, Joughin, Fahnestock, Smith and Scambos2008). Numerical experiments on Helheim Glacier confirmed the critical effect on glacier dynamics of bed geometry near the terminus (Nick and others, Reference Nick, Vieli, Howat and Joughin2009). Sea ice and ice mélange in front of a marine-terminating glacier are thought to have a buttressing effect on the glacier front (Joughin and others, Reference Joughin2008a; Amundson and others, Reference Amundson2010). Therefore, the melting of sea ice and ice mélange has the potential to reduce resistive stress acting on the glacier and cause summer speedup. The meltwater input to the glacier bed and the evolution of a subglacial drainage system also influence seasonal speed variations (Joughin and others, Reference Joughin2008b; Schoof, Reference Schoof2010; Sole and others, Reference Sole2011). In early summer, surface melt increases and water drains to the bed, where a subglacial drainage system is not sufficiently developed. In such a situation, subglacial water pressure increases and basal sliding is enhanced. The drainage system develops further as a greater amount of meltwater drains, resulting in a drop in subglacial water pressure and a decrease in ice speed (e.g. Bartholomew and others, Reference Bartholomew2010; Sundal and others, Reference Sundal2011). These mechanisms have been proposed and tested at some of the glaciers in Greenland, but processes affecting seasonal speed variations are not yet fully understood.
To investigate seasonal speed variations of marine-terminating outlet glaciers, we processed satellite images for high-frequency ice-speed variations in 10 marine-terminating outlet glaciers along the coast of the Prudhoe Land, northwestern Greenland. We compared observed ice-speed changes in years 2014, 2015 and 2016 with the terminus position, sea-ice/ice-mélange conditions and air temperature in order to discover the source of seasonal speed variations. The contribution of summer speedup to annual ice motion was quantified to evaluate its importance to ice discharge from the glaciers. We also computed seasonal variations in frontal ablation, by comparing the ice speed and terminus positions obtained in this study.
2. Study site
We focused on 10 marine-terminating outlet glaciers located along the coast of Prudhoe Land in northwestern Greenland (77.53–77.94°N and 66.04–71.63°W) (Fig. 1). This area has been a central part of previous studies on the dynamics of outlet glaciers (Porter and others, Reference Porter2014; Sugiyama and others, Reference Sugiyama, Sakakibara, Tsutaki, Maruyama and Sawagaki2015; Tsutaki and others, Reference Tsutaki, Sugiyama, Sakakibara and Sawagaki2016; Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018), but not covered by past studies of seasonal speed variations (Moon and others, Reference Moon2014, Reference Moon, Joughin and Smith2015; Vijay and others, Reference Vijay2019). A number of outlet glaciers, having a front width between 1.5 km (Sharp Glacier) and 5.0 km (Tracy Glacier), are located within a relatively small area. Glaciers in the study area have been retreating since the early 2000s and some of the glaciers have shown significant acceleration during this period (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018).
Fig. 1. The locations of the study glaciers along the coast of Prudhoe Land. The background is a true color mosaic Landsat image acquired on 9 July 2014. The inset shows the location of the study area in Greenland.
Heilprin and Tracy Glaciers are the two largest glaciers in this region, having a drainage basin of 1.04 × 104 km2 in total (Rignot and Kanagaratnam, Reference Rignot and Kanagaratnam2006). Tracy Glacier retreated at a rate of 200 m a−1 and accelerated from 810 to 1740 m a−1 during the period from 1980s to 2014 (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018), which was the greatest change recorded among the study glaciers. Farquhar, Melville and Sharp Glaciers flow at a rate more than 400 m a−1 (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018). These five glaciers are located in the eastern half of the study area, draining into a ~20 km wide fjord Inglefield Bredning. Hubbard, Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers are located in the western half of the study area, where each individual fjord connects to western part of Inglefield Bredning or directly to northern Baffin Bay. The ice speed near the termini of these glaciers ranges from 90 to 400 m a−1 (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018), which is smaller than the speed of the glaciers in the eastern half. Heilprin, Farquhar, Bowdoin and Diebitsch Glaciers showed a substantial speedup and retreated rapidly after 2000 (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018).
3. Data and methods
To measure the glacier terminus position, ice speed and sea-ice/ice-mélange conditions from 2014 to 2016, we used satellite images of the Landsat 8 Operational Land Imager (OLI) (level 1T, with 15/30 m resolutions) distributed by the United States Geological Survey. Because the OLI is an optical sensor, the measurements were performed when sunlight was available between March and October.
We obtained hourly air temperatures at Mittarfik Qaanaaq (Qaanaaq Airport; 77.48°N, 69.38°W; 16 m a.s.l.) located within 30–84 km of the glacier termini, from the National Oceanic and Atmospheric Administration. The temperature dataset was utilized to calculate mean air temperature during a period for each ice-speed measurement. The Qaanaaq temperature showed a good correlation (r 2 = 0.94, p < 0.001) with that recorded at Thule Airbase located 108 km south from Qaanaaq. Thus, we assume the data well represent air temperature over the study area.
3.1. Ice speed
The ice speed was measured by automatically tracking glacier surface features on temporally separated satellite-image pairs (Scambos and others, Reference Scambos, Dutkiewicz, Wilson and Bindschadler1992; Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2014, Reference Sakakibara and Sugiyama2018). We used Landsat 8 OLI images obtained from 2014 to 2016 for this purpose. These images were acquired with temporal separations of 16–48 days and study glaciers were covered by multiple tracks, thus enabling high-frequency speed measurements to be made. The frequency of the ice-speed measurement at each glacier was 44–196 times per year, depending on the availability of the scene. Glacier surface displacement was measured with an automatic image matching scheme known as the orientation correlation method in the frequency domain (CCF-O) (Heid and Kääb, Reference Heid and Kääb2012). The correlation window size was 30 pixels, and the Landsat images were co-registered using displacement vectors obtained in off-ice areas. Measurements were taken over areas near the glacier termini, and we sampled speed at the glacier center at a distance of 0.7–2.4 km from the termini in 2014. The distance from the terminus to the sampling site for each glacier was approximately equal to the half of the glacier front width (Moon and others, Reference Moon2014). The magnitude of summer speedup was dependent on the distance from the glacier front, but the timing of the peak speed and seasonal patterns was not influenced by the sampling location. Displacement vectors obtained with a low signal-to-noise ratio (<30) were excluded, and remaining vectors were filtered using directional and median low-pass filters (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2014). Ice-speed data at the sampled locations were smoothed by taking the median values of 7 by 7 neighboring pixels (Howat and others, Reference Howat, Box, Ahn, Herrington and McFadden2010). Ambiguities due to the uncertainty in the cross-correlation peak and co-registration errors ranged from 6 to 300 m a−1 with a mean value of 43 m a−1, depending on image quality and length of temporal separation. After the filtering and smoothing procedures, measurement error was mostly associated with the co-registration of image pairs. Thus, the uncertainty was not influenced by seasonal changes in glacier surface conditions.
3.2. Terminus position and sea-ice/ice-mélange conditions
Glacier terminus positions were mapped from 2014 to 2016, using false color composite images generated by assigning the colors blue, green and red to bands 4, 5 and 6 of the Landsat 8 images, respectively (McNabb and Hock, Reference McNabb and Hock2014). Relatively coarse false color composite images were used for this purpose because the boundary between a glacier terminus and sea ice was more clearly observed than panchromatic images available in a higher resolution. In each image, the glacier frontal margins were delineated manually using the QGIS geographic information system software. Mean displacement of the terminus position was computed by dividing the areal change near the frontal margin by the width of the terminus (Moon and Joughin, Reference Moon and Joughin2008). Relative positioning errors introduced during the process of overlapping the Landsat 8 images was ±8 m (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2018). This value was derived as the co-registration error in the ice-speed measurement, by computing the root-mean-square of displacement obtained in off-ice areas.
Sea ice and ice mélange conditions near the termini were determined using a method proposed by Moon and others (Reference Moon, Joughin and Smith2015). In every image, sea-ice/ice-mélange concentration and potential rigidity within ~5 km from each glacier terminus were classified as likely to be rigid (rigid), potentially rigid (mixed) or unlikely rigid (open) (Moon and others, Reference Moon, Joughin and Smith2015) (example is shown in Fig. 2). The condition was classified as ‘rigid’ when the region in front of the glacier was completely or near-completely covered by sea ice/ice mélange. An ‘open’ classification was given for open water. When the region was covered with a substantial amount of ice but also with open water, the condition was classified as ‘mixed’ (Moon and others, Reference Moon, Joughin and Smith2015).
Fig. 2. False color mosaic Landsat images used for categorization of fjord conditions. The examples show conditions in front of Tracy Glacier, which were categorized as (a) rigid, (b) mixed and (c) open sea ice.
3.3. Frontal ablation
Frontal ablation (sum of calving and subaqueous melting) was calculated by subtracting the terminus displacement from ice motion over a period from March 2014 to September 2016. Mean ice speed across the glacier was estimated as 55–90% (mean: 73%) of the speed obtained at the sampling site at the glacier center. The conversion factor was derived for each glacier by computing a mean speed across the glacier using all available measurements. We assume ice speed was uniformly distributed from the surface to the bottom (perfect sliding). Uncertainty in the calculated frontal ablation was estimated as 15 to 24 m from the error ranges assumed in the terminus position and ice-speed measurements.
4. Results
4.1. Ice speed
The study glaciers showed clear seasonal speed variations (Fig. 3). For example, the ice speed at Tracy Glacier decreased slightly from 2200 to 2100 m a−1 during the period from March to May 2015 (Fig. 3e). The glacier progressively accelerated to 2500 m a−1 from late May to early July, followed by gradual deceleration to 2200 m a−1 from July to September. A similar seasonal speed pattern was observed at Bowdoin Glacier. The ice speed at Bowdoin Glacier was within the range of 300 to 400 m a−1 during the period from March to May 2015 (Fig. 4e). Ice speed increased from 400 to 600 m a−1 from late May to early July, and decreased to 350 m a−1 from July to September. Early summer speedup and late summer slowdown were commonly observed at the other study glaciers (Fig. 4). In general, ice speed showed only minor changes from March to May. Acceleration typically began between late May and early June, followed by the achievement of a peak speed between late June and early July. Ice speed subsequently slowed down from July to September.
Fig. 3. (a–o) Ice speed (red), terminus position (blue) and sea-ice/ice-mélange conditions (magenta = rigid, gray = mixed and cyan = open) for five study glaciers in the eastern part of the study area. Red lines and associated shaded bands indicate daily mean speed and the std dev., respectively. The distance shown next to the glacier names indicates the location of ice-speed sampling as measured from the 2014 termini. (p–r) Daily (orange) and hourly (yellow) air temperatures, and the sum of the positive degree hours (PDH) (blue) at Qaanaaq Airport. The vertical orange lines show the first day of positive hourly temperature in each year.
Fig. 4. (a–o) Ice speed (red), terminus position (blue) and sea-ice/ice-mélange conditions (magenta = rigid, gray = mixed and cyan = open) for five study glaciers in the western part of the study area. Red lines and associated shaded bands indicate daily mean speed and the std dev., respectively. The distance shown next to the glacier names indicates the location of ice-speed sampling as measured from the 2014 termini. (p–r) Daily (orange) and hourly (yellow) air temperatures, and the sum of the PDH (blue) at Qaanaaq Airport. The vertical orange lines show the first day of positive hourly temperature in each year.
The timing of the seasonal speedups and slowdowns is comparable in the study glaciers regardless of size and location. However, the magnitude of summer speedup was significantly different for each glacier. The most significant speedup was observed at Farquhar Glacier, where ice speed increased by more than 300 m a−1. The magnitude of the speedup at other glaciers ranged between 54 and 278 m a−1 corresponding to a range between 8 and 226% increase relative to the speed before the acceleration. Seasonal speed patterns do not always follow the general trend described above. For example, ice speed at Heilprin Glacier showed no significant change from May to July in 2014 (Fig. 3a). Speedup was unclear also at Farquhar Glacier in 2014 (Fig. 3g).
4.2. Terminus position and sea-ice/ice-mélange conditions
Some of the glaciers showed substantial sudden retreat due to large calving events, typically by more than 100 m within the interval of satellite images ranging between 2 and 23 days. For example, Heilprin Glacier retreated by 200 m in June 2014 (Fig. 3a) and in April 2015, and 110 m in May 2015 (Fig. 3b). Similarly, sporadic major calving events were observed at Tracy (Figs 3d to f) and Melville Glaciers (Figs 3j to l). Farquhar (Figs 3g to i) and Bowdoin Glaciers (Figs 3dto f) experienced slightly smaller events (>50 m retreat). These sudden retreats were followed by gradual advances lasting for a few months. Some other glaciers showed seasonal variations in the terminus position. For example, Morris Jesup Glacier advanced by 30 m from May to June and retreated by 130 m from early July to September in 2014 (Fig. 4j). Verhoeff and Diebitsch Glaciers showed similar frontal variations (Figs 4g to i and 4m to o). The change in the trend from advance to retreat occurred typically in June or July.
Sea ice and ice mélange disappeared at most of the glacier fronts in late June to early July (Figs 3, 4). The timing of the sea-ice/ice-mélange disappearance was consistent over the three years studied. Fjords opened earlier at Morris Jesup and Diebitsch Glaciers, both situated in the western part of the study area (Figs 4j to o). Because the fjords of these two glaciers are directly connected to Baffin Bay, sea ice and ice mélange were released earlier in summer that is from mid-April to early June. By contrast, sea ice and ice mélange stayed longer after the onset of melting in Inglefield Bredning. At the fronts of Heilprin, Tracy, Farquhar and Melville Glaciers, mixed conditions lasted for more than two weeks before the glacier front became ice-free. No clear relationship was found between the sea-ice/ice-mélange conditions and the large calving events observed at these glaciers.
5. Discussion
5.1. Mechanism of seasonal speed variations
Our data reveal details of seasonal speed variations along the coast of Prudhoe Land. All the study glaciers show speedup in summer. The speedup initiated between late May and early June, speed gradually increased until early July and then it slowed down in the period from July to September. The magnitude of the seasonal speedup was a few 100 m a−1 which corresponds to 8–226% increase relative to pre-speedup conditions (Figs 3, 4). The observed seasonal speed pattern and magnitude of the speed change are similar to those reported in the Melville Bay region in northwestern Greenland (74.92–76.08°N) (Moon and others, Reference Moon2014). Thus, the observed seasonal speed pattern appears to be a ubiquitous feature across northwestern Greenland, and the ice-speed variations are likely driven by a common mechanism.
To investigate a link between ice speed and atmospheric conditions, the ice-speed time series data were plotted against air temperature averaged over the same period (Fig. 5). This analysis was performed for the data in 2015 because ice speed was measured most frequently for that year out of the study years. The plots demonstrate an influence of air temperature on the speed of the studied glaciers.
Fig. 5. Scatter plots of ice speed and air temperature at Qaanaaq Airport averaged over each ice-speed measurement period. The data covers a period from April to October 2015.
For example, ice speed at Heilprin Glacier was fairly stable within 1300–1600 m a−1 under freezing temperature conditions (arrow 1 in Fig. 5a). This corresponds to a period from March to May (Figs 3b, 3q). About one month after hourly mean temperature exceeded 0°C (vertical line in Fig. 3q), ice speed began to increase, and further accelerated in the course of the temperature rise. Ice speed and temperature showed a positive correlation until speed reached the maximum in early July (arrow 2 in Fig. 5a). Similar speed-temperature relationships were observed at Tracy, Farquhar, Bowdoin and Verhoeff Glaciers (Figs 5b to c and 5g to h). The other glaciers showed some unique characteristics, but timing and trend of the speedup were similar to the relationship discussed above, i.e. insignificant speed change when temperatures are below zero, and acceleration to the summer maximum speed as temperatures increase. The relationship implies that the summer speedup was triggered by meltwater production and enhanced as the melt rate increased. Based on these observations, we assume that subglacial water pressure increased due to meltwater penetration into the glacier bed, and the glacier accelerated due to enhanced basal sliding (e.g. Iken, Reference Iken1972; Andrews and others, Reference Andrews2014; Ryser and others, Reference Ryser2014; Vijay and others, Reference Vijay2019).
After peak speed was attained at Heilprin Glacier, ice motion rapidly decelerated while air temperature remained high (arrow 3 in Fig. 5a). Ice speed remained low until the end of the summer, resulting in a hysteresis loop between ice speed and temperature. A similar drop in ice speed from July to August was observed at all glaciers (Fig. 5). The glaciers decelerated in the middle of summer when meltwater was abundant, which we attribute to the seasonal evolution of the subglacial drainage system. Observations at alpine glaciers have demonstrated that an inefficiently linked cavity system evolves into an efficient conduit drainage system over the course of a melt season (e.g. Iken and Bindschadler, Reference Iken and Bindschadler1986; Fountain and Walder, Reference Fountain and Walder1998). After the efficient drainage system is established, subglacial water pressure drops and glacier speed becomes unaffected by the amount of meltwater input to the bed (e.g. Sugiyama and Gudmundsson, Reference Sugiyama and Gudmundsson2003; Schoof, Reference Schoof2010). This interpretation is consistent with the speedup mechanism described previously.
At Melville Glacier, ice speed abruptly increased to >1300 m a−1 when the air temperature was still below 0°C (Fig. 5d). This speedup occurred in May, when the fjord was still covered with sea ice, immediately after the sudden terminus retreat of ~300 m due to a large calving event (Fig. 3k). Except for this incident, no clear relationship was observed between frontal variations and glacier dynamics. A sudden retreat of >200 m was observed frequently at Heilprin (Figs 3a, b), Tracy (Figs 3d to f) and other glaciers, but ice speed was unaffected by these events. Seasonal variations in the front position were also not relevant to the change in speed. Morris Jesup, Verhoeff and Diebitsch Glaciers began retreat in June or July, whereas summer speedup of these glaciers began already in May. Therefore, reduction in resistive stress due to terminus retreat (Joughin and others, Reference Joughin2008a, Reference Joughin, Smith, Shean and Floricioiu2014) is not the source of the observed seasonal speed variations. We also found no link between sea-ice/ice-mélange disappearance from the fjords and the onset of the speedup. Sea ice and ice mélange disappeared typically in early July when summer speedup ceased (Figs 3, 4).
Based on the discussion above, we conclude that summer speedup observed at the glaciers in Prudhoe Land was driven by subglacial water pressure which increased due to the meltwater input to the bed. The glaciers decelerated in midsummer when an efficient drainage system developed and water pressure dropped. This result is consistent with a recent study on other glaciers in Greenland reporting the primary influence of surface meltwater production and associated subglacial hydraulic conditions on a seasonal ice-speed change (Vijay and others, Reference Vijay2019). The effect of the glacier front position on the ice speed is limited, although sudden retreat due to a large calving event sometimes triggered glacier acceleration.
5.2 Contribution of summer speedup to glacier motion
As demonstrated in recent studies (Moon and others, Reference Moon2014, Reference Moon, Joughin and Smith2015), seasonal speed patterns are crucial for quantifying the annual discharge of marine-terminating outlet glaciers in Greenland. To investigate the relative importance of the observed summer speedup in longer-term ice discharge, we quantified the contribution of summer speedup to annual ice motion. We define the magnitude of summer speedup as the difference between summer speed (June–August) and ‘base speed’ defined by the mean over the rest of the year (September–May of the following year). In other words, summer speedup refers to the additional displacement that occurs on top of the base flow. The magnitude of the summer speedups was compared with annual ice motion, using data from September 2014 to August 2016 (Fig. 6).
Fig. 6. The contribution of summer speedup (June–August) to mean ice speed from September 2014 to August 2016 for (a) Heilprin, Tracy, Farquhar and Melville Glaciers and (b) Sharp, Hubbard, Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers. The numbers indicate the relative contributions of the summer speedup in %. (c) Scatter plots of the summer speedup and its relative contribution vs the mean ice speed.
The contribution of the summer speedup to annual motion of fast-flowing Heilprin, Tracy, Farquhar and Melville Glaciers, which flow at a rate >1000 m a−1, was within the range of 24–62 m a−1 (2–5%) (Fig. 6a). Summer speedup contribution of the other six glaciers ranged from 2 to 23 m a−1 (0.5–13%) (Fig. 6b). Summer increase in total ice discharge from the Greenland ice sheet is reported as ~6% of the annual mean (King and others, Reference King2018), which implies a 2–3% contribution of summer speedup to annual ice discharge. This estimate is consistent with our results obtained for fast-flowing glaciers. Our results demonstrate that the impact of summer speedup on annual ice motion was greater at relatively slow-flowing glaciers, whereas larger magnitude of speedup was observed at fast-flowing glaciers (Fig. 6c). Seasonal speed variation was particularly important at Verhoeff and Morris Jesup Glaciers, where summer speedup accounted for more than 10% of annual ice motion. The numbers obtained for the slow-flowing glacier are similar to those reported for a land-terminating glacier in western Greenland (6–14%) (Bartholomew and others, Reference Bartholomew2010). In the previous studies, the magnitude of summer speedup of marine-terminating glaciers was greater near the terminus than inland (Sole and others, Reference Sole2011; Joughin and others, Reference Joughin, Smith, Shean and Floricioiu2014). Therefore, it is important to understand seasonal speed patterns near the glacier termini to accurately quantify annual ice discharge from marine-terminating glaciers.
5.3. Seasonal variations in the frontal ablation
The high-frequency ice speed and terminus position data enabled us to quantify temporal variations in frontal ablation. Figure 7 shows cumulative forward ice motion near the glacier front, terminus position (positive number shows retreat) and cumulative frontal ablation of Tracy and Bowdoin Glaciers. The frontal ablation of Tracy Glacier occurred at a nearly constant rate throughout the study period, except for several calving events represented by the sudden retreat in May and July 2015, and May 2016 (Fig. 7a). These large calving events occurred irrespective of sea-ice/ice-mélange conditions (Figs 3e, f). Frontal ablation at Heilprin and Farquhar Glaciers showed seasonal patterns similar to that at Tracy Glacier (Figs s2a, b). Large calving events are more frequent during summer periods, but the glacier front steadily loses ice by frontal ablation throughout the year. In contrast to these three glaciers, frontal ablation at Bowdoin Glacier was greater in summer than during the other period (Fig. 7b). Ablation during a period from June to August (25% period of a year) accounted for 59% of the total ablation which occurred from September 2014 to August 2016. This seasonality in the frontal ablation caused seasonal variations in the terminus position, i.e. retreat in summer and advance during the rest of the year (Fig. 7b). Forward ice motion was greater in summer, but enhanced frontal ablation was the dominant factor affecting change in the terminus position. Similar trends in the frontal ablation were observed in Melville, Verhoeff, Morris Jesup and Diebitsch Glaciers (Figs s2c, f, g and h).
Fig. 7. Cumulative ice motion near the glacier front (blue), ice front displacement (positive number shows retreat) (red) and cumulative frontal ablation obtained as the sum of the ice motion and the frontal displacement (green) from 2014 to 2016 for (a) Tracy and (b) Bowdoin Glacier.
The contribution of the summer frontal ablation (June–August) to the total ablation varied for each glacier (Fig. 8). At Heilprin Glacier, summer ablation accounts for 22% (140 m a−1) of the ablation that occurred from September 2014 to August 2016, indicating a near-uniform distribution over the period of a year. The contribution was slightly greater at Tracy, Farquhar, Melville, Sharp and Hubbard Glaciers, ranging from 33 to 41% (60–570 m a−1). Around the half of the annual frontal ablation occurred during summer at Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers (80–180 m a−1) (Figs 8a, b). Summer frontal ablation is generally correlated with annual ablation, but its relative contribution decreases as annual ablation increases (Fig. 8c).
Fig. 8. Frontal ablation rate for (a) Heilprin, Tracy, Farquhar and Melville Glaciers and (b) Sharp, Hubbard, Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers. The numbers indicate the relative contributions of the summer frontal ablation (red) to the annual ablation in %. The data are the mean values over a period between September 2014 and September 2016. (c) Scatter plots of the summer frontal ablation and its relative contribution vs the annual frontal ablation rate.
The first group, glaciers with year-round distributed frontal ablation, was characterized by fast ice flow and large-glacier size. We assume that the influence of summer speedup and sea-ice/ice-mélange melting on calving was less significant at this type of glacier. This is because the relative contribution of summer speedup to annual ice motion is small at these glaciers (Fig. 6). It is also likely that the buttressing effect of sea ice and ice mélange was insufficient to stabilize the fast-flowing glacier front. The greater contribution of summer frontal ablation at other glaciers can be explained by calving enhanced by speedup and the loss of sea ice/ice mélange, which play a more important role in the calving of smaller and slow-flowing glaciers.
We also suspect an influence of submarine melting on the seasonal variations in the frontal ablation. Fjord depth at the glacier front is important for submarine melting because it affects the intrusion of warm Atlantic water (AW) (Porter and others, Reference Porter2014; Rignot and others, Reference Rignot2016). According to the glacier bed elevation measured in 2014 (CReSIS, 2016), the fjords are deeper than 300 m below sea level at the front of Heilprin, Tracy and Farquhar Glaciers. AW intrudes into such deep fjords and provides heat for submarine melting throughout the year (Porter and others, Reference Porter2014). In contrast to these large glaciers, Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers flow into fjords shallower than 200 m (Sugiyama and others, Reference Sugiyama, Sakakibara, Tsutaki, Maruyama and Sawagaki2015; OMG mission, 2016). A submarine melt rate increases in summer because of fjord circulation driven by subglacial meltwater discharge (e.g. Sciascia and others, Reference Sciascia, Straneo, Cenedese and Heimbach2013). We speculate that the contribution of summer submarine melting is more significant at these glaciers because warm AW is not available in winter. Submarine melting increases in summer and enhances calving (Fried and others, Reference Fried2015; Rignot and others, Reference Rignot, Fenty, Xu, Cai and Kemp2015), which results in greater seasonal variations in frontal ablation.
6. Conclusions
We studied ice-speed variations of 10 marine-terminating outlet glaciers along the coast of Prudhoe Land in northwestern Greenland. During the study period from 2014 to 2016, the glaciers showed a common seasonal speed pattern: speedup began between late May and early June, increased gradually until it reached its peak value between late June and early July, and subsequently slowed down from July to September. The magnitude of the speedup was within the range of 54 to 278 m a−1 (8–226%). This seasonal speed pattern was similar to the one reported in the nearby Melville Bay region in northwest Greenland (74.92–76.08°N) (Moon and others, Reference Moon2014), suggesting that the observed seasonal speed pattern is commonly distributed over all of northwestern Greenland.
Observed summer speedups began approximately when hourly air temperatures exceeded 0°C in early summer. Ice speed decelerated in midsummer when air temperature still remained at the highest level. The influence of frontal retreat and sea-ice/ice-mélange conditions was insignificant. Based on these observations, we conclude that summer speedup was driven by subglacial water pressure which increased due to the meltwater input to the glacier bed, and the glacier decelerated when the pressure dropped through the development of an efficient subglacial drainage system. The contribution of summer speedup to annual ice motion ranged from 2 to 62 m a−1 (0.5–13%), and the relative contribution was greater at slow-flowing glaciers.
Frontal ablation of the study glaciers showed seasonal variations as well. Greater ablation was observed in the summer period (June–August) and the seasonality is more enhanced in relatively small and slow-flowing glaciers. For example, more than a half of the annual frontal ablation occurred during the summer period at Bowdoin, Verhoeff, Morris Jesup and Diebitsch Glaciers. These glaciers showed clear seasonal variations in the front position: advance from March to June/July and retreat from July to September. The frontal variations were dominated by the seasonal variations in the frontal ablation rather than by ice-speed variations, suggesting the importance of ocean influence on calving and submarine melting. Observed seasonal front variations were critically controlled by summer frontal ablation particularly at relatively small slow-flowing glaciers. At the front of these glaciers, calving is enhanced by speedup and sea-ice/ice-mélange weakening during the summer period.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/jog.2019.81.
Acknowledgements
We thank Jason Amundson and an anonymous referee for constructive review comments. The paper was handled by the Scientific Editor Shad O'Neel and commented by the Chief Editor Hester Jiskoot. The Landsat images were downloaded from the United States Geological Survey Earth Explorer (http://earthexplorer.usgs.gov/). The dataset of air temperature was downloaded from the National Oceanic and Atmospheric Administration's National Centers for Environmental Information (https://www.ncdc.noaa.gov/). English text was corrected by Susan Braun-Clarke. This study was funded by MEXT Japan (Japanese Ministry of Education, Culture, Sports, Science and Technology) through the Arctic Challenge for Sustainability (ArCS) Project.
Author contribution
DS processed satellite data, analyzed the results and wrote the paper. SS contributed to the interpretation of the results and writing.
Open access
