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.

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.

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.

The spatiotemporal distribution of snow affects hydrological and climatological processes at different scales. Accordingly, quantifying geometric features of snow-cover patterns is important, providing a valuable complement for snow water equivalent (SWE) modelling. This study on satellite-based morphological analysis originally uses two types of geometric indexes: (1) MN, Minkowski numbers (area (MN1), perimeter (MN2), Euler number (MN3)), and (2) CL, average chord length, to describe the morphology of Sentinel-2-derived snow-covered areas (SCAs), within the high-alpine site Zugspitze for a 5 year period. Results indicate that they capture the seasonal variability of snow-cover patterns, particularly during accumulation and ablation. Being to some degree independent from each other, MN2, MN3 and CL provide additional information upon shape, connectivity and length scale of snow cover, compared to most used indexes (e.g. fractional SCA). Correlation values up to +0.7 for MN2, +0.58 for MN3 and +0.46 for CL were observed with selected topographic characteristics, suggesting a close connection between geometric features of snow cover and ground features. Comparing in situ SWE measurements with MN and CL shows a correlation between −0.5 and +0.5. These indexes can hence be applied in combination with in situ data and/or modelling approaches to improve spatially distributed SWE in high-alpine catchments.

Liquid water at the ground–snow interface is thought to play a crucial role in the release of glide-snow avalanches, which can be massive and threaten infrastructure in alpine regions. Several mechanisms have been postulated to explain the formation of this interfacial water. However, these mechanisms remain poorly understood, in part because suitable measurement techniques are lacking. Here, we demonstrate the use of neutron radiography for imaging water transport in soil–snow systems. Columns of sand, gravel and snow were used to simulate the capillary forces of the soil–vegetation–snow layering found in nature. The columns were connected to a water reservoir to maintain a constant-pressure boundary condition and placed in a climatic chamber within the neutron beam. We show that neutron radiography is capable of measuring changes in the optical density distribution (related to liquid water content) within all three layers of the model system. Results suggest that a porous interface between the sand and snow may induce the formation of a water layer in the basal snowpack. Improved understanding of the water transport in soil–snow systems should lead to better prediction of glide-snow avalanche release and could also benefit other fields such as snow hydrology.

Understanding the dynamics of snow avalanches is crucial for predicting their destructive potential and mobility. To gain insight into avalanche dynamics at a particle level, the AvaNode in-flow sensor system was developed. These synthetic particles, equipped with advanced and affordable sensors such as an inertial measurement unit (IMU) and global navigation satellite system (GNSS), travel with the avalanche flow. This study focuses on assessing the feasibility of the in-flow measurement systems. The experiments were conducted during the winter seasons of 2021–2023, both in static snow cover and dynamic avalanche conditions of medium-sized events. Radar measurements were used in conjunction with the particle trajectories and velocities to understand the behaviour of the entire avalanche flow. The dynamic avalanche experiments allowed to identify three distinct particle flow states: (I) initial rapid acceleration, (II) a steady state flow with the highest velocities (9–17 ms−1), and (III) a longer deceleration state accompanied by the largest measured rotation rates. The particles tend to travel towards the tail of the avalanche and reach lower velocities compared to the frontal approach velocities deduced from radar measurements (ranging between 23–28 ms−1). The presented data give a first insight in avalanche particle measurements.

We developed a Handheld Integrating Sphere Snow Grain Sizer (HISSGraS) for field use to measure the specific surface area (SSA) of snow. In addition to snow samples, HISSGraS can directly measure snow surfaces and snow pit walls. The basic measurement principle is the same as that of the IceCube SSA instrument. The retrieval algorithm for SSA from reflectance employs two conversion equations formulated using spherical and nonspherical grain shape models. We observed SSAs using HISSGraS, IceCube and the gas adsorption method in a snowfield in Hokkaido, Japan. Intercomparison of the results confirmed that with HISSGraS direct measurement, SSA profile observations can be completed in just ~1/10 the time required for measurement of snow samples. Our results also suggest that HISSGraS and IceCube have similar accuracy when the same snow samples are measured using the same grain shape model. However, SSAs of near-surface snow layers measured using the three techniques exhibited some biases, possibly due to rapid snow metamorphism or melting during measurement and some technical issues with optical techniques. When excluding SSA data for the surface layer, which metamorphosed remarkably during measurement, IceCube- and HISSGraS-derived SSAs correlated strongly with those obtained by gas adsorption and HISSGraS accuracy is 21–34%.

The relationship between the behavior of water in snow and its microstructure is crucial to improve the prediction of wet snow disasters. X-ray computed tomography (X-ray CT) is frequently used to observe snow microscopically. However, distinguishing between ice and water in the X-ray images is difficult because ice exhibits an X-ray absorption coefficient similar to that of water. In contrast, magnetic resonance imaging (MRI) acquires nuclear magnetic resonance (NMR) signals of protons in a liquid and visualizes the NMR signal intensity, enabling discrimination between water and ice signals. However, snow grains and pore spaces cannot be distinguished in MRI images because they do not generate NMR signals.

To investigate the relationship between the microstructure of snow and the distribution of liquids in snow, we developed a novel method that combines X-ray CT and MRI images to compensate for the disadvantages associated with each method. Using this method, we successfully visualized where liquid (C12H24) occupied pore spaces. We also showed the possibility of using C12H24 instead of water to obtain water retention curve of snow cover, which is a fundamental aspect of hydraulic properties. Although there is room for improvement in the visualization of water in snow, such as shortening the imaging time to escape snow metamorphosis and image superimposition methods, this method is expected to effectively elucidate the behavior of water in snow and clarify the characteristics of wet snow.

Spatial variability of snowpack properties adds uncertainty in the evaluation of avalanche hazard. We propose a combined mechanical–statistical approach to study how spatial variation of slab depth affects the skier-triggering probability and possible release size. First, we generate multiple slab depth maps on a plane fictional slope based on Gaussian Random Fields (GRF) for a specific set of mean, variance and correlation length. For each GRF, we derive analytically the Skier Propagation Index (SPI). We then simulate multiple skier tracks and computed the probability based on the number of skier hits where SPI is below 1. Finally, we use a depth-averaged material point method to evaluate the possible avalanche size for given slab depth variations. The results of this analysis show that large correlation lengths and small variances lead to a lower probability of skier-triggering as it reduces the size and the number of areas with low slab depth. Then, we show the effect of skiing style and skier group size on skier-triggering probability. Spatial variability also affects the possible avalanche size by adding stress fluctuation causing early or late tensile failure. Finally, we demonstrate with our models the well-known relationship between the probability and the size in avalanche forecasting.

In the course of an artificially triggered avalanche, a particle tracking procedure is combined with supplementary measurements, including Global Navigation Satellite System (GNSS) positioning, terrestrial laser scanning and Doppler radar measurements. Specifically, an intertial measurement unit is mounted inside a rigid sphere, which is placed in the avalanche track. The sphere is entrained by the moving snow, recording translational accelerations, angular velocities and the flux density of Earth's magnetic field. Based on the recorded data, we present a threefold analysis: (i) a qualitative data interpretation, identifying different particle motion phases which are associated with corresponding flow regimes, (ii) a quantitative time integration algorithm, determining the corresponding particle trajectory and associated velocities on the basis of standard sensor calibration, and (iii) an improved quantitative evaluation relying on a novel in situ sensor calibration technique, which is motivated by the limitations of the given dataset. The final results, i.e. the evolution of the angular orientation of the sensor unit, translational and rotational velocities and estimates of the sensor trajectory, are assessed with respect to their reliability and relevance for avalanche dynamics as well as for future design of experiments.

Passive-microwave sea ice concentration (SIC) algorithms employ different frequencies and polarisations in their operational implementations. Commonly, these algorithms utilise combinations such as 19/37 GHz, yielding reduced measurement uncertainties but at a coarse spatial resolution. Alternatively, these algorithms can solely use 89 GHz, producing a higher spatial resolution but with increased measurement uncertainties. This study evaluates the application of a resolution-enhancing SIC algorithm (reSICCI3LF), initially developed for the coarser Special Sensor Microwave Imager / Sounder, on the Advanced Microwave Scanning Radiometer. By applying reSICCI3LF, we aim to produce a 5 km SIC for 2013–2020 in the Fram Strait and the Barents and Kara Sea region that gains the benefits of both types of algorithms, high spatial resolution and low measurement uncertainty.

We present the algorithm tuning, spectral analysis of spatial resolutions, and validation against the Round Robin Data Package of 0% and 100% SIC points and SIC derived from Landsat-8. The findings demonstrate that the reSICCI3LF algorithm produces a SIC field with fine details, achieving a balance between high spatial resolution and lower measurement uncertainties compared to a 89 GHz based SIC. Consequently, this resolution-enhanced SIC technique can potentially initialise higher-resolution coupled ocean and sea ice forecasting systems through data assimilation.

For snow-avalanche hazard mapping, one needs efficient tools that nevertheless capture the essential physical processes. The code MoT-PSA (Method of Transport – Powder Snow Avalanche) described here is based on the two-layer depth-averaged formulation for mixed snow avalanches developed by Eglit and co-workers in the 1980s but is extended to 3-D terrain and uses a fast numerical scheme based on the method of transport. Compared to previous works, we introduce novel formulations for the suspension and deposition of snow from the dense core. Snow cover and air entrainment are quantified with physics-based models. A sensitivity study of the model parameters on an idealized topography shows that both the dense core and the parameters of the powder snow cloud (PSC) governing particle suspension and settling significantly affect the dynamics. As expected, we observe that snow cover entrainment favours the formation of large PSCs with long runout. The powder-snow avalanche that occurred in Lom (Norway) on 27 February 2020 is back-calculated using MoT-PSA. With plausible parameter values, the model reproduces the dense core stopping at the gully's base and the dilute PSC travelling across the frozen lake for almost 1 km.

Subglacial seismicity presents the opportunity to monitor inaccessible glacial beds at the epicentral location and time. Glaciers can be underlain by rock or till, a first order control on bed mechanics. Velocity-weakening, necessary for unstable slip, has been shown for each bed type, but is much stronger and evolves over more than an order of magnitude longer distances for till beds. Utilizing a de-stiffened double direct shear apparatus, we found conditions for instability at freezing temperatures and high slip rates for both bed types. During stick–slip stress-drops, we recorded acoustic emissions with piezoelectric transducers frozen into the ice. The two populations of event waveforms appear visually similar and overlap in their statistical features. We implemented a suite of supervised machine learning algorithms to classify the bed type of recorded waveforms and spectra, with prediction accuracy between 65–80%. The Random Forest Classifier is interpretable, showing the importance of initial oscillation peaks and higher frequency energy. Till beds have generally higher friction and resulting stress-drops, with more impulsive first arrivals and more high frequency content compared to rock emissions, but rock beds can produce many till-like events. Seismic signatures could enhance interpretation of bed conditions and mechanics from subglacial seismicity.

We review successes and challenges from five recent subglacial bedrock drilling campaigns intended to find evidence for Antarctic Ice Sheet retreat during warm periods in the geologic past. Insights into times when the polar ice sheets were smaller than present serve as guiding information for modeling efforts that aim to predict the rate and magnitude of future sea level rise that would accompany major retreat of the Antarctic Ice Sheet. One method to provide direct evidence for the timing of deglaciations and minimum extent of prior ice sheets is to extract subglacial bedrock cores for cosmogenic nuclide analysis from beneath the modern ice sheet surface. Here we summarize the lessons learned from five field seasons tasked with obtaining bedrock cores from shallow depths (<120 m beneath ice surface) across West Antarctica since 2016. We focus our findings on drilling efforts and technology and geophysical surveys with ground-penetrating radar. Shallow subglacial drilling provides a high risk, high reward means to test for past instabilities of the Antarctic Ice Sheet, and we highlight key challenges and solutions to increase the likelihood of success for future subglacial drilling efforts in polar regions.

Avalanche modeling is an essential tool to assess snow avalanche hazard. Today, most popular numerical approaches adopt depth-averaged equations. These methods are computationally efficient but limited in capturing processes occurring in the flow depth direction, e.g. erosion or deposition, which are often considered using ad hoc parameterizations or neglected completely. However, processes such as snow erosion, can crucially influence flow dynamics and run-out and are often not negligible. We address these issues by using a new three-dimensional (3-D) model, based on the material point method and finite strain elastoplasticity. To assess the possibilities and challenges associated with these highly detailed but computationally expensive calculations, we simulated the ‘Salezer’ snow avalanche that released in Davos, Switzerland in 2019. To reproduce the event in our simulations, we use the release areas mapped in a photogrammetric drone survey and estimate the snow conditions on the day of the event. We compare macroscopic features, such as flow outline and snow deposition of the simulated avalanche to field observations. An in-depth analysis of transient 3-D flow structures at the avalanche head not only demonstrates the degree of physical detail in the model, but also highlights challenges which still need to be addressed.

This study investigates the impacts of the timing of an extreme cyclone that occurred in August 2012 on the sea-ice volume evolution based on the Arctic Ice Ocean Prediction System (ArcIOPS). By applying a novel cyclone removal algorithm to the atmospheric forcing during 4–12 August 2012, we superimpose the derived cyclone component onto the atmospheric forcing one month later or earlier. This study finds that although the extreme cyclone leads to strong sea-ice volume loss in all runs, large divergence occurs in sea-ice melting mechanism in response to various timing of the cyclone. The extreme cyclone occurred in August, when enhanced ice volume loss is attributed to ice bottom melt primarily and ice surface melt secondarily. If the cyclone occurs one month earlier, ice surface melt dominates ice volume loss, and earlier appearance of open water within the ice zone initiates positive ice-albedo feedback, leading to a long lasting of the cyclone-induced impacts for approximately one month, and eventually a lower September ice volume. In contrast, if the cyclone occurs one month later, ice bottom melt entirely dominates ice volume loss, and the air-open water heat flux in the ice zone tends to offset ice volume loss.