Search

The mixing of water in the ocean is a key influence on our climate and ecosystems, but its importance is often under-recognized. Mixing in the seas around Antarctica—a key process that redistributes heat, carbon, nutrients, plankton, and all other things in the sea, with profound consequences—also affects the stability of the continent's glaciers and ice sheets, with consequences for sea level rise globally. A recent discovery showed that when the fronts of glaciers disintegrate, they "calve" huge chunks of ice that can cause underwater tsunamis in the ocean, which can spread for many miles and cause strong bursts of mixing when they break. This surprising finding is changing the way we think about mixing close to Antarctica, what causes it, and how it matters. [ABSTRACT FROM AUTHOR]

Antarctica, the only continent without permanent human settlements, holds important messages for humanity. Recent studies of carbon dioxide concentrations in Antarctic ice cores reveal the impact of European colonization on the Western Hemisphere, leading to the abandonment of agricultural land and the regrowth of natural vegetation. The breaking off of iceberg A23a highlights the threat of melting ice shelves and the subsequent rise in global sea levels. Record-low levels of sea ice in Antarctica suggest a new state influenced by climate change. The International Association of Antarctic Tourism Operators (IAATO) plays a crucial role in preserving Antarctica's unique species and ecosystem. However, Antarctic tourism has environmental costs that need to be considered. Bird flu has reached Antarctica through migrating seabirds, posing a risk to penguins and marine mammals. Overfishing of krill, the base of the Antarctic food chain, could lead to an ecological catastrophe. Melting ice in Antarctica and Greenland has even affected the Earth's rotation. [Extracted from the article]

Marine‐terminating glacier fjords play a central role in the transport of oceanic heat toward ice sheets, regulating their melt. Mixing processes near glacial termini are key to this circulation but remain poorly understood. We present new summer measurements of circulation and mixing near a marine‐terminating glacier with active sub‐glacial discharge. 65% of the fjord's vertical overturning circulation is driven by the buoyant plume, however we newly report intense vertical and horizontal mixing in the plume's horizontal spreading phase, accounting for the remaining 35%. Buoyant plume theory supports 2%–5% of total glacial melt. Thus, most of the heat associated with vertical overturing short‐circuits the glacial front. We find however that turbulence in the horizontal spreading phase redistributes the short‐circuited heat back into the surface waters of the near‐glacial zone. Our findings highlight the need for further research on the complex mixing processes that occur near the glacier terminus. Plain Language Summary: Melting of glacial ice is the single largest contributor to global sea‐level rise. Many glaciers flow into the ocean where the near‐vertical ice wall is bathed in relatively warm sea water. Freshwater from ice surface melting seeps down through cracks and crevasses in the glacier to be discharged at the base of the ice wall, many tens to hundreds of meters below the sea surface. This fresh water rises from the depths as a highly turbulent plume, drawing in and pushing upwards the surrounding seawater, eventually spreading horizontally as a mixture of fresh discharge and mixed‐in seawater. Few measurements exist in the dangerous zone where iceberg often calve. Using data from a robotic platform we show that the vertical rise and the horizontal spreading of fresh water both play important roles in the total ocean‐induced melting of the glacial face in this type of system. Key Points: Entrainment into the buoyant plume drives ∼65% of fjord vertical overturning circulationIntense turbulent mixing in the horizontal spreading phase of the plume drives the remaining ∼35%95%–98% of the heat in the rising plume short‐circuits the glacier, to be redistributed into glacial proximal waters by horizontal mixing [ABSTRACT FROM AUTHOR]

Projections of sea‐level rise from ice‐sheet shrinkage in a warming world have large uncertainties, linked to limited knowledge of changes at the ocean‐ice sheet interface. This interface most typically is modeled as a grounding line, across which still‐connected ice flows into the ocean to float as an ice shelf, or where icebergs calve from a cliff before the ice begins to float. But, extensive and rapidly increasing evidence shows that this is really a grounding zone, and that processes in this grounding zone omitted from many models could exert major controls on sea‐level rise. Plain Language Summary: Marine‐terminating glaciers flow into the ocean across extensive grounding zones. These kilometers‐long and glacier‐wide zones represent the last broad region of glacier contact with the rock and sediments below before the ice enters the ocean as a floating ice shelf or calved icebergs. Loss of this basal drag along with enhanced basal melting caused by tidally driven seawater intrusion leads to faster outflow and rapid thinning of the overlying ice. As a result of the local thinning, grounding zones retreat inland and sea level rises with more loss of previously grounded ice. Most ice‐sheet models used in sea‐level projections do not include grounding‐zone processes, but rather stop their ice‐ocean interactions at a grounding line. They are thereby omitting important dynamic feedbacks and underestimating future sea‐level contributions from the marine‐based sectors of the Greenland and Antarctic ice sheets. Key Points: Tidally modulated seawater intrusion leads to loss of ice‐bed contact as well as significant (maximal) basal melting within grounding zonesThe future dynamics of marine outlet glaciers are ultimately controlled by coupled processes operating within and through grounding zonesDespite the importance of grounding zones to ice‐sheet dynamics, most ice‐sheet models used in sea‐level projections do not include them [ABSTRACT FROM AUTHOR]

Terminus change is a complex outcome of ice‐ocean boundary processes and poses challenges for ice sheet models due to inadequate calving laws, creating uncertainty in sea level change projections. To address this, we quantify glacier termini sinuosity and convexity, testing the hypothesis that terminus morphology reflects dominant calving processes. Using 10 glaciers with diverse calving styles in Greenland over the period from 1985 to 2021, we establish a supervised classification of calving style by comparing morphology and literature‐derived calving observations. Validation with four of these glaciers and flotation conditions and subglacial discharge routing observations confirms concave, smooth termini indicate buoyant flexure dominated‐calving, while convex, sinuous termini suggest serac failure dominated‐calving. We also identify a mixed style where both calving types may occur. We use these classes to label calving style from 1985 to 2021 for all 10 glaciers and explore how this changes over time as glaciers retreat. Plain Language Summary: The ice‐ocean boundary (terminus) of the Greenland ice sheet changes over space and time due to the melt and calving of ice. There is a range of calving behaviors due to many environmental and geometric factors, including changing ocean temperatures, glacier melting, and the geometry of the glacier. This makes representing how these boundary changes in ice sheet models are difficult, leading to uncertainties in predicting sea level rise from retreating glaciers. The study focuses on the shape (morphology) of 10 glaciers from 1985 to 2021, finding that smooth, concave termini indicate a type of calving related to the flotation of the glacier, while sinuous, convex termini indicate melt‐dominated calving. The research helps classify calving styles and understand how they change over time as glaciers retreat. Key Points: Two calving styles are linked to distinct terminus morphologies but both can exist for a single glacier (over time or space)Mixed styles of calving are more difficult to classifyClassifying terminus morphology over time can aid in identifying how terminus calving mechanisms change over time [ABSTRACT FROM AUTHOR]

The prescription of a simple and robust parameterization for calving is one of the most significant open problems in ice sheet modelling. One common approach to the modelling of crevasse propagation in calving in ice shelves is to view crevasse growth as an example of linear elastic fracture mechanics. Prior work has, however, focused on highly idealized crack geometries, with a single fracture incised into a parallel-sided slab of ice. In this paper, we study how fractures growing from opposite sides of such an ice slab interact with each other, focusing on different simple crack arrangements: we consider either perfectly aligned cracks or periodic arrays of laterally offset cracks. We visualize the dynamics of crack growth using simple tools from dynamical systems theory and find that aligned cracks tend to impede each other's growth due to the torques generated by normal stresses on the crack faces, while periodically offset cracks facilitate simultaneous growth of bottom and top cracks. For periodic cracks, the presence of multiple cracks on one side of the ice slab, however, also generates torques that slow crack growth, with widely spaced cracks favouring calving at lower extensional stresses than closely spaced cracks. [ABSTRACT FROM AUTHOR]

Understanding the drivers of iceberg calving from Antarctic ice shelves is important for future sea level rise projections. Ocean waves promote calving by imposing stresses and strains on the shelves. Previous modeling studies of ice shelf responses to ocean waves have focused on highly idealized geometries with uniform ice thickness and a flat seabed. This study leverages on a recently developed mathematical model that incorporates spatially varying geometries, combined with measured ice shelf thickness and seabed profiles, to conduct a statistical assessment of how 15 Antarctic ice shelves respond to ocean waves over a broad range of relevant wave periods, from swell to infragravity waves to very long period waves. The results show the most extreme responses at a given wave period are generated by features in the ice shelves and/or seabed geometries, depending on the wave regime. Relationships are determined between the median ice shelf response and the median shelf front thickness or the median water cavity depth. The findings provide further evidence of the role of ocean waves in large‐scale calving events for certain ice shelves (particularly the Wilkins) and indicate a possible role of ocean waves in calving events for other shelves (Larsen C and Conger). Further, the relationships determined provide a method to assess the potential for increased calving as ice shelves evolve with climate change, and, hence, contribute to assessments of future sea level rise. Plain Language Summary: Antarctic ice shelves are the floating extensions of the Antarctic Ice Sheet. They play a critical role in many Southern Ocean processes. In particular, they help maintain the stability of the Antarctic Ice Sheet by moderating the flow of grounded ice into the Southern Ocean. Climate change is causing them to thin and retreat, which is a major threat to global sea levels. Iceberg calving accounts for half of ice shelf loss, and ocean waves contribute to the calving process by rhythmically bending ice shelves. The influence of ocean waves on calving is expected to increase as the shelves and their surrounding sea ice barriers become weaker. Therefore, quantifying the responses of ice shelves to ocean waves is needed to project the future of the shelves and, hence, sea level rise. In this study, we use a recently developed mathematical model to analyze the responses of 15 Antarctic ice shelves to ocean waves, ranging from storm waves to tsunamis. We show how features in the geometry can create large responses that could drive calving events, and we derive simple relationships between the responses and the geometry to aid projections of future scenarios. Key Points: Crevasses and seabed protrusions create large ice shelf flexure in response to ocean wavesIce shelves that have experienced large scale calving events had much greater responses to swell than typical shelvesMedian ice shelf responses to swell are strongly correlated to median shelf front thicknesses [ABSTRACT FROM AUTHOR]

Lago Argentino is one of the largest lakes in Patagonia. Some of its arms host calving glaciers that flow from the southern Patagonian Ice Field. The best known of these is the Perito Moreno glacier. Although recent studies have reconstructed part of glacier dynamics in the southern arms of Lago Argentino, there are no studies yet for the northern arms. A high-resolution seismic survey was carried out in the northern area of Punta Bandera, i.e., in the sector between the Brazo Norte and the Canal de los Témpanos to characterize the stratigraphy of the sedimentary fill of the lake and identify the main unconformities related to glacial dynamics. Eight seismic facies grouped into two seismic units were mapped, and ten unconformities were recognized. Based on the acquired data, a bathymetric grid was created for the first time to obtain a morphological representation of the lake floor and access to a previously unmapped area of Lago Argentino. The analysis of the geometries and acoustic properties of the different seismic facies identified, correlated with the available information on land moraines, has allowed us to partially reconstruct the advance and retreat phases of the glacial lobes that formerly were present in the two arms of the lake, although the temporal assignments are still speculative due to the lack of calibrations with boreholes. Nevertheless, this study provides important information for deciphering the glaciation history of a key sector of the southern Patagonian Ice Field during the late Pleistocene/Holocene. The data show that the northern and southern areas have different lake floor morphologies, with the southern, shallower area of the Canal de los Témpanos being glaciated first than the northern, more overdeepened area of the Brazo Norte. [ABSTRACT FROM AUTHOR]

The mass balance of the Greenland Ice Sheet is strongly influenced by the dynamics of its outlet glaciers. Therefore, it is of paramount importance to accurately and continuously monitor these glaciers, especially the variation in their frontal positions. A temporally comprehensive parameterization of glacier calving is essential for understanding dynamic changes and constraining ice sheet modeling. However, many current calving front records are limited in terms of temporal resolution as they rely on manual delineation, which is laborious and not appropriate considering the increasing amount of satellite imagery available. In this contribution, we address this problem by applying an automated method to extract calving fronts from optical satellite imagery. The core of this workflow builds on recent advances in the field of deep learning while taking full advantage of multispectral input information. The performance of the method is evaluated using three independent test datasets. For the three datasets, we calculate mean delineation errors of 61.2, 73.7, and 73.5 m, respectively. Eventually, we apply the technique to Landsat-8 imagery. We generate 9243 calving front positions across 23 outlet glaciers in Greenland for the period 2013–2021. Resulting time series not only resolve long-term and seasonal signals but also resolve subseasonal patterns. We discuss the implications for glaciological studies and present a first application for analyzing the effect of bedrock topography on calving front variations. Our method and derived results represent an important step towards the development of intelligent processing strategies for glacier monitoring, opening up new possibilities for studying and modeling the dynamics of Greenland's outlet glaciers. [ABSTRACT FROM AUTHOR]

The recent calving of Astrolabe Glacier on the Adélie Coast (East Antarctica) in November 2021 presents an opportunity to better understand the processes leading to ice tongue fracturing. To document the fractures and rift evolution that led to the calving, we used the archive of Sentinel-2 optical images to measure the ice motion and strain rates from 2017 to 2021. The long-term evolution of the Astrolabe ice tongue is mapped using airborne and satellite imagery from 1947 to November 2021. These observations are then compared with measurements of sea-ice extent and concentration. We show that calving occurs almost systematically at the onset of or during the melting season. Additionally, we observe a significant change in the periodicity of sea ice surrounding Astrolabe Glacier in the last decade (2011–2021) compared to previous observations (1979–2011), which has resulted in a change in the Astrolabe calving cycle. Indeed, one can observe a decrease in the duration of sea-ice-free conditions during the austral summers after 2011 in the vicinity of the glacier, which seems to have favoured spatial extension of the ice tongue. However, the analysis of strain rate time series revealed that the calving of November 2021 (20 km 2) occurred at the onset of sea-ice melting season but resulted from the glacier dislocation that took place suddenly in June 2021 in the middle of the winter. These observations indicate that while sea ice can protect and promote the spatial extension of a glacier ice tongue, its buttressing is not sufficient to inhibit rifting and ice fracturing. [ABSTRACT FROM AUTHOR]

The calving of glaciers regularly produces tsunami‐like waves that pose a serious threat to coastal environments. Those strong waves are not only able to move ice mélange and redistribute icebergs, growlers, or sea ice across a fjord but also flood and remodel neighbouring cliffs and beaches. Here, we analyze over 90 years (1929–2023) of coastal zone changes that occurred in front of Eqip Sermia. We show that calving waves play a dominant role in transforming the lateral moraine and forming a beach and spit system south of the glacier front. Part of the former moraine has transformed into a boulder‐dominated spit, which closed the lagoon over the years. By multidecadal analysis, we also detected a significant erosion of unconsolidated cliffs located on the opposite side of the bay (~0.53 m per year between 1985 and 2023). In addition, we demonstrate that even a single event (one calving wave) can remodel a beach surface by entrainment of up to 1.8‐m‐diameter boulders and the erosion of the beach surface by washing away sand and gravel from rocky outcrops. Our study constitutes important progress toward modes of paraglacial coastal evolution in regions characterized by rapidly retreating calving glaciers. [ABSTRACT FROM AUTHOR]

The ice-gouging topography of bed of the southwestern Kara Sea is the result of the impact of icebergs and sea ice. During cruise 52 of the R/V Akademik Nikolaj Strakhov using a multibeam echo sounder, we collected representative data on the key parameters of ice scours (location, orientation, depth, width) for the first time, which allows us to draw conclusions on a regional scale. We revealed regularities in the spatial and temporal distribution of ice scours and their density in different parts of the seabed. It was revealed that the maximum dimensions of the ice scours decrease with distance from the sources of iceberg calving from NW to SE. The orientation of the ice scours correlates with the main drift directions of icebergs. Most of the ice scours are located at depths up to 220 m and could have formed both at the modern and lower sea levels (in postglacial time). We identified a high degree of seabed transformation by ice-gouging processes in the southwestern Kara Sea. [ABSTRACT FROM AUTHOR]

During Heinrich events, great armadas of icebergs episodically flooded the North Atlantic Ocean and weakened overturning circulation. The ice discharges of these episodes constrain the sensitivity of overturning circulation to iceberg melting. We reconstructed these ice discharges to be as high as 0.13 sverdrup (Sv) (1 Sv = 1 million cubic meters per second) during Heinrich event 4 and to average 0.029 Sv over all episodes. The present-day Greenland Ice Sheet calving of icebergs is comparable to that of a mid-range Heinrich event. As the future Greenland Ice Sheet recedes from marine-terminating outlets, its iceberg calving likely will not persist long enough for icebergs alone to cause catastrophic disruption to the Atlantic overturning circulation, although the accelerating Greenland runoff and continued global warming remain threats to the circulation stability. [ABSTRACT FROM AUTHOR]

Herein, fast fracture initiation in glacier ice is modeled using a Material Point Method and a simplified constitutive law describing tensile strain softening. Relying on a simple configuration where ice flows over a vertical step, crevasse patterns emerge and are consistent with previous observations reported in the literature. The model's few parameters allows identification of a single dimensionless number controlling fracture spacing and depth. This scaling law delineates two regimes. In the first one, ice thickness does not play a role and only ice tensile strength controls the spacing, giving rise to numerous surface crevasses, as observed in crevasse fields. In this regime, scaling can recover classical values for ice tensile strength from macroscopic field observations. The second regime, governed by ice bending, produces large‐scale, deep fractures resembling serac falls or calving events. Plain Language Summary: In ice sheets and alpine glaciers, fast‐flowing sections are often characterized by crevasse fields that play a significant role in the cryo‐hydrologic system by facilitating meltwater flow, enhancing basal sliding, weakening the ice, and impacting glacier thermodynamics. Modeling these fractures at the glacier scale remains challenging and often necessitates integrating diverse models which hinders the straightforward consideration of physical issues associated with crevasse fields on a large scale. Here, a new numerical framework allows us to conduct field‐scale experiments and paves the way for a scaling law to elucidate the macroscopic factors influencing fracture fields and to easily incorporate crevasse depth and spacing into large‐scale models. A newly discovered scaling law highlights the transition between a mechanical behavior where the regular crevasse spacing is unaffected by geometry to a regime where geometry plays a significant role, particularly in large‐scale fracture processes like glacier calving. While the numerical experiments in this paper focus on glaciers, the model and conceptual framework is versatile and can address the mechanical behavior of fractures in broader geophysical contexts such as snow, rock or ice avalanches, tectonics and landslides. Key Points: Fractures in glacier flow are modeled using material point method with elastoplasticity and tensile strain softeningA dimensional analysis reveals a key dimensionless number characterizing two different regimes of fast fractureOne regime predicts acknowledged ice tensile strength from field observations and characterizes the regular crevasse spacing [ABSTRACT FROM AUTHOR]

This study examines the feature space of seismic waveforms often used in machine learning applications for seismic event detection and classification problems. Our investigation centres on the southern Alaska region, where the seismic record captures diverse seismic activity, notably from the calving of marine-terminating glaciers and tectonic earthquakes along active plate boundaries. While the automated discrimination of earthquakes and glacier quakes is our nominal goal, this data set provides an outstanding opportunity to explore the general feature space of regional seismic phases. That objective has applicability beyond ice quakes and our geographic region of study. We make a noteworthy discovery that features rooted in the spectral content of seismic waveforms consistently outperform statistical and temporal features. Spectral features demonstrate robust performance, exhibiting resilience to class imbalance while being minimally impacted by factors such as epicentral distance and signal-to-noise ratio. We also conduct experiments on the transferability of the model and find that transferability primarily depends on the appearance of the waveforms. Finally, we analyse misclassified events and find examples that are identified incorrectly in the original regional catalogue. [ABSTRACT FROM AUTHOR]

The largest uncertainty in future sea-level rise is loss of ice from the Greenland and Antarctic Ice Sheets. Ice shelves, freely floating platforms of ice that fringe the ice sheets, play a crucial role in restraining discharge of grounded ice into the ocean through buttressing. However, since the 1990s, several ice shelves have thinned, retreated, and collapsed. If this pattern continues, it could expose thick cliffs that become structurally unstable and collapse in a process called marine ice cliff instability (MICI). However, the feedbacks between calving, retreat, and other forcings are not well understood. Here we review observed modes of calving from ice shelves and marine-terminating glaciers, and their relation to environmental forces. We show that the primary driver of calving is long-term internal glaciological stress, but as ice shelves thin they may become more vulnerable to environmental forcing. This vulnerability—and the potential for MICI—comes from a combination of the distribution of preexisting flaws within the ice and regions where the stress is large enough to initiate fracture. Although significant progress has been made modeling these processes, theories must now be tested against a wide range of environmental and glaciological conditions in both modern and paleo conditions. Ice shelves, floating platforms of ice fed by ice sheets, shed mass in a near-instantaneous fashion through iceberg calving. Most ice shelves exhibit a stable cycle of calving front advance and retreat that is insensitive to small changes in environmental conditions. Some ice shelves have retreated or collapsed completely, and in the future this could expose thick cliffs that could become structurally unstable called ice cliff instability. The potential for ice shelf and ice cliff instability is controlled by the presence and evolution of flaws or fractures within the ice. [ABSTRACT FROM AUTHOR]

Forecasts of sea-ice motion and fragmentation are of vital importance for all human interactions with sea ice, ranging from those involving indigenous hunters to shipping in polar regions. Sea-ice models are also important for simulating long-term changes in a warming climate. Here, we apply the Helsinki Discrete Element Model (HiDEM), originally developed for glacier calving, to sea-ice breakup and dynamics. The code is highly optimized to utilize high-end supercomputers to achieve an extreme time and space resolution. Simulated fracture patterns and ice motion are compared with satellite images of the Kvarken region of the Baltic Sea from March 2018. A second application of HiDEM involves ice ridge formation in the Gulf of Riga. With a few tens of graphics processing units (GPUs), the code is capable of reproducing observed ice patterns that in nature may take a few days to form; this is done over an area of ∼100km×100km , with an 8 m resolution, in computations lasting ∼10 h. The simulations largely reproduce observed fracture patterns, ice motion, fast-ice regions, floe size distributions, and ridge patterns. The similarities and differences between observed and computed ice dynamics and their relation to initial conditions, boundary conditions, and applied driving forces are discussed in detail. The results reported here indicate that the HiDEM has the potential to be developed into a detailed high-resolution model for sea-ice dynamics at short timescales, which, when combined with large-scale and long-term continuum models, may form an efficient framework for forecasts of sea-ice dynamics. [ABSTRACT FROM AUTHOR]

Baffin Bay is the travel destination of most icebergs calving from west Greenland. They commonly follow the bay's cyclonic circulation and might end up far south along the coast of Newfoundland and Labrador, where many shipping routes converge. Given the hazard that icebergs pose to marine transportation, understanding their distribution is fundamental. One of the forces driving iceberg drift arises from the presence of sea ice. Observations in the Southern Ocean indicate that icebergs get locked in thick and concentrated sea ice. We present observations that support the occurrence of this sea ice locking mechanism (SIL) in Baffin Bay as well. Most iceberg models, however, represent the sea ice force over an iceberg as a simple drag force. Here, we implement a new parameterization in the iceberg module of the Nucleus for European Modeling of the Ocean (NEMO‐ICB) to represent SIL. We show that, by using this new parameterization, icebergs are more likely to travel outside of the Baffin Island Current during winter, which is supported by satellite observations. There is a slight improvement in the representation of iceberg severity along the coast of Newfoundland and Labrador and a slight shift of iceberg melt toward this region and Lancaster Sound/Hudson Strait. Although the impacts of icebergs on sea ice are still not represented, and targeted observations are needed for model calibration regarding sea ice concentration thresholds from which icebergs get locked, we are confident that this model improvement takes iceberg modeling one step forward toward reality. Plain Language Summary: After they break off Greenland's glaciers, icebergs drift in the ocean in somewhat predictable patterns across Baffin Bay. What makes them predictable is that we know the main natural forces that cause icebergs to move: the winds, the ocean currents, and the sea ice cover. Because it is impossible to monitor all the icebergs coming from Greenland, we use computer models to understand their typical trajectories. In these models, each one of the forces responsible for iceberg drift can be represented mathematically, usually in a very simplistic way. The sea ice force, for example, is proportional to how much friction there is between the iceberg and sea ice. However, observations in nature show that the iceberg‐sea ice interaction is more complicated than that: if the sea ice cover is thick and compact, it is able to trap icebergs, and both will move as one solid block. In this study, we improve a model's representation of the sea ice force by including this trapping mechanism. Results show that iceberg trajectory can be slightly different during winter due to sea ice trapping and that this mechanism potentially improves our ability to predict when icebergs reach the transatlantic shipping routes east of Canada. Key Points: A new sea ice‐iceberg locking parameterization was implemented in Nucleus for European Modeling of the Ocean and tested in a 1/4° resolution regional configurationIcebergs locked in sea ice in Baffin Bay are more likely to travel outside of the Baffin Island Current than "unlocked" onesThe existence of a sea ice locking mechanism in Baffin Bay is supported by observations of co‐varying iceberg and sea ice speeds [ABSTRACT FROM AUTHOR]

The extent, chronology and dynamics of the pre‐Marine Isotope Stage (MIS) 2 last British–Irish Ice Sheet (BIIS) are not well known. Although the BRITICE‐CHRONO Project has detailed the maximum extent and retreat phases of the last BIIS for the period after 30 ka and into the Last Glacial Maximum (LGM), the Project identified several pre‐existing datasets and generated new data that implied glaciation pre‐dating the LGM but which post‐dated the Last Interglacial (Eemian; MIS5e); these data are reviewed here. There are no dated till units but are other indicators clearly indicative of glaciation: deep‐sea ice‐rafted detritus flux into the adjacent NE Atlantic, cosmogenic rock‐exposure age dating from glaciated surfaces in Wales and the island of Lundy (Bristol Channel), and optically stimulated luminescence (OSL) ages of proximal glacifluvial sequences on the Isle of Lewis (Outer Hebrides) and in the Cheshire Basin. Taken together these indicate BIIS inception during MIS5d, growth into MIS4 and evidence for dynamic retreat–advance phases during MIS3. OSL evidence for high relative sea level indicated by raised beaches in southern Ireland during MIS4 and 3 at a time of lowered glacio‐eustatic sea level indicates substantial glacial isostatic loading, explained by the early growth of the BIIS during the last cold stage. High relative sea level during MIS4 and 3 coincident with adjacent calving ice sheet margins provides an explanation for the rafted giant erratics found around the shores of southern Britain and Ireland. [ABSTRACT FROM AUTHOR]

Frontal ablation is responsible for a large fraction of the mass loss from tidewater glaciers. The main contributors to frontal ablation are iceberg calving and submarine melting, with calving often being the largest. However, submarine melting, in addition to its direct contribution to mass loss, also promotes calving through the changes induced in the stress field at the glacier terminus, so both processes should be jointly analysed. Among the factors influencing submarine melting, the formation of a buoyant plume due to the emergence of fresh subglacial water at the glacier grounding line plays a key role. In this study we used Elmer/Ice to develop a 3D glacier dynamics model including calving and subglacial hydrology coupled with a line plume model to calculate the calving front position at every time step. We applied this model to the Hansbreen–Hansbukta glacier–fjord system in southern Spitsbergen, Svalbard, where a large set of data are available for both the glacier and the fjord from September 2008 to March 2011. We found that our 3D model reproduced the expected seasonal cycle of advance–retreat. Besides, the modelled front positions were in good agreement with the observed front positions at the central part of the calving front, with longitudinal differences, on average, below 15 m for the period from December 2009 to March 2011. But there were regions of the front, especially the eastern margin, that presented major differences. [ABSTRACT FROM AUTHOR]

Ice shelves limit the flux of grounded ice into the ocean by buttressing the discharge of land-based ice upstream. Ice shelf weakening and collapse can lead to decreased buttressing and observations increasingly show that some ice shelves have experienced increased melt and increased calving, with recent hypotheses suggesting that increased melt leads to increased fracturing. However, the specific processes that control this correlation are not yet understood, with mechanisms other than melt affecting fracturing. Here we use the topography of the ice shelf base from BedMachine to investigate how basal melting and ice deformation contribute to crevasse and melt channel formation and evolution on the Pine Island Ice Shelf in West Antarctica. We find that high basal melt rates and high first principal strain rates lead to substantial roughening of the ice shelf through a collection of features, including melt channels and crevasses. Critically, melt channels and crevasses are the deepest in all directions at locations where the highest rates of melting and straining occur simultaneously. This suggests that the combination of melt rates and strain rates work in tandem to excavate and seed the deepest melt channels and crevasses on ice shelves. These features then may form lines of weakness that transform into rifts and, ultimately, the detachment boundary for calving events. This implies that melt and fracture play an important role in controlling the dynamics of ice shelves. Plain Language Summary Future sea level rise is tied to floating Antarctic ice shelves which limit the flow of ice from the continent. Ice shelf collapse could trigger an acceleration of land-based ice into the ocean due to the loss of support, leading to future sea level rise. It has been shown that ice shelves with a higher melt have more crevasses, are more fractured, and are rough. However, the specific processes driving the potential connection between melt and fracture are not well understood. Here we use ice shelf topography to study the relationship between roughness and the size and spacing of geometric features on the Pine Island Ice Shelf. We find that high melt rates and high strain rates result in deeper and wider features on the ice shelf. We also find that, in contrast to the results of previous studies, both melt and fracture contribute to the deepest features. This suggests that these processes contribute to enhancing the depth of features. We hypothesize that this amplifying interaction may then promote rifting and iceberg calving, hinting that it is the combination of melt and fracture that is important in controlling the dynamics of ice shelves. [ABSTRACT FROM AUTHOR]

Over the past decades, seismic monitoring has been increasingly used to track glacial activities associated with ice loss. Many seismological studies focus on West Antarctica, whereas glacial seismicity in East Antarctica is much less studied. Here, we apply unsupervised deep learning to a dense nodal seismic array near Dålk Glacier, East Antarctica, operating from 6 December 2019 to 2 January 2020. An autoencoder is used to automatically extract event features, which are then input into a Gaussian mixture model for clustering. We divide the data into 50 clusters and merge them according to their temporal and spectral characteristics. The results reveal five main types of seismic signals: two groups with monochromatic and high frequencies, two groups with broadband frequency and short duration, and a group with mainly low frequency and long duration. By comparing the environmental conditions (wind, temperature and tides), we infer that the two monochromatic groups are wind-induced vibrations of the near-station flag markers and topography; the two broadband groups are likely thermal contractions on the blue ice surface and stick-slip events at the ice base; and the low-frequency events are water-filled basal crevassing and iceberg calving. In particular, we observe one type of low-frequency event preceded by high-frequency onset, which is likely basal crevassing near the grounding line of Dålk Glacier and predominantly occurred during rising tides. Our findings show that deep clustering is effective in identifying a wide range of glacial seismic events and can contribute to the rapid growth of passive glacier seismic monitoring. Plain Language Summary Over the past decades, the acceleration of ice mass loss in Antarctica has contributed to rising sea levels. As the movement, fracturing and collision of ice can produce seismic signatures, seismic monitoring is increasingly used to track glacial activities. Many seismological studies have been conducted in West Antarctica, whereas glacial seismicity in East Antarctica is less understood. Here, we investigate seismic data near Dålk Glacier, East Antarctica, from 6 December 2019 to 2 January 2020, using a dense seismic array and machine learning methods. We identify five categories of seismic events, which are attributed to various processes, including wind-induced vibration, thermal contraction at the ice surface, stickslip at the ice base, water-filled basal crevassing, and iceberg calving. These findings shed light on the diverse glacial activities near Dålk Glacier and their sensitivity to environmental conditions. Our study also shows the potential of unsupervised clustering as a valuable tool for glacial monitoring. [ABSTRACT FROM AUTHOR]

Linear elastic fracture mechanics (LEFM) models have been used to estimate crevasse depths in glaciers and to represent iceberg calving in ice sheet models. However, existing LEFM models assume glacier ice to be homogeneous and utilise the mechanical properties of fully consolidated ice. Using depth-invariant properties is not realistic, as the process of compaction from unconsolidated snow to firn to glacial ice is dependent on several environmental factors, typically leading to a lesser density and Young's modulus in upper surface strata. New analytical solutions for longitudinal stress profiles are derived, using depth-varying properties based on borehole data from the Ronne ice shelf, and used in an LEFM model to determine the maximum penetration depths of an isolated crevasse in grounded glaciers and floating ice shelves. These maximum crevasse depths are compared to those obtained for homogeneous glacial ice, showing the importance of including the effect of the upper unconsolidated firn layers on crevasse propagation. The largest reductions in penetration depth ratio were observed for shallow grounded glaciers, with variations in Young's modulus being more influential than firn density (a maximum difference in crevasse depth of 46 % and 20 % respectively); whereas, firn density changes resulted in an increase in penetration depth for thinner floating ice shelves (95 %–188 % difference in crevasse depth between constant and depth-varying properties). Thus, our study shows that the firn layer can increase the vulnerability of ice shelves to fracture and calving, highlighting the importance of considering depth-dependent firn-layer material properties in LEFM models for estimating crevasse penetration depths and predicting rift propagation. [ABSTRACT FROM AUTHOR]

Ice shelf calving constitutes roughly half of the total mass loss from the Antarctic ice sheet. Although much attention is paid to calving of giant tabular icebergs, these events are relatively rare. More frequent, smaller-scale calving events likely play an important role in the ice shelf frontal dynamics. Here, we investigate the role of bending stresses at the ice shelf front in driving calving on the scale 100 m – 1 km, perpendicular to the ice edge. We focus in particular on how buoyant underwater 'feet' that protrude beyond the above-water ice cliff may cause tensile stresses at the base of the ice and ultimately lead to fracture. Indirect and anecdotal observations of such feet at the Ross Ice Shelf front suggest that this process may be widespread. We consider satellite observations, together with an elastic beam model and a parameterization of frontal wave erosion to estimate the size and frequency of such calving events. Our results suggest that foot-induced mass loss at Ross Ice Shelf may cause up to 25 % of the total frontal ablation. However, stresses induced through this process are likely not sufficient to initiate crevassing but rather act to propagate existing crevasses. In addition, the relatively strong ice thickness dependence of the frontal uplift suggests an important role for internal bending moments due to temperature gradients in the ice. The highly variable environment, irregularity of pre-existing crevasse spacing, and complex rheology of the ice continue to pose challenges in better constraining the drivers behind the observed deformations and resulting calving rates. [ABSTRACT FROM AUTHOR]

Field and remote sensing studies suggest that ice mélange influences glacier-fjord systems by exerting stresses on glacier termini and releasing large amounts of freshwater into fjords. The broader impacts of ice mélange over long time scales are unknown, in part due to a lack of suitable ice mélange flow models. Previous efforts have included modifying existing viscous ice shelf models, despite the fact that ice mélange is fundamentally a granular material, and running computationally expensive discrete element simulations. Here, we draw on laboratory studies of granular materials, which exhibit viscous flow when stresses greatly exceed the yield point, plug flow when the stresses approach the yield point, and stress transfer via force chains. By implementing the nonlocal granular fluidity rheology into a depth- and width-integrated stress balance equation, we produce a numerical model of ice mélange flow that is consistent with our understanding of well-packed granular materials and that is suitable for long time-scale simulations. For parallel-sided fjords, the model exhibits two possible steady state solutions. When there is no calving of new icebergs or melting of previously calved icebergs, the ice mélange is pushed down fjord by the advancing glacier terminus, the velocity is constant along the length of the fjord, and the thickness profile is exponential. When calving and melting are included, the ice mélange evolves to another steady state in which its location is fixed relative to the fjord walls, the thickness profile is relatively steep, and the flow is extensional. For the latter case, the model predicts that the steady-state ice mélange buttressing force depends on the surface and basal melt rates through an inverse power law relationship, decays roughly exponentially with both fjord width and gradient in fjord width, and increases with the iceberg calving flux. The increase in buttressing force with the calving flux, which depends on glacier thickness, appears to occur more rapidly than the force required to prevent the capsize of full-glacier-thickness icebergs, suggesting that glaciers with high calving fluxes may be more strongly influenced by ice mélange than those with small fluxes. [ABSTRACT FROM AUTHOR]

Icebergs modulate the effective location of freshwater input from ice sheets into the ocean and therefore play an important role for the climate, especially during times of increased ice discharge (e.g., Heinrich events). None of the models participating in the Paleo Modeling Intercomparison Project simulations of the Last Glacial Maximum or the last deglaciation included icebergs. Here, we present a newly developed dynamic/thermodynamic iceberg module that was specifically designed to be incorporated in climate models used for long‐term climate simulations with interactive ice sheets. In contrast to the widely used Lagrangian iceberg models, it is formulated in an Eulerian framework. This simplifies coupling to ocean models and enhances computational efficiency for glacial climates. In a set of sensitivity experiments, where the module was implemented into an Earth System Model, we validate the model for present‐day climate conditions and test its sensitivity to key parameters. Further, we investigate the effect of iceberg hosing on the Atlantic meridional overturning circulation (AMOC) as compared to traditional freshwater hosing. Varying the hosing rate slowly in time yields a good approximation of the hysteresis curve of the AMOC. We find that the sensitivity of the AMOC to iceberg hosing is stronger than to freshwater hosing in the same ocean point, but weaker as compared to a latitude belt forcing in the North Atlantic. This emphasizes the 2e necessity to include interactive icebergs in long‐term coupled climate simulations to realistically represent melt patterns and the response of the AMOC to freshwater input from melting ice sheets. Plain Language Summary: Icebergs transport frozen water from calving ice sheets far away from the coasts and release meltwater along their path while slowly melting. The sensitivity of the Atlantic meridional overturning circulation depends on the location of the meltwater input. As the transport and melt of icebergs determine the location of meltwater input into the ocean, it is obvious that icebergs are an important component of the climate system. Unfortunately, they are not included in most climate models, partly as they are computational quite expensive. Existing iceberg modules can follow groups of icebergs, and may be prohibitively computationally expensive when applied on very long timescales. Here we introduce a different approach. We do not follow each group of icebergs but describe the iceberg distribution using equations for concentrations and different size classes. This formulation is adequate for climate purposes, and simplifies the incorporation into ocean models. The difference in the climate response to iceberg calving and the widely used freshwater hosing is relatively large. This emphasizes the importance of including interactive iceberg modules into climate models for a realistic simulation of ice discharge events. Key Points: A new Eulerian iceberg module for climate studies is introducedThe Eulerian approach simplifies the incorporation of icebergs into standard ocean general circulation modelsExperiments with a comprehensive climate model demonstrate that iceberg discharge has different impact on Atlantic meridional overturning circulation compared to freshwater [ABSTRACT FROM AUTHOR]

The mass loss of glaciers outside the polar ice sheets has been accelerating during the past several decades and has been contributing to global sea-level rise. However, many of the mechanisms of this mass loss process are not well understood, especially the calving dynamics of marine-terminating glaciers, in part due to a lack of high-resolution calving front observations. Svalbard is an ideal site to study the climate sensitivity of glaciers as it is a region that has been undergoing amplified climate variability in both space and time compared to the global mean. Here we present a new high-resolution calving front dataset of 149 marine-terminating glaciers in Svalbard, comprising 124 919 glacier calving front positions during the period 1985–2023 (10.5281/zenodo.10407266, Li et al., 2023). This dataset was generated using a novel automated deep-learning framework and multiple optical and SAR satellite images from Landsat, Terra-ASTER, Sentinel-2, and Sentinel-1 satellite missions. The overall calving front mapping uncertainty across Svalbard is 31 m. The newly derived calving front dataset agrees well with recent decadal calving front observations between 2000 and 2020 (Kochtitzky and Copland, 2022) and an annual calving front dataset between 2008 and 2022 (Moholdt et al., 2022). The calving fronts between our product and the latter deviate by 32 ± 65 m on average. The R2 of the glacier calving front change rates between these two products is 0.98, indicating an excellent match. Using this new calving front dataset, we identified widespread calving front retreats during the past four decades, across most regions in Svalbard except for a handful of glaciers draining the ice caps Vestfonna and Austfonna on Nordaustlandet. In addition, we identified complex patterns of glacier surging events overlaid with seasonal calving cycles. These data and findings provide insights into understanding glacier calving mechanisms and drivers. This new dataset can help improve estimates of glacier frontal ablation as a component of the integrated mass balance of marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

The Brunt Ice Shelf, Antarctica, accelerated rapidly from a velocity of 900 to 1500 m a -1 during 6 months, following calving of a 1500 km 2 iceberg on 22 January 2023. The immediate response to calving was observed as a change to the rate of acceleration and not to velocity directly. Acceleration increased by a factor of 10, with a second, smaller calving at the end of June 2023, leading to further tripling of acceleration. The acceleration was caused by the reduction of buttressing at the McDonald Ice Rumples, leading to high localised strain rates, which reduce the strength of the remaining ice shelf. [ABSTRACT FROM AUTHOR]

Ice calved from the Antarctic and Greenland Ice Sheets or tidewater glaciers ultimately melts in the ocean contributing to sea-level rise. Icebergs have also been described as biological hotspots due to their potential roles as platforms for marine mammals and birds, and as micronutrient fertilizing agents. Icebergs may be especially important in the Southern Ocean where availability of the micronutrients iron and manganese extensively limits marine primary production. Whilst icebergs have long been described as a source of iron to the ocean, their nutrient signature is poorly constrained and it is unclear if there are regional differences. Here we show that 589 ice fragments collected from floating ice in contrasting regions spanning the Antarctic Peninsula, Greenland, and smaller tidewater systems in Svalbard, Patagonia and Iceland have similar characteristic (micro)nutrient signatures with limited or no significant differences between regions. Icebergs are a minor or negligible source of macronutrients to the ocean with low concentrations of NOx (NO3 + NO2, median 0.51 µM), PO4 (median 0.04 µM), and dissolved Si (dSi, median 0.02 µM). In contrast, icebergs deliver elevated concentrations of dissolved Fe (dFe; mean 82 nM, median 12 nM) and Mn (dMn; mean 26 nM, median 2.6 nM). A tight correlation between total dissolvable Fe and Mn (R2 = 0.95) and a Mn:Fe ratio of 0.024 suggested a lithogenic origin for the majority of sediment present in ice. Total dissolvable Fe and Mn retained a strong relationship with sediment load (both R2 = 0.43, p<0.001), whereas weaker relationships were observed for dFe, dMn and dSi. Sediment load for Antarctic ice (median 9 mg L-1, n=144) was low compared to prior reported values for the Arctic. A particularly curious incidental finding was that melting samples of ice were observed to rapidly lose their sediment load, even when sediment layers were embedded within the ice and stored in the dark. Our results demonstrated that the nutrient signature of icebergs is consistent with an atmospheric source of NOx and PO4. Conversely, high Fe and Mn, and modest dSi concentrations, are associated with englacial sediment, which experiences limited biogeochemical processing prior to release into the ocean. [ABSTRACT FROM AUTHOR]

The Baltic Ice Stream, a large fast‐flowing sector of the last Fennoscandian Ice Sheet that occupied the present‐day Baltic Sea basin, was first conceptualized in the earliest days of glacial geological research in Scandinavia. Landform and sedimentological evidence from the terrestrial margins support the concept and numerical ice‐sheet models demonstrate its existence and possible evolution. However, with evidence for the Baltic Ice Stream thus far limited to the terrestrial periphery, its true form, scale, function, and role in deglaciation have proven enigmatic. Here we present geomorphological evidence directly from the Baltic seabed that confirms the existence of and sheds light on the behaviour of the Baltic Ice Stream. Based on an extensive collection of high‐, moderate‐ and low‐resolution bathymetric terrain models covering a large proportion of the Baltic Sea floor, and complemented by LiDAR‐data for the Baltic islands, we have identified and mapped >20 000 individual subglacial bedforms, meltwater landforms and grounding line landforms. We reconstruct a six‐stage sequence of ice flow and retreat, finding that streaming was persistent in the Baltic but that pathways were variable in extent, timing and duration: different sectors of the Baltic exhibit asynchronous streaming and out‐of‐phase grounding line changes. During deglaciation, grounding line re‐advances occurred in both the southwestern and the northern Baltic Proper, and, while abundant iceberg ploughmarks attest to calving as a significant ice loss mechanism, lobate margins suggest supply to the Baltic catchment was consistently high. Our reconstruction is limited by a fragmentary geomorphic record. Here we put forward a first hypothesis for how the Baltic Ice Stream evolved, and hope it stimulates new geomorphic, stratigraphical and core data collection to extend the landform record, provide insights into subglacial and grounding line processes, and constrain the chronology for Baltic Ice Stream flow and retreat. [ABSTRACT FROM AUTHOR]

There is currently poor scientific agreement on whether the ice–bed interface is frozen or thawed beneath approximately one third of the Greenland ice sheet. This disagreement in basal thermal state results, at least partly, from differences in the subglacial geothermal heat-flow basal boundary condition used in different ice-flow models. Here, we employ seven widely used Greenland geothermal heat-flow maps in 10 000-year spin-ups of the Community Ice Sheet Model (CISM). We perform two spin-ups: one nudged toward thickness observations and the other unconstrained. Across the seven heat-flow maps, and regardless of unconstrained or nudged spin-up, the spread in basal ice temperatures exceeds 10 ∘ C over large areas of the ice–bed interface. For a given heat-flow map, the thawed-bed ice-sheet area is consistently larger under unconstrained spin-ups than nudged spin-ups. Under the unconstrained spin-up, thawed-bed area ranges from 33.5 % to 60.0 % across the seven heat-flow maps. Perhaps counterintuitively, the highest iceberg calving fluxes are associated with the lowest heat flows (and vice versa) for both unconstrained and nudged spin-ups. These results highlight the direct, and non-trivial, influence of the heat-flow boundary condition on the simulated equilibrium thermal state of the ice sheet. We suggest that future ice-flow model intercomparisons should employ a range of basal heat-flow maps, and limit direct intercomparisons with simulations using a common heat-flow map. [ABSTRACT FROM AUTHOR]

The Greenland Ice Sheet has undergone rapid mass loss over the last four decades, primarily through solid and liquid discharge at marine‐terminating outlet glaciers. The acceleration of these glaciers is in part due to the increase in temperature of ocean water in contact with the glacier terminus. However, quantifying heat transport to the glacier through fjord circulation can be challenging due to iceberg abundance, which threatens instrument survival and fjord accessibility. Here we utilize iceberg movement to infer upper‐layer fjord circulation, as freely floating icebergs (i.e., outside the mélange region) behave as natural drifters. In the summers of 2014 and 2019, we deployed transmitting GPS units on a total of 13 icebergs in Ilulissat Icefjord, an iceberg‐rich and historically data‐poor fjord in west Greenland, to quantify circulation over the upper 0–250 m of the water column. We find that the direction of upper‐layer fjord circulation is strongly impacted by the timing of tributary meltwater runoff, while the speed of this circulation changes in concert with glacier behavior, which includes increases and decreases in glacier speed and meltwater runoff. During periods of increased meltwater runoff entering from tributary fjords, icebergs at these confluences deviated from their down‐fjord trajectory, even reversing up‐fjord, until the runoff pulse subsided days later. This study demonstrates the utility of iceberg monitoring to constrain upper‐layer fjord circulation, and highlights the importance of including tributary fjords in predictive models of heat transport and fjord circulation. Plain Language Summary: The Greenland Ice Sheet has been rapidly losing mass over the last four decades, primarily at its edges through glacier melting and iceberg calving into fjords. Warming ocean water in contact with the glacier terminus can accelerate mass loss. However, quantifying the currents that transport this warm ocean water are challenging to constrain due to the abundance of icebergs in the near‐terminus region. Here, we track freely floating icebergs, natural drifters, to infer surface circulation (0–250 m depth) in an iceberg‐rich fjord. In the summers of 2014 and 2019, we deployed GPS units on 13 icebergs in Ilulissat Icefjord, a historically data‐poor fjord in west Greenland. We find the direction of currents to be strongly impacted by tributary fjord runoff, with changes in iceberg trajectory coinciding with runoff pulses from these tributary fjords. We find the circulation speed to be most closely associated with glacier speed and meltwater runoff from the glacier at the head of Ilulissat Icefjord. This study highlights the utility of using icebergs to infer surface circulation and the importance of including tributary fjords in future circulation models. Key Points: We used 13 on‐iceberg GPS units to constrain upper‐layer (0–250 m) circulation in Ilulissat Icefjord, west GreenlandDeviations in down‐fjord iceberg trajectory coincide with tributary meltwater flux, in both location and timingThe speed of upper‐layer circulation changes in concert with glacier behavior, including glacier speed and meltwater runoff [ABSTRACT FROM AUTHOR]

Frontal ablation at tidewater glaciers, which comprises iceberg calving, submarine and subaerial melting, is a key boundary condition for numerical ice sheet models but remains difficult to measure directly in-situ. Many previous studies have quantified frontal ablation over varying spatio-temporal scales, however most use ice discharge as an approximation for frontal ablation, thereby neglecting the influence of terminus location change. Frontal ablation estimates that do account for terminus location change are spatio-temporally limited by the availability of observational data. Here, we present a processing chain to quantify frontal ablation using open-source observational data. We apply the processing chain to 49 tidewater glaciers in Greenland with reliable near-terminus bathymetry data in the BedMachine V4 dataset. Nearterminus volume change over the time period 1987 - 2020 is determined using a previously published dataset of terminus positions (TermPicks), ice thicknesses from ArcticDEM and AeroDEM, adjusted for surface elevation change over time, and bathymetry data from BedMachine v4. Assuming a vertical terminus geometry and uniform ice density, we estimate frontal ablation as the difference between mass flux towards the terminus (Mankoff et al., 2020) and mass change between consecutive observation. The frontal ablation dataset offers exciting opportunities for developing new insights into ice dynamics, including helping to improve numerical model hindcasting and projections. Lastly, we provide a processing chain that may serve as a community standard for determining frontal ablation from observational data for any tidewater glacier. [ABSTRACT FROM AUTHOR]

Open ocean areas surrounded by sea ice and maintained by ocean heat, or sensible‐heat polynyas, are linked to key ice‐sheet processes, such as ice‐shelf basal melt and ice‐shelf fracture, when they occur near ice‐shelf fronts. However, the lack of detailed multi‐year records of polynya variability prevent assessing coupling between polynya and frontal dynamics. Here, we present the first multi‐decadal polynya area record (2000–2022) at Pine Island Glacier (PIG), West Antarctica, from thermal and optical satellite imagery. We found substantial interannual variability in polynya area, with consistencies in the timing of polynya opening, maximal extent, and closing. Furthermore, the largest polynya in our record (269 km2) occurred at PIG's western margin just 68 days before iceberg B‐27 calved, suggesting that polynya size and position may influence rifting dynamics. Our new data set provides a pathway to assess coevolving polynya and frontal dynamics, demonstrating the importance of building long‐term, year‐round polynya variability records. Plain Language Summary: Persistent sensible‐heat polynyas are areas of open ocean surrounded by sea ice maintained by ocean heat. These surface features occur near ice‐shelf fronts at the coastal margins of Antarctica and therefore have the potential to impact ice‐shelf stability through heat transfer processes. However, our understanding of long‐term polynya variability remains limited due to the lack of multi‐year records documenting polynya evolution. Here, we use satellite imagery to measure polynya area near Pine Island Glacier (PIG) and build the first multi‐decadal record in Antarctica. We observed a large amount of year‐to‐year area variability from 2000 to 2022, with the largest polynya in our record (269 km2) occurring at the western edge of PIG just 68 days before a large iceberg calved from PIG. This correspondence suggests that polynya size and position may influence iceberg calving. Our new data set provides a pathway to assess potentially coupled ice and ocean processes, demonstrating the importance of building long‐term, year‐round polynya variability records. Key Points: We generated a 22 years record of polynya area at Pine Island Glacier from satellite thermal and optical imageryOur data set shows high variability in sensible‐heat polynya area (0–322 km2) from 2000 to 2022Large, marginal, and persistent sensible‐heat polynyas may reduce ice‐shelf buttressing and contribute to rift initiation and propagation [ABSTRACT FROM AUTHOR]

Calving front location plays a crucial role in studying ice–ocean interaction, mapping glacier area change, and constraining ice dynamic models. However, relying solely on visual interpretation to extract annual changes in the calving front of ice shelves is a time-consuming process. In this study, a comparative analysis was conducted on the segmentation obtained from fully convolutional networks (FCN), U-Net, and U2-Net models, revealing that U2-Net exhibited the most effective classification. Notably, U2-Net outperformed the other two models by more than 30 percent in terms of the F1 parameter. Therefore, this paper introduces an automated approach that utilizes the U2-Net model to extract the calving front of ice shelves based on a Landsat image, achieving an extraction accuracy of 58 m. To assess the model's performance on additional ice shelves in the polar region, the calving front of the Totten and Filchner ice shelves were also extracted for the past decade. The findings demonstrated that the ice velocity of the Filchner ice shelf exceeded that of the Totten ice shelf. Between February 2014 and March 2015, the majority of the calving fronts along the Filchner Ice Shelf showed an advancing trend, with the fastest-moving front measuring 3532 ± 58 m/yr. [ABSTRACT FROM AUTHOR]

Many floating ice shelves in Antarctica buttress the ice streams feeding them, thereby reducing the discharge of icebergs into the ocean. The rate at which ice shelves calve icebergs and how fast they flow determine whether they advance, retreat, or remain stable, exerting a first-order control on ice discharge. To parameterize calving within ice sheet models, several empirical and physical calving "laws" have been proposed in the past few decades. Such laws emphasize dissimilar features, including along- and across-flow strain rates (the eigencalving law), a fracture yield criterion (the von Mises law), longitudinal stretching (the crevasse depth law), and a simple ice thickness threshold (the minimum thickness law), among others. Despite the multitude of established calving laws, these laws remain largely unvalidated for the Antarctic Ice Sheet, rendering it difficult to assess the broad applicability of any given law in Antarctica. We address this shortcoming through a set of numerical experiments that evaluate existing calving laws for 10 ice shelves around the Antarctic Ice Sheet. We utilize the Ice-sheet and Sea-level System Model (ISSM) and implement four calving laws under constant external forcing, calibrating the free parameter of each of these calving laws for each ice shelf by assuming that the current position of the ice front is in steady state and finding the set of parameters that best achieves this position over a simulation of 200 years. We find that, in general, the eigencalving and von Mises laws best reproduce observed calving front positions under the steady-state position assumption. These results will streamline future modeling efforts of Antarctic ice shelves by better informing the relevant physics of Antarctic-style calving on a shelf-by-shelf basis. [ABSTRACT FROM AUTHOR]

Seismicity at glaciers, ice sheets, and ice shelves provides observational constraint on a number of glaciological processes. Detecting and locating this seismicity, specifically icequakes, is a necessary first step in studying processes such as basal slip, crevassing, imaging ice fabric, and iceberg calving, for example. Most glacier deployments to date use conventional seismic networks, comprised of seismometers distributed over the entire area of interest. However, smaller-aperture seismic arrays can also be used, which are typically sensitive to seismicity distal from the array footprint and require a smaller number of instruments. Here, we investigate the potential of arrays and array-processing methods to detect and locate subsurface microseismicity at glaciers, benchmarking performance against conventional seismic-network-based methods for an example at an Antarctic ice stream. We also provide an array-processing recipe for body-wave cryoseismology applications. Results from an array and a network deployed at Rutford Ice Stream, Antarctica, show that arrays and networks both have strengths and weaknesses. Arrays can detect icequakes from further distances, whereas networks outperform arrays in more comprehensive studies of a particular process due to greater hypocentral constraint within the network extent. We also gain new insights into seismic behaviour at the Rutford Ice Stream. The array detects basal icequakes in what was previously interpreted to be an aseismic region of the bed, as well as new icequake observations downstream and at the ice stream shear margins, where it would be challenging to deploy instruments. Finally, we make some practical recommendations for future array deployments at glaciers. [ABSTRACT FROM AUTHOR]

Sea level rise (SLR) due to surface melt and to dynamic losses from the ice sheets—that is via accelerated flow of glaciers into the sea—is something that may be potentially mitigated by cooling the ice sheet and oceans via solar geoengineering. We use two ice dynamic models driven by changes in surface mass balance (SMB) from four climate models to estimate the SLR contribution from the Greenland ice sheet under the Intergovernmental Panel on Climate Change (IPCC) Representative Concentration Pathway (RCP) 4.5, and 8.5, and Geoengineering Model Intercomparison Project G4 scenarios. The G4 scenario adds 5 Tg/yr sulfate aerosols to the equatorial lower stratosphere (equivalent of 1/4 the 1991 Mt Pinatubo SO2 eruption) to the IPCC RCP4.5 scenario, which itself approximates the greenhouse gas emission commitments agreed in Paris in 2015. Over the 2020–2090 period, mass loss under G4 is about 31%–38% that under RCP4.5, which is 36%–48% lower than under RCP8.5. Ice lost across the grounding line under both G4 and RCP4.5 is reduced in the future as the termini of many southeast Greenland outlets retreat onto bedrock above sea level. Glaciers with large low‐lying catchments in the west, north, and northeast of Greenland (e.g., Jakobshavn, 79N, Zachariae Isstrøm, and Petermann glaciers) discharge more ice from the ice‐sheet interior under RCP4.5 than under G4. Although calving losses vary much more than the SMB difference between ice dynamic models, both models point to significant ice discharge losses of between 15% and 42% across the scenarios. Plain Language Summary: Sea level rise from Greenland over the next century may be around 10 cm by 2100. But the amount, and the long‐term stability of the ice sheet depend on the degree of summer warming it experiences. Limiting greenhouse gas emissions to the levels pledged by states under the 2015 Paris agreement cuts ice sheet losses by 1/3–1/2 of the losses under the business‐as usual‐emissions scenario. If further cooling is induced by aerosols put into the stratosphere at a rate of about 1/4 of the 1991 eruption of Mt Pinatubo, then the ice sheet loss is reduced by about 30% compared with Paris emissions. This specific aerosol geoengineering scenario maintains both the deep and fast‐flowing glaciers and the smaller mountain glaciers closer to present sizes than they would be under a greenhouse gas emission scenario similar to international pledges in the 2015 Paris agreement. Iceberg calving remains the most difficult to quantify aspect of Greenland ice loss. Key Points: Stratospheric aerosol injection at the rate of 5 Tg/yr (G4) lowers Greenland mass loss relative to RCP4.5 by 31%–38% by 2090Across four Earth System Model and two ice dynamic models (G4–RCP4.5) differences are 34%–40% in surface mass balance and 16%–34% in ice dischargeDynamic mass loss by calving from glaciers is the largest uncertainty between ice dynamics models [ABSTRACT FROM AUTHOR]

The mass loss of glaciers outside the polar ice sheets has been accelerating during the past several decades and has been contributing to global sea-level rise. However, many of the mechanisms of this mass loss process are not well understood, especially the calving dynamics of marine-terminating glaciers, in part due to a lack of high-resolution calving front observations. Svalbard is an ideal site to study the climate sensitivity of glaciers as it is a region that has been undergoing amplified climate variability in both space and time compared to the global mean. Here we present a new high-resolution calving front dataset of 149 marine-terminating glaciers in Svalbard, comprising 124919 glacier calving front positions during the period of 1985–2023 (https://doi.org/10.5281/zenodo.8399899) (Li et al., 2023). This dataset was generated using a novel automated deep learning framework and multiple optical and SAR satellite images from Landsat, Terra-ASTER, Sentinel-2, and Sentinel-1 satellite missions. The overall calving front mapping uncertainty across Svalbard is 46 ± 21 m. The newly derived calving front dataset agrees well with recent decadal calving front observations between 2000 and 2020 (Kochtitzky and Copland, 2022) and an annual calving front dataset between 2008 and 2022 (Moholdt et al., 2022). The R² of the glacier calving front change rates between our product and the latter is 0.98, indicating an excellent match. Using this new calving front dataset, we identified widespread calving front retreats during the past three decades, across most regions in Svalbard except for a handful of glaciers draining the ice caps Vestfonna and Austfonna on Nordaustlandet. In addition, we identified complex patterns of glacier surging events overlaid with seasonal calving cycles. These data and findings provide insights into understanding glacier calving mechanisms and drivers. This new dataset can help improve estimates of glacier frontal ablation as a component of the integrated mass balance of marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

The article focuses on quantifying the annual mass budget of all Antarctic ice shelves from 1997 to 2021. The study reveals that out of 162 ice shelves, 71 lost mass, 29 gained mass, and 62 remained stable and Basal melting was the primary contributor to mass loss, impacting 68 percent of the ice shelves studied, and the research emphasizes the importance of understanding ice shelf mass changes for predicting the impact on grounded ice and global sea levels.

Deep Sea Drilling Program (DSDP) Site 548 was cored in 1984 at a water depth of 1256 m on the Goban Spur, offshore southwest Ireland. Coring retrieved a ~100‐m‐thick Pleistocene contourite sequence. This study uses planktonic foraminiferal assemblage and benthic foraminiferal oxygen isotope analyses to establish an age model for the upper 40 m of this core. This site's multidisciplinary analyses of planktonic foraminiferal assemblages, lithic grains, facies and calcium carbonate concentration reveal a 250 000‐year record of the North Atlantic polar front variability and British–Irish Ice Sheet (BIIS) history. The sequence is characterized by alternations of ice rafted debris (IRD) laden pelagic mud facies with calcium carbonate‐rich silty sand contourite facies that track glacial/interglacial cycles. The polar front migrated southward across the area several times during glacial maxima and stadial periods, while warmer Mediterranean Outflow Water (MOW) flowed northward across the region during interglacial and interstadial periods depositing contourites. Lithic analyses reveal a complex history of IRD deposition associated with iceberg calving from the Laurentide Ice Sheet and northwest European ice sheets, mainly the BIIS. Comparison between the Goban Spur (DSDP Site 548) and the Celtic Margin (MD03‐2692) and central North Atlantic Integrated Ocean Drilling Program (IODP) Site U1308 suggests differences between the 'non‐Laurentide Ice Sheet' Heinrich Events (HE) 6 and 3 at the Goban Spur, with IRD from the BIIS being prominent during HE 6 and IRD from other European ice sheets north of the BIIS likely being more dominant during HE 3. The nature of lithics in IRD‐rich horizons during Terminations 3, 3A, 2 and 1 suggests significant iceberg calving episodes preceding BIIS retreat during the onset of interstadial intervals. [ABSTRACT FROM AUTHOR]

Climate shifts are particularly conspicuous in glaciated areas. Satellite and terrestrial observations show significant increases in the melting and breakup of tidewater glaciers and their influence on sea level rise and ocean mixing. Increasing melt rates are creating an urgency to better understand the link between atmospheric and oceanic conditions and glacier frontal ablation through iceberg calving and melting. Elucidating this link requires a combination of short- and long-timescale measurements of terminus activity. Recent work has demonstrated the potential of using underwater sound to quantify the time and scale of calving events to yield integrated estimates of ice mass loss. Here, we present estimates of subaerial calving flux using underwater sound recorded at Hansbreen, Svalbard, in September 2013 combined with an algorithm for the automatic detection of calving events. The method is compared with ice calving volumes estimated from geodetic measurements of the movement of the glacier terminus and an analysis of satellite images. The total volume of above-water calving during the 26 d of acoustical observation is estimated to be 1.7±0.7×107 m 3 , whereas the subaerial calving flux estimated by traditional methods is 7±2×106 m 3. The results suggest that passive cryoacoustics is a viable technique for long-term monitoring of mass loss from marine-terminating glaciers. [ABSTRACT FROM AUTHOR]

Two composite sedimentary sequences sampled in the ice‐proximal (12CS) and ice‐distal (02CS) areas of Coronation Fjord (Baffin Island, Nunavut, Canada) were investigated in order to reconstruct the effect of climate variability on 600 years of changes in sediment transfer from the eastern Penny Ice Cap (PIC). Detrital proxies, and physical and sedimentological analyses revealed that glacial meltwater discharges led to frequent rapidly deposited layers (RDLs) in ice‐proximal settings. RDLs in ice‐distal settings involved the sudden release of a large quantity of sediment‐laden water during floods probably originating from adjacent fjords with large sandur deltas. Laminated sediments with ice‐rafted debris throughout the Little Ice Age interval in the ice‐proximal environment suggest that colder conditions promoted glacier growth, leading to successive episodes of turbid hyperpycnal meltwater plumes and iceberg calving in Coronation Fjord. Since 1850 ce, the accelerated Coronation retreat in response to modern warming has led to increased sedimentation rates, abrupt mineralogical and grain size proxy variations and more frequent RDLs. Similar trends between the detrital proxies of the ice‐proximal core and Atlantic Multidecadal Oscillation record and Arctic surface air temperature suggest high connectivity between atmospheric and sea surface temperature variations and PIC dynamics over the last 600 years. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The Brunt Ice Shelf, Antarctica, accelerated rapidly from 900 m a-1 to 1500 m a-1 during six months following the calving of a 1500 km2 iceberg on 22nd January 2023. Initially, the rate of acceleration increased by a factor of ten, with a second, smaller calving at the end of June 2023 leading to further tripling of acceleration. The acceleration is caused by reduction of buttressing at the McDonald Ice Rumples due to loss of contact with the sea floor and has led to high strain rates to the south, with potential consequences for the stability of the remaining ice shelf. [ABSTRACT FROM AUTHOR]

Knowledge of frontal ablation from marine-terminating glaciers (i.e., mass lost at the calving face) is critical for constraining glacier mass balance, improving projections of mass change, and identifying the processes that govern frontal mass loss. Here, we discuss the challenges involved in computing frontal ablation and the unique issues pertaining to both glaciers and ice sheets. Frontal ablation estimates require numerous datasets, including glacier terminus area change, thickness, surface velocity, density, and climatic mass balance. Observations and models of these variables have improved over the past decade, but significant gaps and regional discrepancies remain, and better quantification of temporal variability in frontal ablation is needed. Despite major advances in satellite-derived large-scale datasets, large uncertainties remain with respect to ice thickness, depth-averaged velocities, and the bulk density of glacier ice close to calving termini or grounding lines. We suggest ways in which we can move toward globally complete frontal ablation estimates, highlighting areas where we need improved datasets and increased collaboration. [ABSTRACT FROM AUTHOR]

Close-range sensors are employed to observe glaciological processes that operate over short timescales (e.g. iceberg calving, glacial lake outburst floods, diurnal surface melting). However, under poor weather conditions optical instruments fail while the operation of radar systems below 17 GHz do not have sufficient angular resolution to map glacier surfaces in detail. This letter reviews the potential of millimetre-wave radar at 94 GHz to obtain high-resolution 3-D measurements of glaciers under most weather conditions. We discuss the theory of 94 GHz radar for glaciology studies, demonstrate its potential to map a glacier calving front and summarise future research priorities. [ABSTRACT FROM AUTHOR]