Recent advances in remote sensing have produced global glacier surface elevation change data. Parsing these elevation change signals into contributions from the climate (i.e. climatic mass balance) and glacier dynamics (i.e. flux divergence) is critical to enhance our process-based understanding of glacier change. In this study, we evaluate three approaches for direct, continuous measurements of the climatic mass balance, flux divergence and elevation change at a site on Gulkana Glacier in Alaska using low-cost global navigation satellite system (GNSS) sensors, GNSS interferometric reflectometry (GNSS-IR), banded ablation stakes with time-lapse cameras and combinations thereof. Cumulative climatic mass balance over the season was 4.85 m and the three approaches were within 0.08 m through early July before the snowpack melted, and within 0.28 m through mid-August. The flux divergence increased from 0.52 ± 0.03 cm d−1 before June 3 to about 0.73 cm d−1 after June 27. We demonstrate a single GNSS system fixed atop an ablation stake can measure contemporaneous climatic mass balance, flux divergence and elevation change based on the antenna's position and GNSS-IR techniques. The ability of these systems to measure glacier mass balance and flux divergence offers unique opportunities for year-round observations on mountain glaciers in the future.

Tidewater glaciers frequently advance and retreat in ways uncoupled from climate forcing. This complicates the task of forecasting the evolution of individual glaciers and the overall Greenland ice sheet, much of which is drained by tidewater glaciers. Past observational research has identified a set of processes collectively known as the tidewater glacier cycle (TGC) to describe tidewater glacier evolution in four stages: the advancing stage, the extended stage, the retreating stage and the retreated stage. Once glacier retreat is initiated, the TGC is thought to depend largely on the glacier's calving rate, which is controlled by fjord geometry. However, there has been little modeling or systematic observational work on the topic. Measuring calving rates directly is challenging and thus we developed an averaged von Mises stress state at the glacier terminus as a calving rate proxy that can be estimated from surface velocities, ice thickness, a terminus position and subglacial topography. We then analyzed 44 tidewater glaciers in Greenland and assessed the current state in the TGC for them. Of the 44 glaciers, we find that fjord geometry is causing instability in ten cases, vs stability in seven, with 11 in rapid retreat and 16 have been historically stable.

Synthetic Aperture Radar (SAR) has been used extensively to determine the surface ice flow velocity of tidewater glaciers and investigate changes in seasonal or annual ice dynamics at medium spatial resolution (⩾100 m). However, assessing tidewater glacier behaviour at these resolutions risks missing key details of glacier dynamics, which is particularly important for determination of strain rates that relate to crevasse formation, depth, and ice damage. Here we present surface ice velocity and strain maps with a 16 m posting derived from high-resolution (1 m) PAZ Ciencia spotlight mode SAR imagery for Narsap Sermia, SW Greenland, for October 2019 to February 2021. Results reveal fine details in strain rate, including an area of compression proximal to the terminus, with an upstream shift of strains through time. The velocity evolution of Narsap Sermia shows distinct seasonal changes starting in summer 2020, which are largely modulated by the subglacial drainage system. Comparison of our results with medium-resolution velocity products shows that while these can capture general strain and velocity patterns, our high-resolution data reveals considerably larger ranges of strain values. This is likely to have implications for tuning strain rate dependent calving and ice damage parameterisations within numerical models.

Uncrewed aerial vehicles (UAVs) are frequently used in glaciological applications, among other things, for photogrammetric assessments of calving dynamics at glacier termini. However, UAVs are often limited by battery endurance and weight constraints on the scientific payload that can be added. At Sálajiegna, the largest freshwater calving glacier in Sweden, we explored the combined use of a versatile maritime robot (uncrewed surface vehicle, USV) and a UAV to characterise Sálajiegna's short-term and seasonal calving front dynamics and mass loss. For this, a photogrammetric payload suite was integrated into the USV. Consecutive USV surveys of Sálajiegna's front, followed by point cloud based calving detection and surface-reconstruction based volume quantification, allowed for a detailed description of calving-induced terminus changes and is hence suggested as a viable alternative to the differencing of digital elevation models. By combining USV and UAV measurements, we identify sectors of high and low calving activity, a calving front retreat of up to 56 m and a thinning rate in the terminus region of 5.4 cm d−1 during the summer of 2022.

We demonstrate that a ~20 km long valley glacier in the St. Elias Mountains, Yukon, can experience both partial and full surges, likely controlled by the presence of a topographic constriction and the formation and drainage of supraglacial lakes. Based on analysis of air photos, satellite images and field observations since the 1940s, we identify a full surge of ‘Little Kluane Glacier’ from 2013 to 2018, and a partial surge of just the upper north arm between 1963 and 1972. Repeat digital elevation models and velocity profiles indicate that the recent surge initiated from the upper north arm in 2013, which developed into a full surge of the main trunk from 2017 to 2018 with peak velocities of ~3600 m a−1 and frontal advance of ~1.7 km from May to September 2018. In 2016, a mass movement from the north arm to the main trunk generated a surface depression in a region immediately downstream of a topographic constriction, which promoted the formation and rapid drainage of supraglacial lakes to the glacier bed, and likely established the conditions to propel the initial partial surge into a full surge. Our results underscore the complex interplay between glacier geometry, surface hydrology and topography required to drive full surges of this glacier.

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

Calving glaciers respond quickly to atmospheric variability through ice dynamic adjustment. Particularly, single weather extremes may cause changes in ice-flow velocity and terminus position. Occasionally, this can lead to substantial event-driven mass loss at the ice front. We examine changes in terminus position, ice-flow velocity, and calving flux at the grounded lacustrine Schiaparelli Glacier in the Cordillera Darwin using geo-referenced time-lapse camera images and remote sensing data (Sentinel-1) from 2015 to 2022. Lake-level records, lake discharge measurements, and a coupled energy and mass balance model provide insight into the subglacial water discharge. We use downscaled reanalysis data (ERA5) to identify climate extremes and track land-falling atmospheric rivers to investigate the ice-dynamic response on possible atmospheric drivers.

Meltwater controls seasonal variations in ice-flow velocity, with an efficient subglacial drainage system developing during the warm season and propagating up-glacier. Calving accounts for 4.2% of the ice loss. Throughout the year, warm spells, wet spells, and landfalling atmospheric rivers promote calving. The progressive thinning of the ice destabilizes the terminus position, highlighting the positive feedback between glacier thinning, near-terminus ice-flow acceleration, and calving flux.

The future of tidewater glaciers in response to climate warming is one of the largest sources of uncertainty in the contribution of the Greenland ice sheet to global sea-level rise. In this study, we investigate the ability of an ice-sheet model to reproduce the past evolution of the velocity and surface elevation of a tidewater glacier, Upernavik Isstrøm, by prescribing front positions. To achieve this, we run two ensembles of simulations with a Weertman and a regularised-Coulomb friction law. We show that the ice-flow model has to include a reduction in friction in the first 15 km upstream of the ice front in fast-flowing regions to capture the trends observed during the 1985–2019 period. Without this process, the ensemble model overestimates the ice flow before the retreat of the front in 2005 and does not fully reproduce its acceleration during the retreat. This results in an overestimation of the total mass loss between 1985 and 2019 of 50% (300 vs 200 Gt). Using a variance-based sensitivity analysis, we show that uncertainties in the friction law and the ice-flow law have a greater impact on the model results than surface mass balance and initial surface elevation.

The collapse of several ice shelves in the Antarctic Peninsula since the late 20th century has resulted in the upstream acceleration of multiple formerly buttressed outlet glaciers, raising questions about the stability of Antarctica's remaining ice shelves and the effects their demise may have upon inland ice. Here, we use high temporal resolution Sentinel-1A/B synthetic aperture radar-derived observations to assess the velocity response of Larsen C Ice Shelf (LCIS) to the calving of colossal iceberg A-68 in 2017. We find marked oscillations in ice-shelf flow across LCIS in the months following A-68's calving, beginning with a near-ice-shelf-wide slowdown of 11.3 m yr−1 on average. While falling close to the limits of detectability, these ice-flow variations appear to have been presaged by similar oscillations in the years prior to A-68's breakaway, associated primarily with major rifting events, together reflecting potentially hitherto unobserved ice-shelf mechanical processes with important implications for ice-shelf weakening. Such ice-flow oscillations were, however, short-lived, with more recent observations suggesting a deceleration below longer-term rates of ice flow. Collectively, our observations reveal complex spatial-temporal patterns of ice-flow variability at LCIS. Similarly abrupt fluctuations may have important implications for the stability of other ice shelves, necessitating the continued, close observation of Antarctica's coastline in the future.

We investigate unusual discontinuous glacier motion on Thompson Glacier, Umingmat Nunaat, Arctic Canada, using synthetic aperture radar (SAR) images and ice-flow modeling. A novel intensity-rescaling scheme is developed to reduce errors in high-resolution speckle tracking, resulting in a ~25% improvement in accuracy. Interferometric SAR (InSAR) and speckle tracking using high resolution RADARSAT-2 data indicate velocity discontinuities of up to 1 cm d−1 across deep and longitudinally extensive supraglacial channels on Thompson Glacier. We use a cross-sectional finite-element ice-flow model to determine the conditions under which velocity discontinuities of the observed magnitude and signature are possible. The modeling suggests that discontinuous motion across (long and straight) supraglacial channels can occur without ice fracture and under a wide variety of glacier thermal structures, including in fully temperate glaciers. Despite the wide range of conditions conducive to discontinuous motion, the form we observe requires that the associated channels be deep, longitudinally extensive and located in regions of lateral shearing. We speculate that these combined conditions are rare except on polythermal glaciers, where drainage features such as moulins are comparatively scarce and lower deformation rates allow channels to incise consistently and persist over many years.

The debris that covers the ablation areas of high-elevation debris-covered glaciers contributes to the distinctive features and processes occurring both on and within such glaciers. Despite recent advances, knowledge of the subsurface environments of high-elevation debris-covered glaciers is still extremely limited. In particular, targeted field-based data are needed to parameterise and refine the projections of these glaciers in numerical models. Here, we outline the current understanding of the internal properties of high-elevation debris-covered glaciers based on direct field-based methods and suggest future research directions for field-based studies.

Antarctica's ice shelves stabilize the ice sheet and, therefore, understanding processes affecting the mass budgets of ice shelves is important for estimating grounded ice loss. To study the ice shelf dynamics, we analyzed seismological and GNSS data from the Ekström Ice Shelf in Dronning Maud Land. We extracted probabilistic power spectral densities (PPSD) in the frequency band 3.4–6.8 Hz, typical of icequakes, from seismological data and observed pronounced signals in the PPSD with near 3 and 4 cycles per day (cpd) corresponding to tidal overharmonics, in addition to the main tidal constituents near 1 and 2 cpd. GNSS data reveal the same components in ice flow speed but not in vertical displacements. Generally, tide-induced grounding line migration modulates the flow velocity of an entire ice shelf. We find that this velocity modulation causes the increased icequake activity in the tidal overharmonics with 3 and 4 cpd in an ice shear zone where the flow velocity drops to nearly zero.

Spatially variable basal conditions are thought to govern how ice sheets behave at glacial time scales (>1000 years) and responsible for changes in dynamics between the core and peripheral regions of the Laurentide and Fennoscandian ice sheets. Basal motion is accomplished via the deformation of unconsolidated sediments, or via sliding of the ice over an undeformable bed. We present an ice sheet sliding module for the Parallel Ice Sheet Model (PISM) that takes into account changes in sediment cover and incorporates surface meltwater. This model routes meltwater, produced at the surface and base of the ice sheet, toward the margin of the ice sheet. Basal sliding is accomplished through the deformation of water saturated sediments, or sliding at the ice-bed interface. In areas with continuous, water saturated sediments, sliding is almost always accomplished through sediment deformation. In areas with incomplete cover, sliding has a stronger dependence on the supply of water. We find that the addition of surface meltwater to the base is a more important factor for ice sheet evolution than the style of sliding. In a glacial cycle simulation, our model causes a more rapid buildup of the Laurentide Ice Sheet.

Links between proglacial lakes and glacier dynamics are poorly understood but are necessary to predict how mountain glaciers will react to a warmer, wetter climate, where such lakes are expected to increase both in number and volume. Here, we examine a long-term (~120 year) record of terminus retreat, thinning and surface velocities from in-situ and remote sensing observations at the terminus of Kaskawulsh Glacier, Yukon, Canada, and determine the impact of a local proglacial hydrological reorganisation on glacier dynamics. After an initial deceleration during the late 1990s, terminus velocities increased at a rate of 3 m a−2 from 2000–12, while proglacial Slims Lake area increased simultaneously. The rapid drainage of the lake in May 2016 substantially altered the velocity profile, decreasing annual velocities by 48% within 3 km of the terminus between 2015 and 2021, at an average rate of ~ 12.5 m a−2. A key cause of the rapid drop in glacier motion was a reduction in flotation of the lower part of the glacier terminus after lake drainage. This has important implications for glacier dynamics and provides one of the first assessments of the impacts of a rapid proglacial lake drainage event on local terminus velocities.

The magnitude and azimuth of horizontal ice flow at Camp Century, Greenland have been measured several times since 1963. Here, we provide a further two independent measurements over the 2017–21 period. Our consensus estimate of horizontal ice flow from four independent satellite-positioning solutions is 3.65 ± 0.13 m a−1 at an azimuth of 236 ± 2°. A portion of the small, but significant, differences in ice velocity and azimuth reported between studies likely results from spatial gradients in ice flow. This highlights the importance of restricting inter-study comparisons of ice flow estimates to measurements surveyed within a horizontal distance of one ice thickness from each other. We suggest that ice flow at Camp Century is stable on seasonal to multi-decadal timescales. The airborne and satellite laser altimetry record indicates an ice thickening trend of 1.1 ± 0.3 cm a−1 since 1994. This thickening trend is qualitatively consistent with previously inferred ongoing millennial-scale ice thickening at Camp Century. The ice flow divide immediately north of Camp Century may now be migrating southward, although the reasons for this divide migration are poorly understood. The Camp Century flowlines presently terminate in the vicinity of Innaqqissorsuup Oqquani Sermeq (Gade Gletsjer) on the Melville Bay coast.

Future mass loss projections of the Greenland ice sheet require understanding of the processes at a glacier terminus, especially of iceberg calving. We present detailed and high-rate terrestrial radar interferometer observations of Eqip Sermia and Bowdoin Glacier, two outlet glaciers in Greenland with comparable dimensions and investigate iceberg calving, surface elevation, velocity, strain rates and their links to air temperature, tides and topography. The results reveal that the two glaciers exhibit very different flow and calving behaviour on different timescales. Ice flow driven by a steep surface slope with several topographic steps leads to high velocities, areas of extension and intense crevassing, which triggers frequent but small calving events independent of local velocity gradients. In contrast, ice flow under smooth surface slopes leaves the ice relatively intact, such that sporadic large-scale calving events dominate, which initiate in areas with high shearing. Flow acceleration caused by enhanced meltwater input and tidal velocity variations were observed for terminus sections close to floatation. Firmly grounded terminus sections showed no tidal signal and a weak short-term reaction to air temperature. These results demonstrate reaction timescales to external forcings from hours to months, which are, however, strongly dependent on local terminus geometry.

Accurate ice flow velocity data are essential for studying the mass balance of the Antarctic ice sheet. However, there is a lack of ice velocity maps of 1960s–80s in basin-wide regions or the entire ice sheet. In this study, an enhanced hierarchical network densification approach is developed for basin-wide Antarctic velocity mapping using historical ARGON and Landsat images. The produced multiple historical velocity maps from 1963 to 1989 in the region of the Fimbul and Jelbart ice shelves, East Antarctica, achieved an accuracy better than 29 m a−1. They revealed that the ice flow velocity had no significant changes over the period. Combining the surface mass balance estimate with the ice discharge estimated from our historical velocity maps and recently published velocity maps, we estimated a positive mass balance of 8.6 ± 3.9 Gt a−1 in the study area from 1963 and 2015. Our results indicate that the region's positive mass balance, as estimated in recently published studies, has been maintained since the 1960s. It is also in concordance with the low level of mass balance from 1992 to 2017 in East Antarctica. This suggests that the study area has been stable since the 1960s.

Khovsgol Lake, a seasonal ice-covered lake located in boreal Siberia, plays important economic, transportation and agricultural roles for local residents. On 11 January 2021, an Mw 6.7 earthquake struck the center of this lake, providing a unique opportunity to better understand lake ice movement, ridge distribution, and hydraulic and hydrodynamic processes, as a result of the mainshock. We use a pixel-tracking method on both Sentinel-1 synthetic aperture radar and Sentinel-2 optical remote-sensing data covering this lake before, during and after the mainshock to obtain lake-ice displacement. Combining these measurements with different viewing geometries, we determine the locations of pressure ridges in each period and derive their relative 3D displacement maps with a precision of <2.2 m. Our results indicate a large E-W displacement of up to 10.3 m and a vertical displacement of up to 4.6 m during the coseismic period, which has a distinctly different magnitude and pattern of displacement compared to the pre- and post-event periods. Given the ice motion pattern and magnitude, we speculate that the unique ice displacement of Khovsgol Lake may have been caused by a rare seismic seiche after the mainshock.

Since the mid-20th century, ice shelves around the Antarctic Peninsula have declined in extent and thickness, and some have shown signs of structural instability. Here, using satellite imagery from 1999/2000 to 2019/20 (Landsat 7 and 8, Sentinel-2 and ASTER), we measure areal changes, calculate surface flow speeds, and quantify structural changes of Bach, Stange and George VI ice shelves, located in the southwest Antarctic Peninsula. We recorded a total area loss of 797.5 km2 from 2009/10 to 2019/20, though spatial and temporal patterns varied at individual ice fronts. Flow speeds remained largely stable over the observation periods, but notable acceleration was calculated for Bach Ice Shelf, and at the northern and southern extents of George VI Ice Shelf. Open fractures widened and lengthened between 2009/10 and 2019/20 on all three ice shelves. We conclude that Stange Ice Shelf is stable, and not under any immediate threat of enhanced recession. Continued ice-mass loss and consequential speed up of George VI South may cause further fracturing and destabilisation in the coming decades. Of more immediate concern are the glaciological changes noted for Bach Ice Shelf and George VI North; substantial areas of stabilising ice have already, or will soon be removed, that may lead to enhanced recession within the next decade.

This paper introduces the Instructed Glacier Model (IGM) – a model that simulates ice dynamics, mass balance and its coupling to predict the evolution of glaciers, icefields or ice sheets. The novelty of IGM is that it models the ice flow by a Convolutional Neural Network, which is trained from data generated with hybrid SIA + SSA or Stokes ice flow models. By doing so, the most computationally demanding model component is substituted by a cheap emulator. Once trained with representative data, we demonstrate that IGM permits to model mountain glaciers up to 1000 × faster than Stokes ones on Central Processing Units (CPU) with fidelity levels above 90% in terms of ice flow solutions leading to nearly identical transient thickness evolution. Switching to the GPU often permits additional significant speed-ups, especially when emulating Stokes dynamics or/and modelling at high spatial resolution. IGM is an open-source Python code which deals with two-dimensional (2-D) gridded input and output data. Together with a companion library of trained ice flow emulators, IGM permits user-friendly, highly efficient and mechanically state-of-the-art glacier and icefields simulations.