Iceberg Calving Research Articles

Feedbacks Between Fjord Circulation, Mélange Melt, and the Subglacial Discharge Plume at Kangerlussuaq Glacier, East Greenland

AbstractIn recent decades, the Greenland ice sheet has been losing ice and contributing substantially to global sea level rise. Approximately half of this recent loss is due to glacier acceleration, increasing the calving of icebergs into the ocean. This process has been linked with increased ocean heat content on the continental shelf, yet the pathways delivering this heat into Greenland's fjords and its interactions with fjord‐scale processes modulating glacier ice loss are still unclear. In this study, we use a series of numerical ocean model configurations to examine feedbacks between ocean circulation, subglacial discharge, submarine glacier melt, and ice mélange in Kangerlussuaq Fjord—a major fjord system where Greenland's third‐largest glacier terminates. We find that subglacial discharge is a major control on ocean properties, increasing the up‐fjord advection of deep warm water more than 10‐fold over fjords without discharge and modulating ocean temperature on the continental shelf near the fjord mouth. Further, discharge‐driven upwelling increases ice mélange melt 3‐fold, revealing that subglacial discharge is an important control on mélange melt, particularly in the summer when submarine glacier melt and subsequent glacier retreat is highest. These results suggest that subglacial plume activity contributes to the strong correlation between mélange thickness and retreat noted in previous studies and may contribute to extensive future retreat at Kangerlussuaq Glacier.

Wave erosion, frontal bending, and calving at Ross Ice Shelf

Abstract. 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. Here, we investigate the role of frontal melting and stresses at the ice shelf front in driving bending and calving on the scale of ∼100 m, 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. Indirect and anecdotal observations of such feet at the Ross Ice Shelf front suggest that the resulting bending may be widespread and can trigger calving. We consider satellite observations together with an elastic beam model and a parameterization of wave erosion to better understand the dynamics at the ice shelf front. Our results suggest that on average frontal ablation rather consistently accounts for 20±5 m yr−1 of ice loss at Ross Ice Shelf, likely mostly due to wave erosion and smaller-scale, 𝒪(100 m), foot-induced calving. This constitutes only ∼2 % of the total frontal mass loss (since near-front ice velocities are ∼1000 m yr−1). Observational evidence suggests that sporadic larger events can skew this rate (we document one foot-induced calving event of size ∼1 km). Stresses from foot-induced bending are likely not sufficient to initiate crevassing but rather act to propagate existing crevasses. In addition, our results support recent findings by Buck (2024) that additional bending moments, likely due to temperature gradients in the ice, may play a role in driving frontal deflections. 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.

Seasonal changes of mélange thickness coincide with Greenland calving dynamics

Iceberg calving is a major contributor to Greenland’s ice mass loss. Ice mélange, tightly packed sea ice and icebergs, has been hypothesized to buttress the calving fronts. However, quantifying the mélange buttressing force from field observations remains a challenge. Here we show that such quantification can be achieved with a single field measurement: thickness of mélange at the glacier terminus. We develop a three-dimensional discrete element model of mélange along with a simple analytical model to quantify the mélange buttressing using mélange thickness data from ArcticDEM over 32 Greenland glacier termini. We observed a strong seasonality in mélange thickness: thin mélange (averaged thickness 34−15+17 m) in summertime when terminus retreats, and thick mélange (averaged thickness 119−37+31 m) in wintertime when terminus advances. The observed seasonal changes of mélange thickness strongly coincide with observed Greenland calving dynamics and the modeled buttressing effects.

A quasi-one-dimensional ice mélange flow model based on continuum descriptions of granular materials

Abstract. 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 timescales 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 exhibit 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-timescale simulations. For parallel-sided fjords, the model exhibits two possible steady-state solutions. When there is no calving of 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 and treated as constants, the ice mélange evolves into 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 buttressing force appears to increase with calving flux (i.e., glacier thickness) more rapidly than the force required to prevent the capsizing 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.

The influence of firn layer material properties on surface crevasse propagation in glaciers and ice shelves

Abstract. 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 utilize 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 lower 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 are 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 the penetration depth ratio were observed for shallow grounded glaciers, with variations in Young's modulus being more influential than firn density (maximum differences 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.

Accelerated Basal Melt Rates of Ice Shelves in North Greenland From 2013 to 2022 Estimated With the High‐Resolution ArcticDEM

AbstractBasal melting of ice shelves has become one of the main causes of mass loss from the Greenland ice sheet. However, most studies have focused on individual ice shelves, making it difficult to gain a more comprehensive understanding of basal melting across Greenland ice shelves. To address this issue, we utilized timestamped ArcticDEM strip data coregistered with ICESat‐2 data to estimate the basal melt rates of the ice shelves in North Greenland at a resolution of 150 m from 2013 to 2022, employing a mass conservation approach within the Lagrangian framework. Additionally, to investigate the influence of temperature on basal melt rates, a basic analysis correlating the basal melt rates with temperatures was conducted. Overall, the mass loss caused by basal melting of the six ice shelves has amounted to 27.86 ± 35.63 Gt yr−1, accounting for approximately 90% of the non‐calving mass loss, equivalent to a sea level rise of 0.08 ± 0.10 mm yr−1, far exceeding surface mass loss and glacier calving. The two larger ice shelves, Petermann and 79° North (79N), have contributed to 85% of the basal melt mass loss. Regarding the spatiotemporal distribution, the basal melt rates have gradually decreased from near the grounding line to the ice shelf front. Apart from the Ryder ice shelves, the basal melting of the other ice shelves is in a state of accelerated ablation. Moreover, compared to the skin temperature of the ice shelf, the sea water potential temperature has a greater impact on the basal melt rate.

A 3.3‐Million‐Year Record of Antarctic Iceberg Rafted Debris and Ice Sheet Evolution Quantified by Machine Learning

AbstractOver the last 3.3 million years, the Antarctic Ice Sheet (AIS) has undergone phases of ice sheet growth and decay, impacting sea level and climate globally. Presently, the largely marine‐terminating AIS loses mass primarily by iceberg calving and basal melt of ice shelves. Quantifying past rates and timing of AIS melt is vital to understanding future cryosphere and sea level changes. One proxy for past ice sheet instabilities is iceberg rafted debris (IRD) fluxes. However, traditional methods of IRD quantification are labor‐intensive. Here, we present a new method of identifying IRD grains in sediment core X‐ray images using a convolutional neural network machine learning algorithm. We present a 3.3‐million‐year record of AIS IRD melt events using sediment cores from International Ocean Discovery Program Sites U1536, U1537, and U1538 in the Southern Ocean's “Iceberg Alley.” We identify two increases in the IRD fluxes throughout this period, at ∼1.8 and 0.43 Ma. We propose that after 1.8 Ma, the AIS expanded and transitioned from a primarily terrestrial‐terminating to a primarily marine‐terminating ice sheet. Therefore, after 1.8 Ma, glacial terminations and AIS iceberg discharge are associated with variations in global ice volume, presumably through the mechanism of sea level and, therefore, grounding line change. The second AIS regime change occurs during the Mid‐Brunhes Event (∼0.43 Ma). After this time, there are heightened and continuous IRD fluxes at each glacial termination, indicating increased AIS size and instability after this time.

Calving front monitoring at a subseasonal resolution: a deep learning application for Greenland glaciers

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

A Boulder Beach Formed by Waves From a Calving Glacier Revisited: Multidecadal Tsunami–Controlled Coastal Changes in Front of Eqip Sermia, West Greenland

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

Heinrich event ice discharge and the fate of the Atlantic Meridional Overturning Circulation.

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.

A large-scale high-resolution numerical model for sea-ice fragmentation dynamics

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

A 3D glacier dynamics–line plume model to estimate the frontal ablation of Hansbreen, Svalbard

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

Environment‐Modulated Glacial Seismicity Near Dålk Glacier in East Antarctica Revealed by Deep Clustering

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

Uncrewed aerial vehicle with onboard winch system for rapid, cost-effective, and safe oceanographic profiling in hazardous and inaccessible areas

Interactions between coastal waters and marine-terminating glaciers in the Polar Regions play a significant role in global sea level rise fueled by a rapidly warming Arctic.The risk of glacier calving, and the abundance of ice, can make it impossible for surface vessels to access the waters near glacier termini. Alternative methods using manned aircraft are expensive. As a result, oceanographic measurements are limited near glacier termini.We present an uncrewed aerial vehicle (UAV) with an on-board winch system that allows oceanographic profiling in remote, hazardous areas using a commercial conductivity, temperature, and depth (CTD) sensor payload. The UAV is optimized for easy handling and deployment and is capable of high-speed and efficient cruise flight. An autopilot system provides pilot assistance and autonomous flight capabilities. The total weight of the UAV including payload is 6.5 kg with an endurance of 24 min.Testing of the system was conducted in South Greenland during winter conditions in March 2023 with successful profiles collected near a glacier terminus (<5 m) and in small openings in ice mélange (2.2 m). The system proved capable, reliable, and efficient. Further development of the system will allow other sensors for an even more flexible measurement suite.

Discrimination between icequakes and earthquakes in southern Alaska: an exploration of waveform features using Random Forest algorithm

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

Stability of Ice Shelves and Ice Cliffs in a Changing Climate

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.

Sedimentary stratigraphy and provenance off Dronning Maud Land (East Antarctica) during the mid-Pleistocene transition: Implications for paleoclimate and ice dynamics

Glacial-marine sediments from the Antarctic continental margin provide a record of depositional environment, oceanographic variability and ice dynamics that is tapped with scientific ocean drilling. This study focuses on Ocean Drilling Program Core 693A-2R, a 9.7 m sediment core retrieved from near the continental margin of the Archean Grunehogna Craton in Dronning Maud Land (DML), East Antarctica. The results contribute to a better understanding of ice-shelf behavior in DML during the mid-Pleistocene transition (MPT), a well-known transition from 40-kyr to 100-kyr cycle periods. The age model, constructed based on Sr isotope stratigraphy and geomagnetic reversals, indicates that the core spans 1.20 to 0.65 Ma. The dynamic behavior of DML ice shelves with periodic iceberg calving is revealed by the glacial–interglacial variation in sedimentation patterns, with interglacials characterized by higher concentrations of ice-rafted debris (IRD) associated with enhanced paleoproductivity than glacial intervals. The responses of DML ice shelves to warm climates are represented by a prolonged interglacial period at 1.0–1.1 Ma (MIS 31–27) and significant interglacial expressions during MIS 19 and 17. The 40Ar/39Ar ages of individual ice-rafted hornblende grains are compared with the on-land geology of DML and neighboring regions to determine the provenances of IRD. Specifically, 40Ar/39Ar results record primarily late Neoproterozoic to Cambrian ages (600–400 Ma) with a predominant peak of 520–480 Ma. This Pan-African/Ross orogeny signature is very common in East Antarctica but is not found in the most proximal margin of the Grunehogna Craton, and is instead associated with the region of DML several hundred kilometers east of the deposition site. This indicates that significant discharges of icebergs occurred in the remote DML, which were then transported by the westward-flowing Antarctic Coastal Current to deposit IRD at the studied site during the MPT. This study establishes a confirmed MPT sedimentary sequence off DML, against which future MPT proxy records from the Weddell Sea embayment and other sectors in Antarctica can be compared and correlated, and provides a basis for more detailed analyses of the response of DML ice sheet to Pleistocene climate variations.

Array processing in cryoseismology: a comparison to network-based approaches at an Antarctic ice stream

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

Monitoring glacier calving using underwater sound

Abstract. 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 (Glowacki and Deane, 2020). 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 m3, whereas the subaerial calving flux estimated by traditional methods is 7±2×106 m3. The results suggest that passive cryoacoustics is a viable technique for long-term monitoring of mass loss from marine-terminating glaciers.

Annual mass budget of Antarctic ice shelves from 1997 to 2021.

Antarctic ice shelves moderate the contribution of the Antarctic Ice Sheet to global sea level rise; however, ice shelf health remains poorly constrained. Here, we present the annual mass budget of all Antarctic ice shelves from 1997 to 2021. Out of 162 ice shelves, 71 lost mass, 29 gained mass, and 62 did not change mass significantly. Of the shelves that lost mass, 68 had statistically significant negative mass trends, 48 lost more than 30% of their initial mass, and basal melting was the dominant contributor to that mass loss at a majority (68%). At many ice shelves, mass losses due to basal melting or iceberg calving were significantly positively correlated with grounding line discharge anomalies; however, the strength and form of this relationship varied substantially between ice shelves. Our results illustrate the utility of partitioning high-resolution ice shelf mass balance observations into its components to quantify the contributors to ice shelf mass change and the response of grounded ice.