We present the first simulations of Jakobshavn Isbrae (Sermeq Kujalleq), west Greenland, using a 3D Stokes calving model that permits unrestricted advance and retreat. Using the position-based crevasse-depth calving law, the model is applied to simulate the calving dynamics of 2016–2017 season when Jakobshavn Isbrae is assumed to be stable because of the presence of a strong proglacial ice mélange. The calving law needs to be adjusted to avoid an underestimation of calving, but once adjusted the calving model simulates seasonal calving dynamics that reflect observed calving-driven retreat very well. We find that a crevasse penetration threshold of 94.5% best matches observations from satellite imagery. Additional, 2-year transient simulations show that although ice mélange is essential to the glacier's winter readvance, when removed, the glacier only retreats a couple of kilometres before reaching a stable position. While the backstress provided by the ice mélange allows the glacier to advance beyond this point, the retreated terminus position is determined by a combination of bed geometry and glacier dynamics. Ultimately, while the ice mélange allows winter readvance, cessation of the well-documented rapid retreat of Jakobshavn Isbrae will be influenced by the bed geometry.

Surface meltwater can influence subglacial hydrology and ice dynamics if it reaches ice sheet's base. Firn aquifers store meltwater and drain into wide crevasses marking the aquifer's downstream boundary, indicating water from firn aquifers can drive hydrofracture to establish surface-to-bed hydraulic connections at inland locations. Yet, sparse observations limit our understanding of the physical processes controlling firn aquifer drainage. We assess the potential for future inland firn aquifer drainage migration with field observations and linear elastic fracture mechanics (LEFMs) modeling to determine the conditions needed to initiate and sustain hydrofracture on Helheim Glacier, Greenland. We find that local stress conditions alone can drive crevasse tips into the firn aquifer, allowing hydrofracture initiation year-round. We infer inland expansion of crevasses over the firn aquifer from crevasse-nucleated whaleback dune formation and Global Navigation Satellite System-station detected crevasse opening extending 14 and 4 km, respectively, inland from the current, farthest-upstream drainage point. Using our LEFM model, we identify three vulnerable regions with coincidence between dry crevasse depth and water table variability, indicating potential future inland firn aquifer drainage sites. These results suggest the downstream boundary of firn aquifers can migrate inland under future warming scenarios and may already be underway.

The state-of-the-art measurement capabilities of ICESat-2 allow high spatiotemporal resolution of complex ice-dynamic processes that occur during a surge. Detailed and precise mapping of height changes on surging glaciers has previously escaped observations from space due to the limited resolution of space-borne altimeter data and the surface characteristics of glaciers during surge, such as heavy crevassing. This makes geophysical interpretation of deformation and assessment of mass transfer difficult. In this paper, we present an approach that facilitates analysis of the evolution of geophysical processes during a surge, including height changes, crevassing, mass transfer and roughness. We utilize all data from two years of ICESat-2 observations collected during the mature phase of the Negribreen Glacier System surge in 2019 and 2020. The progression of Negribreen's surge has resulted in large-scale elevation changes and wide-spread crevassing, making it an ideal case study to demonstrate ICESat-2 measurement capabilities, which are maximized when coupled with the Density Dimension Algorithm for ice surfaces (DDA-ice). Results show expansion of the surge in upper Negribreen which demonstrates the ability of ICESat-2/DDA-ice to measure a rapidly changing surging glacier and provide the best estimates for cryospheric changes and their contributions to sea-level rise.

Basal crevasses threaten the stability of ice shelves through the potential to form rifts and calve icebergs. Furthermore, it is important to determine the dependence of crevasse stability on temperature due to large vertical temperature variations on ice shelves. In this work, considering the vertical temperature profile through ice viscosity, we compare (1) the theoretical crack depths and (2) the threshold stress causing the transition from basal crevasses to full thickness fractures in several fracture theories. In the Zero Stress approximation, the depth-integrated force at the crevassed and non-crevassed location are unbalanced, violating the volume-integrated Stokes equation. By incorporating a Horizontal Force Balance (HFB) argument, recent work showed analytically that the threshold stress for rift initiation is only half of that predicted by the Zero Stress approximation. We generalize the HFB theory to show that while the temperature profile influences crack depths, the threshold rifting stress is insensitive to temperature. We compare with observations and find that HFB best matches observed rifts. Using HFB instead of Zero Stress for cracks in an ice-sheet model would substantially enlarge the predicted fracture depth, reduce the threshold rifting stress and potentially increase the projected rate of ice shelf mass loss.

Rift propagation, rather than basal melt, drives the destabilization and disintegration of the Thwaites Eastern Ice Shelf. Since 2016, rifts have episodically advanced throughout the central ice-shelf area, with rapid propagation events occurring during austral spring. The ice shelf's speed has increased by ~70% during this period, transitioning from a rate of 1.65 m d−1 in 2019 to 2.85 m d−1 by early 2023 in the central area. The increase in longitudinal strain rates near the grounding zone has led to full-thickness rifts and melange-filled gaps since 2020. A recent sea-ice break out has accelerated retreat at the western calving front, effectively separating the ice shelf from what remained of its northwestern pinning point. Meanwhile, a distributed set of phase-sensitive radar measurements indicates that the basal melting rate is generally small, likely due to a widespread robust ocean stratification beneath the ice–ocean interface that suppresses basal melt despite the presence of substantial oceanic heat at depth. These observations in combination with damage modeling show that, while ocean forcing is responsible for triggering the current West Antarctic ice retreat, the Thwaites Eastern Ice Shelf is experiencing dynamic feedbacks over decadal timescales that are driving ice-shelf disintegration, now independent of basal melt.

Global Navigation Satellite System (GNSS) observations and ground-based timelapse photography obtained over the record-high 2019/2020 melt season are combined to characterise the flexure and fracture behaviour of a previously formed doline on George VI Ice Shelf, Antarctica. The GNSS timeseries shows a downward vertical displacement of the doline centre with respect to the doline rim of ~60 cm in response to loading from a central meltwater lake. The GNSS data also show a tens-of-days episode of rapid-onset, exponentially decaying horizontal displacement, where the horizontal distance between the doline rim and its centre increases by ~70 cm. We interpret this event as the initiation and/or widening of a fracture, aided by stress perturbations associated with meltwater loading in the doline basin. Viscous flexure modelling indicates that the meltwater loading generates tensile surface stresses exceeding 75 kPa. This, together with our timelapse photos of circular fractures around the doline, suggests the first such documentation of meltwater-loading-induced ‘ring fracture’ formation on an ice shelf, equivalent to the fracture type proposed as part of the chain-reaction lake drainage process involved in the 2002 breakup of the Larsen B Ice Shelf.

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.

Crevasses pose severe risks for mountaineers and field glaciologists. Smaller cracks between 0.5 and 2 m are still dangerous, but often not visible in medium resolution satellite imagery. If they are snow covered, they are completely undetectable by optical sensors. We set out to develop an approach to detect potentially crevassed areas by a minimum of geometric data, and to make the method generally applicable to glacier regions. On Vernagtferner, we compared a reference dataset of crevasses observed in high-resolution optical imagery with the curvature of the ice surface and the spatial gradients in driving stress. Both parameters can be derived from a digital surface model and a bedrock model, derived from ice thickness measurements. The correlation patterns show that crevasses preferably form in convex areas and in areas where the driving stress rapidly increases. This corresponds with the theory of crevasse formation. Although the method still misclassifies larger parts, the approach has the potential to define probable non-crevassed areas as well as to aid the planning of safe routes.

In Southeast Greenland, summer melt and high winter snowfall rates give rise to firn aquifers: vast stores of meltwater buried beneath the ice-sheet surface. Previous detailed studies of a single Greenland firn aquifer site suggest that the water drains into crevasses, but this is not known at a regional scale. We develop and use a tool in Ghub, an online gateway of shared datasets, tools and supercomputing resources for glaciology, to identify crevasses from elevation data collected by NASA's Airborne Topographic Mapper across 29000 km2 of Southeast Greenland. We find crevasses within 3 km of the previously mapped downglacier boundary of the firn aquifer at 20 of 25 flightline crossings. Our data suggest that crevasses widen until they reach the downglacier boundary of the firn aquifer, implying that crevasses collect firn-aquifer water, but we did not find this trend with statistical significance. The median crevasse width, 27 meters, implies an aspect ratio consistent with the crevasses reaching the bed. Our results support the idea that most water in Southeast Greenland firn aquifers drains through crevasses. Less common fates are discharge at the ice-sheet surface (3 of 25 sites) and refreezing at the aquifer bottom (1 of 25 sites).

Mass loss from iceberg calving at marine-terminating glaciers is one of the largest and most poorly constrained contributors to sea-level rise. However, our understanding of the processes controlling ice fracturing and crevasse evolution is incomplete. Here, we use Gabor filter banks to automatically map crevasse density and orientation through time on a ~150 km2 terminus region of Narsap Sermia, an outlet glacier of the southwest Greenland ice sheet. We find that Narsap Sermia is dominated by transverse (flow-perpendicular) crevasses near the ice front and longitudinal (flow-aligned) crevasses across its central region. Measured crevasse orientation varies on sub-annual timescales by more than 45$^\circ$ in response to seasonal velocity changes, and also on multi-annual timescales in response to broader dynamic changes and glacier retreat. Our results show a gradual up-glacier propagation of the zone of flow-transverse crevassing coincident with frontal retreat and acceleration occurring in 2020/21, in addition to sub-annual crevasse changes primarily in transition zones between longitudinal to transverse crevasse orientation. This provides new insight into the dynamics of crevassing at large marine-terminating glaciers and a potential approach for the rapid identification of glacier dynamic change from a single pair of satellite images.

Basal crevasses are macroscopic structural discontinuities at the base of ice sheets and glaciers that arise by fracture. Motivated by observations and by the mechanics of elastic fracture, we hypothesise that spatial variations in basal stress (in the presence of basal water pressure) can promote and localise basal crevassing. We quantify this process in the theoretical context of linear elastic fracture mechanics. We develop a model evaluating the effect of shear-stress variation on the growth of basal crevasses. Our results indicate that sticky patches promote the propagation of basal crevasses, increase their length of propagation into the ice and, under some conditions, give them curved trajectories that incline upstream. A detailed exploration of the parameter space is conducted to gain a better understanding of the conditions under which sticky-patch-induced basal crevassing is expected beneath ice sheets and glaciers.

Tasman Glacier, a temperate maritime glacier in the New Zealand Southern Alps, is rapidly receding. Climate warming is resulting in lengthening of the ablation season, meaning crevasses in the accumulation area are becoming exposed at the surface for longer. We combine measurements of air temperature and wind speed from inside crevasses with surface meteorological data, finding that during summer, in-crevasse air temperature is frequently positive, and can at times exceed surface air temperature. Greatest warming occurred in the widest crevasses during clear-sky conditions, but full depth warming of crevasses also occurred at night. Net shortwave radiation contributes to heating of air in the upper regions of crevasses, but turbulent sensible heat transfer was responsible for driving warm air deeper into crevasses. Crevasses orientated to maximise radiation retrieval, and running parallel to wind flow, have the greatest potential for warming and heat storage. We hypothesise a positive feedback loop in the surface energy-balance system, where crevasses entrain and trap heat, which enhances melting, that in turn enlarges the crevasses, enabling greater heat storage and further melting. Energy-balance models that treat accumulation areas of alpine glaciers as homogeneous surfaces will therefore underestimate snow melt and overestimate mass balance.

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.

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

The representation of iceberg calving in numerical models is a key source of uncertainty in century-scale sea-level rise projections. Parameters central to model representations of calving, including the tensile strength of glacier ice, remain poorly constrained. Grain-size and sample-size dependence make it difficult to reconcile laboratory and in situ estimates of ice tensile strength. Further, assumptions of various numerical models obscure comparison of the ‘strength’ parameter with a physically observable value. Here, we address the problem of fracture during calving using an analogous natural laboratory: a viscoelastic analysis of observed surface deformation and associated stresses in the 2015 collapse of eastern Skaftá cauldron, Vatnajökull ice cap, Iceland. We find that the ice within the cauldron could have experienced instantaneous elastic stress on the order of several MPa. We observe surface crevasses at the cauldron edges and center, but find that large areas of ice remain intact despite high stress. Our findings suggest a tensile strength of glacier ice on the order of 1 MPa, consistent with laboratory estimates but exceeding previous glacier-specific estimates.

On 30 September 2017, an Air France Airbus A380-800 suffered a failure of its fourth engine while over Greenland. This failure resulted in the loss of the engine fan hub, fan blades and surrounding structure. An initial search recovered 30 pieces of light debris, but the primary part of interest, a ~220 kg titanium fan hub, was not recovered because it had a different fall trajectory than the light debris, impacted into the ice-sheet's snow surface, and was quickly covered by drifting snow. Here we describe the methods used for the detection of the fan hub and details of the field campaigns. The search area included two crevasse fields of at least 50 snow-covered crevasses 1 to ~30 m wide with similar snow bridge thicknesses. After 21 months and six campaigns, using airborne synthetic aperture radar, ground-penetrating radar, transient electromagnetics and an autonomous vehicle to survey the crevasse fields, the fan hub was found within ~1 m of a crevasse at a depth of ~3.3 to 4 m and was excavated with shovels, chain saws, an electric winch, sleds and a gasoline heater, by workers using fall-arrest systems.

Hydrofracturing can enhance the depth to which crevasses propagate and, in some cases, allow full depth crevasse penetration and iceberg detachment. However, many existing crevasse models either do not fully account for the stress field driving the hydrofracture process and/or treat glacier ice as elastic, neglecting the non-linear viscous rheology. Here, we present a non-local continuum poro-damage mechanics (CPDM) model for hydrofracturing and implement it within a full Stokes finite element formulation. We use the CPDM model to simulate the propagation of water-filled crevasses in idealized grounded glaciers, and compare crevasse depths predicted by this model with those from linear elastic fracture mechanics (LEFM) and zero stress models. We find that the CPDM model is in good agreement with the LEFM model for isolated crevasses and with the zero stress model for closely-spaced crevasses, until the glacier approaches buoyancy. When the glacier approaches buoyancy, we find that the CPDM model does not allow the propagation of water-filled crevasses due to the much smaller size of the tensile stress region concentrated near the crevasse tip. Our study suggests that the combination of non-linear viscous and damage processes in ice near the tip of a water-filled crevasse can alter calving outcomes.

Russian (former Soviet) systematic studies of Antarctica by radio-echo sounding (RES) and ground-penetrating radar technique (GPR) were commenced in 1964. Since that time airborne RES surveys have covered about 5.5 × 106 km2 of the icy continent discovering remarkable geographic objects such as Subglacial Gamburtsev Mountains, and allowed studies of Filchner-Ronne Ice Shelf, Amery Ice Shelf and Lambert Glacier. Ground-based investigations during the 1990s and 2000s revealed the structure of the Lake Vostok area and surveyed along the Mirny to Vostok and Progress to Vostok traverse routes. GPR studies during the 2010s were to select the site for a new snow-runway at Mirny Station, with the resumption of the aviation after a 25 year hiatus.

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.

While large-scale observations of intensified fracture and rifting can be observed through remote-sensing observations, understanding crevasse initiation may best be achieved with small-scale observations in which crevasses can be directly observed. Here we investigate the kinematic drivers of crevasse initiation in the McMurdo Shear Zone (MSZ), Antarctica. We delineated 420 crevasses from ~95 km of 400 MHz frequency ground-penetrating radar data and compared these data with kinematic outputs derived from remotely-sensed ice surface velocities to develop a statistical method to estimate crevasse initiation threshold strain rate values. We found the MSZ to be dominated by simple shear and that surface shear strain rates proved best for predicting crevasse features, with regions of higher shear strain rate more likely to have a greater number of crevasses. In the surveyed portion of our study region, values of shear strain rate and vorticity rate derived from the MEaSUREs2 velocity dataset range between 0.005–0.020 and 0.006–0.022 a−1, respectively, with crevasses located at ≥0.011 and ≥0.013 a−1. While threshold values from this study cannot be directly applied to other glacial environments, the method described here should allow for the study of shear margin evolution and assessment of localized damage and weakening processes in other locations where in situ data are available.