Search

The rapid formation of glacial lakes is one of the most conspicuous landscape changes caused by atmospheric warming in glacierised regions. However, relatively little is known about the history and current state of glacial lakes in the High Arctic. This study aims to address this issue by providing the first inventory of glacial lakes in Svalbard, focusing in particular on the post-Little Ice Age evolution of glacial lakes and their typology. To do so, we used aerial photographs and satellite imagery together with archival topographic data from 1936 to 2020. The inventory comprises the development of 566 glacial lakes (146 km2) that were still in direct contact with glaciers during the period 2008–2012. The results show a consistent increase in the total area of glacial lakes from the 1930s to 2020 and suggest an apparent link between climatic and geological factors, and the formation of specific lake dam types: moraine, ice, or bedrock. We also detected 134 glacial lake drainage events that have occurred since the 1930s. This study shows that Svalbard has one of the highest rates of glacial lake development in the world, which is an indicator of the overall dynamics of landscape change in the archipelago in response to climate change.

This study provides the first comprehensive reconstruction of the dynamics of Iceberg Glacier, located on western Axel Heiberg Island, and reveals detailed observations of a complete surge for the first time in the Canadian Arctic. Historical aerial photographs, declassified intelligence satellite photographs, optical satellite imagery and synthetic aperture radar data were used to quantify changes in terminus position, ice velocity and glacier thickness since the 1950s. A surge initiated at the terminus in 1981 and terminated in 2003, suggesting a 22-year active phase. High surface velocities, reaching ~2300 m a−1 in 1991, were accompanied by a maximum terminus advance of >7 km and a large transfer of mass down-glacier, causing significant median trunk-wide surface elevation changes attaining >3 ± 1 m a−1. We suggest that the retreat from a pinning point, flotation of the terminus, the removal of sea-ice from the ice front, and an increase in subglacial meltwater availability from relatively high air temperatures in 1981 likely contributed to surge initiation. The ensuing quiescent period has seen a continual decrease in surface flow rates to an average centreline velocity of 11.5 m a−1 in 2020–21, a gradual steepening of the glacier surface and a > 2.5 km terminus retreat.

Albedo is a key factor in modulating the absorption of solar radiation on ice surfaces. Satellite measurements have shown a general reduction in albedo across the Greenland ice sheet over the past few decades, particularly along the western margin of the ice sheet, a region known as the Dark Zone (albedo < 0.45). Here we chose a combination of Landsat 4–8 and Sentinel 2 imagery to enable us to derive the longest record of albedo variations in the Dark Zone, running from 1984 to 2020. We developed a simple, pragmatic and efficient sensor transformation to provide a long time series of consistent, harmonized satellite imagery. Narrow to broadband conversion algorithms were developed from regression models of harmonized satellite data and in situ albedo from the Program for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather stations. The albedo derived from the harmonized Landsat and Sentinel 2 data shows that the maximum extent of the Dark Zone expanded rapidly between 2005 and 2007, increasing to ~280% of the average annual maximum extent of 2900 km2 to ~8000 km2 since. The Dark Zone is continuing to darken slowly, with the average annual minimum albedo decreasing at a rate of $\sim \!-0.0006 \pm 0.0004 \, {\rm a}^{-1}$ (p = 0.16, 2001–2020).

Air temperature inversions, a situation in which atmospheric temperature increases with height, are key components of the Arctic planetary boundary layer. The present study investigates the spatial and temporal variations of temperature inversions over different surface types (rock, gravel, snow, ice) along the Mittivakkat valley (southeast Greenland). For this purpose, 113 vertical profiles with high spatio-temporal resolution of air temperature and relative humidity were collected with unoccupied aerial vehicles (UAVs) during a 13-day field campaign in summer 2019. Air temperature inversions were present in 83% of the profiles, of which 24% were surface-based inversions and 76% were elevated inversions. The proglacial area covered with bare rock and gravel induces surface heating and convection during the day and, through interaction with local circulation patterns, leads to the frequent formation of elevated inversions. In contrast, the glacier surface itself acts as a persistent cooling surface and leads to the formation of surface-based inversions. A low-level fog layer that forms under the inversion layer may be causing non-linear vertical ablation gradients on Mittivakkat Gletsjer. Furthermore, we demonstrate that atmospheric measurements using UAVs can better capture small-scale processes than other products like radiosonde or modeled reanalysis data.

Calculation of the calving loss of tidewater glaciers depends on accurate bedrock information. In regional to global-scale projections of future tidewater glacier evolution this dependence is problematic. Bedrock topographies are often unknown and can only be modelled from surface properties. Existing approaches, however, mostly underestimate the ice thickness towards the calving fronts of marine-terminating glaciers. This implies a compromised performance of global-scale projection models which often employ functions of water depth at the calving fronts of tidewater glaciers. Here, we present a sensitivity study that analyses the impact of five different bedrock datasets on projected mass losses from the tidewater glacier Hansbreen in southern Svalbard. Our modelling study calculates the glacier's response to artificial mass-balance forcing. We show that bedrock inaccuracies may lead to a substantially deviating retreat behaviour. The common underestimation of frontal ice thickness/water depth in the modelled bedrock datasets induces an underestimation of sea level-relevant mass losses over the first several decades of modelling. The duration of this period is reduced when assuming warmer climates. Our results thus underline the importance of accurate bedrock topography data for the reliability of glacier evolution projections and for the accuracy of the temporal trajectories of related sea level-relevant mass losses.

Using historical and recent aerial photography and structure from motion (SfM) multiview stereo (MVS) techniques, we reconstruct the 1959 and 2018 ice surface topography and determine the geodetic mass balance of Bowman Glacier, a small mountain glacier on northern Ellesmere Island. This is combined with optical satellite imagery to reconstruct the evolution in extent of the glacier over six decades, and ground-penetrating radar measurements of ice thickness to estimate the remaining ice volume. Between 1959 and 2020, Bowman Glacier lost 78% of its extent (reducing from 2.75 to 0.61 km2), while average annual area loss rates have nearly tripled in the past two decades. Over the 1959–2018 period, glacier-wide ice-thickness change averaged −22.7 ± 4.7 m, corresponding to a mean specific annual mass balance of −347.0 ± 71.4 mm w.e. a−1. Projecting rates of area and volume change into the future indicates that the glacier will likely entirely disappear between 2030 and 2060. This study demonstrates the potential of SfM-MVS processing to generate elevation products from 1950/60s historical aerial photographs, and to extend observations of ice elevation and glacier volume change for the Canadian Arctic, prior to the satellite record.

Mass loss at the Greenland Ice Sheet is influenced by atmospheric processes controlling its surface mass balance, and by submarine melt and calving where glaciers terminate in fjords. There, an ice mélange – a composite matrix of calved ice bergs and sea ice – may provide a buttressing force on a glacier terminus and control terminus dynamics. Kangerlussuaq Glacier is a major outlet of the Greenland Ice Sheet, for which recent major retreat events in 2004/2005 and 2016–2018 coincided with the absence of an ice mélange in Kangerlussuaq Fjord. To better understand the response of Kangerlussuaq Glacier to climatic and oceanic drivers, a 2D flowline model is employed. Results indicate that an ice mélange buttressing force exerts a major control on calving frequency and rapid retreat. When an ice mélange forms in Kangerlussuaq Fjord, it provides stabilising forces and conditions favourable for winter terminus re-advance. When it fails to form during consecutive years, model results indicate that Kangerlussuaq Glacier is primed to retreat into the large overdeepenings in Kangerlussuaq Fjord, and to terminus positions more than 30 km farther inland, implying that excessive mass loss from Kangerlussuaq Glacier by the year 2065 cannot be excluded.

The rapid formation of glacial lakes is one of the most conspicuous landscape changes caused by atmospheric warming in glacierised regions. However, relatively little is known about the history and current state of glacial lakes in the High Arctic. This study aims to address this issue by providing the first inventory of glacial lakes in Svalbard, focusing in particular on the post-Little Ice Age evolution of glacial lakes and their typology. To do so, we used aerial photographs and satellite imagery together with archival topographic data from 1936 to 2020. The inventory comprises the development of 566 glacial lakes (146 km2) that were still in direct contact with glaciers during the period 2008–2012. The results show a consistent increase in the total area of glacial lakes from the 1930s to 2020 and suggest an apparent link between climatic and geological factors, and the formation of specific lake dam types: moraine, ice, or bedrock. We also detected 134 glacial lake drainage events that have occurred since the 1930s. This study shows that Svalbard has one of the highest rates of glacial lake development in the world, which is an indicator of the overall dynamics of landscape change in the archipelago in response to climate change. [ABSTRACT FROM AUTHOR]

This study provides the first comprehensive reconstruction of the dynamics of Iceberg Glacier, located on western Axel Heiberg Island, and reveals detailed observations of a complete surge for the first time in the Canadian Arctic. Historical aerial photographs, declassified intelligence satellite photographs, optical satellite imagery and synthetic aperture radar data were used to quantify changes in terminus position, ice velocity and glacier thickness since the 1950s. A surge initiated at the terminus in 1981 and terminated in 2003, suggesting a 22-year active phase. High surface velocities, reaching ~2300 m a−1 in 1991, were accompanied by a maximum terminus advance of >7 km and a large transfer of mass down-glacier, causing significant median trunk-wide surface elevation changes attaining >3 ± 1 m a−1. We suggest that the retreat from a pinning point, flotation of the terminus, the removal of sea-ice from the ice front, and an increase in subglacial meltwater availability from relatively high air temperatures in 1981 likely contributed to surge initiation. The ensuing quiescent period has seen a continual decrease in surface flow rates to an average centreline velocity of 11.5 m a−1 in 2020–21, a gradual steepening of the glacier surface and a > 2.5 km terminus retreat. [ABSTRACT FROM AUTHOR]

Albedo is a key factor in modulating the absorption of solar radiation on ice surfaces. Satellite measurements have shown a general reduction in albedo across the Greenland ice sheet over the past few decades, particularly along the western margin of the ice sheet, a region known as the Dark Zone (albedo < 0.45). Here we chose a combination of Landsat 4–8 and Sentinel 2 imagery to enable us to derive the longest record of albedo variations in the Dark Zone, running from 1984 to 2020. We developed a simple, pragmatic and efficient sensor transformation to provide a long time series of consistent, harmonized satellite imagery. Narrow to broadband conversion algorithms were developed from regression models of harmonized satellite data and in situ albedo from the Program for Monitoring of the Greenland Ice Sheet (PROMICE) automatic weather stations. The albedo derived from the harmonized Landsat and Sentinel 2 data shows that the maximum extent of the Dark Zone expanded rapidly between 2005 and 2007, increasing to ~280% of the average annual maximum extent of 2900 km2 to ~8000 km2 since. The Dark Zone is continuing to darken slowly, with the average annual minimum albedo decreasing at a rate of $\sim \!-0.0006 \pm 0.0004 \, {\rm a}^{-1}$ (p = 0.16, 2001–2020). [ABSTRACT FROM AUTHOR]

Air temperature inversions, a situation in which atmospheric temperature increases with height, are key components of the Arctic planetary boundary layer. The present study investigates the spatial and temporal variations of temperature inversions over different surface types (rock, gravel, snow, ice) along the Mittivakkat valley (southeast Greenland). For this purpose, 113 vertical profiles with high spatio-temporal resolution of air temperature and relative humidity were collected with unoccupied aerial vehicles (UAVs) during a 13-day field campaign in summer 2019. Air temperature inversions were present in 83% of the profiles, of which 24% were surface-based inversions and 76% were elevated inversions. The proglacial area covered with bare rock and gravel induces surface heating and convection during the day and, through interaction with local circulation patterns, leads to the frequent formation of elevated inversions. In contrast, the glacier surface itself acts as a persistent cooling surface and leads to the formation of surface-based inversions. A low-level fog layer that forms under the inversion layer may be causing non-linear vertical ablation gradients on Mittivakkat Gletsjer. Furthermore, we demonstrate that atmospheric measurements using UAVs can better capture small-scale processes than other products like radiosonde or modeled reanalysis data. [ABSTRACT FROM AUTHOR]

We revise and evaluate frontal ablation fluxes obtained by the Open Global Glacier Model (OGGM) for Greenland's tidewater peripheral glaciers de-coupled from the ice sheet. By making use of new region-wide ice thickness and solid ice discharge data, we re-evaluate model performance and suggest future research directions to improve the ice thickness estimation of glacier models. OGGM is unable to predict individual tidewater glacier dynamics well if it has to rely only on surface mass balance estimates and the assumption of a closed budget to constrain the calving parameterization. Velocity observations are essential to constrain the model and estimate the dynamic mass loss of Greenland's tidewater peripheral glaciers. [ABSTRACT FROM AUTHOR]

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change. [ABSTRACT FROM AUTHOR]

Knowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration. [ABSTRACT FROM AUTHOR]

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes. [ABSTRACT FROM AUTHOR]

Distributed Acoustic Sensing (DAS) is increasingly recognised as a valuable tool for glaciological seismic applications, although analysing the large data volumes generated in acquisitions poses computational challenges. We show the potential of active-source DAS to image and characterise subglacial sediment beneath a fast-flowing Greenlandic outlet glacier, estimating the thickness of sediment layers to be 20–30 m. However, the lack of subglacial velocity constraint limits the accuracy of this estimate. Constraint could be provided by analysing cryoseismic events in a counterpart 3-day record of passive seismicity through, for example, seismic tomography, but locating them within the 9 TB data volume is computationally inefficient. We describe experiments with data compression using the frequency-wavenumber (f-k) transform ahead of training a convolutional neural network, that provides a ~300-fold improvement in efficiency. In combining active and passive-source and our machine learning framework, the potential of large DAS datasets could be unlocked for a range of future applications. [ABSTRACT FROM AUTHOR]

Using historical and recent aerial photography and structure from motion (SfM) multiview stereo (MVS) techniques, we reconstruct the 1959 and 2018 ice surface topography and determine the geodetic mass balance of Bowman Glacier, a small mountain glacier on northern Ellesmere Island. This is combined with optical satellite imagery to reconstruct the evolution in extent of the glacier over six decades, and ground-penetrating radar measurements of ice thickness to estimate the remaining ice volume. Between 1959 and 2020, Bowman Glacier lost 78% of its extent (reducing from 2.75 to 0.61 km2), while average annual area loss rates have nearly tripled in the past two decades. Over the 1959–2018 period, glacier-wide ice-thickness change averaged −22.7 ± 4.7 m, corresponding to a mean specific annual mass balance of −347.0 ± 71.4 mm w.e. a−1. Projecting rates of area and volume change into the future indicates that the glacier will likely entirely disappear between 2030 and 2060. This study demonstrates the potential of SfM-MVS processing to generate elevation products from 1950/60s historical aerial photographs, and to extend observations of ice elevation and glacier volume change for the Canadian Arctic, prior to the satellite record. [ABSTRACT FROM AUTHOR]

Calculation of the calving loss of tidewater glaciers depends on accurate bedrock information. In regional to global-scale projections of future tidewater glacier evolution this dependence is problematic. Bedrock topographies are often unknown and can only be modelled from surface properties. Existing approaches, however, mostly underestimate the ice thickness towards the calving fronts of marine-terminating glaciers. This implies a compromised performance of global-scale projection models which often employ functions of water depth at the calving fronts of tidewater glaciers. Here, we present a sensitivity study that analyses the impact of five different bedrock datasets on projected mass losses from the tidewater glacier Hansbreen in southern Svalbard. Our modelling study calculates the glacier's response to artificial mass-balance forcing. We show that bedrock inaccuracies may lead to a substantially deviating retreat behaviour. The common underestimation of frontal ice thickness/water depth in the modelled bedrock datasets induces an underestimation of sea level-relevant mass losses over the first several decades of modelling. The duration of this period is reduced when assuming warmer climates. Our results thus underline the importance of accurate bedrock topography data for the reliability of glacier evolution projections and for the accuracy of the temporal trajectories of related sea level-relevant mass losses. [ABSTRACT FROM AUTHOR]

Mass loss at the Greenland Ice Sheet is influenced by atmospheric processes controlling its surface mass balance, and by submarine melt and calving where glaciers terminate in fjords. There, an ice mélange – a composite matrix of calved ice bergs and sea ice – may provide a buttressing force on a glacier terminus and control terminus dynamics. Kangerlussuaq Glacier is a major outlet of the Greenland Ice Sheet, for which recent major retreat events in 2004/2005 and 2016–2018 coincided with the absence of an ice mélange in Kangerlussuaq Fjord. To better understand the response of Kangerlussuaq Glacier to climatic and oceanic drivers, a 2D flowline model is employed. Results indicate that an ice mélange buttressing force exerts a major control on calving frequency and rapid retreat. When an ice mélange forms in Kangerlussuaq Fjord, it provides stabilising forces and conditions favourable for winter terminus re-advance. When it fails to form during consecutive years, model results indicate that Kangerlussuaq Glacier is primed to retreat into the large overdeepenings in Kangerlussuaq Fjord, and to terminus positions more than 30 km farther inland, implying that excessive mass loss from Kangerlussuaq Glacier by the year 2065 cannot be excluded. [ABSTRACT FROM AUTHOR]

Arctic ice shelves have declined over the past several decades, one of many indications of a rapidly changing cryosphere. Here we use a collection of off-nadir Landsat 8 images, a 1978 digital orthophotograph and photogrammetrically derived DEM, satellite altimetry and other data to examine the causes of an Arctic ice-shelf retreat in northernmost Greenland, the Hunt Fjord Ice Shelf (HFIS). HFIS has several distinct provenance regions comprised of glacier-derived ice and corrugated multi-decadal fast ice, with varying ice thicknesses (5–64 m). Available imagery shows little change in HFIS between 1978 and 2012, after which several midsummer calving events occurred (2012, 2016 and 2019) that reduced the HFIS by 42.5 km2 (~56%). Shelf area losses began as the number of surface melt days on the adjacent ice sheet more than doubled relative to the 1980s. Recent calving events also occurred during open-water periods at the ice-shelf front. Prior to mid-2012, there were no calving events during similar open-water periods. HFIS tributary glaciers have thinned by 3–20 m near their grounding zones, and may have accelerated since the 1980s, likely due to increased basal melting from contact with warm Atlantic Water.

We assess the evolution of glaciological structures during the 2003–05 surge in the Paulabreen glacier system, Svalbard. Glaciological structures on the glacier surface were mapped using aerial photographs captured in the early stages of the surge (2003) and 5 years after surge termination (2011). Three-dimensional measurements of glaciological structures were collected at the tidewater front in 2013. These datasets document the physical changes during (1) the late quiescent phase; (2) the early phase of the surge as the surge front propagated down Skobreen and advanced into Paulabreen and (3) the final stages of the surge following the surge front reaching the glacier terminus. Crevasse patterns and clusters of arcuate shear planes record zones of compressive and extensional flow associated with the downglacier progression of the surge front. The transfer of surging ice from Skobreen into Paulabreen caused lateral displacement of the medial moraines to the northeast. At the ice front, this movement tilted glaciological structures in the same direction. Structures at the southwest margin record strike–slip faulting and the elevation of debris into the ice in a zone of compression and transpression. We summarise these observations in a schematic reconstruction of structural evolution during the surge.

Iceberg calving is one of the major mechanisms of ice loss from tidewater glaciers and ice sheets, but obtaining accurate estimates of ice discharge that are both continuous and accurate is a challenging task. Recent results have demonstrated the effective application of passive cryoacoustics – the use of naturally generated sounds to study the cryosphere – to quantify subaerial calving fluxes. However, little is known about the acoustic signatures of submarine calving. This study investigates the underwater noise from 656 subaerial and 162 submarine calving events observed at Hansbreen, Svalbard in the summers of 2016 and 2017. Statistical analysis of the acoustic signal shows that the normalized power of the calving noise is log-normally distributed regardless of the calving mode. However, submarine events can be distinguished from subaerial events by using the shape parameter of the log-normal distribution paired with the calving signal duration. The newly developed classification model may potentially be used for two purposes: (1) to study potential causal relationships between these two calving modes and (2) to separate calving fluxes into subaerial and submarine components. The latter will also require knowledge of the relationship between ice mass and sound spectral level for submarine calving events.

The once-contiguous Ellesmere Ice Shelf, first reported in writing by European explorers in 1876, and now almost completely disintegrated, has rolling, wave-like surface topography, the origin of which we investigate using a viscous buckling instability analysis. We show that rolls can develop during a winter season (~ 100 d) if sea-ice pressure (depth-integrated horizontal stress applied to the seaward front of the Ellesmere Ice Shelf) is sufficiently large (1 MPa m) and ice thickness sufficiently low (1–10 m). Roll wavelength initially depends only on sea-ice pressure, but evolves over time depending on amplitude growth rate. This implies that a thinner ice shelf, with its faster amplitude growth rate, will have a shorter wavelength compared to a thicker ice shelf when sea-ice pressure is equal. A drawback of the viscous buckling mechanism is that roll amplitude decays once sea-ice pressure is removed. However, non-Newtonian ice rheology, where effective viscosity, and thus roll change rate, depends on total applied stress may constrain roll decay rate to be much slower than growth rate and allow roll persistence from year to year. Whether the viscous-buckling mechanism we explore here ultimately can be confirmed as the origin of the Ellesmere Ice Shelf rolls remains for future research.

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change.

Knowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration.

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes.

Distributed Acoustic Sensing (DAS) is increasingly recognised as a valuable tool for glaciological seismic applications, although analysing the large data volumes generated in acquisitions poses computational challenges. We show the potential of active-source DAS to image and characterise subglacial sediment beneath a fast-flowing Greenlandic outlet glacier, estimating the thickness of sediment layers to be 20–30 m. However, the lack of subglacial velocity constraint limits the accuracy of this estimate. Constraint could be provided by analysing cryoseismic events in a counterpart 3-day record of passive seismicity through, for example, seismic tomography, but locating them within the 9 TB data volume is computationally inefficient. We describe experiments with data compression using the frequency-wavenumber (f-k) transform ahead of training a convolutional neural network, that provides a ~300-fold improvement in efficiency. In combining active and passive-source and our machine learning framework, the potential of large DAS datasets could be unlocked for a range of future applications.

We revise and evaluate frontal ablation fluxes obtained by the Open Global Glacier Model (OGGM) for Greenland's tidewater peripheral glaciers de-coupled from the ice sheet. By making use of new region-wide ice thickness and solid ice discharge data, we re-evaluate model performance and suggest future research directions to improve the ice thickness estimation of glacier models. OGGM is unable to predict individual tidewater glacier dynamics well if it has to rely only on surface mass balance estimates and the assumption of a closed budget to constrain the calving parameterization. Velocity observations are essential to constrain the model and estimate the dynamic mass loss of Greenland's tidewater peripheral glaciers.

Marine-terminating glaciers lose mass through melting and iceberg calving, and we find that meltwater drainage systems influence calving timing at Helheim Glacier, a tidewater glacier in East Greenland. Meltwater feeds a buoyant subglacial discharge plume at the terminus of Helheim Glacier, which rises along the glacial front and surfaces through the mélange. Here, we use high-resolution satellite and time-lapse imagery to observe the surface expression of this meltwater plume and how plume timing and location compare with that of calving and supraglacial meltwater pooling from 2011 to 2019. The plume consistently appeared at the central terminus even as the glacier advanced and retreated, fed by a well-established channelized drainage system with connections to supraglacial water. All full-thickness calving episodes, both tabular and non-tabular, were separated from the surfacing plume by either time or by space. We hypothesize that variability in subglacial hydrology and basal coupling drive this inverse relationship between subglacial discharge plumes and full-thickness calving. Surfacing plumes likely indicate a low-pressure subglacial drainage system and grounded terminus, while full-thickness calving occurrence reflects a terminus at or close to flotation. Our records of plume appearance and full-thickness calving therefore represent proxies for the grounding state of Helheim Glacier through time.

Basal conditions play an essential role in the dynamics of outlet glaciers, but direct observations at the bed of glaciers are challenging to obtain. Instead, inverse methods can be used to infer basal parameters from surface observations. Here, we use a simple ice-flow model as a forward model in an inversion scheme to retrieve the spatio-temporally variable basal stress parameter for Hagen Bræ, North Greenland, from 1990 to 2020. Hagen Bræ is a surge-type glacier with up to an order of magnitude variability of winter velocities near the grounding line. We find that downstream changes in the basal stress parameter can explain most of the variation of flow velocity, and we further identify a region of high resistance ~20–40 km from the grounding line. We hypothesise that this region of high resistance plays an important role in controlling glacier discharge.

We present the first fully coupled 3D full-Stokes model of a tidewater glacier, incorporating ice flow, subglacial hydrology, plume-induced frontal melting and calving. We apply the model to Store Glacier (Sermeq Kujalleq) in west Greenland to simulate a year of high melt (2012) and one of low melt (2017). In terms of modelled hydrology, we find perennial channels extending 5 km inland from the terminus and up to 41 and 29 km inland in summer 2012 and 2017, respectively. We also report a hydrodynamic feedback that suppresses channel growth under thicker ice inland and allows water to be stored in the distributed system. At the terminus, we find hydrodynamic feedbacks exert a major control on calving through their impact on velocity. We show that 2012 marked a year in which Store Glacier developed a fully channelised drainage system, unlike 2017, where it remained only partially developed. This contrast in modelled behaviour indicates that tidewater glaciers can experience a strong hydrological, as well as oceanic, control, which is consistent with observations showing glaciers switching between types of behaviour. The fully coupled nature of the model allows us to demonstrate the likely lack of any hydrological or ice-dynamic memory at Store Glacier.

We conduct extended versions of the ISMIP6 future climate experiments for the Greenland ice sheet until the year 3000 with the model SICOPOLIS. Beyond 2100, the climate forcing is kept fixed at late-21st-century conditions. For the unabated warming pathway RCP8.5/SSP5-8.5, the ice sheet suffers a severe mass loss, which amounts to ~ 1.8 m SLE (sea-level equivalent) for the 12-experiment mean, and ~ 3.5 m SLE (~ 50% of the entire mass) for the most sensitive experiment. For the reduced emissions pathway RCP2.6/SSP1-2.6, the mass loss is limited to a two-experiment mean of ~ 0.28 m SLE. Climate-change mitigation during the next decades will therefore be an efficient means for limiting the contribution of the Greenland ice sheet to sea-level rise in the long term.

We assess the evolution of glaciological structures during the 2003–05 surge in the Paulabreen glacier system, Svalbard. Glaciological structures on the glacier surface were mapped using aerial photographs captured in the early stages of the surge (2003) and 5 years after surge termination (2011). Three-dimensional measurements of glaciological structures were collected at the tidewater front in 2013. These datasets document the physical changes during (1) the late quiescent phase; (2) the early phase of the surge as the surge front propagated down Skobreen and advanced into Paulabreen and (3) the final stages of the surge following the surge front reaching the glacier terminus. Crevasse patterns and clusters of arcuate shear planes record zones of compressive and extensional flow associated with the downglacier progression of the surge front. The transfer of surging ice from Skobreen into Paulabreen caused lateral displacement of the medial moraines to the northeast. At the ice front, this movement tilted glaciological structures in the same direction. Structures at the southwest margin record strike–slip faulting and the elevation of debris into the ice in a zone of compression and transpression. We summarise these observations in a schematic reconstruction of structural evolution during the surge. [ABSTRACT FROM AUTHOR]

Arctic ice shelves have declined over the past several decades, one of many indications of a rapidly changing cryosphere. Here we use a collection of off-nadir Landsat 8 images, a 1978 digital orthophotograph and photogrammetrically derived DEM, satellite altimetry and other data to examine the causes of an Arctic ice-shelf retreat in northernmost Greenland, the Hunt Fjord Ice Shelf (HFIS). HFIS has several distinct provenance regions comprised of glacier-derived ice and corrugated multi-decadal fast ice, with varying ice thicknesses (5–64 m). Available imagery shows little change in HFIS between 1978 and 2012, after which several midsummer calving events occurred (2012, 2016 and 2019) that reduced the HFIS by 42.5 km2 (~56%). Shelf area losses began as the number of surface melt days on the adjacent ice sheet more than doubled relative to the 1980s. Recent calving events also occurred during open-water periods at the ice-shelf front. Prior to mid-2012, there were no calving events during similar open-water periods. HFIS tributary glaciers have thinned by 3–20 m near their grounding zones, and may have accelerated since the 1980s, likely due to increased basal melting from contact with warm Atlantic Water. [ABSTRACT FROM AUTHOR]

Iceberg calving is one of the major mechanisms of ice loss from tidewater glaciers and ice sheets, but obtaining accurate estimates of ice discharge that are both continuous and accurate is a challenging task. Recent results have demonstrated the effective application of passive cryoacoustics – the use of naturally generated sounds to study the cryosphere – to quantify subaerial calving fluxes. However, little is known about the acoustic signatures of submarine calving. This study investigates the underwater noise from 656 subaerial and 162 submarine calving events observed at Hansbreen, Svalbard in the summers of 2016 and 2017. Statistical analysis of the acoustic signal shows that the normalized power of the calving noise is log-normally distributed regardless of the calving mode. However, submarine events can be distinguished from subaerial events by using the shape parameter of the log-normal distribution paired with the calving signal duration. The newly developed classification model may potentially be used for two purposes: (1) to study potential causal relationships between these two calving modes and (2) to separate calving fluxes into subaerial and submarine components. The latter will also require knowledge of the relationship between ice mass and sound spectral level for submarine calving events. [ABSTRACT FROM AUTHOR]

Microstructures, including crystallographic fabric, within the margin of streaming ice can exert strong control on flow dynamics. To characterize a natural setting, we retrieved three cores, two of which reached bed, from the flank of Jarvis Glacier, eastern Alaska Range, Alaska. The core sites lie ~1 km downstream of the source, with abundant water present in the extracted cores and at the base of the glacier. All cores exhibit dipping layers, a combination of debris bands and bubble-free domains. Grain sizes coarsen on average approaching the lateral margin. Crystallographic orientations are more clustered and with c-axes closer to horizontal nearer the lateral margin. The measured fabric is sufficiently weak to induce little mechanical anisotropy, but the data suggest that despite the challenging conditions of warm ice, abundant water and a short flow distance, many aspects of the microstructure, including measurable crystallographic fabric, evolved in systematic ways.

The once-contiguous Ellesmere Ice Shelf, first reported in writing by European explorers in 1876, and now almost completely disintegrated, has rolling, wave-like surface topography, the origin of which we investigate using a viscous buckling instability analysis. We show that rolls can develop during a winter season (~ 100 d) if sea-ice pressure (depth-integrated horizontal stress applied to the seaward front of the Ellesmere Ice Shelf) is sufficiently large (1 MPa m) and ice thickness sufficiently low (1–10 m). Roll wavelength initially depends only on sea-ice pressure, but evolves over time depending on amplitude growth rate. This implies that a thinner ice shelf, with its faster amplitude growth rate, will have a shorter wavelength compared to a thicker ice shelf when sea-ice pressure is equal. A drawback of the viscous buckling mechanism is that roll amplitude decays once sea-ice pressure is removed. However, non-Newtonian ice rheology, where effective viscosity, and thus roll change rate, depends on total applied stress may constrain roll decay rate to be much slower than growth rate and allow roll persistence from year to year. Whether the viscous-buckling mechanism we explore here ultimately can be confirmed as the origin of the Ellesmere Ice Shelf rolls remains for future research. [ABSTRACT FROM AUTHOR]

Knowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration. [ABSTRACT FROM AUTHOR]

Bed topography is a critical parameter for determining the modern-day and future dynamics of ice sheets and their outlet glaciers. This is because the topography controls the state of stress for glaciers. At glacier termini, topography can influence the timing of terminus retreat by controlling access to warm ocean waters and/or by influencing the ability of a glacier terminus to retreat over bed bumps (moraines). Inland from the terminus, the topography can also influence where glacier retreat and thinning can stabilize. In part, this is because of knickpoints in bed topography created through glacial erosion that may influence the extent to which thinning can diffuse inland for an individual glacier and thus, the timing and magnitude of long-term mass loss. Here we provide a review of the current literature on these topics. While much of the reviewed literature assumes that topography is stable on relevant timescales to humans, new research suggests that topography may change much faster than previously thought and this further complicates our ability to project future outlet glacier change. [ABSTRACT FROM AUTHOR]

Ice cliffs are features along ice sheet margins, along tropical mountain glaciers, at termini of mountain glaciers and on debris-covered glacier tongues, that have received scattered attention in literature. They cover small relative areas of glacier or margin surface respectively, but have been involved in two apparent anomalies. On the one hand, they have been identified as potential hotspots of extreme melt rates on debris-covered tongues contributing to their relatively rapid ablation, compared to the surrounding glacier surface. On the other hand, they appear where the ice margin is stable (or temporarily advancing) even under conditions of negative mass balance. In this manuscript, we recapitulate why ice cliffs remain interesting features to investigate and what we know about them so far. We conclude by suggesting to further investigate their genesis and variable morphology and their potential as windows into past climates and processes. [ABSTRACT FROM AUTHOR]

Distributed Acoustic Sensing (DAS) is increasingly recognised as a valuable tool for glaciological seismic applications, although analysing the large data volumes generated in acquisitions poses computational challenges. We show the potential of active-source DAS to image and characterise subglacial sediment beneath a fast-flowing Greenlandic outlet glacier, estimating the thickness of sediment layers to be 20–30 m. However, the lack of subglacial velocity constraint limits the accuracy of this estimate. Constraint could be provided by analysing cryoseismic events in a counterpart 3-day record of passive seismicity through, for example, seismic tomography, but locating them within the 9 TB data volume is computationally inefficient. We describe experiments with data compression using the frequency-wavenumber (f-k) transform ahead of training a convolutional neural network, that provides a ~300-fold improvement in efficiency. In combining active and passive-source and our machine learning framework, the potential of large DAS datasets could be unlocked for a range of future applications. [ABSTRACT FROM AUTHOR]

We revise and evaluate frontal ablation fluxes obtained by the Open Global Glacier Model (OGGM) for Greenland's tidewater peripheral glaciers de-coupled from the ice sheet. By making use of new region-wide ice thickness and solid ice discharge data, we re-evaluate model performance and suggest future research directions to improve the ice thickness estimation of glacier models. OGGM is unable to predict individual tidewater glacier dynamics well if it has to rely only on surface mass balance estimates and the assumption of a closed budget to constrain the calving parameterization. Velocity observations are essential to constrain the model and estimate the dynamic mass loss of Greenland's tidewater peripheral glaciers. [ABSTRACT FROM AUTHOR]

Greenland's future contribution to sea-level rise is strongly dependent on the extent to which dynamic perturbations, originating at the margin, can drive increased ice flow within the ice-sheet interior. However, reported observations of ice dynamical change at distances >~50 km from the margin have a very low spatial and temporal resolution. Consequently, the likely response of the ice-sheet's interior to future oceanic and atmospheric warming is poorly constrained. Through combining GPS and satellite-image-derived ice velocity measurements, we measure multi-decadal (1993–1997 to 2014–2018) velocity change at 45 inland sites, encompassing all regions of the ice sheet. We observe an almost ubiquitous acceleration inland of tidewater glaciers in west Greenland, consistent with acceleration and retreat at glacier termini, suggesting that terminus perturbations have propagated considerable distances (>100 km) inland. In contrast, outside of Kangerlussuaq, we observe no acceleration inland of tidewater glaciers in east Greenland despite terminus retreat and near-terminus acceleration, and suggest propagation may be limited by the influence of basal topography and ice geometry. This pattern of inland dynamical change indicates that Greenland's future contribution to sea-level will be spatially complex and will depend on the capacity for dynamic changes at individual outlet glacier termini to propagate inland.

Distant glacial areas are interconnected by a complex system of fractures and water channels which run in the glacier interior and characterize the englacial realm. Water can slowly freeze in these channels where the slow freezing excludes air bubbles giving the ice a clear aspect. This ice is uplifted to the surface ablation zone by glacial movements and can therefore be observed in the form of clear surface ice bands. We employed an indirect method to sample englacial water by coring these ice bands. We were able, for the first time, to compare microbial communities sampled from clear (i.e. frozen englacial water bands) and cloudy ice (i.e. meteoric ice) through 16S rRNA gene sequencing. Although microbial communities were primarily shaped and structured by their spatial distribution on the glacier, ice type was a clear secondary factor. One area of the glacier, in particular, presented significant microbial community clear/cloudy ice differences. Although the clear ice and supraglacial communities showed typical cold-adapted glacial communities, the cloudy ice had a less defined glacial community and ubiquitous environmental organisms. These results highlight the role of englacial channels in the microbial dispersion within the glacier and, possibly, in the shaping of glacial microbial communities.

Marine-terminating glaciers lose mass through melting and iceberg calving, and we find that meltwater drainage systems influence calving timing at Helheim Glacier, a tidewater glacier in East Greenland. Meltwater feeds a buoyant subglacial discharge plume at the terminus of Helheim Glacier, which rises along the glacial front and surfaces through the mélange. Here, we use high-resolution satellite and time-lapse imagery to observe the surface expression of this meltwater plume and how plume timing and location compare with that of calving and supraglacial meltwater pooling from 2011 to 2019. The plume consistently appeared at the central terminus even as the glacier advanced and retreated, fed by a well-established channelized drainage system with connections to supraglacial water. All full-thickness calving episodes, both tabular and non-tabular, were separated from the surfacing plume by either time or by space. We hypothesize that variability in subglacial hydrology and basal coupling drive this inverse relationship between subglacial discharge plumes and full-thickness calving. Surfacing plumes likely indicate a low-pressure subglacial drainage system and grounded terminus, while full-thickness calving occurrence reflects a terminus at or close to flotation. Our records of plume appearance and full-thickness calving therefore represent proxies for the grounding state of Helheim Glacier through time. [ABSTRACT FROM AUTHOR]

We present the first fully coupled 3D full-Stokes model of a tidewater glacier, incorporating ice flow, subglacial hydrology, plume-induced frontal melting and calving. We apply the model to Store Glacier (Sermeq Kujalleq) in west Greenland to simulate a year of high melt (2012) and one of low melt (2017). In terms of modelled hydrology, we find perennial channels extending 5 km inland from the terminus and up to 41 and 29 km inland in summer 2012 and 2017, respectively. We also report a hydrodynamic feedback that suppresses channel growth under thicker ice inland and allows water to be stored in the distributed system. At the terminus, we find hydrodynamic feedbacks exert a major control on calving through their impact on velocity. We show that 2012 marked a year in which Store Glacier developed a fully channelised drainage system, unlike 2017, where it remained only partially developed. This contrast in modelled behaviour indicates that tidewater glaciers can experience a strong hydrological, as well as oceanic, control, which is consistent with observations showing glaciers switching between types of behaviour. The fully coupled nature of the model allows us to demonstrate the likely lack of any hydrological or ice-dynamic memory at Store Glacier. [ABSTRACT FROM AUTHOR]

We conduct extended versions of the ISMIP6 future climate experiments for the Greenland ice sheet until the year 3000 with the model SICOPOLIS. Beyond 2100, the climate forcing is kept fixed at late-21st-century conditions. For the unabated warming pathway RCP8.5/SSP5-8.5, the ice sheet suffers a severe mass loss, which amounts to ~ 1.8 m SLE (sea-level equivalent) for the 12-experiment mean, and ~ 3.5 m SLE (~ 50% of the entire mass) for the most sensitive experiment. For the reduced emissions pathway RCP2.6/SSP1-2.6, the mass loss is limited to a two-experiment mean of ~ 0.28 m SLE. Climate-change mitigation during the next decades will therefore be an efficient means for limiting the contribution of the Greenland ice sheet to sea-level rise in the long term. [ABSTRACT FROM AUTHOR]

Basal conditions play an essential role in the dynamics of outlet glaciers, but direct observations at the bed of glaciers are challenging to obtain. Instead, inverse methods can be used to infer basal parameters from surface observations. Here, we use a simple ice-flow model as a forward model in an inversion scheme to retrieve the spatio-temporally variable basal stress parameter for Hagen Bræ, North Greenland, from 1990 to 2020. Hagen Bræ is a surge-type glacier with up to an order of magnitude variability of winter velocities near the grounding line. We find that downstream changes in the basal stress parameter can explain most of the variation of flow velocity, and we further identify a region of high resistance ~20–40 km from the grounding line. We hypothesise that this region of high resistance plays an important role in controlling glacier discharge. [ABSTRACT FROM AUTHOR]

Proglacial braided river systems discharge large volumes of meltwater from ice sheets and transport coarse-grained sediments from the glaciated areas to the oceans. Here, we test the hypothesis if high-energy hydrological events can leave distinctive signatures in the sedimentary record of braided river systems. We characterize the morphology and infer a mode of formation of a 25 km long and 1–3 km wide Early Pleistocene incised valley recently imaged in 3-D seismic data in the Hoop area, SW Barents Sea. The fluvial system, named Bjørnelva River Valley, carved 20 m deep channels into Lower Cretaceous bedrock at a glacial paleo-surface and deposited 28 channel bars along a paleo-slope gradient of ~0.64 m km−1. The landform morphologies and position relative to the paleo-surface support that Bjørnelva River Valley was formed in the proglacial domain of the Barents Sea Ice Sheet. Based on valley width and valley depth, we suggest that Bjørnelva River Valley represents a braided river system fed by violent outburst floods from a glacial lake, with estimated outburst discharges of ~160 000 m3 s−1. The morphological configuration of Bjørnelva River Valley can inform geohazard assessments in areas at risk of outburst flooding today and is an analogue for landscapes evolving in areas currently covered by the Greenland and Antarctic ice sheets.

As a result of climate warming, glacial meltwater discharge has been increasing in Greenland. During the summers of 2015 and 2016, there were rapid increases in discharge from Qaanaaq Glacier in northwestern Greenland. These discharges resulted in floods that destroyed the road linking the settlement of Qaanaaq to Qaanaaq Airport. Field measurements were performed and a numerical model of glacier runoff was developed to quantify these discharges. The high discharge associated with the 2015 flood, estimated at 9.1 m3 s−1 (hourly mean), resulted from intensive glacier melting due to warm air temperature and strong winds, while the high discharge associated with the 2016 flood resulted from heavy rainfall (90 mm d−1) that led to a peak discharge estimated at 19.9 m3 s−1. The developed model, when used to investigate future glacier runoff under warming conditions, revealed a nonlinear increase in glacial melt with increasing temperature. Additionally, the model forecasted a threefold increase in total summer discharge, owing to a 4 °C rise in temperature. Thus, this study quantified the impact of a changing climate on glacier runoff, which gives insight into future risks of flood hazards along the coast of Greenland.