We investigate the physical basis of the crevasse-depth (CD) calving law by analysing relationships between glaciological stresses and calving behaviour at Sermeq Kujalleq (Store Glacier), Greenland. Our observations and model simulations show that the glacier has a stable position defined by a compressive arch between lateral pinning points. Ice advance beyond the arch results in calving back to the stable position; conversely, if melt-undercutting forces the ice front behind the stable position, it readvances because ice velocities exceed subaqueous melt rates. This behaviour is typical of self-organising criticality, in which the stable ice-front position acts as an attractor between unstable super-critical and sub-critical regimes. This perspective provides strong support for a ‘position-law’ approach to modelling calving at Sermeq Kujalleq, because any calving ‘rate’ is simply a by-product of how quickly ice is delivered to the critical point. The CD calving law predicts ice-front position from the penetration of surface and basal crevasse fields, and accurately simulates super-critical calving back to the compressive arch and melt-driven calving into the sub-critical zone. The CD calving law reflects the glaciological controls on calving at Sermeq Kujalleq and exhibits considerable skill in simulating its mean position and seasonal fluctuations.

To better constrain meltwater transport and ice viscosity in temperate glaciers, particularly in ice stream shear margins, we use a custom permeameter to study the untested model relationship between the permeability of temperate ice and its liquid water content. The permeability of lab-made ice of two mean grain diameters (1.8 and 4.2 mm) is measured, and water content is controlled with the ice salinity and measured calorimetrically. Fluorescein dye is added to through-flowing, chilled water to highlight flow pathways through the ice after experiments. As predicted by a simple model, permeability increases with approximately the square of the water content and by about three orders of magnitude across water contents of 0.1–4.4%. However, permeability values are less than those of the model by average factors of 2.6 and 4.1 for the finer and coarser ice, respectively. This discrepancy is likely due to tortuous, truncated or air-clogged veins. The order-of-magnitude agreement between measured and modeled values may indicate that reduced permeability from these factors is nearly compensated by preferential flow in oversized veins that are isolated or arborescent. Both kinds of preferred flow pathways are observed but the latter only in fine-grained ice at water contents > 2%.

This study examines the performance of low-cost, low-power GNSS positioning systems for glacier monitoring in high-latitude environments. We compare the positioning performance of co-located low-cost u-blox ZED-F9P GNSS units (a few hundred USDs) and survey-grade Trimble R10 units (> $10,000 USD) under stationary (on land) and dynamic (on glacier) conditions near Terra Nova Bay, Antarctica. Low-cost and survey-grade systems yield almost identical error magnitudes under short (3 m), medium (34 km) and long (390 km) baseline kinematic-positioning scenarios. We further examined the efficacy of low-cost GNSS for glaciological applications by installing four u-blox and two Trimble receivers on Priestley Glacier to observe tide-modulated ice flexure. All receivers successfully detected subtle tidal oscillations with amplitudes < 3 cm, consistent with the predicted phasing from a tide model. These experiments offer a strong rationale for the widespread use of low-cost receivers to expand and densify GNSS monitoring networks, both in Antarctica and in glaciated regions worldwide.

Widespread glacier losses have been observed in most glaciated regions on Earth during recent decades, with a typical pattern of strong thinning in their lower reaches and limited elevation changes in their accumulation areas. Here, we use Pléiades satellite stereo-images of the Mont-Blanc massif (Alps) to reveal that thinning took place through the entire elevation range during the exceptional 2021/22 mass-balance year. Above 3000 m a.s.l. on Argentière glacier and Mer de Glace, thinning rates exceeded 3.5 m a−1 while almost no change occurred during the previous 9 years. Below 3000 m a.s.l., these anomalous thinning rates are essentially explained by changes in surface mass balance. At higher altitudes, other processes such as firn densification may play a role. Our analysis shows that high altitude glaciers, mostly stable during the last 100 years, are now responding to the impact of climate change.

Both submarine melt and calving are important for the overall mass balance of marine-terminating glaciers, but uncertainty is rife with regards to the magnitude of the processes. Modelling allows for these processes to be investigated without the need to visit inaccessible ice marginal zones. This study looks at the impact of different submarine melt and sea-ice back pressure scenarios on modelled calving activity and dynamics at Kronebreen, Svalbard, by running separate summer and winter simulations with various submarine melt parameterisations and sea-ice characteristics. It is found that submarine melt is an important driver of seasonal variation in modelled glacier dynamics and calving activity, with the choice of sliding law also exerting a significant influence on results.

Svalbard has experienced a dramatic increase in air temperature and glacier retreat since the end of the Little Ice Age. In many cases, this retreat has resulted in glaciers transitioning from being marine-terminating to land-terminating. Nordenskiöldbreen is an excellent contemporary example of this transition. A set of historical observations of glacier front positions was used to assess Nordenskiöldbreen's retreat rate and we found that the southern portion of the glacier front retreated by ~3500 m, since records began in 1896. The general retreat rate corresponds well with the air temperature trend during most of the 20th century. However, the average retreat rate has slowed since the 1990s despite increasing air temperatures. We show that this discrepancy between air temperature and retreat rate marks the transition from marine-terminating towards a land-terminating glacier, as the glacier's bedrock topography started to play an essential role in the glacier margin geometry, ice flow and retreat dynamics.

Semantic segmentation is a critical part of observation-driven research in glaciology. Using remote sensing to quantify how features change (e.g. glacier termini, supraglacial lakes, icebergs, crevasses) is particularly important in polar regions, where glaciological features may be spatially small but reflect important shifts in boundary conditions. In this study, we assess the utility of the Segment Anything Model (SAM), released by Meta AI Research, for cryosphere research. SAM is a foundational AI model that generates segmentation masks without additional training data. This is highly beneficial in polar science because pre-existing training data rarely exist. Widely-used conventional deep learning models such as UNet require tens of thousands of training labels to perform effectively. We show that the Segment Anything Model performs well for different features (icebergs, glacier termini, supra-glacial lakes, crevasses), in different environmental settings (open water, mélange, and sea ice), with different sensors (Sentinel-1, Sentinel-2, Planet, timelapse photographs) and different spatial resolutions. Due to the performance, versatility, and cross-platform adaptability of SAM, we conclude that it is a powerful and robust model for cryosphere research.

Recent research has highlighted the potential for high-resolution, high-density, 3D and 4D ground-penetrating radar (GPR) acquisitions on alpine glaciers. When carried out on foot, such surveys are laborious and time consuming, which limits their application to small domains of limited glaciological interest. Further, crevasses and other hazards make the data acquisition risky. To address these issues, we have developed a drone-based GPR system. The system has a payload weight of 2.2 kg and a data output rate of 14 traces per second. An 80-MHz antenna and a recording time of 2800 ns mean that depths of over 100 m can be reached in temperate ice. Differential GPS positioning assures accurate flight paths. At a speed of 4 m s−1 and height of 5 m above the glacier surface, our system can acquire over 4 line-km of GPR data in 20 min on a single set of drone batteries. After presenting the technical specifications of the system and tests required to optimize its performance, we showcase a recently acquired 3D dataset from the Otemma glacier in Switzerland, where 462 parallel GPR profiles were surveyed at a 1-m line spacing, totaling over 112 line-km of data, in only 4 days.

Observations remain sparse for peripheral glaciers and ice caps in Greenland. Here, we present the results of a multi-frequency radar survey of Lyngmarksbræen Ice Cap in West Greenland conducted in April 2017. Radar measurements show thick ice of up to ~120 m in subglacial valleys associated with the largest outlet glaciers, while relatively thin ice cover the upper plateau ice divides, suggesting future vulnerability to ice cap fragmentation. At the time of the radar survey, Lyngmarksbræen Ice Cap had a total volume of 0.82 ± 0.1 km3. Measurements show a 1.5–2 m thick end-of-winter snowpack, and that firn is largely absent, signifying a prolonged period of negative mass balance for most of the ice cap. The thermal regime of Lyngmarksbræen Ice Cap is investigated through analysis of scattering observed along radar profiles. Results show that the ice cap is largely below the pressure melting point, but that temperate ice exists both in deep basal pockets and in shallow zones that some places extend from ~15 m depth and to the ice base. The distribution of shallow temperate ice appears unrelated to variations in ice thickness; instead we find a strong correlation to the presence of nearby surface crevasses.

Stable water isotope records of six firn cores retrieved from two adjacent plateaus on the northern Antarctic Peninsula between 2014 and 2016 are presented and investigated for their connections with firn-core glacio-chemical data, meteorological records and modelling results. Average annual accumulation rates of 2500 kg m−2 a−1 largely reduce the modification of isotopic signals in the snowpack by post-depositional processes, allowing excellent signal preservation in space and time. Comparison of firn-core and ECHAM6-wiso modelled δ18O and d-excess records reveals a large agreement on annual and sub-annual scales, suggesting firn-core stable water isotopes to be representative of specific synoptic situations. The six firn cores exhibit highly similar isotopic patterns in the overlapping period (2013), which seem to be related to temporal changes in moisture sources rather than local near-surface air temperatures. Backward trajectories calculated with the HYSPLIT model suggest that prominent δ18O minima in 2013 associated with elevated sea salt concentrations are related to long-range moisture transport dominated by westerly winds during positive SAM phases. In contrast, a broad δ18O maximum in the same year accompanied by increased concentrations of black carbon and mineral dust corresponds to the advection of more locally derived moisture with northerly flow components (South America) when the SAM is negative.

Subglacial lake water-volume changes produce ice-elevation anomalies that provide clues about water flow beneath glaciers and ice sheets. Significant challenges remain in the quantitative interpretation of these elevation-change anomalies because the surface expression of subglacial lake activity depends on basal conditions, rate of water-volume change, and ice rheology. To address these challenges, we introduce an inverse method that reconstructs subglacial lake activity from altimetry data while accounting for the effects of viscous ice flow. We use a linearized approximation of a Stokes ice-flow model under the assumption that subglacial lake activity only induces small perturbations relative to a reference ice-flow state. We validate this assumption by accurately reconstructing lake activity from synthetic data that are produced with a fully nonlinear model. We then apply the method to estimate the water-volume changes of several active subglacial lakes in Antarctica by inverting data from NASA's Ice, Cloud, and land Elevation Satellite 2 (ICESat-2) laser altimetry mission. The results show that there can be substantial discrepancies (20% or more) between the inversion and traditional estimation methods due to the effects of viscous ice flow. The inverse method will help refine estimates of subglacial water transport and further constrain the role of subglacial hydrology in ice-sheet evolution.

Ice-crystal fabric can induce mechanical anisotropy that significantly affects flow, but ice-flow models generally do not include fabric development or its effect upon flow. Here, we incorporate a new spectral expansion of fabric, and more complete description of its evolution, into the ice-flow model Elmer/Ice. This approach allows us to model the effect of both lattice rotation and migration recrystallization on large-scale ice flow. The fabric evolution is coupled to flow using an unapproximated non-linear orthotropic rheology that better describes deformation when the stress and fabric states are misaligned. These improvements are most relevant for simulating dynamically interesting areas, where recrystallization can be important, tuning data are scarce and rapid flow can lead to misalignment between stress and fabric. We validate the model by comparing simulated fabric to ice-core and phase-sensitive radar measurements on a transect across Dome C, East Antarctica. With appropriately tuned rates for recrystallization, the model is able to reproduce observations of fabric. However, these tuned rates differ from those previously derived from laboratory experiments, suggesting a need to better understand how recrystallization acts differently in the laboratory compared to natural settings.

The complex topography and size of High Mountain Asia (HMA) result in large differences in glacier mass-balance variability and climate sensitivity. Current understanding of these sensitivities is limited by simplifications in past studies’ model structure. This study overcomes this limitation by using a mass-balance model to investigate the climatic mass-balance variability and climate sensitivity of 16 glaciers covering major mountain ranges in HMA. Generally, glaciers in the southeast have higher mass turnover while glaciers at the margins of HMA show higher interannual mass-balance variability. All glaciers are most sensitive to temperature perturbations in summer. The climatic mass balance of 15 glaciers is most sensitive to precipitation perturbations in summer or spring and summer, even if the seasonal accumulation peak is not in summer. Only one glacier's mass balance (Chhota Shigri Glacier) is most sensitive to precipitation perturbations in winter. Glaciers with high mass turnover and high summer-precipitation ratio are more sensitive to temperature perturbations. Sensitivity experiments reveal that besides the non-linearity of mass-balance temperature sensitivity, mass-balance precipitation sensitivity is non-linear as well. Furthermore, resolving the diurnal cycle of albedo, (re)freezing and the differentiation between liquid and solid precipitation are important to assess climate sensitivity of glaciers in HMA.

Converting measurements of ice-sheet surface elevation change to mass change requires measurements of accumulation and knowledge of the evolution of the density profile in the firn. Most firn-densification models are tuned using measured depth–density profiles, a method which is based on an assumption that the density profile in the firn is invariant through time. Here we present continuous measurements of firn-compaction rates in 12 boreholes near the South Pole over a 2 year period. To our knowledge, these are the first continuous measurements of firn compaction on the Antarctic plateau. We use the data to derive a new firn-densification algorithm framed as a constitutive relationship. We also compare our measurements to compaction rates predicted by several existing firn-densification models. Results indicate that an activation energy of 60 kJ mol−1, a value within the range used by current models, best predicts the seasonal cycle in compaction rates on the Antarctic plateau. Our results suggest models can predict firn-compaction rates with at best 7% uncertainty and cumulative firn compaction on a 2 year timescale with at best 8% uncertainty.

Spatial variability in bed topography, characterized as bed roughness, impacts ice-sheet flow and organization and can be used to infer subglacial conditions and processes, yet is difficult to quantify due to sparse observations. Paleo-subglacial beds of formerly expanded glaciers found across the Antarctic continental shelf are well preserved, have relatively limited post-glacial sediment cover and contain glacial landforms that can be resolved at sub-meter vertical scales. We analyze high-resolution bathymetry offshore of Pine Island and Thwaites glaciers in the Amundsen Sea to explore spatial variability of bed roughness where streamlined subglacial landforms allow for the determination of ice-flow direction. We quantify bed roughness using std dev. and Fast Fourier Transform methods, each employed at local (100 km) and regional (101–2 km) scales and in along- and across-flow orientations to determine roughness expressions across spatial scales. We find that the magnitude of roughness is impacted by the parameters selected – which are often not sufficiently reported in studies – to quantify roughness. Important spatial patterns can be discerned from high-resolution bathymetry, highlighting both its usefulness in identifying patterns of streaming ice flow and underscores the need for a standardized way of characterizing topographic variability.

A major consideration for maritime activity in the Southern Hemisphere is the northern limit of icebergs, or the Southern Ocean Limit Of Known Ice (SOLOKI). This analysis of historical reports of icebergs during Southern Hemisphere voyages from 1687 to 1933 provides a basis for examination of their geographical and chronological occurrence during ~250 years. The analyses use tabulated data from 742 voyages and other reports from many sources, some including first-person descriptions. While these data are dependent on icebergs being reported by mariners, as well as the variable frequency of voyages, they demonstrate distinct periods of exceptional frequency of icebergs occurring in certain localities, particularly the far South Atlantic. Based upon historical records the evidence suggests unprecedented numbers of icebergs were present in southern shipping channels in the 1890s. When these historical observations are combined with modern iceberg drift trajectories, their possible origin can be elucidated. Owing to the numbers of icebergs seen and their geographical spread, our results suggest that this was possibly the largest near-synchronous calvings in the last 300 years, and the northernmost extent of the SOLOKI.

We present seismic measurements of the firn column at Korff Ice Rise, West Antarctica, including measurements of compressional-wave velocity and attenuation. We describe a modified spectral-ratio method of measuring the seismic quality factor (Q) based on analysis of diving waves, which, combined with a stochastic method of error propagation, enables us to characterise the attenuative structure of firn in greater detail than has previously been possible. Q increases from 56 ± 23 in the uppermost 12 m to 570 ± 450 between 55 and 77 m depth. We corroborate our method with consistent measurements obtained via primary reflection, multiple, source ghost, and critically refracted waves. Using the primary reflection and its ghost, we find Q = 53 ± 20 in the uppermost 20 m of firn. From the critical refraction, we find Q = 640 ± 400 at 90 m depth. Our method aids the understanding of the seismic structure of firn and benefits characterisation of deeper glaciological targets, providing an alternative means of correcting seismic reflection amplitudes in cases where conventional methods of Q correction may be impossible.

We recently reported on the radio-frequency attenuation length of cold polar ice at Summit Station, Greenland, based on bi-static radar measurements of radio-frequency bedrock echo strengths taken during the summer of 2021. Those data also allow studies of (a) the relative contributions of coherent (such as discrete internal conducting layers with sub-centimeter transverse scale) vs incoherent (e.g. bulk volumetric) scattering, (b) the magnitude of internal layer reflection coefficients, (c) limits on signal propagation velocity asymmetries (‘birefringence’) and (d) limits on signal dispersion in-ice over a bandwidth of ~100 MHz. We find that (1) attenuation lengths approach 1 km in our band, (2) after averaging 10 000 echo triggers, reflected signals observable over the thermal floor (to depths of ~1500 m) are consistent with being entirely coherent, (3) internal layer reflectivities are ≈–60$\to$–70 dB, (4) birefringent effects for vertically propagating signals are smaller by an order of magnitude relative to South Pole and (5) within our experimental limits, glacial ice is non-dispersive over the frequency band relevant for neutrino detection experiments.

From 2017 to 2020, three significant calving events took place on Pine Island Glacier, West Antarctica. Ice-shelf velocities changed over this period and the calving events have been suggested as possible drivers. However, satellite observations also show significant changes in the areal extent of fracture zones, especially in the marginal areas responsible for providing lateral support to the ice shelf. Here, we conduct a model study to identify and quantify drivers of recent ice-flow changes of the Pine Island Ice Shelf. In agreement with recent studies, we find that the calving events caused significant velocity changes over the ice shelf. However, calving alone cannot explain observed velocity changes. Changes in the structural rigidity, i.e. ice damage, further significantly impacted ice flow. We suggest that ice damage evolution of the ice-shelf margins may have influenced recent calving events, and these two processes are linked.