Despite previous studies, glacier–lake interactions and future lake development in the Poiqu River basin, central Himalaya, are still not well understood. We mapped glacial lakes, glaciers, their frontal positions and ice flow from optical remote sensing data, and calculated glacier surface elevation change from digital terrain models. During 1964–2017, the total glacial-lake area increased by ~110%. Glaciers retreated with an average rate of ~1.4 km2 a−1 between 1975 and 2015. Based on rapid area expansion (>150%), and information from previous studies, eight lakes were considered to be potentially dangerous glacial lakes. Corresponding lake-terminating glaciers showed an overall retreat of 6.0 ± 1.4 to 26.6 ± 1.1 m a−1 and accompanying lake expansion. The regional mean glacier elevation change was −0.39 ± 0.13 m a−1 while the glaciers associated with the eight potentially dangerous lakes lowered by −0.71 ± 0.05 m a−1 from 1974 to 2017. The mean ice flow speed of these glaciers was ~10 m a−1 from 2013 to 2017; about double the mean for the entire study area. Analysis of these data along with climate observations suggests that ice melting and calving processes play the dominant role in driving lake enlargement. Modelling of future lake development shows where new lakes might emerge and existing lakes could expand with projected glacial recession.

The accelerated retreat of glaciers and the reduction of glacier ice reserves caused by climate change can significantly affect regional water resources and hydrological cycles. Changes in glacier thickness are among the key indicators that reflect this process. We analyzed changes observed in the elevation of glacier surfaces in the Gongga Mountains (GGM) using multi-temporal Digital Elevation Models (DEMs) derived from topographic maps, SRTM, ICESat and ZY-3 data. The results showed that the mean rate of change in glacier surface altitude in the GGM was ~−26.7 ± 2.03 m (0.54 ± 0.04 m a−1) between 1966 and 2015. The mean melt rates differed over different time periods, latterly showing an accelerating trend. As a general rule, glaciers appear to be losing more volume at lower than at higher elevations. Further analysis of these results suggests that supraglacial debris coverage in the GGM promotes higher rates of mass loss.

A total of eight floating glacier tongues have shrunk in area by >85% from the Yelverton Bay region of Northern Ellesmere Island since 1959, with unusually large losses since 2005. To better understand the causes of these losses, this study undertakes the first examination of ice tongue changes in this region, including an assessment of changes in surrounding marine ice (i.e. sea ice, sikussak and mélange), and atmospheric and oceanographic forcings. From 1959 to 2017, the total ice tongue area decreased by 49.07 km2, with the majority of this loss occurring from 2005 to 2009 (34.68 km2). The loss of ice tongues since 2005 occurred when open water replaced multi-year landfast sea ice (MLSI) and first-year sea ice in the regions adjacent to the ice tongues. These changes were accompanied by an increase in mean annual mid-depth (i.e. 100 and 200 m) ocean temperatures from −0.29°C from 1999 to 2005 to 0.67°C from 2006 to 2012. Despite the recent return of ocean temperatures to below pre-2006 levels, atmospheric summer temperatures have continued to rise (+0.15°C decade−1 between 1948 and 2016), with open water continuing to occur. Without the sustained presence of MLSI in this region the ice tongues are unable to stabilize, making it unlikely that they will re-form in the current climate.

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.

The long-term warming on the Antarctic Peninsula in the second half of the 20th century prompted rapid retreat of glaciers on the peninsula and surrounding islands. Retreat accelerated until the beginning of the new millennium when the regional warming trend significantly decreased. The response of glaciers to the change in temperature trend has been observed around the northern part of the Antarctic Peninsula but the timing of the shift from the surface lowering to mass gain remains unclear. Using historical aerial photographs, DEMs and satellite altimeter data from ICESat, we estimate areal and surface elevation changes of two small ice caps in the northern part of James Ross Island over the last 39 years. The glacierized area on Lachman Crags decreased from 4.337 ± 0.037 to 3.581 ± 0.014 km2 (−17.4%) between 1979 and 2006 and then increased to 3.597 ± 0.047 km2 (0.4%) until 2016. Surface lowering observed on ice caps after 1979 continued at least until 2008 as indicated by the ICESat data. The change from the lowering trend to increase in glacier surface elevation probably occurred after the ablation season 2008/09, which ranks among the warmest summers in the north-eastern Antarctic Peninsula since the mid-20th century.

To obtain information on changes in glacier mass balance in the central Nyainqentanglha Range, a comprehensive study was carried out based on digital-elevation models derived from the 1968 topographic maps, the Shuttle Radar Topography Mission DEM (2000) and TerraSAR-X/TanDEM-X (2013). Glacier area changes between 1968 and 2016 were derived from topographic maps and Landsat OLI images. This showed the area contained 715 glaciers, with an area of 1713.42 ± 51.82 km2, in 2016. Ice cover has been shrinking by 0.68 ± 0.05% a−1 since 1968. The glacier area covered by debris accounted for 11.9% of the total and decreased in the SE–NW directions. Using digital elevation model differencing and differential synthetic aperture radar interferometry, a significant mass loss of 0.46 ± 0.10 m w.e. a−1 has been recorded since 1968; mass losses accelerated from 0.42 ± 0.20 m w.e. a−1 to 0.60 ± 0.20 m w.e. a−1 between 1968–2000 and 2000–2013, with thinning noticeably greater on the debris-covered ice than the clean ice. Surface-elevation changes can be influenced by ice cliffs, as well as debris cover and land- or lake-terminating glaciers. Changes showed spatial and temporal heterogeneity and a substantial correlation with climate warming and decreased precipitation.

Previous studies reported that Good Friday Glacier had been actively surging in the 1950–60s, 1990s and again in 2000–15. Based on observations of terminus position change from air photos and satellite imagery, we fill the gaps between previous studies and conclude that the glacier has been advancing continuously since 1959. Ice surface velocities extracted from optical and synthetic aperture radar satellite images show higher flow rates than on most other marine-terminating glaciers in the region. This behaviour contrasts with the regional trend of glacier retreat over this period. Possible explanations involve a delayed response to positive mass-balance conditions of the Little Ice Age, or a dynamic instability. There is, however, insufficient evidence to attribute this behaviour to classical glacier surging as suggested in previous studies. Based on present-day ice velocity and glacier geometry patterns in the terminus region, we reconstruct the evolution of ice motion throughout the advance, and suggest that what has previously been interpreted as a surge, may instead have been a localised response to small-scale perturbations in bedrock topography.

Global-scale 21st-century glacier mass change projections from six published global glacier models are systematically compared as part of the Glacier Model Intercomparison Project. In total 214 projections of annual glacier mass and area forced by 25 General Circulation Models (GCMs) and four Representative Concentration Pathways (RCP) emission scenarios and aggregated into 19 glacier regions are considered. Global mass loss of all glaciers (outside the Antarctic and Greenland ice sheets) by 2100 relative to 2015 averaged over all model runs varies from 18 ± 7% (RCP2.6) to 36 ± 11% (RCP8.5) corresponding to 94 ± 25 and 200 ± 44 mm sea-level equivalent (SLE), respectively. Regional relative mass changes by 2100 correlate linearly with relative area changes. For RCP8.5 three models project global rates of mass loss (multi-GCM means) of >3 mm SLE per year towards the end of the century. Projections vary considerably between regions, and also among the glacier models. Global glacier mass changes per degree global air temperature rise tend to increase with more pronounced warming indicating that mass-balance sensitivities to temperature change are not constant. Differences in glacier mass projections among the models are attributed to differences in model physics, calibration and downscaling procedures, initial ice volumes and varying ensembles of forcing GCMs.

Changes in Greenland's marine-terminating outlet glaciers have led to changes in the flux of icebergs into Greenland's coastal waters, yet icebergs remain a relatively understudied component of the ice-ocean system. We developed a simple iceberg delineation algorithm for Landsat imagery. A machine learning-based cloud mask incorporated into the algorithm enables us to extract iceberg size distributions from open water even in partially cloudy scenes. We applied the algorithm to the Landsat archive covering Disko Bay, West Greenland, to derive a time series of iceberg size distributions from 2000–02 and 2013–15. The time series captures a change in iceberg size distributions, which we interpret as a result of changes in the calving regime of the parent glacier, Sermeq Kujalleq (Jakobshavn Isbræ). The change in calving style associated with the disintegration and disappearance of Sermeq Kujalleq's floating ice tongue resulted in the production of more small icebergs. The increased number of small icebergs resulted in increasingly negative power law slopes fit to iceberg size distributions in Disko Bay, suggesting that iceberg size distribution time series provide useful insights into changes in calving dynamics.

Arctic sea ice is going through a dramatic change in its extent and volume at an unprecedented rate. Sea-ice thickness (SIT) is a controlling geophysical variable that needs to be understood with greater accuracy. For the first time, a SIT-retrieval method that exclusively uses only airborne SIT data for training the empirical algorithm to retrieve SIT from Soil Moisture Ocean Salinity (SMOS) brightness temperature (TB) at different polarization is presented. A large amount of airborne SIT data has been used from various field campaigns in the Arctic conducted by different countries during 2011–15. The algorithm attempts to circumvent the issue related to discrimination between TB signatures of thin SIT versus low sea-ice concentration. The computed SIT has a rms error of 0.10 m, which seems reasonably good (as compared to the existing algorithms) for analysis at the used 25 km grid. This new SIT retrieval product is designed for direct operational application in ice prediction/climate models.

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.

Changes in glacier length and extent are indicators of contemporary and archives of past climate changes, but this common climate proxy presents a challenge for inferring a climate signal. Modeling studies suggest that length fluctuations can occur due to interannual climate variability within an unchanging mean climate and that changes in interannual climate variability can also drive changes in average length. This paper quantifies the impacts of interannual climate variability on average glacier length and mass balance, using a flowline model coupled to a simplified mass-balance model. Results illustrate that changes in the magnitude of interannual temperature variability can non-linearly affect the mean glacier length through a mass-balance asymmetry between warm and cold years. This asymmetry is present in models where melt only initiates after a temperature threshold is crossed. Glaciers susceptible to this asymmetry can be identified based on the shape of their mass-balance profiles. The presence of mass-balance asymmetries in glaciological databases is evaluated, but current records are too short for high statistical resolving power. While the asymmetry in this study can affect the average length and mass-balance, its impacts are small, and paleoclimate interpretations from glacier-length changes are likely not notably influenced by this process.

Glacier sliding over small obstacles relies on melting on their upstream sides and refreezing downstream. Previous treatments have appealed to ‘pressure melting’ as the cause of the spatial variations in melting temperature that drive this regelation process. However, we show that typical liquid pressure variations across small obstacles are negligible and therefore variations in ice pressure closely approximate variations in effective stress. For a given change in effective stress, the equilibrium melting temperature changes by an order of magnitude more than when the pressure of ice and liquid both change by an equal amount. In consequence, the temperature gradients that drive heat flow across small obstacles are larger than previously recognized and the rate of regelation is faster. Under typical conditions, the transition wavelength at which ice deformation and regelation contribute equally is of m-scale, several times longer than previous predictions, which have been reported to underestimate field inferences.