Establishing an accurate chronology is crucial for interpretation of ice core-based climatic records. While high snow accumulation rates characterise coastal Antarctica, thus enabling recovery of highly resolved climatic records, summertime melting at such low-elevation sites offers challenges in establishing a reliable chronological framework through traditional approaches using the seasonality of stable water isotope and ionic proxy records. Here, we assess visual stratigraphy (VS) obtained from line-scan images as a proxy for annual layer counting in firn section (top 50 m) of the IND-36/B9 ice core (dated 1919–2016 CE) from the Djupranen Ice Rise in central Dronning Maud Land, East Antarctica. We also used these images to obtain melt history for the site and found that traditional thickness-based quantification of melt proportion results in significant overestimations. Since density has dominant control on VS profile over the firn section, we first used circulant single-spectrum analysis to remove the secular trend and then we extracted the seasonal VS signals attributed to dust and sea-salt inclusions. We find that melt layers do not significantly alter the VS records if masked during pre-processing. The age–depth model based on the reconstructed VS profile revealed an excellent match with identified time-markers within an uncertainty of ±2 years.

Avalanches result from an interaction of weather and terrain, where past weather and internal snow cover processes play important roles. So far, climatology was mainly based on weather data, as regional snow instability information, such as avalanche activity, is scarce on climatological time scales. We present a new approach to create a snow avalanche climatology from simulations of avalanche problem types based on snow cover simulations of reanalysis data and a cluster analysis. Analyzing the winters between 1958 and 2020 in the French Alps, wet-snow situations dominated natural release. Dry-snow situations with non-persistent and persistent weak layers occurred each on at least one third of the days. Four typical patterns of avalanche problem types were identified. They follow the main orography with more new snow situations in the northern regions and more cases of persistent weak layers in inner-Alpine regions. In the front-ranges and in southern regions wet-snow situations occurred early in winter – typical for coastal snow climates. Agreement with the standard snow climate classification and the geography of the French Alps suggests that mountain regions with similar conditions can now be outlined. This method for snow avalanche climatology will inform avalanche forecasting and facilitate climate change impact studies.

This work proposes a holistic quantitative snow avalanche risk assessment that evaluates, at reasonable computational costs and for various types of buildings, the impact of forest cover changes on the probability distribution of runout distances, impact pressures and subsequent risk estimates. A typical case study of the French Alps shows that, from a completely deforested to a completely forested path, avalanche risk for a building located downslope decreases by 53–99%, depending on how forest cover is accounted for in avalanche statistical–dynamical modeling. Local forest cover data inferred from old maps and photographs further demonstrates that a 20–60% risk reduction actually occurred between 1825 and 2017 at the site because of the afforestation dynamics, with significant modulations according to the considered building technology. These results (1) assert the protective role of forests against snow avalanches, (2) highlight the potential of combining nature-based solutions with traditional structural measures to reduce risk to acceptable levels at reasonable costs, (3) suggest a significant decrease in risk to settlements in areas that encountered similar forest cover changes and (4) open the door to the quantification of long-term avalanche risk changes as a function of changes of all its hazard, vulnerability and exposure drivers in various mountain context.

We use a mathematical model to investigate the effect of basal topography and ice surface slope on transport and deposition of sediment within a water-filled subglacial channel. In our model, three zones of different behaviour occur. In the zone furthest upstream, variations in basal topography lead to sediment deposition under a wide range of conditions. In this first zone, even very small and gradually varying basal undulations (~5 m amplitude) can lead to the deposition of sediment within a modelled channel. Deposition is concentrated on the downstream gradient of subglacial ridges, and on the upstream gradient of subglacial troughs. The thickness and steepness of the ice sheet has a substantial impact on deposition rates, with shallow ice profiles strongly promoting both the magnitude and extent of sediment deposition. In a second zone, all sediment is transported downstream. Finally, a third zone close to the ice margin is characterised by high rates of sediment deposition. The existence of these zones has implications for esker formation and the dynamics of the subglacial environment.

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.

Surging and streaming of glaciers are modulated by meltwater availability and pressure which controls mechanical coupling at their beds. Using laboratory-scale experimental modelling and palaeoglaciological mapping, we explore how subglacial drainage landsystems control meltwater drainage efficiency and ice flow velocities for terrestrial-based ice lobes resting on flat horizontal and permeable beds. Two end-members regimes, surging and streaming, appear in our experiments. The surge regime is characterised by a rapid increase of drainage efficiency through development of tunnel valleys and their tributaries, thus reducing the duration of ice flow speed-up events by lowering water pressures and increasing ice-bed coupling. Tunnel valleys connected to ice lobe margins, submarginal thrust moraines, reduced ice lobe extensions and ephemeral shear margins are the most distinctive characteristics of this regime. The stream regime is characterised by disconnected channels of smaller dimensions unable to evacuate all the meltwater: this prolonged drainage inefficiency leads to sustained high ice flow velocity and steady shear margins. Small and rectilinear meltwater channels devoid of tributaries, often disconnected from ice lobe margins, and lineation swarms are diagnostic of this regime.

Climate change has had a significant impact on glacier recession, particularly in the Arctic, where glacier meltwater is an important contributor to global sea-level rise. Therefore, it is important to accurately quantify glacier recession within this sensitive region, using multiple observations of glacier extent. In this study, we mapped 480 glaciers in Novaya Zemlya, Russian Arctic, using object-based image analysis applied to multispectral Landsat satellite imagery in Google Earth Engine and quantify the area changes between 1986–89 and 2019–21. The results show that in 1986–89, the total glacierized area was 22 990 ± 301 km2, in 2000–01 the area was 22 525 ± 308 km2 and by 2019–21 the glacier area reduced to 21 670 ± 292 km2, representing a total of 5.8% reduction in glacier area between 1986–89 and 2019–21. Higher glacier area loss was observed on the Barents Sea coast (7.3%) compared to the Kara (4.2%), reflecting previously observed differences in warming trends. The accuracy of the automatically generated outlines of each layer (1986–89, 2000–01 and 2019–21) was evaluated by comparing with manually corrected outlines (reference data) using random sampling, resulting in an overall accuracy estimate of between 96 and 97% compared to the reference data. This automated approach in Google Earth Engine is a promising tool for rapidly mapping glacier change that reduces the amount of time required to generate accurate glacier outlines.

The British Antarctic Survey's (BAS) Rapid Access Isotope Drill (RAID), designed for rapid drilling to survey prospective ice core sites, has been deployed at multiple Antarctic locations over 6 years. This drilling method creates ice chippings that can be discretely sampled and analysed for their chemical and water isotopic composition. Ice sampling methods have evolved since the first uses of the BAS RAID, enabling a more quantifiable sample resolution. Here, we show that water isotope records obtained from RAID ice are comparable to those of equivalent depth resolution from proximal ice cores. Records of chemical impurities also show good agreement with nearby cores. Our findings suggest that the RAID is suitable for both chemical and isotopic reconnaissance of drilling sites. Residual contamination of certain ions is discussed, with proposed design changes to avoid this issue with future use.

Tasman Glacier, a temperate maritime glacier in the New Zealand Southern Alps, is rapidly receding. Climate warming is resulting in lengthening of the ablation season, meaning crevasses in the accumulation area are becoming exposed at the surface for longer. We combine measurements of air temperature and wind speed from inside crevasses with surface meteorological data, finding that during summer, in-crevasse air temperature is frequently positive, and can at times exceed surface air temperature. Greatest warming occurred in the widest crevasses during clear-sky conditions, but full depth warming of crevasses also occurred at night. Net shortwave radiation contributes to heating of air in the upper regions of crevasses, but turbulent sensible heat transfer was responsible for driving warm air deeper into crevasses. Crevasses orientated to maximise radiation retrieval, and running parallel to wind flow, have the greatest potential for warming and heat storage. We hypothesise a positive feedback loop in the surface energy-balance system, where crevasses entrain and trap heat, which enhances melting, that in turn enlarges the crevasses, enabling greater heat storage and further melting. Energy-balance models that treat accumulation areas of alpine glaciers as homogeneous surfaces will therefore underestimate snow melt and overestimate mass balance.

Changes in iceberg calving fluxes and oceanographic conditions around Antarctica have likely influenced the spatial and temporal distribution of iceberg fresh water fluxes to the surrounding ocean basins. However, Antarctic iceberg melt rate estimates have been limited to very large icebergs in the open ocean. Here we use a remote-sensing approach to estimate iceberg melt rates from 2011 to 2022 for 15 study sites around Antarctica. Melt rates generally increase with iceberg draft and follow large-scale variations in ocean temperature: maximum melt rates for the western peninsula, western ice sheet, eastern ice sheet and eastern peninsula are ~50, ~40, ~5 and ~5 m a−1, respectively. Iceberg melt sensitivity to thermal forcing varies widely, with a best-estimate increase in melting of ~24 m a−1°C−1 and range from near-zero to ~100 m a−1°C−1. Variations in water shear likely contribute to the apparent spread in thermal forcing sensitivity across sites. Although the sensitivity of iceberg melt rates to water shear prevents the use of melt rates as a proxy to infer coastal water mass temperature variability, additional coastal iceberg melt observations will likely improve models of Southern Ocean fresh water fluxes and have potential for subglacial discharge plume mapping.

A precise knowledge of landfast sea-ice (fast-ice) thickness is relevant to many different disciplines. Sea Ice Monitoring Stations (SIMS) are used to measure time series of fast-ice thickness at a location. SIMS measure ice and ocean temperature via thermistor strings with many different methods for extracting sea-ice thickness from temperature existing. This study investigates: if thickness results from temperature recorded by SIMS of different designs, and analysed with different methods are comparable; which methods are recommended for their robustness, precision and accuracy and how they compare to independent thickness measurements; how otherwise unuseable data can be salvaged through specific SIMS designs. We present an analysis of fast-ice thickness calculated from SIMS deployed in McMurdo Sound, Antarctica and in the Chukchi Sea near Utqiaġvik, Alaska, over two decades. We find that median thicknesses derived by different methods agree within 1 ± 1.5 cm for McMurdo Sound and 2 ± 3 cm for Utqiaġvik. Thus, it is possible to confidently compare data collected from different stations and analysed with different methods. The vertical gradient of sea-ice temperature gives the best results for fast-ice thickness during the growth season and including standard resistors in a thermistor string can reduce potential data loss due to noise.

We quantify the surface elevation changes along Yala Glacier in Langtang Valley, Nepal Himalaya, since 1981 using geodetic methods to understand the recent evolution and current state of small debris-free glaciers across the region. We analyse differential global positioning system measurements and aerial stereo imagery that were acquired along Yala Glacier in 2007, 2009, 2012 and 2015 to generate digital elevation models for each calculation period. Continuous surface lowering has mainly been observed across the down-glacier area during the calculation periods, although a large degree of variability exists, with this lowering trend propagating up-glacier in recent years. The area-weighted glacier mass balances range from −0.98 ± 0.27 to −0.26 ± 0.30 m w.e. a−1 for the five calculation periods (1981–2007, 2007–2009, 2009–2012, 2012–2015 and 2007–2015). These calculated mass-balance data reveal that Yala Glacier has undergone accelerated mass loss since the late 2000s, which is consistent with the results of previous in situ measurement and remote-sensing studies.

Mass loss from the Amundsen Sea Embayment of the West Antarctic Ice Sheet is a major contributor to global sea-level rise (SLR) and has been increasing over recent decades. Predictions of future SLR are increasingly modelled using ensembles of simulations within which model parameters and external forcings are varied within credible ranges. Accurately reporting the uncertainty associated with these predictions is crucial in enabling effective planning for, and construction of defences against, rising sea levels. Calibrating model simulations against current observations of ice-sheet behaviour enables the uncertainty to be reduced. Here we calibrate an ensemble of BISICLES ice-sheet model simulations of ice loss from the Amundsen Sea Embayment using remotely sensed observations of surface elevation and ice speed. Each calibration type is shown to be capable of reducing the 90% credibility bounds of predicted contributions to SLR by 34 and 43% respectively.

Marine-terminating glaciers cover more than one-fourth of the total glacierized area in the Northern Hemisphere outside the Greenland ice sheet. It is therefore crucial to ensure an adequate representation of these glaciers when projecting large-scale glacier mass changes. We investigate how the introduction of marine frontal processes in the modeling chain influences the results of mass change projections, compared to projections neglecting such processes. We find that including frontal processes reduces the projected glacier mass loss, since incorporating frontal ablation in the model's mass-balance calibration results in a decrease in marine-terminating glaciers’ sensitivity to atmospheric temperatures. We also find that retrograde bed slopes lead to increased frontal ablation as the atmosphere warms, while frontal ablation decreases if bed slopes are prograde. These opposing effects have the potential to partly cancel each other when considering large glacier ensembles. Although we do not account for potential future changes in oceanic climate yet, any effect of these would be moderated by around half of today's marine-terminating glaciers becoming land-terminating in the course of the 21st century. While we find a significant influence of ice flow parameters on our results, boundary conditions remain the largest source of uncertainty in our projections.

We assess the composition and geometry of four individual rock glaciers in Alaska, Wyoming and Colorado by measuring their radio wave speed and applying these results to ground-penetrating radar depth corrections and dielectric mixing models. Our method includes a correction for subsurface reflector dip angle, which we show can lead to an incorrect determination of wave speeds using common midpoint configurations. By observing the radar properties of the rock glaciers and their supraglacial debris, we find that some of the sites exhibit nearly pure ice cores, and all of the sites indicate volumetric ice fractions >50%. These results have implications for terrestrial glaciology and hydrology because the present ice volume is connected to past ice accumulation and subsurface ice preservation, which may affect the future availability of alpine water resources. An understanding of the processes that govern rock glacier evolution over a wide range of latitudes and elevations will also contribute to the exploration of planetary surfaces such as Mars, which hosts a significant population of debris-covered glaciers. Our subsurface composition and geometry estimates will inform simulations of rock glacier formation and evolution to test hypothesized ice origin mechanisms along with the preservation of climate signals.

A recent study (Millan and others, 2022a, Nature Geoscience 15(2), 124–129) claims that ice volume contained in all glaciers outside the ice sheets and its potential contribution to sea level is 20% less than previously estimated. However, the apparent decrease is largely due to differences in choice of domain, as the study excludes 80% of the glacier area in the Antarctic periphery that was included in previous global glacier volume estimates. The issue highlights the difficulty in separating glaciers from the ice-sheet proper, especially in Antarctica, and the need for both the glacier and ice-sheet communities to develop standards and protocols to avoid double-counting in global ice volume and mass-change assessments and projections. Process-based inversion models have replaced earlier scaling methods, but large uncertainties in global glacier volume estimation remain due to the ill-posed nature of the inversion problem and poorly constrained parameters emphasizing the need for more direct ice thickness observations.

This study provides the first comprehensive reconstruction of the dynamics of Iceberg Glacier, located on western Axel Heiberg Island, and reveals detailed observations of a complete surge for the first time in the Canadian Arctic. Historical aerial photographs, declassified intelligence satellite photographs, optical satellite imagery and synthetic aperture radar data were used to quantify changes in terminus position, ice velocity and glacier thickness since the 1950s. A surge initiated at the terminus in 1981 and terminated in 2003, suggesting a 22-year active phase. High surface velocities, reaching ~2300 m a−1 in 1991, were accompanied by a maximum terminus advance of >7 km and a large transfer of mass down-glacier, causing significant median trunk-wide surface elevation changes attaining >3 ± 1 m a−1. We suggest that the retreat from a pinning point, flotation of the terminus, the removal of sea-ice from the ice front, and an increase in subglacial meltwater availability from relatively high air temperatures in 1981 likely contributed to surge initiation. The ensuing quiescent period has seen a continual decrease in surface flow rates to an average centreline velocity of 11.5 m a−1 in 2020–21, a gradual steepening of the glacier surface and a > 2.5 km terminus retreat.

Once fallen, snow settles due to the combined effects of metamorphism and deformation of the ice matrix under gravity. To understand how these coupled processes affect snow evolution, we performed oedometric compression tests and continuously monitored the snow microstructure with X-ray tomography. Centimetric samples with an initial density between 200 and 300 kg m−3 were followed during an initial sintering phase and under two different loads of 2.1 and 4.7 kPa at $-8^\circ$C for ~1 week. The microstructure captured at a voxel size of 8.5 μm was characterized by density, specific surface area (SSA) and two metrics related to bond network, namely the Euler characteristic and the minimum cut surface. Load-induced creep of the ice matrix was observed only for sufficiently low values of initial density (<290 kg m−3 in our tests), and was shown to be associated to a significant increase of the number of bonds. Application of the load, however, did not affect the individual bond size nor the SSA, which appeared to be mainly controlled by isothermal metamorphism. The uniaxial compression did not induce any creation of anisotropy on the microstructural characteristics. Overall, our results show that, for the considered conditions, the deformation of the ice matrix mainly leads to a reduction of the pore space and an increase of the coordination number, while metamorphism mainly affects the grain and bond sizes.

The Ward Hunt and Milne ice shelves are the present-day remnants of a much larger ice shelf that once fringed the coast of Ellesmere Island, Canada. These ice shelves possess a unique surface morphology consisting of wave-like rolls that run parallel to the shoreline. Setting aside the question of how these rolls originally developed, we consider the impact of this roll morphology on the stability of the ice shelf. In particular, we examine whether periodic variations in ice-shelf thickness and water depth implied by the rolls prevent the excitation of Lamb waves in the ice shelf. Using a hierarchy of numerical models, we find that there are band gaps in the flexural and extensional modes of the ice shelf, implying the existence of frequency ranges that lack wave motion. We show that an ice shelf with rolls is able to reflect waves in these frequency ranges that are incident upon its ice front, thereby mitigating undue stress and calving. We speculate that the roll morphology provides a “fitness” for survival that explains why rolls are observed in the oldest and thickest multiyear sea ice of the Arctic.