Supraglacial debris is significant in many regions and complicates modeling of glacier melt, which is required for predicting glacier change and its influences on hydrology and sea-level rise. Temperature-index models are a popular alternative to energy-balance models when forcing data are limited, but their transferability among glaciers and inherent uncertainty have not been documented in application to debris-covered glaciers. Here, melt factors were compiled directly from published studies or computed from reported melt and MERRA-2 air temperature for 27 debris-covered glaciers around the world. Linear mixed-effects models were fit to predict melt factors from debris thickness and variables including debris lithology and MERRA-2 radiative exchange. The models were tested by leave-one-site-out cross-validation based on predicted melt rates. The best model included debris thickness (fixed effect) and glacier and year (random effects). Predictions were more accurate using MERRA-2 than on-site air temperature data, and pooling MERRA-2-derived and reported melt factors improved cross-validation accuracy more than including additional predictors such as shortwave or longwave radiation. At one glacier where monthly ablation was measured over 4 years, seasonal variation of melt factors suggested that heat storage significantly affected the relation between melt and energy exchange at the debris surface.

Glaciers and the periglacial environment in Argentina have been protected by the Law since 2010. This legislation required the development of the first National Glacier Inventory (NGI), which was officially presented in May 2018 and based on satellite images spanning between 2004 and 2016. Here, we present the methods and results of the NGI, summarize the glaciers’ morphological and spatial characteristics, and compare our results to previous regional and global inventories. The NGI reveals an impressive variety of ice masses including rock glaciers, permanent snowfields, mountain and valley glaciers with varying amounts of debris-cover and large outlet glaciers. The Argentinean Andes contain 16 078 ice masses covering an area of 5769 km2 between 200 and 6900 m a.s.l. Comparison of the combined national inventories of Argentina and Chile (~30 000 glaciers and 28 400 km2) with the Randolph Glacier Inventory 6.0 for the Southern Andes (~16 000 glaciers and 29 400 km2), shows that there are large differences in extent and number of glaciers in some sub-regions. The NGI represents an improvement for a better understanding of Argentina's freshwater reservoirs and provides detailed information for the preservation and study of ice masses along 4000 km of the Southern Andes.

The representation of iceberg calving in numerical models is a key source of uncertainty in century-scale sea-level rise projections. Parameters central to model representations of calving, including the tensile strength of glacier ice, remain poorly constrained. Grain-size and sample-size dependence make it difficult to reconcile laboratory and in situ estimates of ice tensile strength. Further, assumptions of various numerical models obscure comparison of the ‘strength’ parameter with a physically observable value. Here, we address the problem of fracture during calving using an analogous natural laboratory: a viscoelastic analysis of observed surface deformation and associated stresses in the 2015 collapse of eastern Skaftá cauldron, Vatnajökull ice cap, Iceland. We find that the ice within the cauldron could have experienced instantaneous elastic stress on the order of several MPa. We observe surface crevasses at the cauldron edges and center, but find that large areas of ice remain intact despite high stress. Our findings suggest a tensile strength of glacier ice on the order of 1 MPa, consistent with laboratory estimates but exceeding previous glacier-specific estimates.

Under the potential to reconstruct the past climatic and atmospheric conditions from a deep ice core in the coastal Antarctic site (Styx Glacier), an 8.84 m long firn core (73°50.975′ S, 163°41.640′ E; 1623 m a.s.l.) was initially studied to propose a reliable age scale for the local estimation of snow accumulation rate. The seasonal variations of δ18O, methanesulfonic acid (MSA) and non-sea-salt sulfate (nssSO42–) were used for the firn core dating and revealed 25 annual peaks (from 1990 to 2014) with volcanic sulfate signal. The observed declining trend in annual accumulation rate with a mean value of 146 ± 60 kg m–2 a–1 is likely to be linked to the changes of sea-ice extent in the Ross Sea region. Moreover, the temporal variation of the annual mean δ18O, an annual flux of MSA and nssSO42– also likely to be under the influence of ice-covered and open water area. This study suggests a potential to recover past changes in an oceanic environment and will be useful for the interpretation of the long ice core drilled at the same site.

Each summer, surface melting of the margin of the Greenland Ice Sheet exposes a distinctive visible stratigraphy that is related to past variability in subaerial dust deposition across the accumulation zone and subsequent ice flow toward the margin. Here we map this surface stratigraphy along the northern margin of the ice sheet using mosaicked Sentinel-2 multispectral satellite imagery from the end of the 2019 melt season and finer-resolution WorldView-2/3 imagery for smaller regions of interest. We trace three distinct transitions in apparent dust concentration and the top of a darker basal layer. The three dust transitions have been identified previously as representing late-Pleistocene climatic transitions, allowing us to develop a coarse margin chronostratigraphy for northern Greenland. Substantial folding of late-Pleistocene stratigraphy is observed but uncommon. The oldest conformal surface-exposed ice in northern Greenland is likely located adjacent to Warming Land and may be up to ~55 thousand years old. Basal ice is commonly exposed hundreds of metres from the ice margin and may indicate a widespread frozen basal thermal state. We conclude that the ice margin across northern Greenland offers multiple opportunities to recover paleoclimatically distinct ice relative to previously studied regions in southwestern Greenland.

Climate-induced cryospheric changes can have a significant impact on the downstream water availability. In this study, the Open Global Glacier Model (OGGM) and the Glacio-hydrological Degree-day Model (GDM) are integrated to project the response of cryospheric and hydrological systems to climate change until 2100. The study area comprises six sub-basins of glacierized Koshi River basin covering Nepalese and Chinese territories. The output from OGGM is provided as input to GDM along with the spatial and hydro-meteorological data. The average glacier area change in all the sub-basins from 2021 to 2100 is estimated as 65 and 85% decrease and the average glacier volume change is estimated as 76 and 86% decrease for RCP 4.5 and 8.5 scenarios, respectively. The future simulated discharge shows an increasing trend in pre-monsoon and monsoon seasons and a decreasing trend in post-monsoon and winter seasons after 2060 in all the sub-basins, which can lead to wetter wet seasons and drier dry seasons in the far future. A shift in peak flow is observed from August to July in most of the sub-basins. The coupled modelling technique used in this study can largely improve our understanding of glacio-hydrological dynamics in the Himalayan region.

Glacier ice thickness is crucial to quantifying water resources in mountain regions, and is an essential input for ice-flow models. Using a surface velocity inversion method, we combine ice thickness measurements with detailed surface elevation and velocity data, and derive ice thickness and volume estimates for the Monte Tronador glaciers, North Patagonian Andes. We test the dependence of the inversion model on surface slope by resampling glacier slopes using variable smoothing filter sizes of 16–720 m. While total glacier volumes do not differ considerably, ice thickness estimates show higher variability depending on filter size. Smaller (larger) smoothing scales give thinner (thicker) ice and higher (lower) noise in ice thickness distribution. A filter size of 300 m, equivalent to four times the mean ice thickness, produces a noise-free thickness distribution with an accuracy of 35 m. We estimate the volume of the Monte Tronador glaciers at 4.8 ± 2 km3 with a mean ice thickness of 75 m. Comparison of our results with earlier regional and global assessments shows that the quality of glacier inventories is a significant source of discrepancy. We show that including surface slope as an input parameter increases the accuracy of ice thickness distribution estimates.

Using an ensemble of close- and long-range remote sensing, lake bathymetry and regional meteorological data, we present a detailed assessment of the geometric changes of El Morado Glacier in the Central Andes of Chile and its adjacent proglacial lake between 1932 and 2019. Overall, the results revealed a period of marked glacier down wasting, with a mean geodetic glacier mass balance of −0.39 ± 0.15 m w.e.a−1 observed for the entire glacier between 1955 and 2015 with an area loss of 40% between 1955 and 2019. We estimate an ice elevation change of −1.00 ± 0.17 m a−1 for the glacier tongue between 1932 and 2019. The increase in the ice thinning rates and area loss during the last decade is coincident with the severe drought in this region (2010–present), which our minimal surface mass-balance model is able to reproduce. As a result of the glacier changes observed, the proglacial lake increased in area substantially between 1955 and 2019, with bathymetry data suggesting a water volume of 3.6 million m3 in 2017. This study highlights the need for further monitoring of glacierised areas in the Central Andes. Such efforts would facilitate a better understanding of the downstream impacts of glacier downwasting.

A volume-area scaling relation is commonly used to estimate glacier volume or its future changes on a global scale. The presence of an insulating supraglacial debris cover alters the mass-balance profile of a glacier, potentially modifying the scaling relation. Here, the nature of scaling relations for extensively debris-covered glaciers is investigated. Theoretical arguments suggest that the volume-area scaling exponent for these glaciers is ~7% smaller than that for clean glaciers. This is consistent with the results from flowline-model simulations of idealised glaciers, and the available data from the Himalaya. The best-fit scale factor for debris-covered Himalayan glaciers is ~60% larger compared to that for the clean ones, implying a significantly larger stored ice volume in a debris-covered glacier compared to a clean one having the same area. These results may help improve scaling-based estimates of glacier volume and future glacier changes in regions where debris-covered glaciers are abundant.

Effective snow grain radius (re) is mapped at high resolution using near-infrared hyperspectral imaging (NIR-HSI). The NIR-HSI method can be used to quantify re spatial variability, change in re due to metamorphism, and visualize water percolation in the snowpack. Results are presented for three different laboratory-prepared snow samples (homogeneous, ice lens, fine grains over coarse grains), the sidewalls of which were imaged before and after melt induced by a solar lamp. The spectral reflectance in each ~3 mm pixel was inverted for re using the scaled band area of the ice absorption feature centered at 1030 nm, producing re maps consisting of 54 740 pixels. All snow samples exhibited grain coarsening post-melt as the result of wet snow metamorphism, which is quantified by the change in re distributions from pre- and post-melt images. The NIR-HSI method was compared to re retrievals from a field spectrometer and X-ray computed microtomography (micro-CT), resulting in the spectrometer having the same mean re and micro-CT having 23.9% higher mean re than the hyperspectral imager. As compact hyperspectral imagers become more widely available, this method may be a valuable tool for assessing re spatial variability and snow metamorphism in field and laboratory settings.

Velocities within an icefall are typically the fastest within a glacier system and experience complex flow. The combination of convergent and fast flow, and steep slope generate a quickly changing and intensely fractured surface. This complicates velocity extraction from repeat satellite images, especially when common pattern matching procedures are used. In this study, we exploit the high temporal revisit of medium-resolution satellite images using a novel image matching technique, ensemble matching, making it possible to generate a high-resolution (30 m) velocity field from high-repeat image sequences despite challenging image conditions. We demonstrate this technique for the first time in the glaciology domain using repeat Sentinel-2 optical data over the famous Khumbu icefall, situated on the southern slopes of Mount Everest. Estimates of velocity go just over 1 m d−1, which is slower than summer velocities from noisy single pair image matching. This icefall is frequently crossed by high-altitude mountaineers who use a route confined by fixed ropes and ladders set out every season. The mountain climbers typically record their trajectory on their personal satellite navigation device. We use such volunteered geographic information to verify our velocity estimates, confirming our underestimation with ensemble matching. Besides unprecedented remotely sensed surface velocities over the icefall, we also note that the generated velocity field can aid with the planning of a safe passage through this icefall.

Turbulent fluxes make a substantial and growing contribution to the energy balance of ice surfaces globally, but are poorly constrained owing to challenges in estimating the aerodynamic roughness length (z0). Here, we used structure from motion (SfM) photogrammetry and terrestrial laser scanning (TLS) surveys to make plot-scale 2-D and 3-D microtopographic estimations of z0 and upscale these to map z0 across an ablating mountain glacier. At plot scales, we found spatial variability in z0 estimates of over two orders of magnitude with unpredictable z0 trajectories, even when classified into ice surface types. TLS-derived surface roughness exhibited strong relationships with plot-scale SfM z0 estimates. At the glacier scale, a consistent increase in z0 of ~0.1 mm d−1 was observed. Space-for-time substitution based on time since surface ice was exposed by snow melt confirmed this gradual increase in z0 over 60 d. These measurements permit us to propose a scale-dependent temporal z0 evolution model where unpredictable variability at the plot scale gives way to more predictable changes of z0 at the glacier scale. This model provides a critical step towards deriving spatially and temporally distributed representations of z0 that are currently lacking in the parameterisation of distributed glacier surface energy balance models.

Using frequency-modulated continuous wave radar data from the 32nd Chinese Antarctic Research Expedition in 2015/16, subsurface profiles were obtained along an East Antarctic inland traverse from Zhongshan station to Dome A, and four distinct regions were selected to analyze the spatiotemporal variability in historical surface mass balance (SMB). Based on depth, density, and age data from ice cores along the traverse, the radar data were calibrated to yield average SMB data. The zone 49–195 km from the coast has the highest SMB (235 kg m−2 a−1). The 780–892 km zone was most affected by the Medieval Warm Period and the Little Ice Age, and the SMB during ad 1454–1836 (71 kg m−2 a−1) was only one-quarter of that in the 20th century. The SMB in the 1080–1157 km zone fluctuates the most, possibly due to erosion or irregular deposition of snow by katabatic winds in low SMB areas with surface elevation fluctuations. Dome A (1157–1236 km) has the lowest SMB (29 kg m−2 a−1) and did not decrease during Little Ice Age. Understanding the spatiotemporal variability of SMB in a larger space can help us understand the complex climate history of Antarctica.

Ice shelves play a critical role in modulating dynamic loss of ice from the grounded portion of the Antarctic Ice Sheet and its contribution to sea-level rise. Measurements of ice-shelf motion provide insights into processes modifying buttressing. Here we investigate the effect of seasonal variability of basal melting on ice flow of Ross Ice Shelf. Velocities were measured from November 2015 to December 2016 at 12 GPS stations deployed from the ice front to 430 km upstream. The flow-parallel velocity anomaly at each station, relative to the annual mean, was small during early austral summer (November–January), negative during February–April, and positive during austral winter (May–September). The maximum velocity anomaly reached several metres per year at most stations. We used a 2-D ice-sheet model of the RIS and its grounded tributaries to explore the seasonal response of the ice sheet to time-varying basal melt rates. We find that melt-rate response to changes in summer upper-ocean heating near the ice front will affect the future flow of RIS and its tributary glaciers. However, modelled seasonal flow variations from increased summer basal melting near the ice front are much smaller than observed, suggesting that other as-yet-unidentified seasonal processes are currently dominant.

Glaciers in the eastern Pamir have reportedly been gaining mass during recent decades, even though glaciers in most other regions in High Mountain Asia have been in recession. Questions still remain about whether the trend is strengthening or weakening, and how far the positive balances extend into the eastern Pamir. To address these gaps, we use three different digital elevation models to reconstruct glacier surface elevation changes over two periods (2000–09 and 2000–15/16). We characterize the eastern Pamir as a zone of transition from positive to negative mass balance with the boundary lying at the northern end of Kongur Tagh, and find that glaciers situated at higher elevations are those with the most positive balances. Most (67% of 55) glaciers displayed a net mass gain since the 21st century. This led to an increasing regional geodetic glacier mass balance from −0.06 ± 0.16 m w.e. a−1 in 2000–09 to 0.06 ± 0.04 m w.e. a−1 in 2000–15/16. Surge-type glaciers, which are prevalent in the eastern Pamir, showed fluctuations in mass balance on an individual scale during and after surges, but no statistical difference compared to non-surge-type glaciers when aggregated across the region.