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.

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.

Ice-contact lakes modify glacier geometry and dynamics by shifting the majority of mass loss from the ice surface to the terminus. Lake-terminating glaciers are known to experience greater thinning rates and higher velocities than land-terminating glaciers, but the controls on variability in surface elevation change and ice flow between lake-terminating glaciers in different regions remain poorly explored. We combined existing datasets of glacier velocity, surface elevation change and glacial lake area to characterise the evolution of 352 lake-terminating and land-terminating glaciers within three Himalayan sub-regions between 2000 and 2019. These analyses show that the influence of ice-contact lakes propagates up-glacier across only the lowermost 30% of the hypsometric distribution, even where lakes are well established. We find that ice-contact lakes only affect glacier behaviour when the lakes reach an advanced evolutionary stage; most clearly manifested in the Eastern Himalaya by statistically robust differences in glacier-wide surface elevation change between lake-terminating (–0.68 ± 0.05 m a–1) and land-terminating (–0.54 ± 0.04 m a–1) glaciers. These differences are driven by the presence of a greater number of well-developed ice-contact lakes in the Eastern Himalaya compared to in the Western and Central Himalaya, resulting from greater mass loss rates to date.

Frontal ablation, the combination of submarine melting and iceberg calving, changes the geometry of a glacier's terminus, influencing glacier dynamics, the fate of upwelling plumes and the distribution of submarine meltwater input into the ocean. Directly observing frontal ablation and terminus morphology below the waterline is difficult, however, limiting our understanding of these coupled ice–ocean processes. To investigate the evolution of a tidewater glacier's submarine terminus, we combine 3-D multibeam point clouds of the subsurface ice face at LeConte Glacier, Alaska, with concurrent observations of environmental conditions during three field campaigns between 2016 and 2018. We observe terminus morphology that was predominately overcut (52% in August 2016, 63% in May 2017 and 74% in September 2018), accompanied by high multibeam sonar-derived melt rates (4.84 m d−1 in 2016, 1.13 m d−1 in 2017 and 1.85 m d−1 in 2018). We find that periods of high subglacial discharge lead to localized undercut discharge outlets, but adjacent to these outlets the terminus maintains significantly overcut geometry, with an ice ramp that protrudes 75 m into the fjord in 2017 and 125 m in 2018. Our data challenge the assumption that tidewater glacier termini are largely undercut during periods of high submarine melting.

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 dynamics of seasonal snow ablation on six glaciers in central Spitsbergen (Dicksonland) were assessed by examining a set of Sentinel-2 satellite images covering the summer ablation season for the period 2016–19. All glaciers lost 80% or more of their surface snow cover during the studied ablation seasons. This bolsters the recently observed trend of local glacier thinning, even at higher altitudes. Snow ablation dynamics are highly dependent on the glaciers altitudes, their position relative to the prevailing wind direction and the exposure to insolation. The accumulation areas of the studied glaciers were delimited based on the overlap of the minimum extent of snow-covered areas in the four consecutive studied summer seasons. The high temporal and spatial resolutions of available images enabled a detailed description of the seasonal snow ablation dynamics. Moreover, an estimate of the average number of days with below threshold glacier snow cover was made. This study contributes to our understanding of recent processes and might further support the modelling of glacier melt and subsequent runoff.

Marine-terminating glaciers on the Antarctic Peninsula (AP) have retreated, accelerated and thinned in response to climate change in recent decades. Ocean warming has been implicated as a trigger for these changes in glacier dynamics, yet little data exist near glacier termini to assess the role of ocean warming here. We use remotely-sensed iceberg melt rates seaward of two glaciers on the eastern and six glaciers on the western AP from 2013 to 2019 to explore connections between variations in ocean conditions and glacier frontal ablation. We find iceberg melt rates follow regional ocean temperature variations, with the highest melt rates (mean ≈ 10 cm d−1) at Cadman and Widdowson glaciers in the west and the lowest melt rates (mean ≈ 0.5 cm d−1) at Crane Glacier in the east. Near-coincident glacier frontal ablation rates from 2014 to 2018 vary from ~450 m a−1 at Edgeworth and Blanchard glaciers to ~3000 m a−1 at Seller Glacier, former Wordie Ice Shelf tributary. Variations in iceberg melt rates and glacier frontal ablation rates are significantly positively correlated around the AP (Spearman's ρ = 0.71, p-value = 0.003). We interpret this correlation as support for previous research suggesting submarine melting of glacier termini exerts control on glacier frontal dynamics around the AP.

We have developed a two-dimensional coupled glacier–fjord model, which runs automatically using Elmer/Ice and MITgcm software packages, to investigate the magnitude of submarine melting along a vertical glacier front and its potential influence on glacier calving and front position changes. We apply this model to simulate the Hansbreen glacier–Hansbukta proglacial–fjord system, Southwestern Svalbard, during the summer of 2010. The limited size of this system allows us to resolve some of the small-scale processes occurring at the ice–ocean interface in the fjord model, using a 0.5 s time step and a 1 m grid resolution near the glacier front. We use a rich set of field data spanning the period April–August 2010 to constrain, calibrate and validate the model. We adjust circulation patterns in the fjord by tuning subglacial discharge inputs that best match observed temperature while maintaining a compromise with observed salinity, suggesting a convectively driven circulation in Hansbukta. The results of our model simulations suggest that both submarine melting and crevasse hydrofracturing exert important controls on seasonal frontal ablation, with submarine melting alone not being sufficient for reproducing the observed patterns of seasonal retreat. Both submarine melt and calving rates accumulated along the entire simulation period are of the same order of magnitude, ~100 m. The model results also indicate that changes in submarine melting lag meltwater production by 4–5 weeks, which suggests that it may take up to a month for meltwater to traverse the englacial and subglacial drainage network.

Cryoconite holes are water-filled cylindrical holes formed on ablation ice surfaces and commonly observed on glaciers worldwide. Temporal changes of cryoconite holes characteristically <5 cm in diameter were monitored with a time-lapse interval camera over 15 d during the melting season on Qaanaaq Glacier in northwest Greenland. The holes drastically changed their dimensions and synchronously collapsed twice during the study period. When the holes collapsed, the coverage of cryoconite on the ice surface increased from 1.0 to 3.5% in the field of view of the camera, and then decreased again to 0.4% after the holes reformed. Comparison with meteorological data showed that the collapses occurred in cloudy and rainy or windy weather conditions, corresponding to low shortwave solar radiation (68–126 W m−2, 40–55% of the incoming flux). In contrast, holes developed in sunny conditions correspond to high solar radiation (186–278 W m−2, 63–88%). Results suggest that the dimensions of holes drastically changed depending on the weather conditions and that frequent cloudy, warm and windy conditions would cause a decay of holes and weathering crust, inducing an increase in the cryoconite coverage on the ice, consequently darkening the glacier surface.

Meltwater is important to understanding glacier health and dynamics. Since melt measurements are uncommon, ice ablation estimates are often based on models including the positive degree day (PDD) model. The PDD estimate is popular since it only requires air temperature as input, but suffers from the lack of physical motivation of an energy-balance model. We present a physics-based alternative to the PDD model that still only takes air/surface temperature as input. The model resembles the PDD model except accounting for time lags in ablation when cold ice needs to be warmed. The model is expressed as a differential equation with a single extra parameter related to the efficiency of heating a near-surface layer of ice. With zero thickness, the model reduces to the PDD model, providing a physical basis for the PDD model. Applying the model to data from Greenland, it improves modestly upon the PDD model, with the main improvement being better prediction of early season melting. This new model is a useful compromise, with some of the physics of more realistic models and the simplicity of a PDD model. The model should improve estimates of meltwater production and help constrain PDD parameters when empirical calibration is challenging.

Band ogives are a striking and enigmatic feature of Mer de Glace glacier flow. The surface mass balances (SMBs) of these ogives have been thoroughly investigated over a period of 12 years. We find similar cumulative SMBs over this period, ranging between −64.1 and −66.2 m w.e., on the dark and light ogives even though the dark ogive albedo is ~40% lower than that of the light ogives. We, therefore, looked for another process that could compensate for the large difference of absorbed short-wave radiation between dark and light ogives. Based on in situ roughness measurements, our numerical modeling experiments demonstrate that a significant difference in turbulent flux over the dark and light ogives due to different surface roughnesses could compensate for the difference in radiative forcing. Our results discard theories for the genesis of band ogives that are based on the assumption of a strong ice ablation contrast between dark and light ogives. More generally, our study demonstrates that future roughness changes are as important to analyze as the radiative impacts of a potential increase of aerosols or debris at the surface of glaciers.

We present the first dedicated study into the phenomenon of ice sails. These are clean ice structures that protrude from the surface of a small number of debris-covered glaciers and can grow to heights of over 25 m. We draw together what is known about them from the academic/exploration literature and then analyse imagery. We show here that ice sails can develop by one of two mechanisms, both of which require clean ice to become surrounded by debris-covered ice, where the debris layer is shallow enough for the ice beneath it to melt faster than the clean ice. Once formed, ice sails can persist for decades, in an apparently steady state, before debris layer thickening eventually causes a reversal in the relative melt rates and the ice sails decay to merge back with the surrounding glacier surface. We support our image-based analysis with a surface energy-balance model and show that it compares well with available observations from Baltoro Glacier in the Karakoram. A sensitivity analysis of the model is performed and confirms the results from our empirical study that ice sails require a relatively high evaporative heat flux and/or a relatively low sensible heat flux in order to exist.

Glacier albedo controls the surface energy budget, the variability of which affects the glacier surface melt rate and, in turn, impacts the mass balance of the glacier. During 2013 and 2014, spatial and temporal variations of albedo were investigated using 18 Landsat images of Urumqi Glacier No. 1. Factors influencing these spatiotemporal profiles were analyzed. An established retrieval process, including geolocation, radiometric calibration, atmospheric, topographic, and anisotropic correction and narrow- to broadband conversion, was applied for the first time to Landsat-8 images. Differences between Landsat image derived albedo values and albedo values measured using a handheld spectroradiometer ranged from −0.024 to 0.049. Spatial and temporal variations of surface albedo were significant, especially in the ablation area. The variability of the values of ice albedo ranged from 0.06 to 0.44 due to topographic effects and light-absorbing impurities. The results suggest that this retrieval method can be used to investigate the spatial and temporal variability of surface albedo from Landsat-8 images on mountain glaciers. Moreover, as constant albedo values for ice and snow cannot be assumed, the distribution of albedo was not completely dependent on altitude under conditions of more intense ablation, and by reason of light-absorbing impurities during the melt season.

This study provides an estimate of fresh water derived from ice melt for the ablation areas of glaciers in the Central Karakoram National Park (CKNP), Pakistan. In the CKNP there are ~700 glaciers, covering ~4600 km2, with widespread debris cover (518 km2). To assess meltwater volume we applied a distributed model able to describe both debris-covered and debris-free ice ablation. The model was calibrated using data collected in the field in the CKNP area and validated by comparison with ablation data collected in the field, independent of the data used in building the model. During 23 July–9 August 2011, the mean model-estimated ablation in the CKNP was 0.024 m w.e. d–1 in debris-covered areas and 0.037 m w.e. d–1 in debris-free areas. We found a mean error of +0.01 m w.e. (corresponding to 2%) and a root-mean-square error equal to 0.09 m w.e. (17%). According to our model, the ablation areas of all the glaciers in the CKNP produced a water volume of 1.963 km3 during the study period. Finally, we performed several sensitivity tests for assessing the impact of the input data variations.

Glacier surface melting can be described using energy-balance models. We conducted a surface energy budget experiment to quantify surface energy fluxes and to identify factors affecting glacial melt in the ablation zone of Laohugou glacier No. 12, western Qilian mountains. The surface energy budget was calculated based on data from an automatic weather station, and turbulent fluxes calculated using the bulk-aerodynamic approach were corrected using measurements from an eddy-covariance system. Simulated mass balances were validated by stake observations. Net shortwave radiation was the primary component of the surface energy balance (126Wm–2), followed by sensible heat flux. Net longwave radiation (–45Wm–2) and latent heat flux (–12.8 Wm–2) represented heat sinks. The bulk-aerodynamic method underestimated sensible and latent heat fluxes by 3.4 and 1.2 W m–2, respectively. The simulated total mass balance of –1703mmw.e. exceeded the observed total by 90 mm w.e. Daily positive accumulated temperature and low albedo were the main factors accelerating glacier melt. An uncertainty assessment showed that mass balance was very sensitive to albedo and varied by 36% when albedo changed by 0.1.

Most of the central Himalayan glaciers have surface debris layers of variable thickness, which greatly affect the ablation rate. An attempt has been made to relate debris-cover thickness to glacier surface melting. Thirty stakes were used to calculate ablation for debris-covered and clean ice of Dokriani Glacier (7 km2) from 2009/10 to 2012/13. Our study revealed significant altitude-wise difference in the rate of clean and debris-covered ice melting. We found a high correlation (R2 = 0.92) between mean annual clean-ice ablation and altitude, and a very low correlation (R2 = 0.14) between debris-covered ice melting and altitude. Debris-covered ice ablation varies with variation in debris thickness from 1 to 40 cm; ablation was maximum under debris thicknesses of 1–6 cm and minimum under 40 cm. Even a small debris-cover thickness (1–2 cm) reduces ice melting as compared to that of clean ice on an annual basis. Overall, debris-covered ice ablation during the study period was observed to be 37% less than clean-ice ablation. Strong downwasting was also observed in the Dokriani Glacier ablation area, with average annual ablation of 1.82 m w.e. a–1 in a similar period. Our study suggests that a thinning glacier rapidly becomes debris-covered over the ablation area, reducing the rate of ice loss.

The McMurdo Dry Valleys, southern Victoria Land, East Antarctica, are a polar desert, and melt from glacial ice is the primary source of water to streams, lakes and associated ecosystems. Previous work found that to adequately model glacier ablation and subsurface ice temperatures with a surface energy-balance model required including the transmission of solar radiation into the ice. Here we investigate the contribution of subsurface melt to the mass balance of (and runoff from) Dry Valley glaciers by including a drainage process in the model and applying the model to three glacier sites using 13 years of hourly meteorological data. Model results for the smooth glacier surfaces common to many glaciers in the Dry Valleys showed that sublimation was typically the largest component of surface lowering, with rare episodes of surface melting, consistent with anecdotal field observations. Results also showed extensive internal melting 5–15 cm below the ice surface, the drainage of which accounted for ~50% of summer ablation. This is consistent with field observations of subsurface streams and formation of a weathering crust. We identify an annual cycle of weathering crust formation in summer and its removal during the 10 months of winter sublimation.

Continuous surface debris cover strongly reduces the ablation of glaciers, but high melt rates may occur at ice cliffs that are too steep to hold debris. This study assesses the contribution of ice-cliff backwasting to total ablation of Miage glacier, Mont Blanc massif, Italy, in 2010 and 2011, based on field measurements, physical melt models and mapping of ice cliffs using a high-resolution (1 m) digital elevation model (DEM). Short-term model calculations closely match the measured melt rates. A model sensitivity analysis indicates that the effects of cliff slope and albedo are more important for ablation than enhanced longwave incidence from sun-warmed debris or reduced turbulent fluxes at sheltered cliff bases. Analysis of the DEM indicates that ice cliffs account for at most 1.3% of the 1 m pixels in the glacier’s debris-covered zone, but application of a distributed model indicates that ice cliffs account for ~7.4% of total ablation. We conclude that ice cliffs make an important contribution to the ablation of debris-covered glaciers, even when their spatial extent is very small.

Blue-ice areas (BIAs) cover ~1% of the East Antarctic ice sheet and are visual evidence of persistent ablation. In these regions, more snow is sublimated and/or eroded than is accumulated. The physical processes driving the formation of BIAs are poorly understood. Here we combine a firndensification model with high-resolution (5.5 km) maps of surface mass balance and ice velocity to simulate the build-up and removal of a firn layer along an ice flowline passing Byrd Glacier. A BIA is formed once the complete firn layer is removed. Feedback processes, which enhance blue-ice formation through the difference in surface characteristics of snow and ice, are examined using sensitivity simulations. The presence of blue ice on Byrd Glacier is found to be mainly determined by (1) ice velocity, (2) surface mass balance and (3) the characteristics (thickness, mass) of the firn layer prior to entering the ablation area. With a moderate decrease of the surface mass balance, the location and extent of the simulated BIA on Byrd Glacier is found to be in good qualitative agreement with MODIS optical imagery.

Ice cliffs have been identified as a reason for higher ablation rates on debris-covered glaciers than are implied by the insulation effects of the debris. This study aims to improve our understanding of cliff backwasting, and the role of radiative fluxes in particular. An energy-balance model is forced with new data gathered in May and October 2013 on Lirung Glacier, Nepalese Himalaya. Observations show substantial variability in melt between cliffs, between locations on any cliff and between seasons. Using a high-resolution digital elevation model we calculate longwave fluxes incident to the cliff from surrounding terrain and include the effect of local shading on shortwave radiation. This is an advance over previous studies, that made simplified assumptions on cliff geometry and radiative fluxes. Measured melt rates varied between 3.25 and 8.6 cm d−1 in May and 0.18 and 1.34 cm d−1 in October. Model results reproduce the strong variability in space and time, suggesting considerable differences in radiative fluxes over one cliff. In October the model fails to reproduce stake readings, probably due to the lack of a refreezing component. Disregarding local topography can lead to overestimation of melt at the point scale by up to ∼9%.