Digital elevation models (DEMs) from the spaceborne interferometric radar mission TanDEM-X hold a large potential for glacier change assessments. However, a bias is potentially introduced through the penetration of the X-band signal into snow and firn. To improve our understanding of radar penetration on glaciers, we compare DEMs derived from the almost synchronous acquisition of TanDEM-X and Pléiades optical stereo-images of Grosser Aletschgletscher in March 2021. We found that the elevation bias – averaged per elevation bin – can reach up to 4–8 m in the accumulation area, depending on post co-registration corrections. Concurrent in situ measurements (ground-penetrating radar, snow cores, snow pits) reveal that the signal is not obstructed by the last summer horizon but reaches into perennial firn. Because of volume scattering, the TanDEM-X surface is determined by the scattering phase centre and does not coincide with a specific firn layer. We show that the bias corresponds to more than half of the decadal ice loss rate. To minimize the radar penetration bias, we recommend to select DEMs from the same time of the year and over long observation periods. A correction of the radar penetration bias is recommended, especially when combining optical and TanDEM-X DEMs.

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.

Measuring surface mass-balance in the accumulation areas of glaciers is challenging because of the high spatial variability of snow accumulation and the difficulty of conducting annual field glaciological measurements. Here, we propose a method that can solve both these problems for many locations. Ground-penetrating radar measurements and firn cores extracted from a site in the French Alps were first used to reconstruct the topography of a buried end-of-summer snow horizon from a past year. Using these data and surface elevation observations from LiDAR and Global Navigation Satellite System instruments, we calculated the submergence velocities over the period between the buried horizon and more recent surface elevation observations. The differences between the changes in surface elevation and the submergence velocities were then used to calculate the annual surface mass-balances with an accuracy of ±0.34 m w.e. Assuming that the submergence velocities remain stable over several years, the surface mass-balance can be reconstructed for subsequent years from the differences in surface elevation alone. As opposed to the glaciological method that requires substantial fieldwork year after year to provide only point observations, this method, once submergence velocities have been calculated, requires only remote-sensing data to provide spatially distributed annual mass-balances in accumulation areas.

Theory and observation show that glacier-flow regimes characterized by high basal slip enhance the projection of topographic detail to the surface, motivating this investigation into the efficacy of using glacier surges to improve bed estimation. Here we adapt a Bayesian inversion scheme and apply it to real and synthetic data as a proof of concept. Synthetic tests show a reduction in mean RMSE between true and inferred beds by more than half, and an increase in the mean correlation coefficient of ~0.5, when data from slip- versus deformation-dominated regimes are used. Multi-epoch inversions, which partition slip- and deformation-dominated regimes, are shown to outperform inversions that average over these flow regimes thereby squandering information. Tests with real data from a surging glacier in Yukon, Canada, corroborate these results, while highlighting the challenges of limited or inconsistent data. With the growing torrent of satellite-based observations, fast-flow events such as glacier surges offer potential to improve bed estimation for some of the world's most dynamic glaciers.

We present geodetic mass-balance estimates for ten glaciers (22.6 ± 1.1 km2) around Volcán Domuyo between 1962 and 2020 (and 46 glaciers covering 29 ± 1.5 km2 between 1984 and 2020), derived from airborne, ASTER and Pléiades imagery. Overall, we find a slightly negative mass balance (−0.15 ± 0.09 m w.e. a–1) for the entire 1962–2020 time span. A closer inspection of sub-periods reveals, however, an increasingly negative mass balance with time. The Domuyo glaciers shifted from a moderately positive mass balance of +0.28 ± 0.13 m w.e. a–1 between 1962 and 1984, to a strongly negative mass balance of −0.99 ± 0.19 m w.e. a–1 between 2012 and 2020. An increase in summer temperatures and a decrease in winter precipitation during the last four decades are likely drivers of the observed glacier changes. We support this finding by implementing a minimal glacier model, relying solely on monthly precipitation and air temperatures. The mass-balance evolution detected in Volcán Domuyo is consistent with other sites spread across the Central Andes, suggesting rapidly increasing glacier wastage impacts are occurring at a geographically wider scale.

The Desert Andes contain >4500 ice masses, but only a handful are currently being monitored. We present the mass changes of the small mountain glacier Agua Negra (1 km2) and of the rest of glaciers in the Jáchal river basin. Remote-sensing data show Agua Negra glacier lost 23% of its area during 1959–2019. Glaciological measurements during 2014–2021 indicate an average annual mass balance of −0.52 m w.e. a−1, with mean winter and summer balances of 0.80 and −1.33 m w.e. a−1, respectively. The Equilibrium Line Altitude (ELA) is estimated to be 5100 ± 100 m a.s.l., which corresponds to an Accumulation Area Ratio (AAR) of 0.28 ± 0.21. Geodetic data from SRTM X and Pléiades show a doubling of the loss rate from −0.32 ± 0.03 m w.e. a−1 in 2000–2013, to −0.66 ± 0.06 m w.e. a−1 in 2013–2019. Comparatively, the ice losses for the entire Jáchal river basin (25 500 km2) derived from ASTER show less negative values, −0.11 ± 16 m w.e. a−1 for 2000–2012 and −0.23 ± 14 m w.e. a−1 for 2012–2018. The regional warming trend since 1979 and a recent decline in snow accumulation are probably driving the observed glacier mass balance.

The Cook Ice Cap (CIC) on the sub-Antarctic Kerguelen Islands recently experienced extremely negative surface mass balance. Further deglaciation could have important impacts on endemic and invasive fauna and flora. To put this exceptional glacier evolution into a multi-centennial-scale context, we refined the evolution of the CIC over the last millennium, investigated the associated climate conditions and explored its potential evolution by 2100 ce. A glaciological model, constrained by cosmic ray exposure dating of moraines, historical documents and recent direct mass balance observations, was used to simulate the ice-cap extents during different phases of advance and retreat between the last millennium and 2100 ce. Cosmogenic dating suggests glacial advance around the early Little Ice Age (LIA), consistent with findings from other sub-Antarctic studies, and the rather cold and humid conditions brought about by the negative phase of the Southern Annular Mode (SAM). This study contributes to our currently limited understanding of palaeoclimate for the early LIA in the southern Indian Ocean. Glaciological modelling and observations confirm the recent decrease in CIC extent linked to the intensification of the SAM. Although affected by large uncertainties, future simulations suggest a complete disappearance of CIC by the end of the century.

Multi-decadal mass loss estimates are available for few glaciers of Central Asia. On Abramov Glacier (Pamir-Alay, Kyrgyzstan), comprehensive long-term glaciological measurements have been carried out from 1968 to 1999 and re-initiated in 2011. A climatological interpretation of this benchmark glacier in Central Asia requires bridging the gap between historical and renewed measurements. This is achieved here by computing the geodetic mass balance from 1975 to 2015 using previously unreleased Soviet aerial imagery and Pléaides stereo-imagery. During 1975–2015, Abramov Glacier lost 2.2 km2 (8.2%) of its area. The mean annual thickness change was − 0.43 ± 0.14 m a−1 for the period 1975–2015, corresponding to a volume change of − 0.45 ± 0.15 km3. The average specific geodetic mass balance amounts to − 0.38 ± 0.12 m w.e. a−1. The 1975–2015 glacier mass loss lies within the range of glaciological and geodetic mass-balance estimates that were previously published for disparate and shorter time intervals since 1968. This study covers a much longer time period than earlier geodetic estimates and demonstrates the capacity to geodetically constrain glacier change at high spatial resolution in Central Asia using historic aerial imagery and Structure from Motion techniques. Therefore, it could serve as a benchmark for future studies of regional mass change.

The Kaskawulsh Glacier is an iconic outlet draining the icefields of the St. Elias Mountains in Yukon, Canada. We determine and attempt to interpret its catchment-wide mass budget since 2007. Using SPOT5/6/7 data we estimate a 2007–18 geodetic balance of −0.46 ± 0.17 m w.e. a−1. We then compute balance fluxes and observed ice fluxes at nine flux gates to examine the discrepancy between the climatic mass balance and internal mass redistribution by glacier flow. Balance fluxes are computed using a fully distributed mass-balance model driven by downscaled and bias-corrected climate-reanalysis data. Observed fluxes are calculated using NASA ITS_LIVE surface velocities and glacier cross-sectional areas derived from ice-penetrating radar data. We find the glacier is still in the early stages of dynamic adjustment to its mass imbalance. We estimate a committed terminus retreat of ~23 km under the 2007–18 climate and a lower bound of 46 km3 of committed ice loss, equivalent to ~15% of the total glacier volume.

The 2007–19 glaciological mass-balance series of Mera Glacier in the Everest Region, East Nepal, is reanalysed using the geodetic mass balance assessed by differencing two DEMs obtained from Pléiades stereo-images acquired in November 2012 and in October 2018. The glaciological glacier-wide annual mass balance of Mera Glacier has to be systematically decreased by 0.11 m w.e. a−1 to match the geodetic mass balance. We attribute part of the positive bias of the glaciological mass balance to an over-estimation of the accumulation above 5520 m a.s.l., likely due to a measurement network unable to capture its spatial variability. Over the period 2007–19, Mera Glacier has lost mass at a rate of −0.41 ± 0.20 m w.e. a−1, in general agreement with regional averages for the central Himalaya. We observe a succession of negative mass-balance years since 2013.

Karakoram glaciers experienced balanced or slightly positive mass budgets since at least the 1970s. Here, we provide an update on the state of balance of Central and Eastern Karakoram glaciers (12 000 km2) between 2008 and 2016 by differencing DEMs derived from satellite optical images. The mass budget of Central Karakoram glaciers was slightly positive (0.12 ± 0.14 m w.e. a−1) while eastern Karakoram glaciers lost mass (−0.24 ± 0.12 m w.e. a−1). The glacier-wide mass balances of surge-type and nonsurge-type glaciers were not statistically different. Our elevation change data also depict the effect of a > 100 Mm3 rock avalanche on Siachen Glacier ablation area in September 2010. It covered a 4 km2 area with a thick debris layer that unexpectedly, led to locally enhanced glacier mass loss during the following years. Enhanced melt opened a > 100 m deep 2 km2 depression and contributed to 6% of the mass loss of Siachen Glacier from 2010 to 2016 (−0.39 m w.e. a−1). We hypothesize that sub- or englacial melt may be responsible for this intriguing behaviour. This study contributes to a better knowledge of the regional pattern of the Karakoram anomaly and of the influence of rock avalanches on glacier mass changes.

Mass-balance measurements of Icelandic glaciers are sparse through the 20th century. However, the large archive of stereo images available allows estimates of glacier-wide mass balance ($\dot{B}$) in decadal time steps since 1945. Combined with climate records, they provide further insight into glacier–climate relationship. This study presents a workflow to process aerial photographs (1945–1995), spy satellite imagery (1977–1980) and modern satellite stereo images (since 2000) using photogrammetric techniques and robust statistics in a highly automated, open-source pipeline to retrieve seasonally corrected, decadal glacier-wide geodetic mass balances. In our test area, Eyjafjallajökull (S-Iceland, ~70 km2), we obtain a mass balance of $<![CDATA[ $ \dot{\curr B}_{\curr 1945}^{\curr 2014} \curr = -0.27 \pm 0.03\,{\rm \curr m\ w}{\rm. \curr e}{\rm.} {\rm \curr a}^{{\rm \ndash \curr 1}}$, with a maximum and minimum of $\dot{\curr B}_{\curr 1984}^{\curr 1989} \curr = 0.77 \curr \pm 0.19\,{\rm \curr m\ \curr w}{\rm\curr . e}{\rm\curr .} {\rm\curr a}^{{\rm\curr \ndash 1}}$ and $\dot{\curr B}_{\curr 1994}^{\curr 1998}\curr = -1.94 \curr \pm 0.34\,{\rm \curr m\ w}{\rm\curr . e}{\rm\curr .} {\rm \curr a}^{{\rm\curr \ndash 1}}$, respectively, attributed to climatic forcing, and $\dot{\curr B}_{\curr 2009}^{\curr 2010} \curr = -3.39{\rm \;} \curr \pm {\rm \;} \curr 0.43\,{\rm \curr m\ w}{\rm\curr . e}{\rm\curr .} {\rm\curr a}^{{\rm\curr \ndash 1}}$, mostly caused by the April 2010 eruption. The reference-surface mass balances correlate with summer temperature and winter precipitation, and linear regression accounts for 80% of the mass-balance variability, yielding a static sensitivity of mass balance to summer temperature and winter precipitation of − 2.1 ± 0.4 m w.e.a–1K–1 and 0.5 ± 0.3 m w.e.a–1 (10%)–1, respectively. This study serves as a template that can be used to estimate the mass-balance changes and glaciers' response to climate.

In this study, we combine remote sensing, in situ and model-derived datasets from 1966 to 2014 to calculate the mass-balance components of Kronebreen, a fast-flowing tidewater glacier in Svalbard. For the well-surveyed period 2009–2014, we are able to close the glacier mass budget within the prescribed errors. During these 5 years, the glacier geodetic mass balance was −0.69 ± 0.12 m w.e. a−1, while the mass budget method led to a total mass balance of −0.92 ± 0.16 m w.e. a−1, as a consequence of a strong frontal ablation (−0.78 ± 0.11 m w.e. a−1), and a slightly negative climatic mass balance (−0.14 ± 0.11 m w.e. a−1). The trend towards more negative climatic mass balance between 1966–1990 (+0.20 ± 0.05 m w.e. a−1) and 2009–2014 is not reflected in the geodetic mass balance trend. Therefore, we suspect a reduction in ice-discharge in the most recent period. Yet, these multidecadal changes in ice-discharge cannot be measured from the available observations and thus are only estimated with relatively large errors as a residual of the mass continuity equation. Our study presents the multidecadal evolution of the dynamics and mass balance of a tidewater glacier and illustrates the errors introduced by inferring one unmeasured mass-balance component from the others.

We simulate the ice dynamics of the San Rafael Glacier (SRG) in the Northern Patagonia Icefield (46.7°S, 73.5°W), using glacier geometry obtained by airborne gravity measurements. The full-Stokes ice flow model (Elmer/Ice) is initialized using an inverse method to infer the basal friction coefficient from a satellite-derived surface velocity mosaic. The high surface velocities (7.6 km a−1) near the glacier front are explained by low basal shear stresses (<25 kPa). The modelling results suggest that 98% of the surface velocities are due to basal sliding in the fast-flowing glacier tongue (>1 km a−1). We force the model using different surface mass-balance scenarios taken or adapted from previous studies and geodetic elevation changes between 2000 and 2012. Our results suggest that previous estimates of average surface mass balance over the entire glacier (Ḃ) were likely too high, mainly due to an overestimation in the accumulation area. We propose that most of SRG imbalance is due to the large ice discharge (−0.83 ± 0.08 Gt a−1) and a slightly positive Ḃ (0.08 ± 0.06 Gt a−1). The committed mass-loss estimate over the next century is −0.34 ± 0.03 Gt a−1. This study demonstrates that surface mass-balance estimates and glacier wastage projections can be improved using a physically based ice flow model.

We present a comprehensive review of the status and changes in glacier length (since the 1850s), area and mass (since the 1960s) along the Himalayan-Karakoram (HK) region and their climate-change context. A quantitative reliability classification of the field-based mass-balance series is developed. Glaciological mass balances agree better with remotely sensed balances when we make an objective, systematic exclusion of likely flawed mass-balance series. The Himalayan mean glaciological mass budget was similar to the global average until 2000, and likely less negative after 2000. Mass wastage in the Himalaya resulted in increasing debris cover, the growth of glacial lakes and possibly decreasing ice velocities. Geodetic measurements indicate nearly balanced mass budgets for Karakoram glaciers since the 1970s, consistent with the unchanged extent of supraglacial debris-cover. Himalayan glaciers seem to be sensitive to precipitation partly through the albedo feedback on the short-wave radiation balance. Melt contributions from HK glaciers should increase until 2050 and then decrease, though a wide range of present-day area and volume estimates propagates large uncertainties in the future runoff. This review reflects an increasing understanding of HK glaciers and highlights the remaining challenges.

Ice surface altimetry from ICESat-1 and NASA aircraft altimeter overflights spanning 2002–09 indicate that a region of lower Crane Glacier, Antarctic Peninsula, shows an unusual temporal pattern of elevation loss: a period of very rapid drawdown (~91ma–1 between September 2004 and September 2005) bounded by periods of large but more moderate rates (23ma–1 until September 2004; 12ma–1 after September 2005). The region of increased drawdown is ~4.5 km ×2.2 km based on satellite (ASTER and SPOT-5) stereo-image digital elevation model (DEM) differencing spanning the event. In a later differential DEM the anomalous drawdown feature is not seen. Bathymetry in Crane Glacier fjord reveals a series of flat-lying, formerly subglacial deeps interpreted as lake sediment basins. We conclude that the elevation-change feature resulted from drainage of a small, deep subglacial lake. We infer that the drainage event was induced by hydraulic forcing of subglacial water past a downstream obstruction. However, only a fraction of Crane Glacier’s increase in flow speed that occurred near the time of lake drainage (derived from image feature tracking) appears to be directly attributable to the event; instead, retreat of the ice front off a subglacial ridge 6 km downstream of the lake is likely the dominant cause of renewed fast flow and more negative mass balance in the subsequent 4 years.

Automatic feature tracking on two Landsat images (acquired inJanuary 2000 and December 2001) generates a complete and accurate velocity field of Mertz Glacier, East Antarctica. This velocity field shows two main tributaries to the ice stream. Between the tributaries, a likely obstruction feature in the bedrock results in a slow-down of the flow. A third Landsat image, acquired in 1989 and combined with the 2000 image, permits the determination of the glacier mean velocity during the 1990s. Although some parts of the Mertz Glacier system show evidence of slight speed increase, we conclude that the Mertz flow speed is constant within our uncertainty (35 m a−1). Using this complete velocity field, new estimates of the ice discharge flux, 17.8 km3 a−1 (16.4 Gt a−1), and of the basal melting of the tongue, 11 m a−1 of ice, are given. Our results lead to an apparent imbalance of the drainage basin (ice discharge 3.5 km3 a−1 lower than the accumulation). Considering previous studies in the Mertz Glacier area, we then discuss the uncertainty of this imbalance and the problems with accumulation mapping for this region.

By subtracting surface topographies from 1979, 1994, 2000 and 2008, we measured icethinning rates increasing from 1 ma-1 (1979-94) to >4 ma-1 (2000-08) on the tongue of Mer de Glace, French Alps. The relative contributions of changes in surface mass balance and ice fluxes to this acceleration in the thinning are estimated using field and remote-sensing measurements. Between 1979-94 and 2000-08, surface mass balance diminished by 1.2mw.e.a-1, mainly because of atmospheric warming. Mass-balance changes induced by the growing debris-covered area and the evolving glacier hypsometry compensated each other. Meanwhile, Mer de Glace slowed down and the ice fluxes through two cross sections at 2200 and 2050ma.s.l. decreased by 60%. Between 1979-94 and 2000-08, two-thirds of the increase in the thinning rates was caused by reduced ice fluxes and one- third by rising surface ablation. However, these numbers need to be interpreted cautiously given our inability to respect mass conservation below our upper cross section. An important implication is that large errors would occur if one term of the continuity equation (e.g. surface mass balance) were deduced from the two others (e.g. elevation and ice-flux changes).

Little is known about the Himalayan glaciers, although they are of particular interest in terms of future water supply, regional climate change and sea-level rise. In 2002, a long-term monitoring programme was started on Chhota Shigri Glacier (32.2° N, 77.5° E; 15.7 km2, 6263–4050 ma.s.l., 9 km long) located in Lahaul and Spiti Valley, Himachal Pradesh, India. This glacier lies in the monsoon–arid transition zone (western Himalaya) which is alternately influenced by Asian monsoon in summer and the mid-latitude westerlies in winter. Here we present the results of a 4 year study of mass balance and surface velocity. Overall specific mass balances are mostly negative during the study period and vary from a minimum value of –1.4 m w.e. in 2002/03 and 2005/06 (equilibrium-line altitude (ELA) ∼5180 m a.s.l.) to a maximum value of +0.1 m w.e. in 2004/05 (ELA 4855 m a.s.l.). Chhota Shigri Glacier seems similar to mid-latitude glaciers, with an ablation season limited to the summer months and a mean vertical gradient of mass balance in the ablation zone (debris-free part) of 0.7mw.e.(100 m)–1, similar to those reported in the Alps. Mass balance is strongly dependent on debris cover, exposure and the shading effect of surrounding steep slopes.

Mass-balance and dynamic behaviour of Chhota Shigri glacier, western Himalaya, India, has been investigated between 2002 and 2010 and compared to data collected in 1987-89. During the period 2002-10, the glacier experienced a negative glacier-wide mass balance of -0.67 ± 0.40 m w.e. a-1. Between 2003 and 2010, elevation and ice-flow velocities slowly decreased in the ablation area, leading to a 24-37% reduction in ice fluxes, an expected response of the glacier dynamics to its recent negative mass balances. The reduced ice fluxes are still far larger than the balance fluxes calculated from the 2002-10 average surface mass balances. Therefore, further slowdown, thinning and terminus retreat of Chhota Shigri glacier are expected over the next few years. Conversely, the 2003/04 ice fluxes are in good agreement with ice fluxes calculated assuming that the glacier-wide mass balance is zero. Given the limited velocity change between 1987−89 and 2003/04 and the small terminus change between 1988 and 2010, we suggest that the glacier has experienced a period of near-zero or slightly positive mass balance in the 1990s, before shifting to a strong imbalance in the 21st century. This result challenges the generally accepted idea that glaciers in the Western Himalaya have been shrinking rapidly for the last few decades.