The mass balance of glaciers requires more detailed and continuous observations to understand their seasonal change in relation to climate. Here, we designed and installed an automated real-time monitoring platform at 4645 m a.s.l. on the Baishui River Glacier No.1 to collect continuous high-resolution observational data, and analyzed the seasonal dynamic from glacier movement and surface mass balance from glacier melting and snow accumulation. The results showed that the platform moved northeastward ~12.9 m at a rate of 0.06 ± 0.02 m d−1 between September 2021 and April 2022. The surface mass balance showed a varied temporal period. July and August were the main ablation periods, while ablation decreased and ceased in September. The glacier neither melted nor accumulated much between October and December, but began to have rapid snow accumulation in January. The glacier surface temperature varied with the air temperature and showed significant inter-seasonal differences among monsoon, post-monsoon and winter seasons. The surface mass balance also exhibited a strong response to the air temperature changes, with an average decrease of 1°C the point mass balance increased by 0.11 m w.e. from monsoon to post-monsoon and 0.22 m w.e. from post-monsoon to winter. Moreover, we found snowfall caused a decrease in the glacier surface temperature by increasing the surface albedo.

The uncertainty of glacier change projections is largely influenced by glacier models. In this study, we focus on temperature-index mass-balance (MB) models and their calibration. Using the Open Global Glacier Model (OGGM), we examine the influence of different surface-type dependent degree-day factors, temporal climate resolutions (daily, monthly) and downscaling options (temperature lapse rates, temperature and precipitation corrections) for 88 glaciers with in-situ observations. Our findings indicate that higher spatial and temporal resolution observations enhance MB gradient representation due to an improved calibration. The addition of surface-type distinction in the model also improves MB gradients, but the lack of independent observations limits our ability to demonstrate the added value of increased model complexity. Some model choices have systematic effects, for example weaker temperature lapse rates result in smaller projected glaciers. However, we often find counter balancing effects, such as the sensitivity to different degree-day factors for snow, firn and ice, which depends on how the glacier accumulation area ratio changes in the future. Similarly, using daily versus monthly climate data can affect glaciers differently depending on the shifting balance between melt and solid precipitation thresholds. Our study highlights the importance of considering minor model design differences to predict future glacier volumes and runoff accurately.

Recent satellite measurements of glacier mass balance show mountain glaciers all over the world had generally negative mass balances in the first decades of the 21st century. We analyse archived data for surface mass balance and summer temperature for 37 Northern Hemisphere glaciers with data for 1961–2020. We compare mean annual balances for 1961–90 and 1991–2020, and for 25 glaciers explain the changes in annual balance by changes in winter and summer balances. Mean balances 1961–90 were already substantially negative for 19 out of the 25 glaciers. Changes in winter balances from 1961–90 to 1991–2020 average close to zero but changes in summer balance are strongly negative. Mean balance 1991–2020 is strongly correlated with change in summer balance, weakly correlated with winter balance change, and strongly correlated with mean balance 1961–90. We estimate 1991–2020 summer temperature anomalies for the 37 glaciers and confirm that summer temperature anomalies for 1991–2020 were higher for the Alps, by nearly 1.5°C, than for other areas. Substantial variations in the temperature – sensitivity of summer balances for individual glaciers of −0.2 to −1.0 m w.e. a−1 °C−1 deserve further study.

As warming intensifies across the Greenland ice sheet, an increasing number of shallow coring and radar studies are targeting the melt-impacted firn column to investigate meltwater processes. Highly inhomogeneous infiltration and refreezing, however, redistributes mass, distorting age–depth relationships and confounding comparisons between different cores. Here, we utilize a dynamic time warping algorithm for time series alignment to investigate and quantify the heterogeneous impact of melt processes on nine closely spaced (within 50 m) firn core-density profiles. The 10 m cores were collected relatively high in Greenland's percolation zone, where melt alteration is minimal compared to lower elevation. Our analysis demonstrates the effectiveness of dynamic time warping as a tool for assessing heterogeneity between ice core records. We find that the optimal alignment of density profiles in the nine cores requires vertical stretching and compression of individual profiles, ranging from, on average, <1 to ~16% of the core lengths. We identified four depth zones of mass redistribution that appear to coincide with observed ice layers. Further, ~75% of density measurements from each core do not align with an age model-derived density profile that assumes no mass redistribution of meltwater, indicating the pervasive impact of melt processes.

The Greenland Ice Sheet has become an increasingly larger contributor to sea level rise in the past two decades and is projected to continue to lose mass. Climate variability is expected to increase under future warming, but the effect of climate variability on the Greenland Ice Sheet volume is poorly understood and is adding to the uncertainty of the projected mass loss. Here we quantify the influence of inter-annual temperature variability on mass loss from the Greenland Ice Sheet using the PISM model. We construct an ensemble of temperature-forcing fields that accounts for inter-annual variability in temperature using reanalysis data from NOAA-CIRES over the period 1851–2014. We investigate the steady-state and transient response of the Greenland Ice Sheet. We find that the simulated steady-state ice-sheet volume decreases by 1.9 ± 0.4 cm of sea level equivalent when forced with a varying temperature forcing compared to a constant temperature forcing, and by 11.5 ± 1.4 cm when the variability is doubled. The northern basins are particularly sensitive with a change in volume of 0.9–1.1%. Our results emphasize the importance of including temperature variability in projections of future mass loss.

The seasonal ice-free period in the Hudson Bay Complex (HBC) has grown longer in recent decades in response to warming, both from progressively earlier sea-ice retreat in summer and later sea-ice advance in fall. Such changes disrupt the HBC ecosystem and ice-based human activities. In this study, we compare 102 simulations from 37 models participating in phase 6 of the Coupled Model Intercomparison Project to the satellite passive microwave record and atmospheric reanalyses. We show that, throughout the HBC, models simulate an ice-free period that averages 30 d longer than in satellite observations. This occurs because seasonal sea-ice advance is unrealistically late and seasonal sea-ice retreat is unrealistically early. We find that much of the ice-season bias can be linked to a warm bias in the atmosphere that is associated with a southerly wind bias, especially in summer. Many models also exhibit an easterly wind bias during winter and spring, which reduces sea-ice convergence on the east side of Hudson Bay and impacts the spatial patterns of summer sea-ice retreat. These results suggest that, for many models, more realistic simulation of atmospheric circulation would improve their simulation of HBC sea ice.

The snow compression test is a snow stability test where an isolated column of snow is progressively loaded by tapping on it to induce failure in a possible weak layer. The test result provides valuable information about the propensity of failure initiation within the snowpack. However, different persons might tap with different force and thus reduce the reproducibility of the test results. The aim of this work was to quantify the influence of different test persons and different snowpacks on snow compression test results. We therefore let 62 persons tap on a stress measurement plate during a workshop of the European Avalanche Warning Services. Moreover, in the field, we performed stress measurements during 116 snow compression tests with 13 persons at eight different locations in the Alps. Data on persons’ body features and snow properties were also collected. Our results show that the stresses that reach a weak snow layer due to tapping are influenced by both the snowpack as well as different persons. Still, the data's scattering is surprisingly small for lower loading steps and decreases with depth. Therefore, we can deduce that, especially when avalanche conditions are particularly dangerous, snow compression test results are quite reproducible.

This study describes the morphological and dynamic changes of Parkachik Glacier, Suru River valley, Ladakh Himalaya, India. We used medium-resolution satellite images; CORONA KH-4, Landsat and Sentinel-2A from 1971–2021, and field surveys between 2015 and 2021. In addition, we used the laminar flow-based Himalayan Glacier Thickness Mapper and provide results for recent margin fluctuations, surface ice velocity, ice thickness, and identified glacier-bed overdeepenings. The results revealed that overall the glacier retreated by −210.5 ± 80 m with an average rate of 4 ± 1 m a−1 between 1971 and 2021. Whereas a field study suggested that the glacier retreat increased to −123 ± 72 m at an average rate of −20 ± 12 m a−1 between 2015 and 2021. Surface ice velocity was estimated using COSI-Corr on the Landsat data. Surface ice velocity in the lower ablation zone was 45 ± 2 m a−1 in 1999–2000 and 32 ± 1 m a−1 in 2020–2021, thus reduced by 28%. Further, the maximum thickness of the glacier is estimated to be ~441 m in the accumulation zone, while for glacier tongue it is ~44 m. The simulation results suggest that if the glacier continues to retreat at a similar rate, three lakes of different dimensions may form in subglacial overdeepenings.

Glide-snow avalanches release due to a loss of friction at the snow–ground interface, which can result in large avalanches that endanger infrastructure in alpine regions. The underlying processes are still relatively poorly understood, in part due to the limited data available on glide processes. Here, we introduce a pixel-based algorithm to detect glide cracks in time-lapse photographs under changing illumination and shadow conditions. The algorithm was applied to 14 years of time-lapse images at Dorfberg (Davos, Switzerland). We analysed 947 glide-snow events at a high-spatial (0.5 m) and temporal (2–15 min) resolution. Avalanche activity and glide-crack opening dynamics were investigated across timescales ranging from seasonally to hourly. Events were separated into surface (meltwater percolation) and interface events (no meltwater percolation). The results show that glide activity is highly variable between and within seasons. Most avalanches released without crack formation or within 24 h after crack opening, and release was favoured in the afternoon hours. Glide rates often showed a stick–slip behaviour. The acceleration of glide rates and non-constant increases in glide crack aspect ratio were indicators for avalanche release. This comprehensive dataset provides the basis for further investigations into glide-snow avalanche drivers.

Snow cover on sea ice poses a challenge for radar measurements as microwave penetration into snow is not yet fully understood. In this study, the aim is to investigate microwave penetration into snow on Arctic sea ice using commercial C (6 GHz) and K (26 GHz) band tank radars. Nadir-looking radar measurements collected at nine study locations over first-year and multiyear landfast sea ice in the Lincoln Sea in May 2018 are analysed together with detailed measurements of the physical properties of the snow cover to determine the dominant scattering horizons at both frequencies. They are evaluated for the feasibility to determine snow depth. The results show that in 39% of the measurements and only on first-year ice a major fraction of the C band radar backscatter originated closer to the snow–ice interface potentially enabling snow depth retrieval. At K band, 81% of the radar returns originated from the snow surface. Partly confirming the findings of previous studies, however, the analysis was potentially hampered by relatively warm air temperatures (up to $-0.9^\circ$C) during the study period as well as stratigraphic features and inconclusive microwave interaction with the saline basal layers found in the snow cover on first-year ice.

Aldegondabreen is a relatively small (5.3 km2) land-terminating glacier located in Nordeskiöld Land of Svalbard, ~10 km southwest of Barentsburg. Cumulative mass balance during 2002–20 equalled −21.79 m w.e., which corresponds to 37% of the total mass loss. The annual mass balance (Ba) varied from −0.24 to −2.19 m w.e., while the winter mass balance (Bw) ranged between 0.36 and 0.85 m w.e. Ba and Bw were strongly correlated with the positive degree-day sums and solid precipitation amounts, respectively, measured at the Barentsburg weather station. There was also a strong correlation (r = 0.76) between Ba and Bw, which shows that winter accumulation amplifies the consecutive summer ice melt by changing the surface albedo. The trend in both observational series is not detectable because the period from 2005 to 2013 was characterized by relatively high Bw and not very negative Ba values. This was also observed on the other Svalbard glaciers, and was related to prevailing north-westerly flows over Svalbard during the summer. Therefore, the decadal periodicity of the Aldegondabreen mass balance follows general archipelago patterns that are determined by regional-scale factors. Thus, the surface mass-balance time series, which is now the longest one in the central part of the Spitsbergen Island, is representative for the archipelago.

Extending the record of glacier area changes into the past improves our understanding of climate change impacts. Although analogue maps showing historic glacier extents are abundant, digital outlines from before the satellite era are sparse as the digitisation of moraines and trimlines on freely available satellite images is challenging. With the now available very high-resolution images provided by Web Map Services (WMS), new doors are open for the precise digitisation. Here, we used the ESRI WMS to digitise Little Ice Age (LIA) glacier extents and present area changes since the LIA in four selected regions along with a detailed uncertainty analysis. We used modern glacier outlines as a starting point and additionally consulted Sentinel-2 images, the ArcticDEM and historic maps for interpretation. Dating records from the literature allowed calculating area change rates. In total, 493 LIA glaciers (4640 km2, now 891 ice bodies with 3590 km2) were digitised, yielding relative area changes of −20% (−0.14% a−1), −15% (−0.10% a−1), −26% (−0.16% a−1) and −61% (−0.19% a−1) for Alaska, Baffin Island, Novaya Zemlya and the tropics, respectively. The ESRI WMS images are a great asset to precisely map moraines and trimlines, but information about the timing of the related extents requires further sources.

Dry-snow slab avalanches release due to crack propagation in a weak snow layer under a cohesive snow slab. Crack propagation speeds can provide insights into the potential size of avalanches and inform fracture and avalanche release models. Despite their importance, slope-scale crack speed measurements from real avalanches are limited. Further, most existing slope-scale measurements utilize the appearance of slab fractures on the snow surface. However, we have no evidence that the appearance of surface cracking is a good indicator of the weak layer crack propagation tip. Here we present a novel method to estimate crack propagation speed from snow surface movements in avalanche videos. Our technique uses changes in frame pixel intensity, allowing us to detect the location of weak layer cracks well before slab fractures appear on the snow surface. We use field experiments and numerical simulations to validate our method before applying it to five avalanches. Our estimates show that cracks propagate faster up and down the slope than in the cross-slope direction; this suggests that different propagation regimes likely govern crack propagation up/down the slope, cross-slope and in flat terrain.

Arctic rain-on-snow (ROS) events can have significant impacts on Arctic wildlife and socio-economic systems. This study addresses the meteorology of two different Arctic ROS events. The first, occurring near Nuuk, Greenland, generated significant impacts, including slush avalanches. The second, less severe, event occurred within the community of Iqaluit, Nunavut, Canada. This research utilizes atmospheric reanalysis, automated surface observation station data and atmospheric soundings to determine the meteorological conditions driving these events and the differences between each case. In both cases, atmospheric blocking played a leading role in ROS initiation, with atmospheric rivers – narrow bands of high water vapor transport, typically originating from the tropics and subtropics – having both direct and indirect effects. Cyclone-induced low-level jets and resultant ‘warm noses’ of higher air temperatures and moisture transport were other key features in ROS generation. To our knowledge, our study is the first to visualize how the varying strength and manifestation of these coupled features contribute to differences in the severity of Arctic ROS events. The meteorological drivers identified here find support from other studies on Arctic ROS events and are similar to weather features associated with Arctic precipitation events of extreme magnitude.

Using operational sea-ice maps, we provide first insight into the seasonal evolution of fast ice in the East Siberian Sea for the period between 1999 and 2021. The fast ice season tends to start later by 4.7 d per decade and to end earlier by 9.7 d per decade. As a result, there is a trend towards a shorter length of fast ice season by 2 weeks per decade. The analysis of air temperatures indicates that onset and end of the fast ice season are largely driven by thermodynamic processes. Two spatial modes (large, L-mode and small, S-mode) of East Siberian fast ice cover which have significant areal differences were distinguished. The occurrence of L- and S-modes was linked to the polarity of the Arctic Oscillation (AO) index. Negative AO phase leads to increased sea-ice convergence in the region, which in turn favours sea-ice grounding and promotes the development of large fast ice extent (L-mode). Lower deformation rates in the region during positive AO phase does not allow the formation of grounded features which results in small fast ice extent (S-mode). An analysis of sea-ice divergence confirms that L-mode seasons are characterised by higher on-shore convergence compared with S-mode seasons.

Over the last two decades, anomalous warming events have been observed in coastal Antarctic regions. While these events have been documented in the Ross Sea sector, the Antarctic interior is believed to have been buffered from warming. In this work, we present data from lakes located near Mt. Heekin and Thanksgiving Valley (~85° S) along the Shackleton Glacier, which are believed to be the southern-most Antarctic dry valley lakes. In 2018, the lakes were characterized, repeat satellite images were examined, and lake water chemistry was measured. Our analysis shows that lake areas recently increased, and the water-soluble ion chemistry indicates a flushing of salts from periglacial soils, likely from increased glacial melt as illustrated by water isotope data. Our results show that high southern latitude ice-free areas have likely been affected by warm pulses over the past 60 years and these pulses may be quasi-synchronous throughout the Transantarctic Mountains.

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs <4 km length and 60 days for larger icebergs. Complete shatter producing many icebergs <1 km length is rare. A side attrition rate of 0.23 m d−1 combined with drift speed of 6 km d−1, or any proportional change in both numbers fits the observed changes in iceberg distribution. The largest injection into the Southern Ocean of fresh water and any iceberg-transported material takes place in a ~2.3 × 10⁶ km2 zone extending east-northeast from the Antarctic Peninsula to the Greenwich meridian. The iceberg contribution to salinities and temperatures, with maximum contribution north of the Weddell Sea, differs in some regions, from those indicated by tracking large icebergs.

Due to adverse snow and cloud conditions, only a few inventories are available for the maritime glaciers in New Zealand. These are difficult to compare as different approaches and baseline data have been used to create them. In consequence, glacier fluctuations in New Zealand over the past two decades are only known for a few glaciers based on field observations. Here we present the results of a new inventory for the ‘year 2000’ (some scenes are from 2001 and 2002) that is based on glacier outlines from a recently published inventory for the year 2016 and allowed consistent change assessment for nearly 3000 glaciers over this period. The year 2000 inventory was created by manual on-screen digitizing using Landsat ETM+ satellite imagery (15 m panchromatic band) in the background and the year 2016 outlines as a starting point. Major challenges faced were late and early seasonal snow, clouds and shadow, the geo-location mismatch between Landsat and Sentinel-2 as well as the correct interpretation of ice patches and ice under debris cover. In total, we re-mapped 2967 glaciers covering an area of 885.5 km2 in 2000, which is 91.7 km2 (or 10.4%) more than the 793.8 km2 mapped in 2016. Area change rates (mean rate −0.65% a−1) increase towards smaller glaciers. Strongest area loss from 2000 to 2016 occurred at elevations ~1900 m but the highest relative loss was found below 800 m a.s.l. In total, 109 glaciers split into two or more entities and 264 had wasted away by 2016.