Outlet glaciers in Greenland are undergoing retreat and diffusive thinning in response to external forcings, but the rates and magnitudes of these responses differ from glacier to glacier for unclear reasons. We test how changes in ice overburden pressure and basal lubrication affect diffusive thinning rates and their spatial patterns by conducting numerical experiments over various idealized Greenland-like glacier domains. We find that ~10 km frontal retreat over a decade can produce sustained thinning rates as large as 16 m a−1 due to ice overburden pressure changes, at outlet glaciers with high basal drag (>60 kPa) and lateral resistive stress (>70 kPa). Localized basal lubrication perturbations induce upstream thinning and downstream thickening up to 12 m a−1; the duration of the lubrication forcing generally has a greater effect than its intensity on induced thickness changes. Lastly, episodic grounding line retreats over a rough bed produce a stepped time series of thinning broadly consistent with observations of dynamic elevation change on multiple Greenland glaciers. Our findings highlight the critical role of the total grounding zone – not ice front position – through the resistive stress change in relation to total glacier thinning.

A surge of Mittie Glacier, a 50 km-long tidewater outlet glacier of Manson Icefield, occurred between approximately 1992 and 2007. Velocities increased slowly at first, but then increased rapidly to reach a peak of 4800 m a−1 in early 1996, the highest ever reported for a glacier in the Canadian Arctic. The surge initiated at the terminus and propagated up-glacier, with a maximum terminus advance of 7.3 km between 1994 and 1999. The surge was almost entirely restricted to the lower ~30 km of the glacier, in a region which radio-echo sounding shows to be grounded below sea level. A 3 km-wide crevasse with a 150 m opening occurred at the separation between faster moving ice downstream and slower moving ice upstream. Surge initiation appears to have been triggered by flotation of the lower terminus, caused by long-term thinning of this region during quiescence.

Glaciers in High-Mountain Asia are experiencing varying rates and patterns of mass loss due to a complex interplay between glacier surface processes, local conditions and climate forcing. Spatially distributed surface mass balance (SMB) estimates can provide valuable insight into these drivers, but observations are currently limited in both space and time. We used very-high-resolution optical stereo images acquired by commercial satellites to prepare time series of digital elevation models (DEMs), and derived contemporaneous surface velocity and elevation change products for six debris-covered glaciers in Nepal. We developed new methods to produce flow-corrected Lagrangian SMB maps to isolate local surface ablation signals with enough detail to study individual ice cliffs. Our results show reduced ablation under thick debris cover and enhanced ablation over ice cliffs. Ablating ice cliffs were responsible for $10\!-\!38\%$ of the total ablation over debris-covered areas, even though they covered $\leq \!11\%$ of the total area. Seasonal SMB products reveal the timing and patterns of summer accumulation and ablation, underscoring the importance of snow avalanches for low-elevation debris-covered glaciers in the region. Our approach can be applied to other glaciers with repeat high-resolution DEM coverage and extended for regional analyses of SMB on seasonal to interannual timescales.

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.

The flow of glacial ice is typically approximated as a nonNewtonian viscous fluid, with the momentum balance described by (an approximation to) the Stokes equations, and the nonlinear rheology described by a flow law. The most commonly used rheological law for glacial ice, Glen's flow law, yields infinite viscosity in the case of zero deformation, which can be the case at the ice surface. This poses a problem when solving the momentum balance numerically. We show that two commonly-used discretisation schemes for the boundary conditions at the ice surface and base, which yield proper numerical convergence when applied to simpler problems, produce poor numerical convergence and large errors, when used to solve the momentum balance with Glen's flow law. We show that a discretisation scheme based on the concept of ghost nodes, which substitutes the boundary conditions directly into the momentum balance equations, yields second-order numerical convergence and errors that can be up to four orders of magnitude smaller. These results are robust across different momentum balance approximations. We show that the improved boundary conditions are particularly useful for solving the 3-D higher-order Blatter-Pattyn Approximation (BPA). In general, this work underlines the importance of thoroughly verifying the numerical solvers used in ice-sheet models, before applying them to future projections of ice-sheet mass loss.

The observed retreat of the grounding line of the present-day ice sheets and the simulated grounding-line retreat of ice sheets under changing climate conditions are often interpreted as indications of marine ice-sheet instability (MISI), driven by a positive feedback between the ice discharge and conditions at the grounding line. However, the arguments that support this feedback are valid only for steady-state conditions. Here, we assess how unconfined marine ice sheets may behave if atmospheric conditions and basal conditions evolve with time. We find that the behavior of grounding lines can exhibit a range from unstoppable advance and retreat to irregular oscillation irrespective of the stability of the corresponding steady-state configurations obtained with time-invariant conditions. Our results show that numerical simulations with a parameterization of the ice flux through the grounding line used in large-scale ice-sheet models produce markedly different results from simulations without the parameterization. Our analysis demonstrates that the grounding-line migration can be driven by the temporal variability in the atmospheric and basal conditions and not by MISI, which assumes unchanging conditions. Instead, the grounding-line advance or retreat is determined by interactions between ice flow, basal processes and environmental conditions throughout the length of a marine ice sheet in addition to the circumstances at its grounding line.

Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.

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.

We present a new map of bed topography and ice thickness together with a corresponding ice volume estimate representative of the years ~2010 for all Scandinavian ice caps and glaciers. Starting from surface observations, we invert for ice thickness by iteratively running an innovative ice dynamics model on a distributed grid and updating bed topography until modelled and observed glacier dynamics as represented by their rate of surface elevation change (dh/dt) fields align. The ice flow model used is the instructed glacier model (Jouvet and Cordonnier, 2023, Journal of Glaciology 1–15), a generic physics-informed deep-learning emulator that models higher-order ice flow with high-computational efficiency. We calibrate the modelled thicknesses against >11 000 ice thickness observations, resulting in a final ice volume estimate of 302.7 km3 for Norway, 18.4 km3 for Sweden and 321.1 km3 for the whole of Scandinavia with an error estimate of ~$\pm 11\%$. The validation statistics computed indicate good agreement between modelled and observed thicknesses (RMSE = 55 m, Pearson's r = 0.87, bias = 0.8 m), outperforming all other ice thickness maps available for the region. The modelled bed shapes thus provide unprecedented detail in the subglacial topography, especially for ice caps where we produce the first maps that show ice-dynamically realistic flow features.

Ice tongues are unconfined by land on their lateral margins and are sensitive to external forcing from the ocean. They are found sporadically around the Antarctic coast but are common in the western Ross Sea. Lateral flexure creates bending stresses within these ice tongues which is likely to contribute to their fragility and may restrict their spatial distribution. A three-year time series (2017–2020) of synthetic aperture radar interferometry (InSAR) and differential interferometry (DInSAR) is used to observe the lateral flexure of the 10 km long Erebus Ice Tongue as a result of ocean currents. The fast ice area around the ice tongue was mapped during the same period. It was found that when fast ice was absent (34.7% of the time), the average lateral movement of the ice tongue was twice as high (0.44 m) as when it was embedded in fast ice (0.19 m). A significant correlation (0.45) between flexure and tidal currents was found when fast ice was absent. An analytical model tuned to observations suggests that even without sea ice for stabilisation, the lateral bending stresses induced by the ocean are insufficient to cause calving without additional amplifying factors.

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.

We use satellite image processing techniques to measure surface elevation and velocity changes on a temperate surging glacier, Sít’ Kusá, throughout its entire 2013–2021 surge cycle. We present detailed records of its dynamic changes during quiescence (2013–2019) and its surge progression (2020–2021). Throughout quiescence, we observe order-of-magnitude speedups that propagate down-glacier seasonally from the glacier's upper northern tributary, above a steep icefall, into the reservoir zone for the surging portion of the glacier. The speedups initiate in fall and gradually accelerate through winter until they peak in late spring, ~1 − 2 months after the onset of melt. Propagation distance of the speedups controls the distribution of mass accumulation in the reservoir zone prior to the surge. Furthermore, the intensity and propagation distance of each year's speedup is correlated with the positive degree day sum from the preceding melt season, suggesting that winter melt storage drives the seasonal speedups. We demonstrate that the speedups are kinematically similar to the 2020–2021 surge, differing mainly in that the surge propagates past the dynamic balance line at the lower limit of the reservoir zone, likely triggered by the exceedance of a tipping point in mass accumulation and basal enthalpy in the reservoir zone.

Perennial ice can be studied for many purposes, including paleoclimate records or rheological properties. For most of those purposes, the ice microstructure must be studied, often through optical microscopy. The aim of this work is to assess the viability of immersion microscopy for the study of ice microstructures. It consists of using an oil between the objective lens and the specimen, to increase image resolution. Immersion microscopy is a technique well-developed for the investigation of diverse materials, but it has so far not been explored for ice research. Here we investigate the challenges and advantages of that technique. The main challenge is related to the selection of the immersion oil itself, which must satisfy a number of criteria, ranging from refractive index and viscosity to toxicity and reactivity. We identify pure silicone oil (dimethicone) as a simple and safe option for immersion microscopy of inner ice structures. Among its advantages, it provides higher resolution (compared to standard ‘dry’ microscopy) and it can be simultaneously used as a long-term coating to prevent undesired sublimation of the ice-sample surfaces. For the observation of surface structures, however, another type of oil with higher refractive index should be used.

This study examines the performance of low-cost, low-power GNSS positioning systems for glacier monitoring in high-latitude environments. We compare the positioning performance of co-located low-cost u-blox ZED-F9P GNSS units (a few hundred USDs) and survey-grade Trimble R10 units (> $10,000 USD) under stationary (on land) and dynamic (on glacier) conditions near Terra Nova Bay, Antarctica. Low-cost and survey-grade systems yield almost identical error magnitudes under short (3 m), medium (34 km) and long (390 km) baseline kinematic-positioning scenarios. We further examined the efficacy of low-cost GNSS for glaciological applications by installing four u-blox and two Trimble receivers on Priestley Glacier to observe tide-modulated ice flexure. All receivers successfully detected subtle tidal oscillations with amplitudes < 3 cm, consistent with the predicted phasing from a tide model. These experiments offer a strong rationale for the widespread use of low-cost receivers to expand and densify GNSS monitoring networks, both in Antarctica and in glaciated regions worldwide.

Both submarine melt and calving are important for the overall mass balance of marine-terminating glaciers, but uncertainty is rife with regards to the magnitude of the processes. Modelling allows for these processes to be investigated without the need to visit inaccessible ice marginal zones. This study looks at the impact of different submarine melt and sea-ice back pressure scenarios on modelled calving activity and dynamics at Kronebreen, Svalbard, by running separate summer and winter simulations with various submarine melt parameterisations and sea-ice characteristics. It is found that submarine melt is an important driver of seasonal variation in modelled glacier dynamics and calving activity, with the choice of sliding law also exerting a significant influence on results.

Monte Perdido Glacier, located in the central Pyrenees, is one of the southernmost glaciers in Europe. Due to climate change, this glacier is suffering an accelerated mass loss, especially in the last decades. If the current trends persist, this glacier is expected to disappear in the next 50 years. As part of the efforts of the scientific community to increase the knowledge about this glacier, this research presents the first microstructural characterization of the Monte Perdido Glacier, focused on a high-impurity concentration segment that belongs to an ice core drilled in 2017. The results reveal the ice has a layering defined by air bubbles and non-soluble impurities. The bubble-defined layering exhibits features of both a primary (sedimentary) and a secondary (strain-induced) origin. We found a clear inverse correspondence between the particle concentration and the grains' size and roundness index. A preliminary micro-Raman characterization of the particles shows the occurrence of atacamite, anatase (likely related to ancient mining activities in the vicinity of the glacier) and quartz. The latter could be an indicator of mineral dust, probably suggesting the arrival of dust-laden air masses from the north of the African continent.

From 2017 to 2020, three significant calving events took place on Pine Island Glacier, West Antarctica. Ice-shelf velocities changed over this period and the calving events have been suggested as possible drivers. However, satellite observations also show significant changes in the areal extent of fracture zones, especially in the marginal areas responsible for providing lateral support to the ice shelf. Here, we conduct a model study to identify and quantify drivers of recent ice-flow changes of the Pine Island Ice Shelf. In agreement with recent studies, we find that the calving events caused significant velocity changes over the ice shelf. However, calving alone cannot explain observed velocity changes. Changes in the structural rigidity, i.e. ice damage, further significantly impacted ice flow. We suggest that ice damage evolution of the ice-shelf margins may have influenced recent calving events, and these two processes are linked.

The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.