We developed automated ablation stakes to measure colocated in situ changes in relative glacier-surface elevation and climatological drivers of ablation. The designs, refined over 10 years of development and deployments, implement open-source hardware and common building materials. The ablation stakes record distance to the snow/ice surface, air temperature and relative humidity every 1–15 min. Using these high-frequency data, we demonstrate that melt factors calculated using standard melt-rate vs temperature regressions converge over averaging windows of approximately 12 h or greater. Furthermore, we evaluate an integral approach to estimating temperature-index melt factors for ablation. In a test case on Glaciar Perito Moreno, Argentina, this integral approach reveals an overall positive-degree-day melt factor of 7.5 mm w.e. $^\circ$C−1 d−1. We describe four deployments with iteratively improved designs and provide a list of materials required to construct an automated ablation stake.

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.

The number of glacial lakes has grown globally concurrently with the retreat of glaciers in the last few decades, increasing the risk of potentially hazardous glacial lake outburst floods (GLOFs) and posing a threat to downstream communities. Norway has several known ice-dammed lakes that produce repeated GLOFs but as we show here, the existing GLOF database is incomplete and needs to be improved through continuous monitoring of glaciers and glacial lakes. This study examines the case of an ice cap in central Norway hosting at least four drainage-prone lakes. We reconstruct the sequence of lake drainage events through a combination of remote sensing, ground-truthing and citizen science while scrutinising the applicability of the PlanetScope imagery vs Sentinel-2 and Landsat-8 OLI products. As opposed to the Landsat imagery that often fails to resolve even the largest glacial lakes of Tystigbreen, both PlanetScope and Sentinel-2 are helpful in identifying previously unrecognised glacial lakes and undocumented drainage events. Our analysis suggests that a fusion of the two satellite products may be beneficial for automated tracking of glacial lake changes. We also demonstrate that local knowledge and systematic involvement of citizens in data collection have a potential to enrich GLOF databases.

Extreme heat events in the summer of 2022 were observed in Eurasia, North America and China. Glaciers are a unique indicator of climate change, and the European Alps experienced substantial glacier mass loss as a result of the conditions in 2022, which prompted a wide range of community concerns. However, relevant findings for glaciers in China have not been currently reported. Here, we document the response of Urumqi Glacier No. 1 in the eastern Tien Shan to the extreme heat observed in 2022 based on in situ measurements that span more than 60 years. In 2022, Urumqi Glacier No. 1 exhibited the second largest annual mass loss on record, and the summer mass balance was the most negative on record. The hottest summer on record and relatively lower solid precipitation ratio contributed to the exceptional mass losses at Urumqi Glacier No. 1 in 2022, demonstrating the significant influence of heatwaves on extreme glacier melt in China.

We developed a multi-frequency, multi-Global Navigation Satellite System (GNSS) positioning instrument optimized for autonomous applications in the cryosphere. At lower power requirements and a fraction of the cost and weight compared to commercially available options, this instrument simplifies field usage and associated logistics. In this paper, we assess several baseline aspects of performance in a polar environment relative to geodetic receivers commonly used for glaciological applications. Evaluations of precision and relative accuracy of positioning show millimeter to centimeter-level (‘geodetic-grade’) quality of this instrument, making it a competitive alternative for GNSS glaciological and geophysical applications such as monitoring surface elevation change and ice flow. An array of these instruments, tested in the field on the Greenland Ice Sheet, also demonstrated robustness throughout the polar winter and met power and reliability requirements.

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.

Widespread glacier losses have been observed in most glaciated regions on Earth during recent decades, with a typical pattern of strong thinning in their lower reaches and limited elevation changes in their accumulation areas. Here, we use Pléiades satellite stereo-images of the Mont-Blanc massif (Alps) to reveal that thinning took place through the entire elevation range during the exceptional 2021/22 mass-balance year. Above 3000 m a.s.l. on Argentière glacier and Mer de Glace, thinning rates exceeded 3.5 m a−1 while almost no change occurred during the previous 9 years. Below 3000 m a.s.l., these anomalous thinning rates are essentially explained by changes in surface mass balance. At higher altitudes, other processes such as firn densification may play a role. Our analysis shows that high altitude glaciers, mostly stable during the last 100 years, are now responding to the impact of climate change.

The geodetic mass balance of a glacier corresponds to glacier-wide volume changes, converted to mass changes using density assumptions. It is typically calculated by differencing multi-temporal digital elevation models. In this study, we show how the annual geodetic mass balance of a glacier can be derived from uncrewed aerial vehicle (UAV) data. The presented workflow is applied to two small- to medium-sized glaciers in the Kyrgyz Tien Shan (Central Asia): Bordu glacier and Sary-Tor glacier. The obtained geodetic mass balance is compared with the glaciological mass balance derived from a network of ablation stakes and snow pits. A previously calibrated mass-balance model is used to correct for the difference in acquisition dates. The results show that the determined geodetic mass balance matches closely with the glaciological mass balance. Besides, for both glaciers the geodetic mass balance does not seem to be particularly sensitive to the assumptions regarding volume-to-mass conversion. Therefore, our results demonstrate that UAVs can serve as a valuable instrument to quantify the annual geodetic mass balance and to validate the glaciological mass balance. The conventional glaciological mass-balance estimation often relies on interpolation and extrapolation methods, whereas UAVs offer the potential for direct data acquisition over the entire glacier surface.

Glaciological mass-balance measurements have been the backbone of internationally coordinated glacier monitoring. The resulting annual observations have been used to understand glacier reactions to climate change, and to assess both regional and global glacier mass changes and related contributions to runoff and sea-level rise. However, the comparability of annual observations is hampered by differences in survey periods and mass-balance amplitudes between glaciers, regions and hemispheres. This study presents a simple approach to temporally downscale glaciological mass balance using seasonal observations and sine functions. The proposed analytical model allows reconstruction of the seasonal course of glacier mass balance at weekly to monthly resolution from only annual or seasonal observations. Strengths and limitations of this analytical approach are discussed and compared with results from numerical mass-balance modelling. Potential applications include seasonal corrections of glaciological and geodetic observations and comparisons to monthly results from spaceborne gravimetry and altimetry.

Aldegondabreen is a relatively small (5.3 km2) land-terminating glacier located in Nordeskiöld Land of Svalbard, ~10 km southwest of Barentsburg. Cumulative mass balance during 2002–20 equalled −21.79 m w.e., which corresponds to 37% of the total mass loss. The annual mass balance (Ba) varied from −0.24 to −2.19 m w.e., while the winter mass balance (Bw) ranged between 0.36 and 0.85 m w.e. Ba and Bw were strongly correlated with the positive degree-day sums and solid precipitation amounts, respectively, measured at the Barentsburg weather station. There was also a strong correlation (r = 0.76) between Ba and Bw, which shows that winter accumulation amplifies the consecutive summer ice melt by changing the surface albedo. The trend in both observational series is not detectable because the period from 2005 to 2013 was characterized by relatively high Bw and not very negative Ba values. This was also observed on the other Svalbard glaciers, and was related to prevailing north-westerly flows over Svalbard during the summer. Therefore, the decadal periodicity of the Aldegondabreen mass balance follows general archipelago patterns that are determined by regional-scale factors. Thus, the surface mass-balance time series, which is now the longest one in the central part of the Spitsbergen Island, is representative for the archipelago.

Close-range sensors are employed to observe glaciological processes that operate over short timescales (e.g. iceberg calving, glacial lake outburst floods, diurnal surface melting). However, under poor weather conditions optical instruments fail while the operation of radar systems below 17 GHz do not have sufficient angular resolution to map glacier surfaces in detail. This letter reviews the potential of millimetre-wave radar at 94 GHz to obtain high-resolution 3-D measurements of glaciers under most weather conditions. We discuss the theory of 94 GHz radar for glaciology studies, demonstrate its potential to map a glacier calving front and summarise future research priorities.

Climate change has had a significant impact on glacier recession, particularly in the Arctic, where glacier meltwater is an important contributor to global sea-level rise. Therefore, it is important to accurately quantify glacier recession within this sensitive region, using multiple observations of glacier extent. In this study, we mapped 480 glaciers in Novaya Zemlya, Russian Arctic, using object-based image analysis applied to multispectral Landsat satellite imagery in Google Earth Engine and quantify the area changes between 1986–89 and 2019–21. The results show that in 1986–89, the total glacierized area was 22 990 ± 301 km2, in 2000–01 the area was 22 525 ± 308 km2 and by 2019–21 the glacier area reduced to 21 670 ± 292 km2, representing a total of 5.8% reduction in glacier area between 1986–89 and 2019–21. Higher glacier area loss was observed on the Barents Sea coast (7.3%) compared to the Kara (4.2%), reflecting previously observed differences in warming trends. The accuracy of the automatically generated outlines of each layer (1986–89, 2000–01 and 2019–21) was evaluated by comparing with manually corrected outlines (reference data) using random sampling, resulting in an overall accuracy estimate of between 96 and 97% compared to the reference data. This automated approach in Google Earth Engine is a promising tool for rapidly mapping glacier change that reduces the amount of time required to generate accurate glacier outlines.

Seismic Rayleigh wave ellipticity measurements are the horizontal-to-vertical ratio of the Rayleigh wave particle motion, and are sensitive to the subsurface structure beneath a seismic station. H/V ratios measured from the ambient vibrations of the Earth are being increasingly used in glaciological applications to determine glacier and ice sheet thickness, seismic velocities and firn properties. Using the newly developed degree-of-polarisation (DOP-E) method which exploits the polarisation properties of seismic noise, we identify and extract Rayleigh waves from seismic stations in Greenland, and relate them to sea ice processes and the geology of the upper crust. Finally, we provide some suggestions for future applications of DOP-E method to gain greater insight into seasonal and long-term variability of sea ice formation and breakup as well as the monitoring of ice sheet thickness, subglacial environment and firn layers in the poles.

We use the apparent horizontal shift of an englacial tephra layer outcrop to calculate local glacier mass balance on Sléttjökull, a lobe of Mýrdalsjökull in Southern Iceland. For this approach, the dipping angle of the englacial tephra layer in the glacier upstream of the outcrop and the flow velocity of the ice need to be known. An earlier investigation was expanded by the application of ground-penetrating radar, detecting the depth of the tephra along tracks with a total length of 10 km. Interpolation between the tracks enables us to derive the dipping angle of the layer along several flow lines. Together with glacier surface velocities, determined from feature tracking, we are able to estimate the local surface mass balance from the horizontal displacement of the tephra outcrop using freely available satellite imagery without additional fieldwork. The earlier local balance series was extended to the period 2014/15 to 2019/20. Although the results for the individual profiles differ slightly from each other, they show the same temporal pattern and clear variations from year to year. The results are compared to traditional mass-balance data from Hofsjökull. The two series show a good agreement in their interannual variability.

We quantify the surface elevation changes along Yala Glacier in Langtang Valley, Nepal Himalaya, since 1981 using geodetic methods to understand the recent evolution and current state of small debris-free glaciers across the region. We analyse differential global positioning system measurements and aerial stereo imagery that were acquired along Yala Glacier in 2007, 2009, 2012 and 2015 to generate digital elevation models for each calculation period. Continuous surface lowering has mainly been observed across the down-glacier area during the calculation periods, although a large degree of variability exists, with this lowering trend propagating up-glacier in recent years. The area-weighted glacier mass balances range from −0.98 ± 0.27 to −0.26 ± 0.30 m w.e. a−1 for the five calculation periods (1981–2007, 2007–2009, 2009–2012, 2012–2015 and 2007–2015). These calculated mass-balance data reveal that Yala Glacier has undergone accelerated mass loss since the late 2000s, which is consistent with the results of previous in situ measurement and remote-sensing studies.

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.

Glacier change in Kenai Fjords National Park in southcentral Alaska affects local terrestrial, fresh water and marine ecosystems and will likely impact ecotourism. We used Landsat 4–8 imagery from 1984 through 2021 to manually map lower glacier ice margins for 19 maritime glaciers in Kenai Fjords National Park. Of these glaciers, six are tidewater, three are lake-terminating, six are land-terminating and four terminated in more than one environment throughout the study period. We used the mapped ice margins to quantify seasonal terminus position and areal change, including distinguishing between ice loss at glacier termini and along glacier margins. Overall, 13 glaciers substantially retreated (more than 2σ), 14 lost substantial area and only two underwent both net advance and area gain. The glaciers that had insubstantial length and area changes were predominantly tidewater. Cumulatively, the lower reaches of these 19 glaciers lost 42 km2 of ice, which was nearly evenly distributed between the terminus and the lateral margins. The rapid rate of glacier change and subsequent land cover changes are highly visible to visitors and locals at Kenai Fjords National Park, and this study quantifies those changes in terms of glacier length and area.

Glacier monitoring has been internationally coordinated for more than 125 years. Despite this long history, there is no authoritative answer to the popular question: ‘Which glaciers are the largest in the world?’ Here, we present the first systematic assessment of this question and identify the largest glaciers in the world – distinct from the two ice sheets in Greenland and Antarctica but including the glaciers on the Antarctic Peninsula. We identify the largest glaciers in two domains: on each of the seven geographical continents and in the 19 first-order glacier regions defined by the Global Terrestrial Network for Glaciers. Ranking glaciers by area is non-trivial. It depends on how a glacier is defined and mapped and also requires differentiating between a glacier and a glacier complex, i.e. glaciers that meet at ice divides such as ice caps and icefields. It also depends on the availability of a homogenized global glacier inventory. Using separate rankings for glaciers and glacier complexes, we find that the largest glacier complexes have areas on the order of tens of thousands of square kilometers whereas the largest glaciers are several thousands of square kilometers. The world's largest glaciers and glacier complexes are located in the Antarctic, Arctic and Patagonia.

Regional assessments of ice elevation change provide insight into the processes controlling an ice sheet's geometric response to climate forcing. In Southwest Greenland's land terminating sector (SWLTS), it is presumed that ice surface elevation changes result solely from changing surface mass balance (SMB). Here we test this assumption by developing a multi-decadal (1985–2017) record of elevation change from digital elevation models (DEMs) and comparing it to regional climate model output and available records of ice speed. The SWLTS thinned by >12 m on average over the full 32-year period, but the change was highly variable in time and space. Thinning was amplified in the central region of the SWLTS, relative to the north and south. During 1985–2007, the north and south regions demonstrated net thickening while the central region thinned. Regional differences in elevation change are inconsistent with SMB anomalies, indicating that enhanced ice flow in the north and south contributed to thickening during this early time interval. While clear validation in the south is prevented by incomplete velocity data, historical surface speeds in the north were elevated. These findings support the interpretation that changing ice flow can influence ice surface elevation in the slow-moving SWLTS.

Glacier motion responds dynamically to changing meltwater inputs, but the multi-decadal response of basal sliding to climate remains poorly constrained due to its sensitivity across multiple timescales. Observational records of glacier motion provide critical benchmarks to decode processes influencing glacier dynamics, but multi-decadal records that precede satellite observation and modern warming are rare. Here we present a record of motion in the ablation zone of Saskatchewan Glacier that spans seven decades. We combine in situ and remote-sensing observations to inform a first-order glacier flow model used to estimate the relative contributions of sliding and internal deformation on dynamics. We find a significant increase in basal sliding rates between melt-seasons in the 1950s and those in the 1990s and 2010s and explore three process-based explanations for this anomalous behavior: (i) the glacier surface steepened over seven decades, maintaining flow-driving stresses despite sustained thinning; (ii) the formation of a proglacial lake after 1955 may support elevated basal water pressures; and (iii) subglacial topography may cause dynamic responses specific to Saskatchewan Glacier. Although further constraints are necessary to ascertain which processes are of greatest importance for Saskatchewan Glacier's dynamic evolution, this record provides a benchmark for studies of multi-decadal glacier dynamics.