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

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.

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.

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.

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.

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.

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

Glaciological mass-balance measurements have been the backbone of internationally coordinated glacier monitoring. The resulting annual observations have been used to understand glacier reactions to climate change, and to assess both regional and global glacier mass changes and related contributions to runoff and sea-level rise. However, the comparability of annual observations is hampered by differences in survey periods and mass-balance amplitudes between glaciers, regions and hemispheres. This study presents a simple approach to temporally downscale glaciological mass balance using seasonal observations and sine functions. The proposed analytical model allows reconstruction of the seasonal course of glacier mass balance at weekly to monthly resolution from only annual or seasonal observations. Strengths and limitations of this analytical approach are discussed and compared with results from numerical mass-balance modelling. Potential applications include seasonal corrections of glaciological and geodetic observations and comparisons to monthly results from spaceborne gravimetry and altimetry.

Innovations in geodesy enable widespread analysis of glacier surface elevation change and geodetic mass balance. However, coincident glacier area data are less widely available, causing inconsistent handling of glacier area change. Here we quantify the bias introduced into meters water equivalent (m w.e.) specific geodetic mass balance results when using a fixed, maximum glacier area, and illustrate the bias for five North American glaciers. Sites span latitudes from the northern U.S. Rocky Mountains (48°N) to the Central Alaska Range (63°N) between 1948 and 2021. Results show that fixed (maximum) area treatment subdues the m w.e. mass change signal, underestimating mass balance by up to 19% in our test cases. This bias scales with relative glacier area change and the mass balance magnitude. Thus, the bias for specific geodetic mass balances will be most pronounced across rapidly deglaciating regions. Our analysis underscores the need for temporally resolved glacier area in geodetic mass balance studies.