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.

The marginal areas of the Greenland ice sheet develop streams and lakes each summer, documenting that surface runoff of meltwater is a major component of ice-sheet mass balance. Here we map the slush limit, a proxy for the extent of surface runoff, using daily MODIS data for the years 2000–2021. We develop an automated algorithm capable of detecting daily slush limits, provided sufficient image quality. The algorithm is applied to the ice sheet's western flank (61.7 $^{\circ }$N to 76.5 $^{\circ }$N). We find significant increasing trends in maximum slush limits until the year 2012, but not thereafter. We show that the slush limit typically rises quickly early in the ablation season but stabilizes before melting ceases. The data provide evidence that upward migration of surface runoff in summer 2012 stopped early at the upper margin of the ice slabs. These thick and continuous ice layers are located close to the surface, in the firn, and impede percolation of melt into deeper pore space. Had the ice slabs extended higher, the summer 2012 provided sufficient energy to raise the slush limit by another $\sim$300 m in elevation.

Refreezing of meltwater in firn is a major component of Greenland ice-sheet's mass budget, but in situ observations are rare. Here, we compare the firn density and total ice layer thickness in the upper 15 m of 19 new and 27 previously published firn cores drilled at 15 locations in southwest Greenland (1850–2360 m a.s.l.) between 1989 and 2019. At all sites, ice layer thickness covaries with density over time and space. At the two sites with the earliest observations (1989 and 1998), bulk density increased by 15–18%, in the top 15 m over 28 and 21 years, respectively. However, following the extreme melt in 2012, elevation-detrended density using 30 cores from all sites decreased by 15 kg m−3 a−1 in the top 3.75 m between 2013 and 2019. In contrast, the lowest elevation site's density shows no trend. Thus, temporary build-up in firn pore space and meltwater infiltration capacity is possible despite the long-term increase in Greenland ice-sheet melting.

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.

Glaciers located in western High Mountain Asia (HMA) have shown mass gain or limited mass losses compared to other mountain regions since ~2000. Increases in accumulation may be responsible. Although no contemporary measurements exist to explore this hypothesis, extensive historic measurements including firn density, stratigraphy and accumulation rates at ~4400 m a.s.l. on Abramov Glacier, Pamir Alay, provide valuable indications of accumulation changes. Abramov Glacier is located at the northern margin of western HMA. In this study, we assess unique historical firn data of Abramov Glacier from the 1970s to evaluate past firn conditions in this data sparse region. The current firn state is investigated based on a series of in situ measurements including firn cores and ground-penetrating radar measurements performed in 2018. We compare the legacy data with contemporary firn characteristics. Our results indicate a high year-to-year variability, but generally increasing net accumulation during the last 60 years on Abramov Glacier, with 0.84 ± 0.22 m w.e. for 2011–18 compared to 0.68 ± 0.32 m w.e. for 1965–72 and 0.59 ± 0.22 m w.e. for 1970–97. These results from in situ data provide ground truth for the data-sparse western HMA.

We show a strong difference in surface mass and energy balance of a mountain glacier and two sites on the ice sheet at 64°N in West Greenland using stake and automated weather station observations. Net surface mass balance is on average 2.2 m w.e. less negative at the coast compared with the ice sheet in the same elevation. We find a larger energy turnover at the ice sheet margin on Qamanarssup Sermia than measured on the coastal mountain glacier Qassigiannguit with both energy input and output being of larger absolute value. More cloudiness and a thicker snow cover at the relatively humid coastal glacier result in smaller gains in net-shortwave radiation and smaller losses in net-longwave radiation and a less negative mass balance. Lower wind speeds at the coastal glacier result in weaker turbulent heat exchange between atmosphere and ice surface. On annual average, 17 W m−2 more energy is available for melt at the ice-sheet margin compared with the coastal glacier in the same elevation.

As a consequence of ongoing atmospheric temperature rise, tropical glaciers belong to the unique and threatened ecosystems on Earth, as defined by the Intergovernmental Panel on Climate Change (Houghton and others, 2001). Worldwide glacier monitoring, especially as part of the Global Climate Observing System (GCOS), includes the systematic collection of data on such perennial surface ice masses. Several peaks in the sierras of Colombia have lost their glacier cover during recent decades. Today, high-altitude glaciers still exist in Sierra Nevada de Santa Marta, in Sierra Nevada del Cocuy and on the volcanoes of Nevados del Ruiz, de Santa Isabel, del Tolima and del Huila. Comparison of reconstructions of maximum glacier area extent during the Little Ice Age with more recent information from aerial photographs and satellite images clearly documents a fast-shrinking tendency and potential disappearance of the remaining glaciers within the next few decades. In the past 50 years, Colombian glaciers have lost 5 0% or more of their area. Glacier shrinkage has continued to be strong in the last 15 years, with a loss of 10−50% of the glacier area. The relationship between fast glacier retreat and local, regional and global climate change is now being investigated. Preliminary analyses indicate that the temperature rise of roughly 1°C in the last 30 years recorded at high-altitude meteorological stations exerts a primary control on glacier retreat. The investigations on the Colombian glaciers thus corroborate earlier findings concerning the high sensitivity of glaciers in the wet inner tropics to temperature rise. To improve understanding of fast glacier retreat in Colombia, a modern monitoring network has been established according to the multilevel strategy of the Global Terrestrial Network for Glaciers (GTN-G) within GCOS. The observations are also contributions to continued assessments of hazards from the glacier-covered volcanoes and to integrated global change research in mountain biosphere reserves.

Modern concepts of worldwide glacier monitoring include numerical models for (1) interconnecting the different levels of observations (local mass balance, representative length change, glacier inventories for global coverage) and (2) extrapolations in space (coupling with climate models) and time (backward and forward). In this context, one important new tool is distributed mass-balance modelling in complex mountain topography. This approach builds on simplified energy-balance models and can be applied for investigating the spatio-temporal representativity of the few mass-balance measurements, for estimating balance values at the tongue of unmeasured glaciers in order to derive long-term average balance values from a great number of glaciers with known length change, and for assessing special effects such as the influence of Sahara dust falls on the albedo and mass balance or autocorrelation effects due to surface darkening of glaciers with strongly negative balances. Experience from first model runs in the Swiss Alps and from applications to the extreme conditions in summer 2003 provides evidence about the usefulness of this approach for glacier monitoring and analysis of glacier changes in high-mountain regions. The main difficulties concern the spatial variability of the input parameters (e.g. precipitation, snow cover and surface albedo) and the uncertainties in the parameterizations of the components of the energy balance. Field measurements remain essential to tie the models to real ground conditions.

The direct comparison of modeled glacier mass-balance distribution with field measurements could be problematic, as the methodology of determination and the processes considered at the point or catchment scale could differ strongly. Moreover, direct measurements cover only small parts of a glacier, and model performance is thus difficult to assess outside these regions. Remaining opportunities for model validation include comparison of snowlines, as derived from remote-sensing data, and maps of the mass-balance distribution, as interpreted by observers with local knowledge. This study compares such hand-drawn maps of the Vernagtferner (Oetzal Alps) mass balance with the modeled pattern as obtained from a distributed energy-/mass-balance model of intermediate complexity. The model is driven by measured daily values of temperature, global radiation and precipitation from a nearby climate station and grids of mean daily potential global radiation and climatologic annual precipitation sums. Compared with the direct measurements, the calculated mean mass balance and equilibrium-line altitude agreed very well in both balance years, although the spatial pattern of the mass-balance distribution displayed larger deviations in regions influenced by snowdrift or local topographic undulations. However, compared with the remaining snow as visible on a satellite image from 1999, the modeled snowline pattern for a larger sample of glaciers is in very good agreement.

Glacier mass-balance parameters such as the equilibrium-line altitude (ELA) play an important role when working with large glacier samples. While the number of observational mass-balance series to derive such parameters is limited, more and more modeled data are becoming available.Here we explore the possibilities of analyzing such 'synthetic' mass-balance data with respect to mass-balance parameters. A simplified energy-balance model is driven by bias-corrected regional climate model output to model mass-balance distributions for 94 glaciers in the Swiss Alps over 15 years. The modeling results in realistic interannual variability and mean cumulative mass balance. Subsequently model output is analyzed with respect to 18 topographic and mass-balance parameters and a correlation analysis is performed. Well-known correlations such as for ELA and median elevation are confirmed from the synthetic data. Furthermore, previously unreported parameter relationships are found such as a correlation of the balance rate at the tongue with the accumulation-area ratio (AAR) and of the glacier elevation range with the AAR. Analyzing modeled data complements in situ observations and highlights their importance: the small number of accurate mass-balance observations available for validation is a major challenge for the presented approach.

Ongoing atmospheric warming causes rapid shrinking of glaciers in the European Alps, with a high chance of their near-complete disappearance by the end of the 21st century. Here we present a comparison of three independent approaches to model the possible evolution of the glaciers in the Swiss Alps over the 21st century. The models have different levels of complexity, work at a regional scale and are forced with three scenarios of temperature increase (low, moderate, high). The moderate climate scenario gives an increase in air temperature of ∼2°C and ∼4°C for the two scenario periods 2021-50 and 2070-99, respectively, resulting in an area loss of 60-80% by 2100. In reality, the shrinkage could be even faster, as the observed mean annual thickness loss is already stronger than the modelled one. The three approaches lead to rather similar results with respect to the overall long-term evolution. The choice of climate scenarios produces the largest spread (∼40%) in the final area loss, while the uncertainty in present-day ice-thickness estimation causes about half this spread.

Spatial heterogeneity of snow and firn properties on glaciers introduces uncertainty in interpretation of point and profile observations and complicates modelling of meltwater percolation and runoff. Here we present a study of the temporal and spatial dynamics of firn density and stratigraphy at the plot-scale (≈10 m × 10 m × 10 m) repeated annually during 2012–14 at the Lomonosovfonna icefield, Svalbard. Results from cores, video inspections in boreholes and radar grid surveys are compared. Ice layers 0.1–50 cm thick comprised ≈8% of the borehole length. Most of them are 1–3 cm thick and could not be traced between boreholes separated by 3 m. Large lateral variability of firn structure affects representativeness of observations in single holes and calls for repeated studies in multiple points to derive a representative stratigraphy signal. Radar reflections are poorly correlated with ice layers in individual boreholes. However, the match between the high amplitude peaks in the grid-averaged radar signal and horizons of preferential ice layer formation revealed by averaging the video surveys over multiple boreholes is higher. These horizons are interpreted as buried firn layers previously exposed to melt-freeze or wind-driven densification and several of them are consistently recovered throughout three field campaigns.

Glacier surface mass-balance measurements on Greenland started more than a century ago, but no compilation exists of the observations from the ablation area of the ice sheet and local glaciers. Such data could be used in the evaluation of modelled surface mass balance, or to document changes in glacier melt independently from model output. Here, we present a comprehensive database of Greenland glacier surface mass-balance observations from the ablation area of the ice sheet and local glaciers. The database spans the 123 a from 1892 to 2015, contains a total of ~3000 measurements from 46 sites, and is openly accessible through the PROMICE web portal (http://www.promice.dk). For each measurement we provide X, Y and Z coordinates, starting and ending dates as well as quality flags. We give sources for each entry and for all metadata. Two thirds of the data were collected from grey literature and unpublished archive documents. Roughly 60% of the measurements were performed by the Geological Survey of Denmark and Greenland (GEUS, previously GGU). The data cover all regions of Greenland except for the southernmost part of the east coast, but also emphasize the importance of long-term time series of which there are only two exceeding 20 a. We use the data to analyse uncertainties in point measurements of surface mass balance, as well as to estimate surface mass-balance profiles for most regions of Greenland.

We present in situ firn temperatures from the extreme 2012 melt season in the southwestern lower accumulation area of the Greenland ice sheet. The upper 2.5 m of snow and firn was temperate during the melt season, when vertical meltwater percolation was inefficient due to a ~5.5 m thick ice layer underlying the temperate firn. Meltwater percolation and refreezing beneath 2.5 m depth only occurred after the melt season. Deviations from temperatures predicted by pure conductivity suggest that meltwater refroze in discrete bands at depths of 2.0–2.5, 5.0–6.0 and 8.0–9.0 m. While we find no indication of meltwater percolation below 9 m depth or complete filling of pore volume above, firn at 10 and 15 m depth was respectively 4.2–4.5°C and 1.7°C higher than in a conductivity-only simulation. Even though meltwater percolation in 2012 was inefficient, firn between 2 and 15 m depth the following winter was on average 4.7°C warmer due to meltwater refreezing. Our observations also suggest that the 2012 firn conditions were preconditioned by two warm summers and ice layer formation in 2010 and 2011. Overall, firn temperatures during the years 2009–13 increased by 0.6°C.