Search

At least three primary iceberg calving styles have been identified in Greenland: serac collapse, which produces falling icebergs tens of meters in length; slab capsize, which produces rotating icebergs hundreds of meters in length; and tabular rifting, which produces kilometer‐scale icebergs. However, calving styles are mostly undocumented across Greenland. Here, we develop a method to disentangle the sizes of individual calving events and map the dominant calving style at glaciers, using the characteristic properties of step retreats in satellite‐derived terminus positions. At glaciers known to frequently produce calving teleseisms, step retreats greater than 200 m account for >80% ${ >} 80\%$ of net calved length since 2018. In contrast, at glaciers known to calve by serac failure, 200 m step retreats account <20% ${< } 20\%$ of net calving. Thus, terminus change timeseries can offer promising insight into the dominant calving styles at marine‐terminating glaciers. Plain Language Summary: Iceberg calving is ubiquitous at the ocean‐bounded edges of ice sheets. However, not all iceberg calving is the same: at least three distinct mechanisms govern the discharge of solid ice from the Greenland Ice Sheet. We show these are associated with a characteristic range of iceberg sizes. Because calving constitutes a major component of Greenland's continued mass loss, it is important that these distinct mechanisms be understood. We develop a method to identify which calving mechanisms are important at various glaciers. We do so by using satellite observations to estimate the lengths of terminus retreats producing icebergs at individual glaciers. Key Points: Greenland outlet glaciers with different calving styles have different characteristic retreat magnitudes in consecutive satellite imagesAt glaciers which frequently produce calving teleseisms, step retreats in terminus position are dominated by retreats >200 mThe distinct processes underlying different calving styles will likely necessitate several distinct calving laws [ABSTRACT FROM AUTHOR]

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs <4 km length and 60 days for larger icebergs. Complete shatter producing many icebergs <1 km length is rare. A side attrition rate of 0.23 m d−1 combined with drift speed of 6 km d−1, or any proportional change in both numbers fits the observed changes in iceberg distribution. The largest injection into the Southern Ocean of fresh water and any iceberg-transported material takes place in a ~2.3 × 10⁶ km2 zone extending east-northeast from the Antarctic Peninsula to the Greenwich meridian. The iceberg contribution to salinities and temperatures, with maximum contribution north of the Weddell Sea, differs in some regions, from those indicated by tracking large icebergs. [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]

Abstract At least three primary iceberg calving styles have been identified in Greenland: serac collapse, which produces falling icebergs tens of meters in length; slab capsize, which produces rotating icebergs hundreds of meters in length; and tabular rifting, which produces kilometer‐scale icebergs. However, calving styles are mostly undocumented across Greenland. Here, we develop a method to disentangle the sizes of individual calving events and map the dominant calving style at glaciers, using the characteristic properties of step retreats in satellite‐derived terminus positions. At glaciers known to frequently produce calving teleseisms, step retreats greater than 200 m account for >80% of net calved length since 2018. In contrast, at glaciers known to calve by serac failure, 200 m step retreats account

Glacier contributions to streamflow moderate summer water temperatures and help keep rivers cool during periods of hot, dry weather. As glaciers retreat, the length of channel exposed to solar radiation above a given point of interest would increase, which should ultimately result in higher water temperatures. However, formation of proglacial lakes may complicate river temperature response to glacier retreat. This study examined river and lake temperature response to the retreat of lake‐terminating Bridge Glacier, British Columbia, using a combination of manual spot and continuous river temperatures at a stream gauge, lake surface temperatures derived from Landsat imagery, and in situ lake temperatures made during a field campaign during summer 2013. Temperature data spanned the period from 1984 to 2022. Lake and river temperatures were relatively constant prior to and during a period of rapid retreat, but exhibited an upward trend following 2012, when iceberg density declined due to an apparent decrease in calving rate, resulting in an increase in lake surface area and river temperature of about 5°C even though the glacier remained lake terminating. In addition to providing the first longitudinal study of water temperature response to the retreat of a lake‐terminating glacier, this study adds to a body of literature that demonstrates that routine spot river temperatures, if collected in a consistent manner over a number of years, can provide valid information about river and stream temperature trends. [ABSTRACT FROM AUTHOR]

This study employs the two-dimensional nonlinear Stokes ice sheet model and a Galerkin least-squares (GLS) finite element method to investigate iceberg calving at the terminus of tidewater glaciers. We propose an approach based on pressure and normal stress solutions to adjust the grounding line position and present effective principal stress contours and profiles for grounded and notch glaciers with basal crevasses at the grounding line. Our results indicate that the openings of these basal crevasses are significantly affected by water pressure. In addition, stress profiles in ungrounded tidewater glaciers vary from those in fully grounded tidewater glaciers, which could affect iceberg calving. We also conduct numerical experiments to analyze the effects of slip length, notch length, and surface slope and examine the effectiveness of the GLS method in numerical solutions. Our results are in agreement with prior findings in the literature that basal crevasses are significantly affected by water pressure, and stress profiles are significantly different in grounded and ungrounded tidewater glaciers. [ABSTRACT FROM AUTHOR]

A major consideration for maritime activity in the Southern Hemisphere is the northern limit of icebergs, or the Southern Ocean Limit Of Known Ice (SOLOKI). This analysis of historical reports of icebergs during Southern Hemisphere voyages from 1687 to 1933 provides a basis for examination of their geographical and chronological occurrence during ~250 years. The analyses use tabulated data from 742 voyages and other reports from many sources, some including first-person descriptions. While these data are dependent on icebergs being reported by mariners, as well as the variable frequency of voyages, they demonstrate distinct periods of exceptional frequency of icebergs occurring in certain localities, particularly the far South Atlantic. Based upon historical records the evidence suggests unprecedented numbers of icebergs were present in southern shipping channels in the 1890s. When these historical observations are combined with modern iceberg drift trajectories, their possible origin can be elucidated. Owing to the numbers of icebergs seen and their geographical spread, our results suggest that this was possibly the largest near-synchronous calvings in the last 300 years, and the northernmost extent of the SOLOKI.

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

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

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs

Abstract The Antarctic ice sheet is buttressed by floating ice shelves that calve icebergs along large fractures called rifts. Despite the significant influence exerted by rifting on ice shelf geometry and buttressing, the scarcity of in situ observations of rift propagation contributes considerable uncertainty to understanding rift dynamics. Here, we report the first‐ever seismic recording of a multiple‐kilometer rift propagation event. Remote sensing and seismic recordings reveal that a rift in the Pine Island Glacier Ice Shelf extended 10.53 km at a speed of 35.1 m/s, the fastest known ice fracture at this scale. We simulate ocean‐coupled rift propagation and find that the dynamics of water flow within the rift limit the propagation rate, resulting in rupture two orders of magnitude slower than typically predicted for brittle fracture. Using seismic recordings of the elastic waves generated during rift propagation, we estimate that ocean water flows into the rift at a rate of at least 2,300 m3/s during rift propagation and causes mixing in the subshelf cavity. Our observations support the hypotheses that large ice shelf rift propagation events are brittle, hydrodynamically limited, and exhibit sensitive coupling with the surrounding ocean.

A major consideration for maritime activity in the Southern Hemisphere is the northern limit of icebergs, or the Southern Ocean Limit Of Known Ice (SOLOKI). This analysis of historical reports of icebergs during Southern Hemisphere voyages from 1687 to 1933 provides a basis for examination of their geographical and chronological occurrence during ~250 years. The analyses use tabulated data from 742 voyages and other reports from many sources, some including first-person descriptions. While these data are dependent on icebergs being reported by mariners, as well as the variable frequency of voyages, they demonstrate distinct periods of exceptional frequency of icebergs occurring in certain localities, particularly the far South Atlantic. Based upon historical records the evidence suggests unprecedented numbers of icebergs were present in southern shipping channels in the 1890s. When these historical observations are combined with modern iceberg drift trajectories, their possible origin can be elucidated. Owing to the numbers of icebergs seen and their geographical spread, our results suggest that this was possibly the largest near-synchronous calvings in the last 300 years, and the northernmost extent of the SOLOKI. [ABSTRACT FROM AUTHOR]

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

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

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

At any one time 130 000 icebergs are afloat in the Southern Ocean; 97% of these are too small to be registered in current satellite-based databases, yet the melting of these small icebergs provides a major input to the Southern Ocean. We use a unique set of visual size observations of 53 000 icebergs in the South Atlantic Ocean, the SCAR International Iceberg Database, to derive average iceberg dissolution rates. Fracture into two parts is the dominant dissolution process for tabular icebergs, with an average half-life of 30 days for icebergs <4 km length and 60 days for larger icebergs. Complete shatter producing many icebergs <1 km length is rare. A side attrition rate of 0.23 m d−1 combined with drift speed of 6 km d−1, or any proportional change in both numbers fits the observed changes in iceberg distribution. The largest injection into the Southern Ocean of fresh water and any iceberg-transported material takes place in a ~2.3 × 10⁶ km2 zone extending east-northeast from the Antarctic Peninsula to the Greenwich meridian. The iceberg contribution to salinities and temperatures, with maximum contribution north of the Weddell Sea, differs in some regions, from those indicated by tracking large icebergs. [ABSTRACT FROM AUTHOR]

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

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

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

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

Along the lateral margins of floating ice shelves in Greenland and Antarctica, ice flow past confining margins and pinning points is often accompanied by extensive rifting. Rifts in zones of marginal decoupling ("detachment zones") typically propagate inward from the margins and result in many of Earth's largest calving events. Velocity maps of detachment zones indicate that flow through these regions is spatially transitioning from confined to unconfined ice shelf flow. We employ the software package icepack to demonstrate that longitudinally decreasing marginal resistance reproduces observed transitions in flow regime, and we show that these spatial transitions are accompanied by near‐margin tension sufficient to explain full‐thickness rifts. Thus, we suggest that zones of progressive decoupling are a primary control on ice shelf calving. The steadiness of detachment zone positions may be a good indicator of ice shelf vulnerability, with migratory or thinning detachment zones indicating shelves at risk of dynamic speedup and increased fracture. Plain Language Summary: The massive icebergs released from the ice shelves of Greenland and Antarctica all originate from fractures in the ice. Fracture occurs where extensional stresses are great enough to break the ice. We observe that zones of extensive fracture tend to coincide with regions where shelf ice flows just beyond lateral confinements such as fjord walls or islands, and that the fractures accrued through such zones often result in icebergs. In this study, we show that the high stresses within these zones occur as a result of a loss of contact with features at ice shelf edges. We propose that the regions where sidewall or island contact is lost may be an important indicator of future iceberg calving and changes in larger scale ice shelf stability. Key Points: Many calving events originate as near‐margin rifts where ice shelves flow beyond lateral obstructions or boundariesPhysical detachment from lateral obstructions is a major source of near‐margin tension and fractureWithin shelf‐marginal detachment zones, thinning/changes to fracture patterns may presage ice shelf destabilization [ABSTRACT FROM AUTHOR]

Uncertainty about sea-level rise is dominated by uncertainty about iceberg calving, mass loss from glaciers or ice sheets by fracturing. Review of the rapidly growing calving literature leads to a few overarching hypotheses. Almost all calving occurs near or just downglacier of a location where ice flows into an environment more favorable for calving, so the calving rate is controlled primarily by flow to the ice margin rather than by fracturing. Calving can be classified into five regimes, which tend to be persistent, predictable, and insensitive to small perturbations in flow velocity, ice characteristics, or environmental forcing; these regimes can be studied instrumentally. Sufficiently large perturbations may cause sometimes-rapid transitions between regimes or between calving and noncalving behavior, during which fracturing may control the rate of calving. Regime transitions underlie the largest uncertainties in sea-level rise projections, but with few, important exceptions, have not been observed instrumentally. This is especially true of the most important regime transitions for sea-level rise. Process-based models informed by studies of ongoing calving, and assimilation of deep-time paleoclimatic data, may help reduce uncertainties about regime transitions. Failure to include calving accurately in predictive models could lead to large underestimates of warming-induced sea-level rise. Iceberg calving, the breakage of ice from glaciers and ice sheets, affects sea level and many other environmental issues. Modern rates of iceberg calving usually are controlled by the rate of ice flow past restraining points, not by the brittle calving processes. Calving can be classified into five regimes, which are persistent, predictable, and insensitive to small perturbations. Transitions between calving regimes are especially important, and with warming might cause faster sea-level rise than generally projected. [ABSTRACT FROM AUTHOR]

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

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

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

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

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

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

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

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

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

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

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

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

Mass loss from Antarctica's three largest ice shelves is dominated by calving, primarily of large tabular icebergs every few decades. Smaller, more frequent calving events also occur, but it is more difficult to detect them and quantify their contribution to total ice‐shelf mass loss. We used surface elevation data from NASA's ICESat‐2 laser altimeter to examine the structure of the Ross Ice Shelf front between October 2018 and July 2020. Profiles frequently show a depression a few meters deep about 200–800 m upstream of the front, with higher values on the eastern portion of the ice shelf. This structure results from bending due to buoyancy of a submerged ice bench generated by ice‐front melting near the waterline when warm water is present in summer. These bending stresses may cause small‐scale calving events whose frequency would change as summer sea ice and atmosphere–ocean heat exchanges vary over time. Plain Language Summary: Mass loss from Antarctica's floating ice shelves, which form as the ice sheet extends into the Southern Ocean, influences how quickly grounded ice flows into the ocean. Estimating future sea‐level change from grounded‐ice loss therefore requires understanding, and developing models for, the processes that affect ice shelves. We used measurements of surface height from NASA's recently launched ICESat‐2 mission to explore one such process, the calving of small icebergs due to upper‐ocean melting of the ice front. We focus on the large Ross Ice Shelf. This local melting leads to bending of the ice shelf that can be seen in ICESat‐2 profiles that cross the ice front. The bending may also fracture the ice shelf to create small icebergs. We found that these surface structures are generally larger on the eastern portion of Ross Ice Shelf than on the western portion. We suggest that this pattern is due to differences in ice, ocean, and sea ice conditions that promote or impede the melting responsible for the ice‐shelf bending. ICESat‐2 will allow us to monitor changes in these small‐scale structures and any associated calving events, which will provide clues about how ice shelves will change in the future. Key Points: We used ICESat‐2 laser altimetry to map prevalent rampart‐moat (R‐M) structures along the Ross Ice Shelf frontR‐M structures indicate a submerged buoyant bench seaward of the aerial ice front and ice stresses that may cause small‐scale calvingAlong‐front variation of R‐M height differences provides insight into sensitivity to oceanic and glaciological conditions [ABSTRACT FROM AUTHOR]

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

Tabular calving events occur from Antarctica's large ice shelves only every few decades, and are preceded by rift propagation. We used high‐resolution imagery and ICESat‐2 data to determine the propagation rates for the three active rifts on Amery Ice Shelf (AIS; T1, T2, and E3) and observe the calving of D‐28 on September 25, 2019 along T1. AIS front advance accelerated downstream of T1 in the years before calving, possibly increasing stress at the rift tip. T1 experienced significant acceleration for 12 days before calving, coinciding with a jump in propagation of E3. ICESat‐2's high resolution and repeat acquisitions every 91 days allowed for analysis of the ice front before and after calving, and rift detection where it was not visible in imagery as a ∼1 m surface depression, suggesting that it propagates as a basal fracture. Our results show that ICESat‐2 can provide process‐scale information about iceberg calving. Plain Language Summary: Antarctica's ice shelves surround the ice sheet where it meets the ocean, and modulate the flow of grounded ice. Ice shelves lose the mass they gain through snowfall via two main processes: calving of icebergs at the front, and melting underneath. If losses exceed gains, they shrink, and the glaciers feeding them will feel less resistance to flow and speed up, and sea level will rise. Overall, Antarctica's ice shelves are experiencing excess mass loss; however, the ice sheet is losing mass only in key sectors, and the large ice shelves are mostly in equilibrium. For those ice shelves, large calving events are normal "background" mass loss processes. However, since they occur infrequently (only every few decades) we have not had many events to study with modern, high spatial resolution satellite data. Antarctica's third largest ice shelf, Amery Ice Shelf, calved a tabular iceberg in September 2019, providing an opportunity for a detailed process study of iceberg calving. We use satellite imagery and ICESat‐2 altimetry data to look at the precursors to calving and the effect of the calving event on the ice shelf. This will inform models about the calving process. Key Points: We use satellite imagery and ICESat‐2 repeat‐track data to examine a large tabular calving event from Amery Ice Shelf in September 2019ICESat‐2 repeat tracks reveal rift structures with uplifted flanking walls that appear to modulate future ice front shapeICESat‐2 reveals the subtle topographic signatures of active rifts, even beyond the limit where they are visible in satellite imagery [ABSTRACT FROM AUTHOR]

Projections of future ice sheet mass loss and thus sea level rise rely on the parametrization of iceberg calving in ice sheet models. The interconnection between submarine melt-induced undercutting and calving is still poorly understood, which makes predicted contributions of tidewater glaciers to sea level rise uncertain. Here, we compare detailed 3-D simulations of fracture initiation obtained with the Helsinki Discrete Element Model (HiDEM) to observations, prior to a major calving event at Bowdoin Glacier, Northwest Greenland. Observations of a plume surfacing at the calving location suggest that local melt-undercutting influenced the size of the major calving event. Therefore, several experiments are conducted with various local and distributed (front-wide) undercut geometries. Although the number of undercut experiments is limited by computational requirements, one of the conjectured undercut geometries reproduces the crevasse leading to the observed major calving event in great detail. Our simulations show that undercutting leads to initiation of wider fractures more than 100 m upstream of the terminus, well-beyond the directly undercut region. When combining a moderate distributed undercut with local amplified undercuts at the two observed plumes, fracture initiation also increases in between the local undercuts. Thus, our results agree with previous studies suggesting the existence of a “calving amplifier” effect by submarine melt, both upglacier and across-glacier. Consequently, the simulations show the potentially large impact of submarine melt-induced undercutting on iceberg size.

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

The rate of land ice loss due to iceberg calving is a key source of variability among model projections of the 21st century sea level rise. It is especially challenging to account for mass loss due to iceberg calving in Greenland, where ice drains to the ocean through hundreds of outlet glaciers, many smaller than typical model grid scale. Here, we apply a numerically efficient network flowline model (SERMeQ) forced by surface mass balance to simulate an upper bound on decadal calving retreat of 155 grounded outlet glaciers of the Greenland Ice Sheet—resolving five times as many outlets as was previously possible. We show that the upper bound holds for 91% of glaciers examined and that simulated changes in terminus position correlate with observed changes. SERMeQ can provide a physically consistent constraint on forward projections of the dynamic mass loss from the Greenland Ice Sheet associated with different climate projections. Key Points: We test an upper‐bound model of calving retreat of 155 ocean‐terminating outlet glaciers that drain the Greenland Ice SheetOur physics‐based method produces terminus positions that correlate with observed positions for 103 glaciers without model tuningOur model bounds retreat rates on 91% of glaciers tested, providing a constraint for future sea level projections [ABSTRACT FROM AUTHOR]

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

The Amery Ice Shelf (AIS) dynamics and mass balance caused by iceberg calving and basal melting are significant in the ocean climate system. Using satellite imagery from Sentinel-1 SAR, we monitored the temporal and spatial variability of the frontal positions on the Amery Ice Shelf, Antarctica, from 2015 to 2021. In this paper, we propose an automatic algorithm based on the SO-CFAR strategy and a profile cumulative method for frontal line extraction. To improve the accuracy of the extracted frontal lines, we developed a framework combining the Constant False Alarm Rate (CFAR) and morphological image-processing strategies. A visual comparison between the proposed algorithm and state-of-the-art algorithm shows that our algorithm is effective in these cases including rifts, icebergs, and crevasses as well as ice-shelf surface structures. We present a detailed analysis of the temporal and spatial variability of fronts on AIS that we find, an advance of the AIS frontal line before the D28 calving event, and a continuous advance after the event. The study reveals that the AIS extent has been advanced at the rate of 1015 m/year. Studies have shown that the frontal location of AIS has continuously expanded. From March 2015 to May 2021, the frontal location of AIS expanded by 6.5 km; while the length of the AIS frontal line is relatively different after the D28 event, the length of the frontal line increased by about 7.5% during 2015 and 2021 (255.03 km increased to 273.5 km). We found a substantial increase in summer advance rates and a decrease in winter advance rates with the seasonal characteristics. We found this variability of the AIS frontal line to be in good agreement with the ice flow velocity.

Glacial earthquakes are caused by large iceberg calving events, which are an important mechanism for mass loss from the Greenland ice sheet. The number of glacial earthquakes in Greenland has increased sixfold over the past two decades. We use teleseismic surface waves to analyze the 145 glacial earthquakes that occurred in Greenland from 2011 through 2013, and successfully determine source parameters for 139 events at 13 marine-terminating glaciers. Our analysis increases the number of events in the glacial-earthquake catalog by nearly 50% and extends it to 21 years. The period 2011–13 was the most prolific 3-year period of glacial earthquakes on record, with most of the increase over earlier years occurring at glaciers on Greenland's west coast. We investigate changes in earthquake productivity and geometry at several individual glaciers and link patterns in glacial-earthquake production and cessation to the absence or presence of a floating ice tongue. We attribute changes in earthquake force orientations to changes in calving-front geometry, some of which occur on timescales of days to months. Our results illustrate the utility of glacial earthquakes as a remote-sensing tool to identify the type of calving event, the grounded state of a glacier, and the orientation of an active calving front.

The simple calving laws currently used in ice-sheet models do not adequately reflect the complexity and diversity of calving processes. To be effective, calving laws must be grounded in a sound understanding of how calving actually works. Here, we develop a new strategy for formulating calving laws, using (a) the Helsinki Discrete Element Model (HiDEM) to explicitly model fracture and calving processes, and (b) the continuum model Elmer/Ice to identify critical stress states associated with HiDEM calving events. A range of observed calving processes emerges spontaneously from HiDEM in response to variations in ice-front buoyancy and the size of subaqueous undercuts. Calving driven by buoyancy and melt under-cutting is under-predicted by existing calving laws, but we show that the location and magnitude of HiDEM calving events can be predicted in Elmer/Ice from characteristic stress patterns. Our results open the way to developing calving laws that properly reflect the diversity of calving processes, and provide a framework for a unified theory of the calving process continuum.

On 2 March 2016, several small en échelon tabular icebergs calved from the seaward front of the McMurdo Ice Shelf, and a previously inactive rift widened and propagated by ~3 km, ~25% of its previous length, setting the stage for the future calving of a ~14 km2 iceberg. Within 24 h of these events, all remaining land-fast sea ice that had been stabilizing the ice shelf broke-up. The events were witnessed by time-lapse cameras at nearby Scott Base, and put into context using nearby seismic and automatic weather station data, satellite imagery and subsequent ground observation. Although the exact trigger of calving and rifting cannot be identified definitively, seismic records reveal superimposed sets of both long-period (>10 s) sea swell propagating into McMurdo Sound from storm sources beyond Antarctica, and high-energy, locally-sourced, short-period (

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

Many large calving events at Greenland's marine‐terminating glaciers generate globally detectable glacial earthquakes. We perform a cross‐correlation analysis using regional seismic data to identify events below the teleseismic detection threshold, focusing on the 24 hr surrounding known glacial earthquakes at Greenland's three largest glaciers. We detect additional seismic events in the minutes prior to more than half of the glacial earthquakes we study and following one third of them. Waveform modeling shows source mechanisms like those of previously known glacial earthquakes, a result consistent with available imagery. The seismic events thus do not represent a failure of the high subaerial ice cliff like that expected to trigger large‐scale calving and a marine ice‐cliff instability but, rather, rotational, buoyancy‐driven calving events, likely of the full glacier thickness. A limited investigation of the prevalence of smaller seismic events at times outside glacial‐earthquake windows identifies several additional events. However, we find that calving at the three glaciers we study—Jakobshavn Isbræ, Helheim Glacier, and Kangerdlugssuaq Glacier—often occurs as sequences of discrete buoyancy‐driven events in which multiple icebergs ranging in size over as much as three orders of magnitude are all lost within ∼30 min. We demonstrate a correlation between glacial‐earthquake magnitude and iceberg size for events with well‐constrained iceberg‐area estimates. Our results suggest that at least 10–30% more dynamic mass loss occurs through buoyancy‐driven calving at Greenland's glaciers than previously appreciated. Key Points: Seismic data show ice motion prior to half of large‐scale calving events we study and following one third of eventsSeismogenic ice motion prior to large calving events represents not subaerial ice‐cliff failure but rotational, buoyancy‐driven calvingWe present the first empirical demonstration of a relationship between the seismic magnitude of a glacial earthquake and iceberg size [ABSTRACT FROM AUTHOR]