Abstract Deciphering the localisation of solid and dissolved impurities on the micron-scale in glacial ice remains a challenge, but is critical to understand the integrity of ice core records and internal deformation. Here we report on the state-of-the-art in microstructural impurity research by highlighting recent progress in bringing together cryo-Raman spectroscopy and laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS). We show the potential of both methods and discuss possibilities to improve inter-method approaches aiming for a more holistic understanding of the evolution of impurity localisation throughout the ice column, including post-depositional processes. In this framework, we elaborate on future research priorities such as LA-ICP-MS imaging on firn samples and integrating a large cryo-cell with imaging capabilities.

Abstract 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.

Abstract 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.

Abstract 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.

Abstract Icebergs in proglacial fjords serve as pupping, resting and molting habitat for some of the largest seasonal aggregations of harbor seals (Phoca vitulina richardii) in Alaska. One of the largest aggregations in Southeast Alaska occurs in Johns Hopkins Inlet, Glacier Bay National Park, where up to 2000 seals use icebergs produced by Johns Hopkins Glacier. Like other advancing tidewater glaciers, the advance of Johns Hopkins Glacier over the past century has been facilitated by the growth and continual redistribution of a submarine end moraine, which has limited mass losses from iceberg calving and submarine melting and enabled glacier thickening by providing flow resistance. A 15-year record of aerial surveys reveals (i) a decline in iceberg concentrations concurrent with moraine growth and (ii) that the iceberg size distributions can be approximated as power law distributions, with relatively little variability and no clear trends in the power law exponent despite large changes in ice fluxes over seasonal and interannual timescales. Together, these observations suggest that sustained tidewater glacier advance should typically be associated with reductions in the number of large, habitable icebergs, which may have implications for harbor seals relying on iceberg habitat for critical life-history events.

Abstract Snow ice (SI) forms from freezing wet snow, known as slush, and contributes to the thickness of level and brash ice. However, the mechanism of snow-slush-snow ice transformation has not been extensively investigated to date, despite the difference in the freezing rate of slush in comparison with water is important for estimating the ice thickness. In this study, we examined the growth of initial congelation ice (CI) and snow ice (SI) in a fresh water tank exposed to outdoor weather conditions in Luleå, northern Sweden. The tank of size 1.8 × 0.65 × 1.2 m in length, width and height was divided into two compartments to facilitate the simultaneous growth of CI and SI. A total of 12 experiments were conducted in the years 2021 and 2022. The transformation from slush to snow ice was achieved by submerging various amounts of snow in the compartments. It was observed that approximately 35% of the initial snow transformed into SI. Snow ice grew 4 mm°C−0.5 d−0.5 faster than congelation ice. The CI growth under SI was 1 mm°C−0.5 d−0.5 slower than the CI growth under CI. This study provides valuable insights for modelling snow-slush-snow ice transformation and designing future laboratory-scale experiments.

Abstract 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.

Abstract 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.

Abstract 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.