Hostname: page-component-848d4c4894-nr4z6 Total loading time: 0 Render date: 2024-06-06T06:50:09.030Z Has data issue: false hasContentIssue false
  • English
  • Français

Terminus thinning drives recent acceleration of a Greenlandic lake-terminating outlet glacier

Published online by Cambridge University Press:  01 April 2024

Ed Holt Open the ORCID record for Ed Holt [Opens in a new window] ,
Peter Nienow  and
Encarni Medina-Lopez
Ed Holt*
Affiliation:
School of Engineering, University of Edinburgh, Edinburgh, UK
Peter Nienow
Affiliation:
School of Geosciences, University of Edinburgh, Edinburgh, UK
Encarni Medina-Lopez
Affiliation:
School of Engineering, University of Edinburgh, Edinburgh, UK
*
Corresponding author: Ed Holt; Email: ed.holt@ed.ac.uk
Rights & Permissions [Opens in a new window]

Abstract

Ice-contact proglacial lakes affect ice dynamics and the transition of glacier termini from land- to lake-terminating has been shown to cause ice flow acceleration. In recent decades, the number and size of Greenlandic ice-marginal lakes has increased, highlighting the need to further understand these lake-terminating ice-margins as their influence on ice sheet mass balance increases. Here, time series of satellite-derived observations of ice velocity, surface elevation, and terminus position were generated at a lake-terminating outlet glacier, Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia in south-west Greenland. At Isortuarsuup Sermia, annual surface velocity at the terminus increased by a factor of 2.5 to 214 ± 4 m yr−1 (2013–2021), with the magnitude of this acceleration declining with distance up-glacier. Meanwhile, near-terminus surface elevation changed at a rate of −2.3 ± 1.1 m yr−1 (2012–2021). Conversely, velocity change at Kangaasarsuup Sermia was minimal, while surface elevation change was approximately half at comparable elevations (−1.2 ± 0.3 m yr−1). We attribute these dynamic differences to thinning at Isortuarsuup Sermia and subsequent retreat from a stabilising sublacustrine moraine, and emphasise the potential of proglacial lakes to enhance future rates of mass loss from the Greenland Ice Sheet.

Keywords

glacier monitoringice dynamicsice velocityremote sensing
Type
Article
Information
Journal of Glaciology , First View , pp. 1 - 13
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
Copyright © The Author(s), 2024. Published by Cambridge University Press on behalf of International Glaciological Society

1. Introduction

Rates of mass loss from the Greenland Ice Sheet (GrIS) increased six-fold between the 1980s and 2018 (Mouginot and others, Reference Mouginot2019), raising sea levels by 10.8 ± 0.9 mm (1992–2018) (The IMBIE Team, 2020). The GrIS is projected to continue losing mass, and estimates of GrIS sea level rise (SLR) contributions vary with emissions scenario. By 2100 an additional 70 ± 40 mm of SLR is projected under RCP2.6, increasing to between 80 and 270  mm under RCP8.5 (Fox-Kemper and others, Reference Fox-Kemper2021). Refining these SLR estimates requires greater understanding of the controls on ice sheet mass loss.

The GrIS is currently losing mass via both surface and dynamic processes (The IMBIE Team, 2020), with dynamics being, in part, dependent on surface processes. For example, long-term negative surface mass balance induced thinning (negative surface elevation change) may lead to acceleration at lake or marine-terminating margins if thinning causes greater reductions in resistive stresses than driving stresses (Pfeffer, Reference Pfeffer2007). Furthermore, at lake- and marine-terminating outlets, acceleration and surface lowering may be self-sustaining if accompanied by retreat into deeper water and further dynamic thinning (Weertman, Reference Weertman1974; Meier and Post, Reference Meier and Post1987; O'Neel and others, Reference O'Neel, Pfeffer, Krimmel and Meier2005; Pfeffer, Reference Pfeffer2007).

At present, however, the influence of proglacial lakes on ice mass loss and SLR are either absent or poorly represented in ice sheet models (Carrivick and others, Reference Carrivick, Tweed, Sutherland and Mallalieu2020). Since proglacial lakes modify ice dynamics, by altering terminus profile, subglacial hydrology and local force balance (e.g. Warren and Kirkbride, Reference Warren and Kirkbride2003; Sugiyama and others, Reference Sugiyama2011; Baurley and others, Reference Baurley, Robson and Hart2020), there is a need to determine the extent to which proglacial lakes will impact ice sheet mass loss over the coming century and beyond (Carrivick and others, Reference Carrivick2022).

There are many subglacial bedrock overdeepenings beneath the GrIS (Patton and others, Reference Patton, Swift, Clark, Livingstone and Cook2016; Morlighem and others, Reference Morlighem2017) which fill with meltwater runoff during margin recession, forming ice-marginal proglacial lakes (Costa and Schuster, Reference Costa and Schuster1988; Carrivick and Tweed, Reference Carrivick and Tweed2013). In recent decades, the area and number of ice-marginal lakes has increased in south-west Greenland, as well as globally (Carrivick and Quincey, Reference Carrivick and Quincey2014; Shugar and others, Reference Shugar2020; How and others, Reference How2021; Rick and others, Reference Rick, McGrath, Armstrong and McCoy2022). Such changes in ice margin configuration are significant because outlet glaciers that terminate in lakes typically have greater rates of mass loss and terminus retreat than their land-terminating counterparts (Kirkbride, Reference Kirkbride1993; Warren and Kirkbride, Reference Warren and Kirkbride2003; Schomacker, Reference Schomacker2010; Tsutaki and others, Reference Tsutaki, Nishimura, Yoshizawa and Sugiyama2011; Carr and others, Reference Carr, Bell, Killick and Holt2017; King and others, Reference King, Dehecq, Quincey and Carrivick2018; Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Raby2021). This change in dynamics reflects differences in the boundary conditions at (a) the bed, and (b) the terminus (Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). The presence of a lake leads to a reduction in the effective pressure at the terminus, and up-glacier, enabling greater rates of basal sliding (Bindschadler, Reference Bindschadler1983; Benn and others, Reference Benn, Warren and Mottram2007; Sugiyama and others, Reference Sugiyama2011). Effective pressure is the difference between ice-overburden pressure and basal water pressure, hence thinner ice at the terminus and or deeper lake water will promote greater ice velocities (Kirkbride and Warren, Reference Kirkbride and Warren1997; Tsutaki and others, Reference Tsutaki, Sugiyama, Nishimura and Funk2013). Lake depth at the terminus sets the base level and thus minimum basal water pressure for the near-terminus subglacial hydraulic system, and the influence of the pressure head due to the lake declines with distance from the terminus at a rate dependent on the bed topography (Meier and Post, Reference Meier and Post1987; Benn and others, Reference Benn, Warren and Mottram2007). Water at the terminus also initiates a suite of complementary processes: sub-aqueous melt, thermal-notch erosion, and calving (both sub-aerial and sub-aqueous) (Röhl, Reference Röhl2006; Sugiyama and others, Reference Sugiyama2011, Reference Sugiyama, Minowa and Schaefer2019; Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Smith2020). Water-depth impacts terminus buoyancy with calving rates increasing as water depth increases (Benn and others, Reference Benn, Warren and Mottram2007; Boyce and others, Reference Boyce, Motyka and Truffer2007; Dykes and others, Reference Dykes, Brook, Robertson and Fuller2011). These processes may contribute to terminus retreat, surface steepening, and further increases in velocity and longitudinal strain rates (e.g. Warren and Kirkbride, Reference Warren and Kirkbride2003; Tsutaki and others, Reference Tsutaki, Sugiyama, Nishimura and Funk2013; King and others, Reference King, Dehecq, Quincey and Carrivick2018). Furthermore, modelling results suggest that glacier response to the development of proglacial lakes may be partially decoupled from short-term changes in climate and contribute to rapid and sustained retreat (Sutherland and others, Reference Sutherland2020).

The effect of proglacial lakes, and their development, on ice dynamics has been observed in many glaciated regions. For example, at Breiðamerkurjökull, Iceland (which has both lake- and land-terminating distributaries), between 1991–2015 there was no change in ice velocity adjacent to the land-terminating margins, whereas velocity increased by a factor of three to 3.5  m day−1 proximate to the terminus of the lake-terminating arm (Baurley and others, Reference Baurley, Robson and Hart2020). This change was linked to increases in surface air temperature initiating terminus retreat into a 200–300 m deep subglacial trough, triggering a positive feedback mechanism (Pfeffer, Reference Pfeffer2007; Nick and others, Reference Nick, Vieli, Howat and Joughin2009). Retreat into deeper water enabled ice flow acceleration, dynamic thinning (in addition to thinning from changes in surface mass balance), greater rates of calving, and further retreat into deeper water.

The contrasting patterns in ice dynamics between lake- and land-terminating glaciers have also been observed in the Himalaya (e.g. King and others, Reference King, Bhattacharya, Bhambri and Bolch2019; Tsutaki and others, Reference Tsutaki2019; Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). For example, greater rates of mass loss and terminus retreat have been observed at lake-terminating outlets (King and others, Reference King, Bhattacharya, Bhambri and Bolch2019), and their centreline velocities are typically double those measured at land-terminating glaciers (18.83 vs. 8.24 m yr−1, for the period 2017–2019) (Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). One notable difference between lake- and land-terminating glaciers is the centreline velocity profile, whereby lake-, like marine-terminating outlet glaciers, accelerate toward the terminus (Tsutaki and others, Reference Tsutaki2019; Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). This extensional flow contributes to surface thinning, whereas the compressive flow regime at land-terminating glaciers can lead to thickening, which may offset surface mass balance induced surface lowering (Tsutaki and others, Reference Tsutaki, Sugiyama, Nishimura and Funk2013, Reference Tsutaki2019). Additionally, glacier geometry influences both terminus stability and velocity, with valley constrictions or submerged sills and their associated impacts on lateral and back-stress causing reductions in velocity and enhanced terminus stability (van der Veen and Whillans, Reference van der Veen and Whillans1989; O'Neel and others, Reference O'Neel, Pfeffer, Krimmel and Meier2005; Benn and others, Reference Benn, Warren and Mottram2007). These observations are supported by sensitivity modelling which suggests that thicker ice, a wider terminus and steeper surface slopes lead to elevated near-terminus velocities at lake-terminating glaciers (Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021).

Lake-terminating glacier dynamics can also be affected by surface meltwater input to the glacier bed (Sugiyama and others, Reference Sugiyama2011). For example, at Glacier Perito Moreno, a lake-terminating glacier in Patagonia, Argentina, hourly variations in measured basal water pressures 4 km from the terminus correspond closely with changes in surface temperature and ice velocity measured in-situ using dGPS (Sugiyama and others, Reference Sugiyama2011). These observations suggest surface meltwater can reach the bed rapidly, and that ice velocity is sensitive to small changes in basal water pressures. This sensitivity is evidenced by the observed difference in relative change: over a ten day period ice velocity varied by 37% about its mean (1.43 m day−1), whereas basal water pressures only varied by 5%. This corresponded to an increase in velocity of 0.053 m day−1 per 1 $^\circ$C (Sugiyama and others, Reference Sugiyama2011).

In addition to the above geometric and climatic controls, the characteristics of individual proglacial lakes will influence the behaviour of lake-terminating glaciers (e.g. Sugiyama and others, Reference Sugiyama2016; Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Smith2020; Watson and others, Reference Watson2020; Dye and others, Reference Dye, Bryant, Dodd, Falcini and Rippin2021). For example, a reduction in calving frequency and volume corresponds to the timing of lake-ice freeze up (Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Smith2020). Additionally, sub-aqueous melt and associated thermal-notch erosion are both functions of the thermal structure of the lake (e.g. Haresign and Warren, Reference Haresign and Warren2005; Röhl, Reference Röhl2006; Minowa and others, Reference Minowa, Sugiyama, Sakakibara and Skvarca2017). Observations in Patagonia (Sugiyama and others, Reference Sugiyama2016, Reference Sugiyama, Minowa and Schaefer2019) revealed a layer of cold turbid water derived from subglacial discharge underlying warmer surface waters. This stratification can prohibit the upwelling of meltwater that is seen in glacial fjords, and allows for the formation of ice terraces below the waterline (Kirkbride and Warren, Reference Kirkbride and Warren1997; Sugiyama and others, Reference Sugiyama, Minowa and Schaefer2019). The thermal state of a proglacial lake is strongly coupled to climate, and is dependent on incident shortwave radiation, surface air temperatures, winds, precipitation and runoff (e.g. Schomacker, Reference Schomacker2010; Richards and others, Reference Richards, Moore and Forrest2012).

The observed impacts of proglacial lakes on ice dynamics (Kirkbride, Reference Kirkbride1993; Warren and Kirkbride, Reference Warren and Kirkbride2003; Sugiyama and others, Reference Sugiyama2011; Baurley and others, Reference Baurley, Robson and Hart2020; Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021) are furthermore supported by modelling studies (e.g. Sutherland and others, Reference Sutherland2020). Collectively, these works suggest the presence of proglacial lakes leads to greater rates of terminus retreat and mass loss than those at land-terminating glaciers, thereby contributing to accelerated rates of deglaciation.

Given the clear potential of lakes to perturb ice dynamics (e.g. Kirkbride, Reference Kirkbride1993), and the projected prevalence of ice-marginal lakes in Greenland (Carrivick and others, Reference Carrivick2022), it is important to evaluate how, against a backdrop of a warming climate, lakes are impacting ice motion and terminus positions. In south-west Greenland mean annual changes in margin position have varied since the 1980s, with notable differences between the lake- and land-terminating sectors. Average annual rates of margin recession increased by an order of magnitude from 1.1  m yr−1 (1987–1992) to 11.5  m yr−1 (2010–2015) along lacustrine margins, whereas the magnitude of changes at terrestrial margins was more modest: from advance of 1.2  m yr−1 to recession of 2.8  m yr−1 (Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Raby2021). Furthermore, observations from across Greenland indicate ice-marginal lakes enhance the flow of adjacent ice by ~25% (Carrivick and others, Reference Carrivick2022).

This study aims to investigate recent (2013–2021) changes in ice velocity, surface elevation, and terminus retreat at two proximate but contrasting outlet glaciers in south-west Greenland; the lake-terminating Isortuarsuup Sermia, and the nearby land-terminating Kangaasarsuup Sermia. We additionally consider the varying processes that impact dynamics at this lake-terminating system, and the potential significance of these at the ice-sheet scale.

2. Data and methodology

2.1 Study site

We generated ice velocities, rates of surface elevation change, and terminus position from two outlet glaciers in south-west Greenland: Isortuarsuup Sermia (IS) (63$^\circ 50'$ N, 49$^\circ 59'$ W) and Kangaasarsuup Sermia (KS) ($64^\circ 07'$ N, $49^\circ 54'$ W) (Fig. 1). Isortuarsuup Sermia is a lake-terminating glacier that drains into one of the largest proglacial lakes in south-west Greenland, the ~60 km2 Isortuarsuup Tasia, and KS is a nearby land-terminating glacier. These outlet glaciers were selected due to their close proximity (~30 km) and their similar morphological characteristics including terminus elevation (~315 m at KS and ~500 m at IS), surface slope, and valley width. Based on meltwater routing via subglacial hydraulic potential (Shreve, Reference Shreve1972), KS drains 660  km2, whereas IS drains 122 km2 (Mankoff, Reference Mankoff2020) (Fig. 1a).

Figure 1. (a) The lake-terminating Isortuarsuup Sermia (IS) and land-terminating Kangaasarsuup Sermia (KS) with their respective centrelines (solid lines) and runoff catchments from Mankoff and others (Reference Mankoff2020) (dashed lines). Ice surface contours at 250 m intervals generated from BedMachine v5 (Morlighem and others, Reference Morlighem2022); background image: Sentinel-2 acquisition from 19 September 2022 (ESA Copernicus, 2022); Location of study area, south-west Greenland, shown by pink box in inset (upper right). Dashed red box at terminus of IS denotes extent of (b); (b) terminus region of IS illustrating the trim-line (white arrows), terminal moraine (black arrow), grounded icebergs (orange arrow). Sentinel-2 acquisition from 18 September 2019 (ESA Copernicus, 2022).

There is clear evidence of a distinct terminal moraine at IS (black arrow in Fig. 1b), which likely formed during a period of prolonged stability, as evidenced by the pronounced trim-line (Fig. 1b) inferred to be of Little Ice Age origin (during the 19th century) by Weidick and others (Reference Weidick, Bennike, Citterio and Nørgaard-Pedersen2012). Furthermore, icebergs are often grounded in the lake approximately 400 m from the glacier terminus suggesting a sublacustrine extension of the visible terminal moraine (Fig. 1b).

2.2 Ice velocity

Ice velocities were obtained from the NASA MEaSUREs ITS_LIVE version 2 data-cubes (Gardner and others, Reference Gardner2018, Reference Gardner, Fahnestock and Scambos2022). This data set is derived from optical (Landsat 8 and Sentinel-2A/B) and radar acquisitions (Sentinel-1A/B). Velocities are determined using the autonomous repeat image feature tracking algorithm applied to pairs of overlapping images from a given sensor, separated by time (date_dt) (Gardner and others, Reference Gardner2018; Lei and others, Reference Lei, Gardner and Agram2021). Velocity fields obtained from image pairs with small time separations (date_dt ≤ 30 days) reveal short-term changes in velocity while long time separations (date_dt ≥ 300 days) provide better estimates of annual averages. Calculated velocities are posted onto a uniform 120 m grid, with a spatially variable effective resolution of 240–1920 m (Lei and others, Reference Lei, Gardner and Agram2021). The resulting data-cube has dimensions easting (x), northing (y), and time (t), where t is the mid-date between the two satellite acquisitions used to generate the velocity field. Each grid cell contains the velocity component in the x (v x) and y (v y) directions, and image-pair time dependent error estimates ($\sigma _{v_t}$) are also supplied. To minimise the effects of point sampling, and to allow for the spatially variable effective resolution, which is a function of the search window size used when feature tracking (Lei and others, Reference Lei, Gardner and Agram2022), the v x and v y fields were first spatially averaged using a 3 × 3 window, and subsequently sampled every 250  m along the glacier centrelines, and the resultant velocity calculated (Eqn (1)):

(1)$$v_t = \sqrt{\bar{v_{x_t}}^2 + \bar{v_{y_t}}^2}$$

Error-weighted average velocities ($\bar {v}$) were determined (Eq. (2)):

(2)$$\overline{v} = {\sum_{t}{{\bar{v_t}}/ \sigma_{v_t}^2}\over \sum_{t}{1/\sigma_{v_t}^2}}$$

With the uncertainty in $\bar {v}$ calculated as (Eq. (3)):

(3)$$\sigma_{v} = \sqrt{{1\over \sum_{t}{1/\sigma_{v_t}^2}}}$$

The error-weighted average annual velocities were calculated using all velocity fields derived from image-pairs separated by between 300 and 430 days, with assignment to a given year on the basis of the mid-date of the image-pair. Annual average velocity profiles were constructed for 2013–2021, and the velocity trends computed from these averages for each point along the centreline using linear regression. To further assess differences in velocity regime, seasonal averages were constructed from velocity fields where date_dt ≤ 30 days, and rates of acceleration along the centreline determined for each season. Seasons were defined as: winter (December, January, February), spring (March, April, May), summer (June, July, August), and autumn (September, October, November). Due to the limited availability of velocity fields with small image separations prior to 2016, the seasonal analysis presented here is confined to the period 2016–2021. For both the annual and seasonal trend analysis, the null hypothesis that there is no trend (i.e. a regression slope coefficient of zero), was evaluated using a two-tailed Wald test, as implemented in the scipy python package. Following convention, when the returned p-value was ≤ 0.05, the trends are taken to be significant.

2.3 Surface elevation change

Rates of surface elevation change were derived from ArcticDEM 4.1 (Porter and others, Reference Porter2022). ArcticDEM is a collection of high resolution (2 m) time dependent (2007–2021) digital elevation models (DEMs). These are constructed from stereo auto-correlation methods (Noh and Howat, Reference Noh and Howat2015), applied to sub-metre resolution optical satellite imagery from the Maxar constellation. The absolute accuracy of individual ArcticDEM strips is approximately 4 m both horizontally and vertically (Porter and others, Reference Porter2022). The volume of data processed in the construction of ArcticDEM allows the accidental inclusion of errors in the data set (Błaszczyk and others, Reference Błaszczyk2019). Despite this, once DEMs have been co-registered, accuracy has been shown to be, in many cases, better than 4 m and with precision of approximately 1  m, making ArcticDEM suitable for measuring changes ≥ 1 m (Błaszczyk and others, Reference Błaszczyk2019).

To co-register ArcticDEMs, a 5 km buffer was constructed for each glacier centreline, and all ArcticDEM strips that intersected this region were selected. To avoid detecting changes in elevation as a function of seasonal snow cover, the list of DEMs was filtered to include only those constructed from images acquired in June, July, August and September. After this filtering, there were 20 DEMs spanning the period September 2012–June 2021 over the study region at IS, and 51 (June 2011– September 2021) at KS. The supplied bit-mask was applied to each DEM to preserve only those pixels marked as good data; cloud, water and edge pixels were masked. All DEMs were visually inspected and the DEM with the best coverage of each glacier was selected to be the reference DEM, to which all others were co-registered. At IS the reference DEM was from 15/06/2016; at KS 04/08/2014. A cloud-free Sentinel-2 scene (18/08/2022 at IS, and 24/07/2022 at KS) was used to identify regions of stable terrain (SI. 1). During the co-registration process, differences in elevation over these regions of stable terrain are minimised, enabling changes in ice surface elevation to be observed accurately. Each DEM was co-registered to the reference using the method proposed by Nuth and Kääb (Reference Nuth and Kääb2011), as implemented in the python package XDEM (Xdem Contributors, 2021). Once co-registered each DEM was down-sampled from 2 m to 20 m using bi-linear interpolation, to reduce both file size and the effects of point sampling. To provide a measure of width-averaged rates of thinning, the co-registered DEMs were sampled every 50 m along a series of parallel lines spaced every 250 m. Rates of change were computed from these elevation samples using linear regression. As per the velocity trend analysis, significant trends were evaluated using a two-tailed Wald test. For each co-registered DEM a quadratic surface (Eq. (4)) was fitted to elevation values (z) within a 21 × 21 (420 × 420 m) window using least squares regression, and the surface slope (S) calculated from the coefficients (Eq. (5)) (e.g. Hurst and others, Reference Hurst, Mudd, Walcott, Attal and Yoo2012). These slope surfaces were also sampled every 100 m along offset parallel lines, and the width-averaged (median) change in slope between the earliest and latest DEM calculated.

(4)$$z = Ax^2 + Bx + 2 + Cxy + Dx + Ey + F$$
(5)$$S = \sqrt{d^2 + e^2}$$

Uncertainty in the estimated rates of surface elevation change was minimised through the co-registration process. Prior to co-registration, the median difference in elevations over stable terrain ranged between −2.2 and 5.0 m at IS and −8.6 and 2.0 m at KS. After co-registration, the median differences in elevation over stable terrain were much closer to zero (−0.02 − 0.02 m, and −0.03 − 0.01 m, at IS and KS, respectively). The normalised median absolute deviation (NMAD), which measures a sample's dispersion, of elevation differences was similarly reduced in the co-registration process (SI. 2) At IS, NMADs were reduced from between 0.41–1.71 m to 0.24–0.69 m, with similar improvements at KS (0.34–2.41 m, to 0.23–1.97 m). We are therefore confident in our ability to detect changes in elevation ≥ 1 m, which is equivalent to annual rates of change of ~0.1 m yr−1 over the period which ArcticDEM is available.

2.4 Terminus positions

Terminus positions were manually digitised from optical imagery captured by Landsat 8 and Sentinel-2 using the Google Earth Engine Digitization Tool (Lea, Reference Lea2018) from 2013 to 2022. Relative changes in terminus position were determined using the rectilinear box method (Moon and Joughin, Reference Moon and Joughin2008) which allows for uneven retreat across the terminus. Uncertainty in relative terminus positions arise from image co-registration errors and manual digitisation errors (e.g. Carr and others, Reference Carr, Stokes and Vieli2014). Image co-registration errors are a function of poor spatial alignment between satellite image acquisitions. The Landsat scenes used in this study had a median image registration accuracy of 4.6 m, and the Sentinel-2 technical reference indicates geolocation uncertainties are ≤ 11 m (S2 MSI ESL Team, 2022). To quantify the digitisation precision, termini were re-digitised five times (Paul and others, Reference Paul2013), and the standard error was determined to be 16.7 m, which yields an uncertainty in rate of terminus position change of ±3.7 m yr−1 over the study period, and is a lower bound on the measurable rate of terminus position change. At IS, digitising errors arose due to difficulties in discriminating between the terminus and either recently calved icebergs or lake ice cover. Digitising errors at KS were principally due to snow cover at the beginning of the melt season, debris cover at the end of the melt season, and shadows cast by the ridge to the south. As such, for KS, the results presented here are average relative terminus positions over each summer.

2.5 Runoff

The runoff data set used (Mankoff and others, Reference Mankoff2020) comprises liquid water discharge estimates at hydrological outlets derived from two regional climate models: MAR (Modele Atmospherique Regional) and RACMO (Regional Atmospheric Climate Model) (Noël and others, Reference Noël2016; Fettweis and others, Reference Fettweis2017). The subglacial stream network is determined from ice surface elevations (Porter and others, Reference Porter2018) and ice thickness estimates from BedMachine (Morlighem and others, Reference Morlighem2017) using a model for the subglacial pressure head (Shreve, Reference Shreve1972). This data set accounts for surface melt, rainfall, meltwater retention and refreezing. Supraglacial flow is discounted and meltwater is assumed to generate and reach the bed within the same model grid cell (Mankoff and others, Reference Mankoff2020). At the study site, the basin output closest to each glacier termini was selected and cumulative runoff calculated for the years 2011–2021. Linear regression was used to evaluate trends in cumulative annual runoff, for both MAR and RACMO, and tested for significance using a two-tailed Wald test. Similarly, to assess the relationship between runoff and ice velocity, average annual ice velocity was regressed against cumulative annual runoff (2013–2021) (derived from the mean average of MAR and RACMO).

This data set is not supplied with uncertainty estimates, however it contains three principal sources of uncertainty. (1) Temporal uncertainty which is a function of how the routing model handles the time lag between meltwater generation and runoff within each grid cell; (2) basin uncertainty, which arises from the data used in computing the subglacial stream network and catchments; and (3) uncertainties in the regional climate models from which the liquid water discharge estimates are derived (Mankoff and others, Reference Mankoff2020). Here, temporal uncertainties are mitigated by computing cumulative annual totals.

3. Results

3.1 Ice velocity

Ice velocities show clear and contrasting patterns in behaviour at the two outlet glaciers between 2013 and 2021 (Fig. 2). Annual average ice velocity within 500 m of the terminus more than doubled from ~80 to 220 m yr−1 at IS, while there was minimal change at KS (from ~25 to 20 m yr−1) (Fig. 2). The magnitude of the acceleration at IS decreases from 15.0 ± 2.4 m yr−2 at the terminus to 1.4 ± 0.5 m yr−2 15 km up-glacier (Fig. 2e); furthermore, the increase in velocity extends across the full width of the terminus (Fig. 3a). Notably, between 2018 and 2019, there was a substantial (~30%) increase in near terminus surface velocity at IS from 147 to 193  m yr−1 (Fig. 2a & c). By contrast, barring a region of low magnitude (1.0–1.9 m yr−2) acceleration 2.5–8.5 km from the terminus, there is no discernible trend in ice velocity at KS (Figs. 2f and 4a). The differences in velocity magnitude and velocity trend between these two glaciers declines with distance from the terminus, such that at a distance of ~15 km from their respective termini, differences in the average annual velocity change over the study period at IS (from 127 to 141 m yr−1) and KS (from 117 to 125  m yr−1) are minimal (Figs. 2 and 3).

Figure 2. Annual average ice velocity (2013–2021) profiles along centrelines shown in Figure 1 at (a) Isortuarsuup Sermia, and (b) Kangaasarsuup Sermia. (c) and (d) show percentage change relative to 2013. (e) and (f) show linear trends where regression slope coefficients are significant at p $\leq 0.05$; error bars denote 95% confidence interval.

Figure 3. Isortuarsuup Sermia: (a) change in average annual velocity between 2013–2021; (b) rate of surface elevation change (September 2012– June 2021) from ArcticDEM (negative denotes thinning); manually digitised ice margin shown in black in (a) and (b); (c) terminus positions 2014–2021. Red box in (a) denotes extent of (c). White arrows in (c) indicates the Little Ice Age trim-line and the black arrow points to the associated terminal moraine with icebergs grounded on its sublacustrine extension indicated by the orange arrow.

Figure 4. Kangaasarsuup Sermia: (a) change in average annual velocity between 2013–2021; (b) rate of surface elevation change (June 2011– September 2021) from ArcticDEM (negative denotes thinning); manually digitised ice margin shown in black in (a) and (b); (c) terminus positions 2014–2021. Red box in (a) denotes extent of (c).

Seasonal cycles were observed at both glaciers, however, there were key differences between IS and KS with regards to seasonal velocity trends over the five year period 2016–2021 (Figs. 5 and 6). At IS, along the entire lower 10 km of the glacier, winter ice surface velocity increased, with the magnitude of acceleration declining from 20.2 (±12) m yr−2 at the terminus to 4.4 (±4.3)  m yr−2 at 10 km (Fig. 6). There were also significant positive trends in autumn and summer within the 2 km closest to terminus, although these were of slightly lower magnitude than those in winter. Further up-glacier (6–9  km), summer ice surface velocity also accelerated (~10 m yr−2). Conversely, at KS, where velocities decline toward the terminus (Fig. 2b), there have been minimal changes in seasonal velocity with just a few locations along the centreline exhibiting statistically significant trends. Specifically, in autumn at ~4 km from the terminus, velocities were decreasing at approximately 6.2 m yr−2 (Fig. 6). The absence of clear changes in seasonal velocities at KS is consistent with the minimal change in annual velocity (Figs. 2d and 4a).

Figure 5. Time series of ice surface velocity at (a) 1  km and (b) 10 km from the terminus at the lake-terminating Isortuarsuup Sermia (blue) and the land-terminating Kangaasarsuup Sermia (orange, dashed). These velocities are computed from image-pairs separated by ≤ 30 days. Lines show the rolling 28 day median.

Figure 6. Seasonal velocity trends (2016–2021) along glacier centrelines at Isortuarsuup Sermia (blue) and Kangaasarsuup Sermia (orange) for winter (DJF, circles), summer (JJA, triangles) and autumn (SON, crosses). Seasonal trends derived from velocity fields where date_dt ≤ 30 days. Trends are linear fits, and only significant (p ≤ 0.05) trends are shown; error bars denote 95% confidence interval. Points at the same distance from the terminus have been offset from one another to aid readability.

3.2 Surface elevation change

There is a clear thinning signal at both IS and KS since 2011 (Figs. 3 and 4) with the rate and magnitude of thinning increasing towards both termini (Fig. 7). At ~15 km from the terminus, KS was thinning at a rate of 0.8 ± 0.3 m yr−1, and IS at 0.3 ± 0.2 m yr−1. These increased to their greatest width-averaged rates of thinning of 3.1 ± 0.7 m yr−1 at the terminus of KS, and 2.1 ± 0.6 m yr−1 900 m from the terminus at IS, where the elevation is 550 m. However, at the equivalent altitude at KS, the thinning rate is ~1.1  m yr−1 and accounting for differences in elevation, and the associated changes in lapse rate and thus surface melt processes, rates of thinning at IS are typically between 0.33–0.65 m yr−1 (interquartile range; median difference: 0.5 m yr−1) greater than at KS (Fig. 7c). The net differences in elevation change (up to −21 m and −13 m, at KS and IS, respectively) are substantially more than the differences measured over stable terrain after co-registration, giving us confidence in these observations.

Figure 7. Rates of surface elevation change at (a) Isortuarsuup Sermia (September 2012–June 2021) and (b) Kangaasarsuup Sermia (June 2011– September 2021) from ArcticDEM. Rates determined from linear regression. Only significant trends (p $\leq 0.05$) are shown. Error bars denote 95% confidence interval. ArcticDEM was sampled every 100 m along 9(7) parallel lines spaced every 250 m across the glacier at IS(KS). Coloured lines show 500 m rolling width-averaged median. Purple line (right hand axis) denotes number of ArcticDEMs with valid elevation measurements at each point along centreline. Shading represents standard deviation of number of ArcticDEMs at each point along centreline to account for the parallel offsets. (c) Rate of surface elevation change from ArcticDEM, against surface elevation at Isortuarsuup Sermia (blue) and Kangaasarsuup Sermia (orange); coloured lines show rolling median over 50 m bins. In (a), (b) and (c) negative values indicate surface thinning. (d) Net change in surface slope (positive indicates surface steepening) between first and last DEM; line represents median change across parallel lines, shading denotes median absolute deviation.

Due to greater rates of thinning at lower elevations, changes in surface gradient are generally positive at both IS and KS. Whilst these increases are of low magnitude (Fig. 7d), they do indicate some surface steepening. At IS, the median change in gradient along the centreline was $0.03\pm 0.2^\circ$, and at KS was $0.08\pm 0.3^\circ$. Immediately proximate to the terminus at IS, there is evidence of a lessening of the surface gradient toward the terminus, whereas at KS surface slopes increased steadily over the lower ~3 km by approximately 0.5$^\circ$ (Fig. 7d). This change in gradient at IS corresponds with a clear decrease in rate of surface elevation change over the 1 km closest to the terminus (Fig. 7a)

3.3 Terminus position

Terminus retreat differed between the two outlets over the course of the study period (Fig. 8). Between August 2014 and September 2021 KS retreated 210 ± 46 m, at an average rate of 30 ± 4 m yr−1, which is greater than the uncertainty in our method. By contrast, IS retreated more slowly at a rate of 9 ± 4 m yr−1.

Figure 8. Relative terminus position at Isortuarsuup Sermia (blue) and KS (orange) July– September 2014–2022. Blue circles denote relative position measured using box method; dashed blue lines representative of winter advance. Orange line illustrates average annual relative terminus position, shading denotes 95% confidence interval.

A seasonal cycle is observed at IS with a median winter advance of 30 m (median absolute deviation (MAD): 18 m), and median summer retreat of 19 m (MAD: 12 m). Furthermore, the terminus at IS showed distinct across-glacier spatial variability with the most pronounced localised retreat (~200  m) occurring in 2019 when a more advanced section of the northern terminus retreated from the sublacustrine terminal moraine (orange arrow in Fig. 3c). Following retreat from this moraine, the glacier did not re-advance back on to it during the remainder of our observation period. By contrast, retreat at KS appears progressive and sustained.

3.4 Runoff

Modelled runoff showed liquid water discharge (2011–2021) at KS was greater than that at IS by a factor of ~3 (Fig. 9), which is consistent with its larger catchment. At IS, there is good agreement between the two climate models whereas at KS, cumulative annual runoff is 5–25% greater in RACMO. Annual peaks in average daily runoff were typically between 60 and 90 m3 s−1 at IS, and 230–320 m3 s−1 at KS. Linear regression of cumulative annual runoff, from both RACMO and MAR, against time showed no significant trend at either IS or KS (SI. 3). Over the study period there was no consistent change in the timing of runoff onset or cessation. Runoff typically started between the end of April and mid-May at KS, and a week later at IS, as expected given its higher elevation, and had generally ceased by early- and mid-October at IS and KS, respectively.

Figure 9. Cumulative (a & b) and average daily runoff rate (5-day rolling average) (c & d) at Isortuarsuup Sermia (a & c) and Kangaasarsuup Sermia (b & d). Colours denote regional climate model with MAR in turquoise, and RACMO in orange.

4. Discussion

The near-terminus increase in ice surface velocity at IS (Fig. 2a) is similar to the accelerations seen at many other lake- and marine-terminating outlet glaciers in recent years (e.g. Joughin and others, Reference Joughin, Smith and Howat2018; Baurley and others, Reference Baurley, Robson and Hart2020). Our findings suggest that the presence of water at the terminus of IS, and the associated effects on near-terminus force balance, has enabled the observed changes in ice dynamics. This is reflected in the shape of the surface velocity profile, which for any given year shows velocity increasing towards the terminus from ~10  km up-glacier (Fig. 2a), and the significant increase in velocity over time along the entire 15 km centreline. We consider below evidence for potential drivers of the dynamic changes observed at IS in contrast to the behaviour at the neighbouring land-terminating KS.

In agreement with previous work (e.g. Tedstone and others, Reference Tedstone2015), these data show no statistically significant relationship between cumulative runoff and average annual velocity at either glacier (SI. 3), and no significant melt trend was observed, either with respect to melt volume or timing. This result does not preclude a relationship between runoff and ice velocity over shorter timescales than are frequently measured from remotely sensed observations of velocity and gridded runoff estimates derived from regional climate models. Indeed, clear seasonal velocity cycles are evident at both glaciers, reflecting the commonly observed coupling between seasonal runoff, the hydraulic evolution of the subglacial drainage system and ice-dynamics (Davison and others, Reference Davison, Sole, Livingstone, Cowton and Nienow2019). Nevertheless, the absence of a relationship between annual runoff and ice velocity at the lake-terminating IS is important, and suggests that the observed acceleration is not directly attributable to enhanced sliding due to increased meltwater input to the bed.

The observed acceleration is therefore likely the result of a dynamic feedback (Benn and others, Reference Benn, Warren and Mottram2007) driven by sustained negative surface mass balance induced thinning (The IMBIE Team, 2020), which is seen at both IS and KS (Fig. 7). Thinning will have brought the terminus closer to flotation and enhanced rates of basal sliding (Pfeffer, Reference Pfeffer2007; Tsutaki and others, Reference Tsutaki2019). This suggestion is supported by both the glacier-wide acceleration (Fig. 3a) and by the acceleration in all seasons (Fig. 6) at IS. Furthermore, due to minimal changes in lake water level, basal water pressures proximate to the terminus will be held approximately constant while ice-overburden pressure will decrease due to surface thinning (Figs. 3 and 7), leading to a long-term decrease in effective pressure. While a sustained decrease in effective pressure promotes the observed increase in ice motion at IS (Fig. 2a), a pronounced acceleration at the terminus occurs between 2018 and 2019 (Fig. 2b). This results from ongoing thinning at the glacier terminus promoting flotation and the subsequent retreat of the northern part of the terminus away from the Little Ice Age sublacustrine moraine (Fig. 3c). The timing of this flotation and retreat is further evidenced by the substantial change in rate of surface elevation change post-2018 (SI. 4) in conjunction with the pronounced increase in velocity, presumably in response to the associated removal of buttressing. The change in surface elevation suggests a clear hinge (above the grounding line) ~1 km up-glacier from the terminus, akin to what has been observed at Helheim Glacier (James and others, Reference James, Murray, Selmes, Scharrer and O'Leary2014). Furthermore, acceleration following retreat from a sublacustrine moraine replicates the behaviour observed at the lake-terminating Yakutat Glacier, Alaska (Trüssel and others, Reference Trüssel, Motyka, Truffer and Larsen2013).

The acceleration and associated extensional flow has led to enhanced rates of surface lowering (2.1 ± 0.6 m yr−1) across the whole width of the terminus region (Figs. 3b and 7c). These spatially variable rates of thinning change the glacier surface slope, which have generally steepened up-glacier (Fig. 7d), likely leading to an increase in driving stresses, which may promote further acceleration and thinning (Howat and others, Reference Howat, Joughin, Tulaczyk and Gogineni2005). This is consistent with modelling work that demonstrated near-terminus velocities at lake-terminating glaciers increased with surface slope (Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). The absence of bed topography estimates near the terminus inhibit efforts to quantify stress-coupling lengths (e.g. Enderlin and others, Reference Enderlin2016); however, it is suggested here that a high degree of longitudinal coupling allows these thinning-driven near-terminus accelerations to propagate ~15  km up-glacier (Fig. 2e).

The observed rates of change in terminus position at these two outlets (Fig. 8) are consistent with those previously documented (Warren, Reference Warren1991). Between 1943 and 1983 KS retreated at an average rate of 38 m yr−1, whilst the terminus at IS was shown to be stable (1949–1985). Additionally, there has only been ~400 m of retreat at IS since the Little Ice Age during which time KS has retreated approximately 3.4 km (Weidick and others, Reference Weidick, Bennike, Citterio and Nørgaard-Pedersen2012). This observation of greater retreat at the land-terminating KS, as opposed to the lake-terminating IS, differs from the regional average pattern along the entire south-west margin where recent rates of margin recession were typically greater along lacustrine sections (mean annual rates of margin change: −11.5 m yr−1, 2010–2015) than land-terminating margins (−2.8 m yr−1) (Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Raby2021). Nevertheless, this anomalous result is not unexpected given that the range of annual rates of margin recession at lacustrine (n=374) and terrestrial margins (n=3325) are approximately equivalent (Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Raby2021, Figure 4b), and we are only presenting results from two such points. Furthermore, sustained long term stability at IS is evidenced by the minimal retreat and ongoing proximity of the terminus to the Little Ice Age maximum (Weidick and others, Reference Weidick, Bennike, Citterio and Nørgaard-Pedersen2012) (leftmost white arrow in Fig. 3c) and the large sublacustrine moraine. This suggests that the topographic configuration at IS has enabled this stable terminus position in a manner similar to the long-term stability observed at other lake- (Trüssel and others, Reference Trüssel, Motyka, Truffer and Larsen2013) and marine- (Catania and others, Reference Catania218) terminating glaciers. For example, the bed topography at Store Glacier (e.g. Catania and others, Reference Catania218; Box and Decker, Reference Box and Decker2011) has promoted stability at its terminus, contrary to the regional trend, and in spite of being in sustained negative balance. Additionally, with respect to terminus position at IS, the recent dramatic increase in ice surface velocity (Figs. 2a and 3a), may in part offset any retreat from frontal ablation.

At Breiðamerkurjökull, Iceland, between 1982 and 2018 there was terminus retreat and a corresponding increase in proglacial lake area at Jökulsárlón and Breiðárlón (Baurley and others, Reference Baurley, Robson and Hart2020). However, the net retreat at these two proximate lake-terminating margins differed by a factor of three. Furthermore, and in line with our observations, cumulative retreat at Breiðárlón over this period was less than at an adjacent land-terminating section. These authors suggest that the relative stability of the terminus at Breiðárlón is attributable to the shallow subglacial trough (60 m vs. 300 m at Jökulsárlón). At present, there are no available lake bathymetry estimates at Isortuarsuup Tasia, and ice thickness estimates from BedMachine v5 (Morlighem and others, Reference Morlighem2022) are not available proximate to either terminus. Additionally, the long-term stability of the terminus at IS (Weidick and others, Reference Weidick, Bennike, Citterio and Nørgaard-Pedersen2012) suggests that bed topography, lateral support from the valley sides and the development of the large terminal moraine, have all exerted a strong control on terminus position (Warren, Reference Warren1991).

The differences in velocity between IS and KS support recent work demonstrating that ice-contact lakes amplify near-terminus velocities (Baurley and others, Reference Baurley, Robson and Hart2020; Pronk and others, Reference Pronk, Bolch, King, Wouters and Benn2021). In Greenland, this velocity uplift has been estimated to be ~25% (Carrivick and others, Reference Carrivick2022). However, this value of 25% contrasts ice-motion adjacent to all styles of ice-marginal lake, including ice-dammed lakes along valley sides tangential to the main ice-flow, with that of all land margins, regardless of whether the margin is terminal or not. However, it is likely that this velocity difference is greater when considering only terminal ice-contact proglacial lakes (i.e. those where the main flow unit within an outlet glacier is flowing directly into the lake). While ice- and moraine-dammed lakes are susceptible to periodic draining and catastrophic outburst floods (Costa and Schuster, Reference Costa and Schuster1988; Carrivick and Tweed, Reference Carrivick and Tweed2013), bedrock-dammed lakes are inherently stable and maintain a greater level of hydraulic connectivity to the up-glacier system (Sugiyama and others, Reference Sugiyama2011; Carrivick and Tweed, Reference Carrivick and Tweed2013). This is illustrated by the ratio of ice surface velocities at IS and KS (2 km from their respective termini) differing by a factor ~1.5 in 2013 and ~3.4 in 2021. Furthermore, this ratio increases toward the terminus, with velocities at IS an order of magnitude greater than those at KS at the end of the study period. Additionally, bedrock overdeepenings are typically sited in regions where ice flow is laterally constrained by topography, and their association with outlet glacier confluences means that they are often within large ice catchments (Patton and others, Reference Patton, Swift, Clark, Livingstone and Cook2016). Consequently, bedrock-dammed proglacial lakes occupying glacially eroded valley bottom overdeepenings, are likely to be of greater importance in controlling ice sheet mass balance than ice-dammed marginal lakes, and future investigations should prioritise these.

In summary, the near terminus thinning and acceleration observed at IS highlight the potential importance of proglacial lakes, as the combined effects on ice dynamics reach inland and can lead to greater rates of mass loss. The behaviour replicates the expected positive feedback effects associated with sustained terminus thinning at calving glaciers (Benn and others, Reference Benn, Warren and Mottram2007) and we argue that as both ice-marginal melt-rate (The IMBIE Team, 2020) and lake number (Carrivick and Quincey, Reference Carrivick and Quincey2014; Shugar and others, Reference Shugar2020) increase, so too will the significance of ice-marginal lake processes for GrIS mass loss. The importance of this behaviour can currently be seen in Alaska (Trüssel and others, Reference Trüssel, Motyka, Truffer and Larsen2013; Larsen and others, Reference Larsen2015), Novaya Zemlya (Carr and others, Reference Carr, Bell, Killick and Holt2017), and the Patagonian ice fields where ice mass loss is strongly controlled by fast-flowing lake-terminating outlets (Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2014).

5. Conclusion

As the margin of the Greenland Ice Sheet recedes, ice-marginal lakes are expected to increase in both number and area in the coming decades, with an attendant increase in their influence on the wider ice sheet (Carrivick and others, Reference Carrivick2022). Our findings suggest that the distinct dynamic differences between the land- and lake-terminating outlets in this study are largely attributable to the presence of the lake. More specifically, we argue that the doubling of near-terminus ice velocity at Isortuarsuup Sermia is likely driven by ongoing negative surface mass balance and glacier thinning. This has reduced ice-overburden pressures near the terminus, where the lake maintains high basal water pressures year-round, and facilitated ice acceleration in all seasons. Furthermore, ongoing thinning and subsequent flotation off a sublacustrine moraine has instigated retreat, thereby promoting enhanced acceleration across the terminus region through the removal of buttressing. In contrast, reductions in ice thickness at the land-terminating KS have not led to flow acceleration due to the profound differences in terminus processes and subglacial hydrological setting. The acceleration and attendant extensional flow at Isortuarsuup Sermia has also led to enhanced rates of thinning near-terminus of between 0.33 and 0.65 m yr−1.

Our observations show the effect of recent mass balance change on the ice-dynamics of a lake-terminating glacier reaches ~15 km up-glacier, highlighting the ability of proglacial lakes to perturb inland ice. This supports earlier observations (Kirkbride, Reference Kirkbride1993; Warren and Kirkbride, Reference Warren and Kirkbride2003; Sakakibara and Sugiyama, Reference Sakakibara and Sugiyama2014; Tsutaki and others, Reference Tsutaki2019; Mallalieu and others, Reference Mallalieu, Carrivick, Quincey and Raby2021) and modelling work (Sutherland and others, Reference Sutherland2020) that stress the importance of proglacial lakes on glacier and ice sheet mass loss. We suggest future work should discriminate between dam type and lake setting (as per Rick and others, Reference Rick, McGrath, Armstrong and McCoy2022) when evaluating ice-marginal lake impacts on ice dynamics, as we contend that the relative importance of proglacial bedrock-dammed lakes on ice sheet mass loss is likely greater than ice-dammed marginal lakes. Additionally, there is a need to establish whether the recent pattern of behaviour seen at Isortuarsuup Sermia is typical for other Greenlandic lake-terminating outlets. Accurately quantifying the effect of ice-marginal lakes on these glaciers demands greater knowledge of ice-marginal lake characteristics, including bathymetry. This work is timely, as climate warming is seeing the ice margin retreat towards the many glacially eroded overdeepenings beneath the Greenland Ice Sheet. An increased incidence of lake-terminating glaciers would likely enhance the dynamic mass loss from Greenland due to accelerated glacier flow, in line with expected positive feedbacks associated with melt induced thinning of these glacier termini (Benn and others, Reference Benn, Warren and Mottram2007), and as witnessed already across numerous glaciated regions including Alaska (Trüssel and others, Reference Trüssel, Motyka, Truffer and Larsen2013), Iceland (Baurley and others, Reference Baurley, Robson and Hart2020), and Patagonia (Sugiyama and others, Reference Sugiyama2011, Reference Sugiyama, Minowa and Schaefer2019).

Supplementary material

The supplementary material for this article can be found at https://doi.org/10.1017/jog.2024.30

Data

The code necessary to reproduce the figures in this study are available at doi.org/10.5281/zenodo.7824988. All secondary data used in this paper are freely available and cited in the reference list.

Acknowledgements

We thank the authors of the several open data sets used in this study, specifically Gardner and others (Reference Gardner, Fahnestock and Scambos2022) for the velocity data; Porter and others (Reference Porter2018) for the surface elevation data. We also thank Mankoff and others (Reference Mankoff2020) for generating runoff estimates; Lea (Reference Lea2018) for creating the Google Earth Engine Digitization Tool; and the three anonymous reviewers whose comments and suggestions considerably improved the manuscript. Additionally, we acknowledge the open-source Python 3 packages used to process these data: Geopandas, Xarray, Pandas, Shapely, NumPy, Scipy, Dask, Seaborn and Matplotlib. Funding for this research was provided by NERC through a SENSE CDT studentship to EH (NE/T00939X/1).

Author contributions

E.H. and P.N. conceptualised the study and led the interpretation. E.H. conducted the data analysis with guidance from E.M.L., and E.H. prepared the manuscript with contributions from all co-authors.

Competing interest

None.

References

Baurley, NR, Robson, BA and Hart, JK (2020) Long-term impact of the proglacial lake Jökulsárlón on the flow velocity and stability of Breiðamerkurjökull glacier, Iceland. Earth Surface Processes and Landforms 45(11), 2647–2663. doi: 10.1002/esp.4920CrossRefGoogle Scholar
Benn, DI, Warren, CR and Mottram, RH (2007) Calving processes and the dynamics of calving glaciers. Earth-Science Reviews 82(3-4), 143–179. doi: 10.1016/j.earscirev.2007.02.002CrossRefGoogle Scholar
Bindschadler, R (1983) The importance of pressurized subglacial water in separation and sliding at the glacier bed. Journal of Glaciology 29(101), 3–19. doi: 10.3189/S0022143000005104CrossRefGoogle Scholar
Błaszczyk, M and 10 others (2019) Quality assessment and glaciological applications of digital elevation models derived from Space-Borne and Aerial images over two Tidewater Glaciers of Southern Spitsbergen. Remote Sensing 11(9), 1121. doi: 10.3390/rs11091121CrossRefGoogle Scholar
Box, JE and Decker, DT (2011) Greenland marine-terminating glacier area changes: 2000–2010. Annals of Glaciology 52(59), 91–98. doi: 10.3189/172756411799096312CrossRefGoogle Scholar
Boyce, ES, Motyka, RJ and Truffer, M (2007) Flotation and retreat of a lake-calving terminus, Mendenhall Glacier, southeast Alaska, USA. Annals of Glaciology 53(181), 211–224. doi: 10.3189/172756507782202928CrossRefGoogle Scholar
Carr, JR, Stokes, C and Vieli, A (2014) Recent retreat of major outlet glaciers on Novaya Zemlya, Russian Arctic, influenced by fjord geometry and sea-ice conditions. Journal of Glaciology 60(219), 155–170. doi: 10.3189/2014JoG13J122CrossRefGoogle Scholar
Carr, R, Bell, H, Killick, R and Holt, T (2017) Exceptional retreat of Novaya Zemlya's marine-terminating outlet glaciers between 2000 and 2013. Cryosphere 11, 2149–2174. doi: 10.5194/tc-11-2149-2017CrossRefGoogle Scholar
Carrivick, JL and Quincey, DJ (2014) Progressive increase in number and volume of ice-marginal lakes on the western margin of the Greenland Ice Sheet. Global and Planetary Change 116, 156–163. doi: 10.1016/j.gloplacha.2014.02.009CrossRefGoogle Scholar
Carrivick, JL and Tweed, FS (2013) Proglacial lakes: character, behaviour and geological importance. Quaternary Science Reviews 78, 34–52. doi: 10.1016/j.quascirev.2013.07.028CrossRefGoogle Scholar
Carrivick, JL, Tweed, FS, Sutherland, JL and Mallalieu, J (2020) Toward numerical modeling of interactions between ice-marginal proglacial lakes and glaciers. Frontiers in Earth Science 8, 577068. doi: 10.3389/feart.2020.577068CrossRefGoogle Scholar
Carrivick, JL and 8 others (2022) Ice-marginal proglacial lakes across greenland: present status and a possible future. Geophysical Research Letters 49(12), e2022GL099276. doi: 10.1029/2022GL099276CrossRefGoogle Scholar
Catania, GA and 7 others (2018) Geometric controls on tidewater glacier retreat in central western greenland. Journal of Geophysical Research. Earth Surface 123(8), 2024–2038. doi: 10.1029/2017JF004499CrossRefGoogle Scholar
Costa, JE and Schuster, RL (1988) The formation and failure of natural dams. Geological Society of America Bulletin 100(7), 1054–1068. doi: 10.1130/0016-7606(1988)100<1054:TFAFON>2.3.CO;22.3.CO;2>CrossRefGoogle Scholar
Davison, BJ, Sole, AJ, Livingstone, SJ, Cowton, TR and Nienow, PW (2019) The influence of hydrology on the dynamics of land-terminating sectors of the Greenland ice sheet. Frontiers in Earth Science 7, 10. doi: 10.3389/feart.2019.00010CrossRefGoogle Scholar
Dye, A, Bryant, R, Dodd, E, Falcini, F and Rippin, DM (2021) Warm arctic proglacial lakes in the aster surface temperature product. Remote Sensing 13(15), 2987. doi: 10.3390/rs13152987CrossRefGoogle Scholar
Dykes, RC, Brook, MS, Robertson, CM and Fuller, IC (2011) Twenty-first century calving retreat of tasman glacier, Southern Alps, New Zealand. Arctic, Antarctic, and Alpine Research 43(1), 1–10. doi: 10.1657/1938-4246-43.1.1CrossRefGoogle Scholar
Enderlin, EM and 5 others (2016) An empirical approach for estimating stress-coupling lengths for marine-terminating glaciers. Frontiers in Earth Science 4, 00104.CrossRefGoogle Scholar
ESA Copernicus (2022) Sentinel-2 Data Products.Google Scholar
Fettweis, X and 8 others (2017) Reconstructions of the 1900–2015 Greenland ice sheet surface mass balance using the regional climate MAR model. The Cryosphere 11(2), 1015–1033. doi: 10.5194/tc-11-1015-2017CrossRefGoogle Scholar
Fox-Kemper, B and 17 others (2021) Ocean, Cryosphere and Sea Level Change. In Masson-Delmotte V, and 18 others (eds), Climate Change 2021: The Physical Science Basis. Contribution of Working Group I to the Sixth Assessment Report of the Intergovernmental Panel on Climate Change, Cambridge University Press, Cambridge, United Kingdom and New York, NY, USA, pp. 1211–1362.Google Scholar
Gardner, AS and 7 others (2018) Increased West Antarctic and unchanged East Antarctic ice discharge over the last 7 years. The Cryosphere 12(2), 521–547. doi: 10.5194/tc-12-521-2018CrossRefGoogle Scholar
Gardner, AS, Fahnestock, MA and Scambos, TA (2022) ITS_LIVE Regional Glacier and Ice Sheet Surface Velocities.Google Scholar
Haresign, E and Warren, CR (2005) Melt rates at calving termini: a study at Glaciar León, Chilean Patagonia. Geological Society Special Publication 242, 99–109. doi: 10.1144/GSL.SP.2005.242.01.09CrossRefGoogle Scholar
How, P and 10 others (2021) Greenland-wide inventory of ice marginal lakes using a multi-method approach. Scientific Reports 11(4481), 1–13. doi: 10.1038/s41598-021-83509-1CrossRefGoogle ScholarPubMed
Howat, IM, Joughin, I, Tulaczyk, S and Gogineni, S (2005) Rapid retreat and acceleration of Helheim Glacier, East Greenland. Geophysical Research Letters 32(22), L22502. doi: 10.1029/2005GL024737CrossRefGoogle Scholar
Hurst, MD, Mudd, SM, Walcott, R, Attal, M and Yoo, K (2012) Using hilltop curvature to derive the spatial distribution of erosion rates. Journal of Geophysical Research 117(F2), 2011JF002057. doi: 10.1029/2011JF002057CrossRefGoogle Scholar
James, TD, Murray, T, Selmes, N, Scharrer, K and O'Leary, M (2014) Buoyant flexure and basal crevassing in dynamic mass loss at Helheim Glacier. Nature Geoscience 7(8), 593–596. doi: 10.1038/ngeo2204CrossRefGoogle Scholar
Joughin, I, Smith, BE and Howat, I (2018) Greenland ice mapping project: ice flow velocity variation at sub-monthly to decadal timescales. The Cryosphere 12(7), 2211–2227. doi: 10.5194/tc-12-2211-2018CrossRefGoogle Scholar
King, O, Dehecq, A, Quincey, D and Carrivick, J (2018) Contrasting geometric and dynamic evolution of lake and land-terminating glaciers in the central Himalaya. Global and Planetary Change 167(2017), 46–60. doi: 10.1016/j.gloplacha.2018.05.006CrossRefGoogle Scholar
King, O, Bhattacharya, A, Bhambri, R and Bolch, T (2019) Glacial lakes exacerbate Himalayan glacier mass loss. Scientific Reports 9(1), 1–9. doi: 10.1038/s41598-019-53733-xCrossRefGoogle ScholarPubMed
Kirkbride, MP (1993) The temporal significance of transitions from melting to calving termini at glaciers in the central Southern Alps of New Zealand. Holocene 3(3), 232–240. doi: 10.1177/095968369300300305CrossRefGoogle Scholar
Kirkbride, MP and Warren, CR (1997) Calving processes at a grounded ice cliff. Annals of Glaciology 24, 116–121. doi: 10.3189/s0260305500012039CrossRefGoogle Scholar
Larsen, CF and 5 others (2015) Surface melt dominates Alaska glacier mass balance. Geophysical Research Letters 42(14), 5902–5908. doi: 10.1002/2015GL064349CrossRefGoogle Scholar
Lea, JM (2018) The Google Earth Engine Digitisation Tool (GEEDiT) and the Margin change Quantification Tool (MaQiT) – simple tools for the rapid mapping and quantification of changing Earth surface margins. Earth Surface Dynamics 6(3), 551–561. doi: 10.5194/esurf-6-551-2018CrossRefGoogle Scholar
Lei, Y, Gardner, A and Agram, P (2021) Autonomous repeat image feature tracking (autoRIFT) and its application for tracking ice displacement. Remote Sensing 13(4), 749. doi: 10.3390/rs13040749CrossRefGoogle Scholar
Lei, Y, Gardner, AS and Agram, P (2022) Processing methodology for the ITS_LIVE Sentinel-1 ice velocity products. Earth System Science Data 14(11), 5111–5137. doi: 10.5194/essd-14-5111-2022CrossRefGoogle Scholar
Mallalieu, J, Carrivick, JL, Quincey, DJ and Smith, MW (2020) Calving seasonality associated with melt-undercutting and lake ice cover. Geophysical Research Letters 47(8), 1–11. doi: 10.1029/2019GL086561CrossRefGoogle Scholar
Mallalieu, J, Carrivick, JL, Quincey, DJ and Raby, CL (2021) Ice-marginal lakes associated with enhanced recession of the Greenland Ice Sheet. Global and Planetary Change 202(8), 103503. doi: 10.1016/j.gloplacha.2021.103503CrossRefGoogle Scholar
Mankoff, K (2020) Streams, Outlets, Basins, and Discharge [k = 1.0]. GEUS Dataverse, V5. doi: 10.22008/FK2/XKQVL7CrossRefGoogle Scholar
Mankoff, KD and 9 others (2020) Greenland liquid water discharge from 1958 through 2019. Earth System Science Data 12(4), 2811–2841. doi: 10.5194/essd-12-2811-2020CrossRefGoogle Scholar
Meier, MF and Post, A (1987) Fast tidewater glaciers. Journal of Geophysical Research: Solid Earth 92(B9), 9051–9058. doi: 10.1029/JB092iB09p09051CrossRefGoogle Scholar
Minowa, M, Sugiyama, S, Sakakibara, D and Skvarca, P (2017) Seasonal variations in ice-front position controlled by frontal ablation at glaciar perito moreno, the Southern Patagonia icefield. Frontiers in Earth Science 5, 00001. doi: 10.3389/feart.2017.00001CrossRefGoogle Scholar
Moon, T and Joughin, I (2008) Changes in ice front position on Greenland's outlet glaciers from 1992 to 2007. Journal of Geophysical Research: Earth Surface 113, F02022. doi: 10.1029/2007JF000927CrossRefGoogle Scholar
Morlighem, M and 31 others (2017) BedMachine v3: complete bed topography and ocean bathymetry mapping of Greenland from multibeam echo sounding combined with mass conservation. Geophysical Research Letters 44(21), 11,051–11,061. doi: 10.1002/2017GL074954CrossRefGoogle ScholarPubMed
Morlighem, M and 31 others (2022) IceBridge BedMachine Greenland, Version 5. Boulder, Colorado USA. NASA National Snow and Ice Data Center Distributed Active Archive Center. doi: 10.5067/GMEVBWFLWA7XCrossRefGoogle Scholar
Mouginot, J and 8 others (2019) Forty-six years of Greenland Ice Sheet mass balance from 1972 to 2018. Proceedings of the National Academy of Sciences of the United States of America 116(19), 9239–9244. doi: 10.1073/pnas.1904242116CrossRefGoogle ScholarPubMed
Nick, FM, Vieli, A, Howat, IM and Joughin, I (2009) Large-scale changes in Greenland outlet glacier dynamics triggered at the terminus. Nature Geoscience 2(2), 110–114. doi: 10.1038/ngeo394CrossRefGoogle Scholar
Noël, B and 6 others (2016) A daily, 1 km resolution data set of downscaled Greenland ice sheet surface mass balance (1958–2015). The Cryosphere 10(5), 2361–2377. doi: 10.5194/tc-10-2361-2016CrossRefGoogle Scholar
Noh, MJ and Howat, IM (2015) Automated stereo-photogrammetric DEM generation at high latitudes: surface extraction with TIN-based Search-space Minimization (SETSM) validation and demonstration over glaciated regions. GIScience & Remote Sensing 52(2), 198–217. doi: 10.1080/15481603.2015.1008621CrossRefGoogle Scholar
Nuth, C and Kääb, A (2011) Co-registration and bias corrections of satellite elevation data sets for quantifying glacier thickness change. The Cryosphere 5(1), 271–290. doi: 10.5194/tc-5-271-2011CrossRefGoogle Scholar
O'Neel, S, Pfeffer, WT, Krimmel, R and Meier, M (2005) Evolving force balance at Columbia Glacier, Alaska, during its rapid retreat. Journal of Geophysical Research 110(F3), F03012. doi: 10.1029/2005JF000292CrossRefGoogle Scholar
Patton, H, Swift, DA, Clark, CD, Livingstone, SJ and Cook, SJ (2016) Distribution and characteristics of overdeepenings beneath the Greenland and Antarctic ice sheets: Implications for overdeepening origin and evolution. Quaternary Science Reviews 148, 128–145. doi: 10.1016/j.quascirev.2016.07.012CrossRefGoogle Scholar
Paul, F and 19 others (2013) On the accuracy of glacier outlines derived from remote-sensing data. Annals of Glaciology 54(63), 171–182. doi: 10.3189/2013AoG63A296CrossRefGoogle Scholar
Pfeffer, WT (2007) A simple mechanism for irreversible tidewater glacier retreat. Journal of Geophysical Research: Earth Surface 112, F03S25. doi: 10.1029/2006JF000590CrossRefGoogle Scholar
Porter, C and 28 others (2018) ArcticDEM, Version 3.Google Scholar
Porter, C and 17 others (2022) ArcticDEM - Strips, Version 4.1.Google Scholar
Pronk, JB, Bolch, T, King, O, Wouters, B and Benn, DI (2021) Contrasting surface velocities between lake- and land-terminating glaciers in the Himalayan region. The Cryosphere 15(12), 5577–5599. doi: 10.5194/tc-15-5577-2021CrossRefGoogle Scholar
Richards, J, Moore, RD and Forrest, AL (2012) Late-summer thermal regime of a small proglacial lake. Hydrological Processes 26(18), 2687–2695. doi: 10.1002/hyp.8360CrossRefGoogle Scholar
Rick, B, McGrath, D, Armstrong, W and McCoy, SW (2022) Dam type and lake location characterize ice-marginal lake area change in Alaska and NW Canada between 1984 and 2019. The Cryosphere 16(1), 297–314. doi: 10.5194/tc-16-297-2022CrossRefGoogle Scholar
Röhl, K (2006) Thermo-erosional notch development at fresh-water-calving Tasman Glacier, New Zealand. Journal of Glaciology 52(177), 203–213. doi: 10.3189/172756506781828773CrossRefGoogle Scholar
S2 MSI ESL Team (2022) Data Quality Report. Sentinel-2 L1C MSI - September 2022. Technical Report OMPC.CS.DQR.01.08-2022, ESA, ESA, published in Sentinel Online.Google Scholar
Sakakibara, D and Sugiyama, S (2014) Ice-front variations and speed changes of calving glaciers in the Southern Patagonia Icefield from 1984 to 2011: calving glaciers in southern Patagonia. Journal of Geophysical Research: Earth Surface 119(11), 2541–2554. doi: 10.1002/2014JF003148CrossRefGoogle Scholar
Schomacker, A (2010) Expansion of ice-marginal lakes at the Vatnajökull ice cap, Iceland, from 1999 to 2009. Geomorphology 119(3-4), 232–236. doi: 10.1016/j.geomorph.2010.03.022CrossRefGoogle Scholar
Shreve, RL (1972) Movement of Water in Glaciers. Journal of Glaciology 11(62), 205–214. doi: 10.3189/S002214300002219XCrossRefGoogle Scholar
Shugar, DH and 9 others (2020) Rapid worldwide growth of glacial lakes since 1990. Nature Climate Change 10(10), 939–945. doi: 10.1038/s41558-020-0855-4CrossRefGoogle Scholar
Sugiyama, S and 7 others (2011) Ice speed of a calving glacier modulated by small fluctuations in basal water pressure. Nature Geoscience 4(9), 597–600. doi: 10.1038/ngeo1218CrossRefGoogle Scholar
Sugiyama, S and 7 others (2016) Thermal structure of proglacial lakes in Patagonia. Journal of Geophysical Research: Earth Surface 121(12), 2270–2286. doi: 10.1002/2016JF004084CrossRefGoogle Scholar
Sugiyama, S, Minowa, M and Schaefer, M (2019) Underwater ice terrace observed at the front of Glaciar Grey, a freshwater calving glacier in Patagonia. Geophysical Research Letters 46(5), 2602–2609. doi: 10.1029/2018GL081441CrossRefGoogle Scholar
Sutherland, JL and 5 others (2020) Proglacial Lakes Control Glacier Geometry and Behavior During Recession. Geophysical Research Letters 47, e2020GL088865. doi: 10.1029/2020GL088865CrossRefGoogle Scholar
Tedstone, AJ and 5 others (2015) Decadal slowdown of a land-terminating sector of the Greenland Ice Sheet despite warming. Nature 526(7575), 692–695. doi: 10.1038/nature15722CrossRefGoogle ScholarPubMed
The IMBIE Team (2020) Mass balance of the Greenland Ice Sheet from 1992 to 2018. Nature 579(7798), 233–239. doi: 10.1038/s41586-019-1855-2CrossRefGoogle Scholar
Trüssel, BL, Motyka, RJ, Truffer, M and Larsen, CF (2013) Rapid thinning of lake-calving Yakutat Glacier and the collapse of the Yakutat Icefield, southeast Alaska, USA. Journal of Glaciology 59(213), 149–161. doi: 10.3189/2013J0G12J081CrossRefGoogle Scholar
Tsutaki, S, Nishimura, D, Yoshizawa, T and Sugiyama, S (2011) Changes in glacier dynamics under the influence of proglacial lake formation in Rhonegletscher, Switzerland. Annals of Glaciology 52(58), 31–36. doi: 10.3189/172756411797252194CrossRefGoogle Scholar
Tsutaki, S, Sugiyama, S, Nishimura, D and Funk, M (2013) Acceleration and flotation of a glacier terminus during formation of a proglacial lake in Rhonegletscher, Switzerland. Journal of Glaciology 59(215), 559–570. doi: 10.3189/2013JoG12J107CrossRefGoogle Scholar
Tsutaki, S and 6 others (2019) Contrasting thinning patterns between lake- and land-terminating glaciers in the Bhutanese Himalaya. The Cryosphere 13(10), 2733–2750. doi: 10.5194/tc-13-2733-2019CrossRefGoogle Scholar
van der Veen, C and Whillans, I (1989) Force budget: I. theory and numerical methods. Journal of Glaciology 35(119), 53–60. doi: 10.3189/002214389793701581CrossRefGoogle Scholar
Warren, CR (1991) Terminal environment, topographic control and fluctuations of West Greenland glaciers. Boreas 20(1), 1–15. doi: 10.1111/j.1502-3885.1991.tb00453.xCrossRefGoogle Scholar
Warren, CR and Kirkbride, MP (2003) Calving speed and climatic sensitivity of New Zealand lake-calving glaciers. Annals of Glaciology 36, 173–178. doi: 10.3189/172756403781816446CrossRefGoogle Scholar
Watson, CS and 5 others (2020) Mass loss from calving in Himalayan Proglacial Lakes. Frontiers in Earth Science 7(1), 1–19. doi: 10.3389/feart.2019.00342CrossRefGoogle Scholar
Weertman, J (1974) Stability of the junction of an ice sheet and an ice shelf. Journal of Glaciology 13(67), 3–11. doi: 10.3189/S0022143000023327CrossRefGoogle Scholar
Weidick, A, Bennike, O, Citterio, M and Nørgaard-Pedersen, N (2012) Neoglacial and historical glacier changes around Kangersuneq fjord in southern West Greenland. The Geological Survey of Denmark and Greenland 27, 1–68.CrossRefGoogle Scholar
Xdem Contributors (2021) Xdem. Zenodo. doi: 10.5281/zenodo.4809697CrossRefGoogle Scholar
Figure 0

Figure 1. (a) The lake-terminating Isortuarsuup Sermia (IS) and land-terminating Kangaasarsuup Sermia (KS) with their respective centrelines (solid lines) and runoff catchments from Mankoff and others (2020) (dashed lines). Ice surface contours at 250 m intervals generated from BedMachine v5 (Morlighem and others, 2022); background image: Sentinel-2 acquisition from 19 September 2022 (ESA Copernicus, 2022); Location of study area, south-west Greenland, shown by pink box in inset (upper right). Dashed red box at terminus of IS denotes extent of (b); (b) terminus region of IS illustrating the trim-line (white arrows), terminal moraine (black arrow), grounded icebergs (orange arrow). Sentinel-2 acquisition from 18 September 2019 (ESA Copernicus, 2022).

Figure 1

Figure 2. Annual average ice velocity (2013–2021) profiles along centrelines shown in Figure 1 at (a) Isortuarsuup Sermia, and (b) Kangaasarsuup Sermia. (c) and (d) show percentage change relative to 2013. (e) and (f) show linear trends where regression slope coefficients are significant at p$\leq 0.05$; error bars denote 95% confidence interval.

Figure 2

Figure 3. Isortuarsuup Sermia: (a) change in average annual velocity between 2013–2021; (b) rate of surface elevation change (September 2012– June 2021) from ArcticDEM (negative denotes thinning); manually digitised ice margin shown in black in (a) and (b); (c) terminus positions 2014–2021. Red box in (a) denotes extent of (c). White arrows in (c) indicates the Little Ice Age trim-line and the black arrow points to the associated terminal moraine with icebergs grounded on its sublacustrine extension indicated by the orange arrow.

Figure 3

Figure 4. Kangaasarsuup Sermia: (a) change in average annual velocity between 2013–2021; (b) rate of surface elevation change (June 2011– September 2021) from ArcticDEM (negative denotes thinning); manually digitised ice margin shown in black in (a) and (b); (c) terminus positions 2014–2021. Red box in (a) denotes extent of (c).

Figure 4

Figure 5. Time series of ice surface velocity at (a) 1  km and (b) 10 km from the terminus at the lake-terminating Isortuarsuup Sermia (blue) and the land-terminating Kangaasarsuup Sermia (orange, dashed). These velocities are computed from image-pairs separated by ≤ 30 days. Lines show the rolling 28 day median.

Figure 5

Figure 6. Seasonal velocity trends (2016–2021) along glacier centrelines at Isortuarsuup Sermia (blue) and Kangaasarsuup Sermia (orange) for winter (DJF, circles), summer (JJA, triangles) and autumn (SON, crosses). Seasonal trends derived from velocity fields where date_dt ≤ 30 days. Trends are linear fits, and only significant (p ≤ 0.05) trends are shown; error bars denote 95% confidence interval. Points at the same distance from the terminus have been offset from one another to aid readability.

Figure 6

Figure 7. Rates of surface elevation change at (a) Isortuarsuup Sermia (September 2012–June 2021) and (b) Kangaasarsuup Sermia (June 2011– September 2021) from ArcticDEM. Rates determined from linear regression. Only significant trends (p$\leq 0.05$) are shown. Error bars denote 95% confidence interval. ArcticDEM was sampled every 100 m along 9(7) parallel lines spaced every 250 m across the glacier at IS(KS). Coloured lines show 500 m rolling width-averaged median. Purple line (right hand axis) denotes number of ArcticDEMs with valid elevation measurements at each point along centreline. Shading represents standard deviation of number of ArcticDEMs at each point along centreline to account for the parallel offsets. (c) Rate of surface elevation change from ArcticDEM, against surface elevation at Isortuarsuup Sermia (blue) and Kangaasarsuup Sermia (orange); coloured lines show rolling median over 50 m bins. In (a), (b) and (c) negative values indicate surface thinning. (d) Net change in surface slope (positive indicates surface steepening) between first and last DEM; line represents median change across parallel lines, shading denotes median absolute deviation.

Figure 7

Figure 8. Relative terminus position at Isortuarsuup Sermia (blue) and KS (orange) July– September 2014–2022. Blue circles denote relative position measured using box method; dashed blue lines representative of winter advance. Orange line illustrates average annual relative terminus position, shading denotes 95% confidence interval.

Figure 8

Figure 9. Cumulative (a & b) and average daily runoff rate (5-day rolling average) (c & d) at Isortuarsuup Sermia (a & c) and Kangaasarsuup Sermia (b & d). Colours denote regional climate model with MAR in turquoise, and RACMO in orange.

Supplementary material: File

Holt et al. supplementary material

Holt et al. supplementary material
Download Holt et al. supplementary material(File)
File 4.3 MB