Search

Mass loss from iceberg calving at marine-terminating glaciers is one of the largest and most poorly constrained contributors to sea-level rise. However, our understanding of the processes controlling ice fracturing and crevasse evolution is incomplete. Here, we use Gabor filter banks to automatically map crevasse density and orientation through time on a ~150 km2 terminus region of Narsap Sermia, an outlet glacier of the southwest Greenland ice sheet. We find that Narsap Sermia is dominated by transverse (flow-perpendicular) crevasses near the ice front and longitudinal (flow-aligned) crevasses across its central region. Measured crevasse orientation varies on sub-annual timescales by more than 45 $^\circ$ in response to seasonal velocity changes, and also on multi-annual timescales in response to broader dynamic changes and glacier retreat. Our results show a gradual up-glacier propagation of the zone of flow-transverse crevassing coincident with frontal retreat and acceleration occurring in 2020/21, in addition to sub-annual crevasse changes primarily in transition zones between longitudinal to transverse crevasse orientation. This provides new insight into the dynamics of crevassing at large marine-terminating glaciers and a potential approach for the rapid identification of glacier dynamic change from a single pair of satellite images. [ABSTRACT FROM AUTHOR]

Iceberg calving accounts for up to half of mass loss from the Greenland Ice Sheet (GrIS), with their size distributions providing insights into glacier calving dynamics and impacting fjord environments through their melting and subsequent freshwater release. Iceberg area and volume data for the GrIS are currently limited to a handful of fjord locations, while existing approaches to iceberg detection are often time-consuming and are not always suited for long time series analysis over large spatial scales. This study presents a highly automated workflow that detects icebergs and appends their associated metadata within Google Earth Engine using high spatial resolution timestamped ArcticDEM (Arctic Digital Elevation Model) strip data. This is applied to three glaciers that exhibit a range of different iceberg concentrations and size distributions: Sermeq Kujalleq (Jakobshavn Isbræ), Umiammakku Isbræ and Kangiata Nunaata Sermia. A total of 39 ArcticDEM scenes are analysed, detecting a total of 163 738 icebergs with execution times of 6 min to 2 h for each glacier depending on the number of DEMs available and total area analysed, comparing well with the mapping of manually digitised outlines. Results reveal two distinct iceberg distributions at Sermeq Kujalleq and Kangiata Nunaata Sermia where iceberg density is high, and one distribution at Umiammakku Isbræ where iceberg density is low. Small icebergs (< 1000 m 2) are found to account for over 80 % of each glacier's icebergs; however, they only contribute to 10 %–37 % of total iceberg volume suggesting that large icebergs are proportionally more important for glacier mass loss and as fjord freshwater reservoirs. The overall dataset is used to construct new area-to-volume conversions (with associated uncertainties) that can be applied elsewhere to two-dimensional iceberg outlines derived from optical or synthetic aperture radar imagery. When data are expressed in terms of total iceberg count and volume, insight is provided into iceberg distributions that have potential applicability to observations and modelling of iceberg calving behaviour and fjord freshwater fluxes. Due to the speed and automated nature of our approach, this workflow offers the potential to interrogate iceberg data on a pan-Arctic scale where ArcticDEM strip data coverage allows. [ABSTRACT FROM AUTHOR]

Our ability to improve prognostic modeling of the Greenland Ice Sheet relies on understanding the long-term relationships between climate and mass flux (via iceberg calving) from marine-terminating tidewater glaciers (TWGs). Observations of recent TWG behavior are widely available, but long-term records of TWG advance are currently lacking. We present glacial geomorphological, sedimentological, archaeological, and modeling data to reconstruct the -20 km advance of Kangiata Nunaata Sermia (KNS; the largest tidewater glacier in southwest Greenland) during the first half of the past millennium. The data show that KNS advanced -15 km during the 12th and 13th centuries CE at a rate of -115 m a-1, contemporaneous with regional climate cooling toward the Little Ice Age and comparable to rates of TWG retreat witnessed over the past -200 years. Presence of Norse farmsteads proximal to KNS demonstrates their resilience to climate change, manifest as a rapidly advancing TWG in a cooling climate. The results place limits on the magnitude of ice-margin advance and demonstrate TWG sensitivity to climate cooling as well as warming. These data combined with our grounding-line stability analysis provide a long-term record that validates approaches to numerical modeling aiming to link calving to climate. [ABSTRACT FROM AUTHOR]

Gaining knowledge of tidewater glacier (TWG) margin evolution, solid ice flux and their responses to climate over large spatio-temporal scales provides valuable context for the projection of future Greenland ice sheet (GrIS) change. Although studies of sector-wide responses of TWGs exist, studies at an ice-sheet-wide scale have only just become feasible. Here, we present a dataset of 224 annual TWG margins for 1984–2017 (n = 3801), showing that averaged over regional scales, normalised terminus change is linear. Regionally linear retreat trends were identified across most sectors of the GrIS starting in the mid-1990s, although in contrast to previous studies, the northeastern sector is shown to have experienced sustained retreat since the mid-1980s. Through cointegration analyses, individual glaciers are shown to have differing sensitivities to potential climate drivers, though on a sector-wide scale the northwest and southeast are shown to be especially sensitive to annual sea surface temperature and June–July–August air temperature, respectively. Although 92% of the analysed glaciers experience retreat across the GrIS, observed increases in absolute flux for the entire ice sheet can be explained by changes in just 11 of these TWGs. [ABSTRACT FROM AUTHOR]

Changes in margins derived from satellite imagery are quantitative indicators of the environmental processes and drivers acting on the Earth's surface, for example retreating ice margins or coastal changes with rising sea level. However, the large-scale rapid visualisation and analysis of the satellite record is often impractical due to factors such as computer processing power, software availability, internet connection speed and/or user expertise in remote sensing. Here are presented three new, freely accessible tools that together can be used to process, visualise and review data from the full Landsat 4'8 and Sentinel 1'2 satellite records in seconds, enabling efficient mapping (through manual digitisation) and automated quantification of margin changes. These tools are highly accessible for users from a range of remote-sensing expertise (from academics to high school students), with minimal computational, licensing and knowledge-based barriers to access. The Google Earth Engine Digitisation Tool (GEEDiT) allows users to define a point anywhere on the planet and filter data from each satellite for user-defined time frames, maximum acceptable cloud cover extent, and options of predefined or custom image band combinations via a simple graphical user interface (GUI). GEEDiT allows georeferenced vectors to be easily and rapidly mapped from each image, with image metadata and user notes automatically appended to each vector, which can then be exported for subsequent analysis. The GEEDiT Reviewer tool allows users to quality control their own/others' data and also filter existing datasets based on the spatial/temporal requirements for their particular research question. The Margin change Quantification Tool (MaQiT) is complementary to GEEDiT and GEEDiT Reviewer, allowing the rapid quantification of these margin changes by utilising two well-established methods that have previously been used to measure glacier margin change and two new methods via a similarly simple GUI. A case study of the lake-terminating glacier Breiðamerkurjökull, Iceland, is used to demonstrate the complementary functionality of GEEDiT, GEEDiT Reviewer and MaQiT, though it should be noted that MaQiT is also suitable for the (re-)analysis of existing datasets not generated by GEEDiT. MaQiT has been developed with the original aim of quantifying tidewater glacier terminus change, though the methods included within the tool have potential for wide applications in multiple areas of Earth surface science (e.g. coastal and vegetation extent change). It is hoped that these tools will allow a wide range of researchers and students across the geosciences to efficiently map, analyse and access volumes of data that would have previously proven prohibitive. [ABSTRACT FROM AUTHOR]

Glacial varves can give significant insights into recession and melting rates of decaying ice sheets. Moreover, varve chronologies can provide an independent means of comparison to other annually resolved climatic archives, which ultimately help to assess the timing and response of an ice sheet to changes across rapid climate transitions. Here we report a composite 1257-year-long varve chronology from southeastern Sweden spanning the regional late Allerød-late Younger Dryas pollen zone. The chronology was correlated to the Greenland Ice-Core Chronology 2005 using the time-synchronous Vedde Ash volcanic marker, which can be found in both successions. For the first time, this enables secure placement of the Lateglacial Swedish varve chronology in absolute time. Geochemical analysis from new varve successions indicate a marked change in sedimentation regime accompanied by an interruption of ice-rafted debris deposition synchronous with the onset of Greenland Stadial 1 ( GS-1; 12 846 years before AD 1950). With the support of a simple ice-flow/calving model, we suggest that slowdown of sediment transfer can be explained by ice-sheet margin stabilization/advance in response to a significant drop of the Baltic Ice Lake level. A reassessment of chronological evidence from central-western and southern Sweden further supports the hypothesis of synchronicity between the first (penultimate) catastrophic drainage of the Baltic Ice Lake and the start of GS-1 in Greenland ice-cores. Our results may therefore provide the first chronologically robust evidence linking continental meltwater forcing to rapid atmosphere-ocean circulation changes in the North Atlantic. [ABSTRACT FROM AUTHOR]

Despite extensive study and debate regarding the significance of turbate (also known as 'rotational') microstructures in glacially deformed sediments, characteristics regarding the dimensions of these features remain unresolved. This study presents the first explicitly quantitative measurement and analysis of turbate microstructure dimensions, and their relation to till texture through thin section analysis. Samples were taken from coarse-resolution horizontal and vertical transects of a macroscopically homogenous subglacial till, with subset areas of each thin section (30 mm2) analysed. The frequency and apparent a-axes and b-axes of both coreless and cored turbate structures (and their corestones) were measured, and simple univariate statistical methods used to establish the (in-)variability of these dimensions through the till profile. Summarizing findings, (i) the dimensions of both cored and coreless turbate populations display log-normal distributions when all measurements are analysed together, although not all individual sample populations possess these same distributions; (ii) turbate dimension populations are inconsistent within a sample block, precluding evaluation of turbate variability through a profile from single thin sections; (iii) analysis of turbate morphology and variability provisionally indicate that the three-dimensional structure of turbates are likely to be cylindrical or flared, while weak relationships are also observed amongst till texture, turbate dimensions and frequency. [ABSTRACT FROM AUTHOR]

Ice‐marginal lakes can affect glacier dynamics but are ignored in studies of the evolution of the Greenland ice sheet (GrIS) and of peripheral mountain glaciers and ice caps (PGICs). Here we show that lakes occupy 10% of the GrIS ice margin and occur on 5% of PGICs. Ice velocity at the GrIS margin is enhanced by ∼ 25% at lakes versus on land. Mean ice discharge into lakes is ∼4.9 Gt.yr, which is ∼1% of ice discharged through marine termini. We locate thousands of subglacial overdeepenings within which 7,404 km2 of future lakes could form, all of which will be ice‐marginal at some time. Future lakes in the west and east will be restricted to the margin of the GrIS and within alpine valleys, respectively. This status and possible future leads us to contend that lakes should be incorporated into projections of Greenland ice loss. Plain Language Summary: Lakes forming on the edges of glaciers can affect ice dynamics but are not included within numerical models of the Greenland ice sheet (GrIS) or of peripheral mountain glaciers and ice caps (PGICs). We show that lakes presently occupy 10% of the GrIS ice margin and occur on 5% of PGICs. We determine that ice is typically moving 25% faster in the vicinity of lakes than at land margins. We estimate that ice discharge into lakes is ∼4.9 Gt.yr, which is 1% of that through marine glaciers. We identify thousands of depressions in the landscape that could hold future lakes as glaciers retreat from them. Overall, understanding of lake evolution should be incorporated into projections of Greenland ice loss and runoff. Key Points: Lakes occupy 10% of the ice margin and vary between 2% and 26% by regionIce velocity is ∼25% greater into lakes compared to onto land and total ice discharge into lakes is ∼4.9 Gt.yr−1Thousands of overdeepenings have potential to host more and larger lakes in the future but not all will be ice‐marginal simultaneously [ABSTRACT FROM AUTHOR]

Greenlandic tidewater glaciers (TWGs) have been undergoing widespread retreat since the mid-1990s and contribute up to 50% of mass loss from the Greenland Ice Sheet (GrIS), though a comprehensive annual record of their retreat during the satellite age is currently lacking. Here we present a Greenland-wide dataset of annual terminus positions for 220 TWGs for the period 1984-2017 based on Landsat 4-8 and Sentinel 2 imagery (n = 3833). These were manually digitised using the cloud-computing based Google Earth Engine Digitisation Tool (GEEDiT), and their changes quantified using the curvilinear box method within the Margin Change Quantification Tool (MaQiT; Lea, 2018). Results were analysed alongside regional climate data (including air and sea surface temperatures), and with the supervised machine learning tree ensemble method Random Forests to determine the existence of threshold-type behaviour that may influence terminus stability. After normalising terminus behaviour of all glaciers, our analysis highlights distinct linear trends in the regional response of TWG termini. The south-east, south-west and north-west regions are all found to behave comparably (advance/stability until the mid-1990s followed by sustained retreat until 2017). However, in the north-east sustained retreat occurred since the mid-1980s, which then accelerated in 2008/2009.The multi-decadal timescale data set generated by this study enables the identification of a regional linear trend of TWG behaviour for the first time, and has allowed the application of Random Forests to determine the relative influence of climate forcings on termini positions. We anticipate that the results of this study hold a potential for creating new, empirically based predictive models of TWG change. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet (GrIS) is losing mass as the climate warms through both increased meltwater runoff and ice discharge at marine-terminating sectors. At the ice sheet surface, meltwater runoff forms a dynamic supraglacial hydrological system which includes stream and river networks and large supraglacial lakes (SGLs). Streams and rivers can route water into crevasses or into supraglacial lakes with crevasses underneath, both of which can then hydrofracture to the ice sheet base, providing a mechanism for the surface meltwater to access the bed. Understanding where, when, and how much meltwater is transferred to the bed is important because variability in meltwater supply to the bed can increase ice flow speeds, potentially impacting the hypsometry of the ice sheet in grounded sectors, and iceberg discharge to the ocean. Here we present a new, physically based, supraglacial hydrology model for the GrIS that is able to simulate (a) surface meltwater routing and SGL filling; (b) rapid meltwater drainage to the ice sheet bed via the hydrofracture of surface crevasses both in and outside of SGLs; (c) slow SGL drainage via overflow in supraglacial meltwater channels; and, by offline coupling with a second model, (d) the freezing and unfreezing of SGLs from autumn to spring. We call the model the Supraglacial Hydrology Evolution and Drainage (or SHED) model. We apply the model to three study regions in southwest Greenland between 2015 and 2019 (inclusive) and evaluate its performance with respect to observed supraglacial lake extents and proglacial discharge measurements. We show that the model reproduces 80 % of observed lake locations and provides good agreement with observations in terms of the temporal evolution of lake extent. Modelled moulin density values are in keeping with those previously published, and seasonal and inter-annual variability in proglacial discharge agrees well with that which is observed, though the observations lag the model by a few days since they include transit time through the subglacial system, while the model does not. Our simulations suggest that lake drainage behaviours may be more complex than traditional models suggest, with lakes in our model draining through a combination of both overflow and hydrofracture and with some lakes draining only partially and then refreezing. This suggests that, in order to simulate the evolution of Greenland's surface hydrological system with fidelity, a model that includes all of these processes needs to be used. In future work, we will couple our model to a subglacial model and an ice flow model and thus use our estimates of where, when, and how much meltwater gets to the bed to understand the consequences for ice flow. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet is a major contributor to current and projected sea level rise in the warming climate. However, uncertainties in Greenland's contribution to future sea level rise remain, partly due to challenges in constraining the role of ice dynamics. Transient ice accelerations, or ice speed-up events, lasting from 1 d to 1 week, have the potential to indirectly affect the mass budget of the ice sheet. They are triggered by an overload of the subglacial drainage system due to an increase in water supply. In this study, we identify melt-induced ice speed-up events at the Russell Glacier, southwest Greenland, in order to analyse synoptic patterns driving these events. The short-term speed-up events are identified from daily ice velocity time series collected from six GPS stations along the glacier for each summer (May–October) from 2009 to 2012. In total, 45 ice speed-up events are identified, of which we focus on the 36 melt-induced events, where melt is derived from two in situ observational datasets and one regional climate model forced by ERA5 reanalysis. We identify two additional potential water sources, namely lake drainages and extreme rainfall, which occur during 14 and 4 out of the 36 melt-induced events, respectively. The 36 melt-induced speed-up events occur during synoptic patterns that can be grouped into three main clusters: (1) patterns that resemble atmospheric rivers with a landfall in southwest Greenland, (2) patterns with anticyclonic blocking centred over southwest Greenland, and (3) patterns that show low-pressure systems centred either south or southeast of Greenland. Out of these clusters, the one resembling atmospheric river patterns is linked to the strongest speed-up events induced by 2 to 3 d continuously increasing surface melt driven by anomalously high sensible heat flux and incoming longwave radiation. In the other two clusters, the net shortwave radiation dominates the contribution to the melt energy. As the frequency and intensity of these weather patterns may change in the warming climate, so may the frequency and intensity of ice speed-up events, ultimately altering the mass loss of the ice sheet. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet (GrIS) is losing mass as the climate warms through both increased meltwater runoff and ice discharge at marine terminating sectors. At the ice sheet surface, meltwater runoff forms a dynamic supraglacial hydrological system which includes stream/river networks and large supraglacial lakes (SGLs). Streams/rivers can route water into crevasses, or into supraglacial lakes with crevasses underneath, both of which can then hydrofracture to the ice sheet base, providing a mechanism for the surface meltwater to access the bed. Understanding where, when and how much meltwater is transferred to the bed is important because variability in meltwater supply to the bed can increase ice flow speeds, potentially impacting the hypsometry of the ice sheet in grounded sectors, and iceberg discharge to the ocean. Here we present a new, physically-based, supraglacial hydrology model for the GrIS that is able to simulate a) surface meltwater routing and SGL filling, b) rapid meltwater drainage to the ice-sheet bed via the hydrofracture of surface crevasses both in, and outside of, SGLs, c) slow SGL drainage via overflow in supraglacial meltwater channels and, by offline coupling with a second model, d) the freezing and unfreezing of SGLs from autumn to spring. We call the model Supraglacial Hydrology Evolution and Drainage (or SHED). We apply the model to three study regions in South West Greenland between 2015 and 2019 inclusive and evaluate its performance with respect to observed supraglacial lake extents, and proglacial discharge measurements. We show that the model reproduces 80% of observed lake locations, and provides good agreement with observations in terms of the temporal evolution of lake extent. Modelled moulin density values are in keeping with those previously published and seasonal and inter-annual variability in proglacial discharge agrees well with that observed, though the observations lag the model by a few days since they include transit time through the subglacial system and the model does not. Our simulations suggest that lake drainage behaviours may be more complex than traditional models suggest, with lakes in our model draining through a combination of both overflow and hydrofracture, and some lakes draining only partially and then refreezing. This suggests that in order to simulate the evolution of Greenland's surface hydrological system with fidelity, then a model that includes all of these processes needs to be used. In future work we will couple our model to a subglacial model and an ice flow model, and thus use our estimates of where, when and how much meltwater gets to the bed to understand the consequences for ice flow. [ABSTRACT FROM AUTHOR]

Marine-terminating outlet glacier terminus traces, mapped from satellite and aerial imagery, have been used extensively in understanding how outlet glaciers adjust to climate change variability over a range of timescales. Numerous studies have digitized termini manually, but this process is labor intensive, and no consistent approach exists. A lack of coordination leads to duplication of efforts, particularly for Greenland, which is a major scientific research focus. At the same time, machine learning techniques are rapidly making progress in their ability to automate accurate extraction of glacier termini, with promising developments across a number of optical and synthetic aperture radar (SAR) satellite sensors. These techniques rely on high-quality, manually digitized terminus traces to be used as training data for robust automatic traces. Here we present a database of manually digitized terminus traces for machine learning and scientific applications. These data have been collected, cleaned, assigned with appropriate metadata including image scenes, and compiled so they can be easily accessed by scientists. The TermPicks data set includes 39 060 individual terminus traces for 278 glaciers with a mean of 136 ± 190 and median of 93 of traces per glacier. Across all glaciers, 32 567 dates have been digitized, of which 4467 have traces from more than one author, and there is a duplication rate of 17 %. We find a median error of ∼ 100 m among manually traced termini. Most traces are obtained after 1999, when Landsat 7 was launched. We also provide an overview of an updated version of the Google Earth Engine Digitization Tool (GEEDiT), which has been developed specifically for future manual picking of the Greenland Ice Sheet. [ABSTRACT FROM AUTHOR]

The increasing volume and spatio-temporal resolution of satellite-derived ice velocity data have created new exploratory opportunities for the quantitative analysis of glacier dynamics. One potential technique, proper orthogonal decomposition (POD), also known as empirical orthogonal functions, has proven to be a powerful and flexible technique for revealing coherent structures in a wide variety of environmental flows. In this study we investigate the applicability of POD to an openly available TanDEM-X/TerraSAR-X-derived ice velocity dataset from Sermeq Kujalleq (Jakobshavn Isbræ), Greenland. We find three dominant modes with annual periodicity that we argue are explained by glaciological processes. The primary dominant mode is interpreted as relating to the stress reconfiguration at the glacier terminus, known to be an important control on the glacier's dynamics. The second and third largest modes together relate to the development of the spatially heterogenous glacier hydrological system and are primarily driven by the pressurisation and efficiency of the subglacial hydrological system. During the melt season, variations in the velocity shown in these two subsidiary modes are explained by the drainage of nearby supraglacial melt ponds, as identified with a Google Earth Engine Moderate Resolution Imaging Spectroradiometer (MODIS) dynamic thresholding technique. By isolating statistical structures within velocity datasets and through their comparison to glaciological theory and complementary datasets, POD indicates which glaciological processes are responsible for the changing bulk velocity signal, as observed from space. With the proliferation of optical- and radar-derived velocity products (e.g. MEaSUREs, ESA CCI, PROMICE), we suggest POD, and potentially other modal decomposition techniques, will become increasingly useful in future studies of ice dynamics. [ABSTRACT FROM AUTHOR]

Surface meltwater is widespread around the Antarctic Ice Sheet margin and has the potential to influence ice shelf stability, ice flow and ice–albedo feedbacks. Our understanding of the seasonal and multi-year evolution of Antarctic surface meltwater is limited. Attempts to generate robust meltwater cover time series have largely been constrained by computational expense or limited ice surface visibility associated with mapping from optical satellite imagery. Here, we add a novel method for calculating visibility metrics to an existing meltwater detection method within Google Earth Engine. This enables us to quantify uncertainty induced by cloud cover and variable image data coverage, allowing time series of surface meltwater area to be automatically generated over large spatial and temporal scales. We demonstrate our method on the Amery Ice Shelf region of East Antarctica, analysing 4164 Landsat 7 and 8 optical images between 2005 and 2020. Results show high interannual variability in surface meltwater cover, with mapped cumulative lake area totals ranging from 384 to 3898 km 2 per melt season. By incorporating image visibility assessments, however, we estimate that cumulative total lake areas are on average 42 % higher than minimum mapped values. We show that modelled melt predictions from a regional climate model provide a good indication of lake cover in the Amery region and that annual lake coverage is typically highest in years with a negative austral summer SAM index. Our results demonstrate that our method could be scaled up to generate a multi-year time series record of surface water extent from optical imagery at a continent-wide scale. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet (GrIS) is known to have experienced widespread retreat over the last century. Information on outlet glacier dynamics, prior to this, are limited due to both a lack of observations and a paucity of mapped or mappable deglacial evidence which restricts our understanding of centennial to millennial timescale dynamics of the GrIS. Here we present glacial geomorphological mapping, for upper Godthåbsfjord, covering 5800 km2 at a scale of 1:92,000, using a combination of ASTER GDEM V2, a medium-resolution DEM (error <10 m horizontal and <6 m vertical accuracy), panchromatic orthophotographs and ground truthing. This work provides a detailed geomorphological assessment for the area, compiled as a single map, comprising of moraines, meltwater channels, streamlined bedrock, sediment lineations, ice-dammed lakes, trimlines, terraces, gullied sediment and marine limits. Whilst some of the landforms have been previously identified, the new information presented here improves our understanding of ice margin behaviour and can be used for future numerical modelling and landform dating programmes. Data also form the basis for palaeoglaciological reconstructions and contribute towards understanding of the centennial to millennial timescale record of this sector of the GrIS. [ABSTRACT FROM AUTHOR]

High Mountain Asia (HMA) hosts the largest glacier concentration outside of polar regions. It is also distinct glaciologically as it forms one of two major surge clusters globally, and many glaciers there contradict the globally observed glacier recession trend. Surging glaciers are critical to our understanding of HMA glacier dynamics, threshold behaviour and flow instability, and hence have been the subject of extensive research, yet many dynamical uncertainties remain. Using the cloud-based geospatial data platform, Google Earth Engine (GEE) and GEE-developed tool, GEEDiT, to identify and quantify trends in the distribution and phenomenological characteristics of surging glaciers synoptically across HMA, we identified 137 glaciers as surging between 1987–2019. Of these, 55 were newly identified, 15 glaciers underwent repeat surges, and 18 were identified with enhanced glaciological hazard potential, most notably from Glacier Lake Outburst Floods (GLOFs). Terminus position time series analysis from 1987–2019 facilitated the development of a six-part phenomenological classification of glacier behaviour, as well as quantification of surge variables including active phase duration, terminus advance distance and rate, and surge periodicity. This research demonstrates the application of remote sensing techniques and the GEE platform to develop our understanding of surging glacier distribution and terminus phenomenology across large areas, as well as their ability to highlight potential geohazard locations, which can subsequently be used to focus monitoring efforts. [ABSTRACT FROM AUTHOR]

Kangerdlugssuaq Glacier is one of Greenland's largest tidewater outlet glaciers, accounting for approximately 5 % of all ice discharge from the Greenland Ice Sheet. In 2018 the Kangerlussuaq ice front reached its most retreated position, since observations began in 1932. We determine the relationship between retreat and: (i) ice velocity; and (ii) surface elevation change, to assess the impact of the retreat on the glacier trunk. Between 2016 and 2018 the glacier retreated ~5 km and brought the Kangerlussuaq ice front into a major (~15 km long) overdeepening. Coincident with this retreat, the glacier thinned as a result of near-terminus acceleration in ice flow. The subglacial topography means that 2016-18 terminus recession is likely to trigger a series of feedbacks between retreat, thinning and glacier acceleration, leading to a rapid and high-magnitude increase in discharge and sea level rise contribution. Dynamic thinning may continue until the glacier reaches the upward sloping bed ~10 km inland of its current position. Given the complexity and scale of the processes involved, such changes will not be represented in prognostic models of the Greenland Ice Sheet to 2100 and beyond. [ABSTRACT FROM AUTHOR]