Icebergs are part of the glacial mass balance and they interact with the ocean and with sea ice. Optical satellite remote sensing is often used to retrieve the above-waterline area of icebergs. However, varying solar angles introduce an error to the iceberg area retrieval that had not been quantified. Herein, we approximate the iceberg area error for top-of-atmosphere Sentinel-2 near-infrared data at a range of solar zenith angles. First, we calibrate an iceberg threshold at a $56^\circ$ solar zenith angle with reference to higher resolution airborne imagery at Storfjorden, Svalbard. A reflectance threshold of 0.12 yields the lowest relative error of 0.19% ± 15.74% and the lowest interquartile spread. Second, we apply the 0.12 reflectance threshold to Sentinel-2 data at 14 solar zenith angles between $45^\circ$ and $81^\circ$ in the Kangerlussuaq Fjord, south-east Greenland. Here we quantify the error variation with the solar zenith angle for a consistent set of large icebergs. The error variation is then standardized to the error obtained in Svalbard. Up to a solar zenith angle of $65^\circ$, the mean standardized iceberg area error remains between 5.9% and −5.67%. Above $65^\circ$, iceberg areas are underestimated and inconsistent, caused by a segregation into shadows and sun-facing slopes.

Synthetic Aperture Radar (SAR) has been used extensively to determine the surface ice flow velocity of tidewater glaciers and investigate changes in seasonal or annual ice dynamics at medium spatial resolution (⩾100 m). However, assessing tidewater glacier behaviour at these resolutions risks missing key details of glacier dynamics, which is particularly important for determination of strain rates that relate to crevasse formation, depth, and ice damage. Here we present surface ice velocity and strain maps with a 16 m posting derived from high-resolution (1 m) PAZ Ciencia spotlight mode SAR imagery for Narsap Sermia, SW Greenland, for October 2019 to February 2021. Results reveal fine details in strain rate, including an area of compression proximal to the terminus, with an upstream shift of strains through time. The velocity evolution of Narsap Sermia shows distinct seasonal changes starting in summer 2020, which are largely modulated by the subglacial drainage system. Comparison of our results with medium-resolution velocity products shows that while these can capture general strain and velocity patterns, our high-resolution data reveals considerably larger ranges of strain values. This is likely to have implications for tuning strain rate dependent calving and ice damage parameterisations within numerical models.

The collapse of several ice shelves in the Antarctic Peninsula since the late 20th century has resulted in the upstream acceleration of multiple formerly buttressed outlet glaciers, raising questions about the stability of Antarctica's remaining ice shelves and the effects their demise may have upon inland ice. Here, we use high temporal resolution Sentinel-1A/B synthetic aperture radar-derived observations to assess the velocity response of Larsen C Ice Shelf (LCIS) to the calving of colossal iceberg A-68 in 2017. We find marked oscillations in ice-shelf flow across LCIS in the months following A-68's calving, beginning with a near-ice-shelf-wide slowdown of 11.3 m yr−1 on average. While falling close to the limits of detectability, these ice-flow variations appear to have been presaged by similar oscillations in the years prior to A-68's breakaway, associated primarily with major rifting events, together reflecting potentially hitherto unobserved ice-shelf mechanical processes with important implications for ice-shelf weakening. Such ice-flow oscillations were, however, short-lived, with more recent observations suggesting a deceleration below longer-term rates of ice flow. Collectively, our observations reveal complex spatial-temporal patterns of ice-flow variability at LCIS. Similarly abrupt fluctuations may have important implications for the stability of other ice shelves, necessitating the continued, close observation of Antarctica's coastline in the future.

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.

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.

The hypothesis that ice-sheet evolution is only controlled by the long-term non-Newtonian viscous behavior of ice has been challenged by observations indicating that effects like brittle failure, stick-slip sliding, tides and wave action may affect ice-sheet evolution on sub-daily timescales. Over these timescales, the quasi-static-creep approximation is no longer appropriate and elastic effects become important. Simulating elastic effects in ice-sheet models over relevant timescales, however, remains challenging. Here, we show that by including a visco-elastic rheology and reintroducing the oft neglected acceleration term back into the ice-sheet stress balance, we can create a visco-elastic system where the velocity is locally determined and information propagates at the elastic wave speed. Crucially, the elastic wave speed can be treated like an adjustable parameter and set to any value to reproduce a range of phenomena, provided the wave speed is large compared to the viscous velocity. We illustrate the system using three examples. The first two examples demonstrate that the system converges to the steady-state viscous and elastic limits. The third example examines ice-shelf rifting and iceberg calving. This final example hints at the utility of the visco-elastic formulation in treating both long-term evolution and short-term environmental effects.

Rifts are full-thickness fractures that propagate laterally across an ice shelf. They cause ice-shelf weakening and calving of tabular icebergs, and control the initial size of calved icebergs. Here, we present a joint inverse and forward computational modeling framework to capture rifting by combining the vertically integrated momentum balance and anisotropic continuum damage mechanics formulations. We incorporate rift–flank boundary processes to investigate how the rift path is influenced by the pressure on rift–flank walls from seawater, contact between flanks, and ice mélange that may also transmit stress between flanks. To illustrate the viability of the framework, we simulate the final 2 years of rift propagation associated with the calving of tabular iceberg A68 in 2017. We find that the rift path can change with varying ice mélange conditions and the extent of contact between rift flanks. Combinations of parameters associated with slower rift widening rates yield simulated rift paths that best match observations. Our modeling framework lays the foundation for robust simulation of rifting and tabular calving processes, which can enable future studies on ice-sheet–climate interactions, and the effects of ice-shelf buttressing on land ice flow.

Frontal ablation, the combination of submarine melting and iceberg calving, changes the geometry of a glacier's terminus, influencing glacier dynamics, the fate of upwelling plumes and the distribution of submarine meltwater input into the ocean. Directly observing frontal ablation and terminus morphology below the waterline is difficult, however, limiting our understanding of these coupled ice–ocean processes. To investigate the evolution of a tidewater glacier's submarine terminus, we combine 3-D multibeam point clouds of the subsurface ice face at LeConte Glacier, Alaska, with concurrent observations of environmental conditions during three field campaigns between 2016 and 2018. We observe terminus morphology that was predominately overcut (52% in August 2016, 63% in May 2017 and 74% in September 2018), accompanied by high multibeam sonar-derived melt rates (4.84 m d−1 in 2016, 1.13 m d−1 in 2017 and 1.85 m d−1 in 2018). We find that periods of high subglacial discharge lead to localized undercut discharge outlets, but adjacent to these outlets the terminus maintains significantly overcut geometry, with an ice ramp that protrudes 75 m into the fjord in 2017 and 125 m in 2018. Our data challenge the assumption that tidewater glacier termini are largely undercut during periods of high submarine melting.

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.

Iceberg calving is one of the major mechanisms of ice loss from tidewater glaciers and ice sheets, but obtaining accurate estimates of ice discharge that are both continuous and accurate is a challenging task. Recent results have demonstrated the effective application of passive cryoacoustics – the use of naturally generated sounds to study the cryosphere – to quantify subaerial calving fluxes. However, little is known about the acoustic signatures of submarine calving. This study investigates the underwater noise from 656 subaerial and 162 submarine calving events observed at Hansbreen, Svalbard in the summers of 2016 and 2017. Statistical analysis of the acoustic signal shows that the normalized power of the calving noise is log-normally distributed regardless of the calving mode. However, submarine events can be distinguished from subaerial events by using the shape parameter of the log-normal distribution paired with the calving signal duration. The newly developed classification model may potentially be used for two purposes: (1) to study potential causal relationships between these two calving modes and (2) to separate calving fluxes into subaerial and submarine components. The latter will also require knowledge of the relationship between ice mass and sound spectral level for submarine calving events.

Tidewater glaciers have been observed to experience instantaneous, stepwise increases in velocity during iceberg-calving events due to a loss of resistive stresses. These changes in stress can potentially impact tidewater glacier stability by promoting additional calving and affecting the viscous delivery of ice to the terminus. Using flow models and perturbation theory, we demonstrate that calving events and subsequent terminus readvance produce quasi-periodic, sawtooth oscillations in stress that originate at the terminus and propagate upstream. The stress perturbations travel at speeds much greater than the glacier velocities and, for laterally resisted glaciers, rapidly decay within a few ice thickness of the terminus. Consequently, because terminus fluctuations due to individual calving events tend to be much higher frequency than climate variations, individual calving events have little direct impact on the viscous delivery of ice to the terminus. This suggests that the primary mechanism by which calving events can trigger instability is by causing fluctuations in stress that weaken the ice and lead to additional calving and sustained terminus retreat. Our results further demonstrate a stronger response to calving events in simulations that include the full stress tensor, highlighting the importance of accounting for higher order stresses when developing calving parameterizations.

Inaccurate representations of iceberg calving from ice shelves are a large source of uncertainty in mass-loss projections from the Antarctic ice sheet. Here, we address this limitation by implementing and testing a continuum damage-mechanics model in a continental scale ice-sheet model. The damage-mechanics formulation, based on a linear stability analysis and subsequent long-wavelength approximation of crevasses that evolve in a viscous medium, links damage evolution to climate forcing and the large-scale stresses within an ice shelf. We incorporate this model into the BISICLES ice-sheet model and test it by applying it to idealized (1) ice tongues, for which we present analytical solutions and (2) buttressed ice-shelf geometries. Our simulations show that the model reproduces the large disparity in lengths of ice shelves with geometries and melt rates broadly similar to those of four Antarctic ice shelves: Erebus Glacier Tongue (length ~ 13 km), the unembayed portion of Drygalski Ice Tongue (~ 65 km), the Amery Ice Shelf (~ 350 km) and the Ross Ice Shelf (~ 500 km). These results demonstrate that our simple continuum model holds promise for constraining realistic ice-shelf extents in large-scale ice-sheet models in a computationally tractable manner.

Iceberg calving, the process where icebergs detach from glaciers, remains poorly understood. Moreover, few parameterizations of the calving process can easily be integrated into numerical models to accurately capture observations, resulting in large uncertainties in projected sea level rise. Recent efforts have focused on estimating crevasse depths assuming tensile failure occurs when crevasses fully penetrate the glacier thickness. However, these approaches often ignore the role of advecting crevasses on calculations of crevasse depth. Here, we examine a more general crevasse depth calving model that includes crevasse advection. We apply this model to idealized prograde and retrograde bed geometries as well as a prograde geometry with a sill. Neglecting crevasse advection results in steady glacier advance and ice tongue formation for all ice temperatures, sliding law coefficients and constant slope bed geometries considered. In contrast, crevasse advection suppresses ice tongue formation and increases calving rates, leading to glacier retreat. Furthermore, crevasse advection allows a grounded calving front to stabilize on top of sills. These results suggest that crevasse advection can radically alter calving rates and hint that future parameterizations of fracture and failure need to more carefully consider the lifecycle of crevasses and the role this plays in the calving process.

Marine-terminating glaciers lose mass through melting and iceberg calving, and we find that meltwater drainage systems influence calving timing at Helheim Glacier, a tidewater glacier in East Greenland. Meltwater feeds a buoyant subglacial discharge plume at the terminus of Helheim Glacier, which rises along the glacial front and surfaces through the mélange. Here, we use high-resolution satellite and time-lapse imagery to observe the surface expression of this meltwater plume and how plume timing and location compare with that of calving and supraglacial meltwater pooling from 2011 to 2019. The plume consistently appeared at the central terminus even as the glacier advanced and retreated, fed by a well-established channelized drainage system with connections to supraglacial water. All full-thickness calving episodes, both tabular and non-tabular, were separated from the surfacing plume by either time or by space. We hypothesize that variability in subglacial hydrology and basal coupling drive this inverse relationship between subglacial discharge plumes and full-thickness calving. Surfacing plumes likely indicate a low-pressure subglacial drainage system and grounded terminus, while full-thickness calving occurrence reflects a terminus at or close to flotation. Our records of plume appearance and full-thickness calving therefore represent proxies for the grounding state of Helheim Glacier through time.

Iceberg calving strongly controls glacier mass loss, but the fracture processes leading to iceberg formation are poorly understood due to the stochastic nature of calving. The size distributions of icebergs produced during the calving process can yield information on the processes driving calving and also affect the timing, magnitude, and spatial distribution of ocean fresh water fluxes near glaciers and ice sheets. In this study, we apply fragmentation theory to describe key calving behaviours, based on observational and modelling data from Greenland and Antarctica. In both regions, iceberg calving is dominated by elastic-brittle fracture processes, where distributions contain both exponential and power law components describing large-scale uncorrelated fracture and correlated branching fracture, respectively. Other size distributions can also be observed. For Antarctic icebergs, distributions change from elastic-brittle type during ‘stable’ calving to one dominated by grinding or crushing during ice shelf disintegration events. In Greenland, we find that iceberg fragment size distributions evolve from an initial elastic-brittle type distribution near the calving front, into a steeper grinding/crushing-type power law along-fjord. These results provide an entirely new framework for understanding controls on iceberg calving and how calving may react to climate forcing.

Frontal ablation from tidewater glaciers is a major component of the total mass loss from the Greenland ice sheet. It remains unclear, however, how changes in atmospheric and oceanic temperatures translate into changes in frontal ablation, in part due to sparse observations at sufficiently high spatial and temporal resolution. We present high-frequency time-lapse imagery (photos every 30 min) of iceberg calving and meltwater plumes at Kangiata Nunaata Sermia (KNS), southwest Greenland, during June–October 2017, alongside satellite-derived ice velocities and modelled subglacial discharge. Early in the melt season, we infer a subglacial hydrological network with multiple outlets that would theoretically distribute discharge and enhance undercutting by submarine melt, an inference supported by our observations of terminus-wide calving during this period. During the melt season, we infer hydraulic evolution to a relatively more channelised subglacial drainage configuration, based on meltwater plume visibility indicating focused emergence of subglacial water; these observations coincide with a reduction in terminus-wide calving and transition to an incised planform terminus geometry. We suggest that temporal variations in subglacial discharge and near-terminus subglacial hydraulic efficiency exert considerable influence on calving and frontal ablation at KNS.

Retreat of calving glaciers worldwide has contributed substantially to sea-level rise in recent decades. Mass loss by calving contributes significantly to the uncertainty of sea-level rise projections. At Bowdoin Glacier, Northwest Greenland, most calving occurs by a few large events resulting from kilometre-scale fractures forming parallel to the calving front. High-resolution terrestrial radar interferometry data of such an event reveal that crevasse opening is fastest at low tide and accelerates during the final 36 h before calving. Using the ice flow model Elmer/Ice, we identify the crevasse water level as a key driver of modelled opening rates. Sea water-level variations in the range of local tidal amplitude (1 m) can reproduce observed opening rate fluctuations, provided crevasse water level is at least 4 m above the low-tide sea level. The accelerated opening rates within the final 36 h before calving can be modelled by additional meltwater input into the crevasse, enhanced ice cliff undercutting by submarine melt, ice damage increase due to tidal cyclic fatigue, crevasse deepening or a combination of these processes. Our results highlight the influence of surface meltwater and tides on crevasse opening leading to major calving events at grounded tidewater glaciers such as Bowdoin.

The Barents Sea offers a suitable location for documenting advection of warm and saline Atlantic Water (AW) into the Arctic and its impact on deglaciation and postglacial conditions. We investigate the timing, succession, and mechanisms of the transition from proximal glaciomarine to marine environment in the northwestern Barents Sea. Two studied sediment cores demonstrate diachronous retreat of the grounded ice sheet from the Kvitøya Trough (core S2528) to Erik Eriksen Trough (core S2519). Oxygen isotope records from core S2528 depict a two-step pattern, with lower δ18O values prior to the Younger Dryas (YD), and higher values afterward because of advection of the more saline, 18O-enriched AW. At this location, subsurface AW penetration increased during the Allerød and YD/Preboreal transition. In the study area, foraminiferal and dinocyst data from the YD interval suggest cold conditions, extensive sea-ice cover, and brine formation, along with the flow of chilled AW at subsurface and the development of a high-productivity polynya in the Erik Eriksen Trough. Dense winter sea-ice cover with seasonal productivity persisted in the Kvitøya Trough area during the early Holocene, whereas surface warming seems to have occurred during the middle Holocene interval.

Since the 2000s, Greenland ice sheet mass loss has been accelerating, followed by increasing numbers of glacial earthquakes (GEs) at near-grounded glaciers. GEs are caused by calving of km-scale icebergs which capsize against the terminus. Seismic record inversion allows a reconstruction of the history of GE sources which captures capsize dynamics through iceberg-to-terminus contact. When compared with a catalog of contact forces from an iceberg capsize model, seismic force history accurately computes calving volumes while the earthquake magnitude fails to uniquely characterize iceberg size, giving errors up to 1 km3. Calving determined from GEs recorded ateight glaciers in 1993–2013 accounts for up to 21% of the associated discharge and 6% of the Greenland mass loss. The proportion of discharge attributed to capsizing calving may be underestimated by at least 10% as numerous events could not be identified by standard seismic detections (Olsen and Nettles, 2018). While calving production tends to stabilize in East Greenland, Western glaciers have released more and larger icebergs since 2010 and have become major contributors to Greenland dynamic discharge. Production of GEs and calving behavior are controlled by glacier geometry with bigger icebergs being produced when the terminus advances in deepening water. We illustrate how GEs can help in partitioning and monitoring Greenland mass loss and characterizing capsize dynamics.

We investigate the appropriateness of calving or crevasse models from the literature using linear elastic fracture mechanics (LEFM). To this end, we compare LEFM model-predicted stress intensity factors (SIFs) against numerically computed SIFs using the displacement correlation method in conjunction with the finite element method. We present several benchmark simulations wherein we calculate the SIF at the tips of water-filled surface and basal crevasses penetrating through rectangular ice slabs under different boundary conditions, including grounded and floating conditions. Our simulation results indicate that the basal boundary condition significantly influences the SIF at the crevasse tips. We find that the existing calving models using LEFM are not generally accurate for evaluating SIFs in grounded glaciers or floating ice shelves. We also illustrate that using the ‘single edge crack’ weight function in the LEFM formulations may be appropriate for predicting calving from floating ice shelves, owing to the low fracture toughness of ice; whereas, using the ‘double edge crack’ or ‘central through crack’ weight functions is more appropriate for predicting calving from grounded glaciers. To conclude, we recommend using the displacement correlation method for SIF evaluation in real glaciers and ice shelves with complex geometries and boundary conditions.