Search

The collapse of ice shelves could expose tall ice cliffs at ice sheet margins. The marine ice cliff instability (MICI) is a hypothesis that predicts that, if these cliffs are tall enough, ice may fail structurally leading to self-sustained retreat. To date, projections that include MICI have been performed with a single model based on a simple parameterization. Here, we implement a physically motivated parameterization in three ice sheet models and simulate the response of the Amundsen Sea Embayment after a hypothetical collapse of floating ice. All models show that Thwaites Glacier would not retreat further in the 21st century. In another set of simulations, we force the grounding line to retreat into Thwaites' deeper basin to expose a taller cliff. In these simulations, rapid thinning and velocity increase reduce the calving rate, stabilizing the cliff. These experiments show that Thwaites may be less vulnerable to MICI than previously thought, and model projections that include this process should be re-evaluated. [ABSTRACT FROM AUTHOR]

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

The stress state of ice shelves using numerical modeling is studied. An ice shelf is modeled by an elastic plate floating on water and attached to the ice cover at a grounding point. Stress calculations are carried out using the FLAC3D program. The resulting stresses significantly exceed the tensile strength of the glacier during bending. The maximum tensile stresses on the lower surface of the glacier near the grounding point reach 5 × 107 Pa and significantly exceed the maximum strength values of the glacier. The mechanism of rapid transformation of a solid ice shelf into mélange is shown. During the deformation process, vertical displacements (ice ridges) can reach tens of meters. [ABSTRACT FROM AUTHOR]

The intrusion of relatively warm water onto the continental shelf is widely recognized as a threat to Antarctic ice shelves and glaciers grounded below sea level, as enhanced ocean heat increases their basal melt. While the circulation of warm water has been documented on the East Antarctic continental shelf, the modes of warm water transport from the deep ocean onto the shelf are still uncertain. This makes predicting the future responses of major East Antarctic marine-grounded glaciers, such as Totten and Ninnis glaciers, particularly challenging. Here, we outline the key role of submarine canyons to convey southward flowing currents that transport warm Circumpolar Deep Water toward the East Antarctic shelf break, thus facilitating warm water intrusion on the continental shelf. Sediment drifts on the eastern flank of the canyons provide evidence for sustained southward-directed flows. These morpho-sedimentary features thus highlight areas potentially prone to enhanced ocean heat transport toward the continental shelf, with repercussions for past, present, and future glacial melting and consequent sea level rise. Sediment drifts record long-lasting southward flowing currents within submarine canyons, which convey warm Circumpolar Deep Water toward the East Antarctic shelf break, thus facilitating warm water intrusion on the continental shelf. [ABSTRACT FROM AUTHOR]

Understanding grounding line dynamics is critical for projecting glacier evolution and sea level rise. Observations from satellite radar interferometry reveal rapid grounding line migration forced by oceanic tides that are several kilometers larger than predicted by hydrostatic equilibrium, indicating the transition from grounded to floating ice is more complex than previously thought. Recent studies suggest seawater intrusion beneath grounded ice may play a role in driving rapid ice loss. Here, we investigate its impact on the evolution of Petermann Glacier, Greenland, using an ice sheet model. We compare model results with observed changes in grounding line position, velocity, and ice elevation between 2010 and 2022. We match the observed retreat, speed up, and thinning using 3‐km‐long seawater intrusion that drive peak ice melt rates of 50 m/yr; but we cannot obtain the same agreement without seawater intrusion. Including seawater intrusion in glacier modeling will increase the sensitivity to ocean warming. Plain Language Summary: Relatively warm seawater melts marine‐terminating glaciers from below. Recent observations suggest that seawater flows below grounded ice at high tide. The presence of seawater at this boundary, referred to as seawater intrusion, has the potential to increase glacier mass loss. We test this hypothesis on Petermann Glacier, Greenland, using an ice sheet flow model. We run the model to reconstruct the glacier's behavior from 2010 to 2022 with and without seawater intrusion. We compare the results with satellite observations of velocity, grounding line position, and ice thinning. When we use enhanced ice melt rates from kilometer‐scale seawater intrusion, we match the observed retreat, speed up, and thinning. When we do not, the model fails to replicate the observations. Seawater intrusion may play a critical role in glacier evolution. Adding this process to ice flow models will increase their sensitivity to ocean warming and projections of ice mass loss and sea level rise. Key Points: Ice melt caused by seawater intrusion in the grounding zone explains the observed grounding line retreat of Petermann GlacierWithout seawater intrusion in the grounding zone, we do not replicate the full extent of observed retreatIncluding seawater intrusion in the grounding zone increases glacier mass loss [ABSTRACT FROM AUTHOR]

Over the past decade, a wealth of research has been devoted to the detection of crevasses in glaciers and ice sheets via remote sensing and machine learning techniques. It is often argued that remotely sensed damage maps can function as early warning signals for shifts in ice shelf conditions from intact to damaged states and can serve as an important tool for ice sheet modellers to improve future sea level rise predictions. Here, we provide evidence for the Filchner–Ronne and Pine Island ice shelves that remotely sensed damage maps are only weakly related to the ice rate factor field A derived by an ice flow model when inverting for surface velocities. This technique is a common procedure in ice flow models, as it guarantees that any inferred changes in A relate to changes in ice flow measured through observations. The weak relationship found is improved when investigating heavily damaged shear margins, as observed on the Pine Island Ice Shelf; however, even in this setting, this association remains modest. Our findings suggest that many features identified as damage through remote sensing methods are not of direct relevance to present-day ice shelf flow. While damage can clearly play an important role in ice shelf processes and thus be relevant for ice sheet behaviour and sea level rise projections, our results imply that mapping ice damage directly from satellite observations may not directly help improve the representation of these processes in ice flow models. [ABSTRACT FROM AUTHOR]

The ice streams flowing into the Amundsen Sea, West Antarctica, are losing mass due to changes in oceanic basal melting of their floating ice shelves. Rapid ice‐shelf melting is sustained by the delivery of warm Circumpolar Deep Water to the ice‐shelf cavities, which is first supplied to the continental shelf by an undercurrent that flows eastward along the shelf break. Temporal variability of this undercurrent controls ice‐shelf basal melt variability. Recent work shows that on decadal timescales the undercurrent variability opposes surface wind variability. Using a regional model, we show that undercurrent variability is induced by sea‐ice freshwater fluxes, particularly those north of the shelf break, which affect the cross‐shelf break density gradient. This sea‐ice variability is linked to tropical Pacific variability impacting atmospheric conditions over the Amundsen Sea. Ice‐shelf melting also feeds back onto the undercurrent by affecting the on‐shelf density, thereby influencing shelf‐break density gradient anomalies. Plain Language Summary: The glaciers that flow toward the Amundsen Sea, West Antarctica, are losing ice faster than most others about the continent. Once these glaciers reach the coast, they extend out onto the ocean surface, forming ice shelves. The rapid loss of ice is caused by changes in melting by relatively warm ocean waters beneath the floating ice shelves. In the Amundsen Sea, a deep ocean current is responsible for delivering warm water from the deep ocean to the ice shelves. We present model results that show that this deep current varies on decadal timescales as a consequence of systematic sea‐ice melt and formation patterns. A faster current drives more rapid ice shelf melting which, via a feedback process, further accelerates the current. Climate variability originating in the tropical Pacific Ocean is responsible for the variability in the sea‐ice, and is therefore also responsible for the effects on melting of the ice shelves. Key Points: In the Amundsen Sea decadal variability of an undercurrent flowing along the shelf break drives decadal variability in ice‐shelf basal meltSea‐ice freshwater fluxes and positive feedbacks from ice‐shelf basal melt drive the undercurrent variabilityTropical Pacific teleconnections induce atmospheric anomalies over the Amundsen Sea which drive the sea‐ice freshwater flux variability [ABSTRACT FROM AUTHOR]

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

The Pine Island and Thwaites glaciers are the two largest contributors to sea level rise from Antarctica. Here we examine the influence of basal friction and ice shelf basal melt in determining projected losses. We examine both Weertman and Coulomb friction laws with explicit weakening as the ice thins to flotation, which many friction laws include implicitly via the effective pressure. We find relatively small differences with the choice of friction law (Weertman or Coulomb) but find losses to be highly sensitive to the rate at which the basal traction is reduced as the area upstream of the grounding line thins. Consistent with earlier work on Pine Island Glacier, we find sea level contributions from both glaciers to vary linearly with the melt volume averaged over time and space, with little influence from the spatial or temporal distribution of melt. Based on recent estimates of melt from other studies, our simulations suggest that the combined melt-driven and sea level rise contribution from both glaciers may not exceed 10 cm by 2200, although the uncertainty in model parameters allows for larger increases. We do not include other factors, such as ice shelf breakup, that might increase loss, or factors such as increased accumulation and isostatic uplift that may mitigate loss. [ABSTRACT FROM AUTHOR]

In recent decades, the Greenland ice sheet has been losing mass through glacier retreat and ice flow acceleration. This mass loss is linked with variations in submarine melt, yet existing ocean models are either coarse global simulations focused on decadal‐scale variability or fine‐scale simulations for process‐based investigations. Here, we unite these scales with a framework to downscale from a global state estimate (15 km) into a regional model (3 km) that resolves circulation on the continental shelf. We further downscale into a fjord‐scale model (500 m) that resolves circulation inside fjords and quantifies melt. We demonstrate this approach in Scoresby Sund, East Greenland, and find that interannual variations in submarine melt at Daugaard‐Jensen glacier induced by ocean temperature variability are consistent with the decadal changes in glacier ice dynamics. This study provides a framework by which coarse‐resolution models can be refined to quantify glacier submarine melt for future ice sheet projections. Plain Language Summary: Over the past several decades, the Greenland ice sheet has been losing ice and contributing to sea‐level rise. About half of this ice loss is induced by melt that occurs where glaciers meet the ocean. Using coarse‐scale ocean models that simulate circulation around the globe, previous studies have noted a strong link between ocean temperature and enhanced glacier ice loss. However, due to the small scale of Greenland's fjords, coarse models are unable to directly quantify circulation in these fjords and melt on submerged glaciers. In this study, we develop a new framework to "zoom in" on a fjord, using high‐resolution models driven by larger coarse‐resolution models. In this approach, we simulate melt on one of Greenland's biggest glaciers and find that periods of higher melt coincide with more ice loss as observed from satellites. Since this framework is adaptable to other regions, it could also be used to simulate melt on other glaciers and support estimates of future sea‐level rise. Key Points: Subsurface temperature variability is simulated in a narrow fjord network using regional models downscaled from a global state estimateModeled increases in ocean melt at Daugaard‐Jensen glacier coincide with the onset of acceleration in 2005 and retreat and thinning in 2011Model variations in shelf‐to‐fjord ocean properties match with observations, providing a basis to estimate ocean forcing in ice projections [ABSTRACT FROM AUTHOR]

The grounding line, delineating the boundary where a grounded glacier goes afloat in ocean water, shifts in response to tidal cycles. Here we analyze COSMO-SkyMed Differential Interferometric Synthetic Aperture Radar data acquired in 2020 and 2021 over Totten, Moscow University, and Rennick glaciers in East Antarctica, detecting tide-induced grounding line position variations from 0.5 to 12.5 km along prograde slopes ranging from ~0 to 5 %. Considering a glacier as a non-Newtonian fluid, we provide two-dimensional formulations of the viscous and viscoelastic short-term behavior of a glacier in partial frictional contact with the bedrock, and partially floating on sea water. Since the models' equations are not amenable to analytical treatment, numerical solutions are obtained using FEniCS, an open-source Python package. We establish the dependence of the grounding zone width on glacier thickness, bed slope, and glacier flow speed. The predictions of the viscoelastic model match ~93 % of all the DInSAR grounding zone measurements and are 71 % more accurate than those of the viscous model. The results of this study underscore the critical role played by ice elasticity in continuum mechanics-based glacier models, and being validated with the DInSAR measurements, can be used in other studies on glaciers. [ABSTRACT FROM AUTHOR]

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

Given the high number and diversity of events in a typical cryoseismic dataset, in particular those recorded on ice sheet margins, it is desirable to use a semi-automated method of grouping similar events for reconnaissance and ongoing analysis. We present a workflow for employing semi-unsupervised cluster analysis to inform investigations of the processes occurring in glaciers and ice sheets. In this demonstration study, we make use of a seismic event catalogue previously compiled for the Whillans Ice Stream, for the 2010–2011 austral summer outlined in Part 1,. We address the challenges of seismic event analysis for a complex wave field by clustering similar seismic events into groups using characteristic temporal, spectral, and polarization attributes of seismic time series with the k -means++ algorithm. This provides the basis for a reconnaissance analysis of a seismic wave field that contains local events (from the ice stream) set in an ambient wave field that itself contains a diversity of signals (mostly from the Ross Ice Shelf). As one result, we find that two clusters include stick-slip events that diverge in terms of length and initiation locality (i.e., central sticky spot and/or the grounding line). We also identify a swarm of high-frequency signals on 16–17 January 2011 that are potentially associated with a surface melt event from the Ross Ice Shelf. Used together with the event detection presented in Part 1, the semi-automated workflow could readily be generalized to other locations and, as a possible benchmark procedure, could enable the monitoring of remote glaciers over time and comparisons between locations. [ABSTRACT FROM AUTHOR]

In late March 2011, landfast sea ice (hereafter, "fast ice") formed in the northern Larsen B embayment and persisted continuously as multi-year fast ice until January 2022. In the 11 years of fast-ice presence, the northern Larsen B glaciers slowed significantly, thickened in their lower reaches, and developed extensive mélange areas, leading to the formation of ice tongues that extended up to 16 km from the 2011 ice fronts. In situ measurements of ice speed on adjacent ice shelf areas spanning 2011 to 2017 show that the fast ice provided significant resistive stress to ice flow. Fast-ice breakout began in late January 2022 and was closely followed by retreat and breakup of both the fast-ice mélange and the glacier ice tongues. We investigate the probable triggers for the loss of fast ice and document the initial upstream glacier responses. The fast-ice breakup is linked to the arrival of a strong ocean swell event (>1.5 m amplitude; wave period waves >5 s) originating from the northeast. Wave propagation to the ice front was facilitated by a 12-year low in sea ice concentration in the northwestern Weddell Sea, creating a near-ice-free corridor to the open ocean. Remote sensing data in the months following the fast-ice breakout reveals an initial ice flow speed increase (>2 -fold), elevation loss (9 to 11 m), and rapid calving of floating and grounded ice for the three main embayment glaciers Crane (11 km), Hektoria (25 km), and Green (18 km). [ABSTRACT FROM AUTHOR]

Glaciers along the Amundsen Sea coastline in West Antarctica are dynamically adjusting to a change in ice-shelf mass balance that triggered their retreat and speed-up prior to the satellite era. In recent decades, the ice shelves have continued to thin, albeit at a decelerating rate, whilst ice discharge across the grounding lines has been observed to have increased by up to 100 % since the early 1990s. Here, the ongoing evolution of ice-shelf mass balance components is assessed in a high-resolution coupled ice–ocean model that includes the Pine Island, Thwaites, Crosson, and Dotson ice shelves. For a range of idealized ocean-forcing scenarios, the combined evolution of ice-shelf geometry and basal-melt rates is simulated over a 200-year period. For all ice-shelf cavities, a reconfiguration of the 3D ocean circulation in response to changes in cavity geometry is found to cause significant and sustained changes in basal-melt rate, ranging from a 75 % decrease up to a 75 % increase near the grounding lines, irrespective of the far-field forcing. These previously unexplored feedbacks between changes in ice-shelf geometry, ocean circulation, and basal melting have a demonstrable impact on the net ice-shelf mass balance, including grounding-line discharge, at multi-decadal timescales. They should be considered in future projections of Antarctic mass loss alongside changes in ice-shelf melt due to anthropogenic trends in the ocean temperature and salinity. [ABSTRACT FROM AUTHOR]

Understanding drivers of marine-terminating ice sheet behavior is important for constraining ice contributions to global sea level rise. In part, the stability of marine-terminating ice is influenced by solid Earth conditions at the grounded-ice margin. While the Cordilleran Ice Sheet (CIS) contributed significantly to global mean sea level during its final post-Last-Glacial-Maximum (LGM) collapse, the drivers and patterns of retreat are not well constrained. Coastal outcrops in the deglaciated Puget Lowland of Washington State – largely below sea level during glacial maxima, then uplifted above sea level via glacial isostatic adjustment (GIA) – record the late Pleistocene history of the CIS. The preservation of LGM glacial and post-LGM deglacial sediments provides a unique opportunity to assess the variability in marine ice sheet behavior of the southernmost CIS. Based on paired stratigraphic and geochronological work, with a newly developed marine reservoir correction for this region, we identify that the late-stage CIS experienced stepwise retreat into a marine environment between 15 000 and 14 000 years before present, consistent with timing of marine incursion into the region reported in earlier works. Standstill of marine-terminating ice for at least 500 years, paired with rapid vertical landscape evolution, was followed by continued retreat of ice in a subaerial environment. These results suggest rapid rates of solid Earth uplift and topographic support (e.g., grounding zone wedges) stabilized the ice margin, supporting final subaerial ice retreat. This work leads to a better understanding of shallow-marine and coastal-ice-sheet retreat and is relevant to sectors of the contemporary Antarctic and Greenland ice sheets and marine-terminating outlet glaciers. [ABSTRACT FROM AUTHOR]

Ice surface speed increases dramatically from upstream to downstream in many ice streams and glaciers. This speed-up is thought to be associated with a transition from internal distributed deformation to highly localized deformation or sliding at the ice-bedrock interface. The physical processes governing this transition remain unclear. Here, we argue that highly localized deformation does not necessarily initiate at the ice-bedrock interface, but could also take the form of an internal shear band inside the ice flow that connects topographic highs. The power-law exponent n in the ice rheology amplifies the feedback between shear heating and shear localization, leading to the spontaneous formation of an internal shear band that can create flow separation within the ice. We model the thermomechanical ice flow over a sinusoidal basal topography by building on the high-resolution Stokes solver FastICE v1.0. We compile a regime diagram summarizing cases in which a sinusoidal topography with a given amplitude and wavelength leads to shear band formation for a given rheology. We compare our model results to borehole measurements from Greenland and find evidence to support the existence of an internal shear band. Our study highlights the importance of re-evaluating the degree to which internal deformation contributes to total deformation in the ice column and to the flow-to-sliding transition. Plain Language Summary On its way toward the ocean, ice speeds up dramatically from less than 1 m/year upstream to more than a kilometer per year downstream. In this paper, we investigate the physical processes controlling this speed-up. Specifically, we focus on the role of the bedrock topography and rheology in facilitating the transition from this slow to rapid motion. We use a two-dimensional numerical model to simulate the flow field within a slab of ice flowing down a ramp over a simplified topography. We find that including the bedrock topography can lead to a zone of highly localized deformation within the ice above topographic highs. We also find that a non-linear rheology amplifies this localization. We compare our model results to borehole measurements from Greenland and find evidence that supports the existence of this highly localized deformation zone. This study indicates that the localized deformation induced by bedrock topography and amplified by non-linear rheology could be one physical mechanism that governs the speed-up of the ice flow. [ABSTRACT FROM AUTHOR]

The heat transfer between the warm oceanic water and the floating portion of the Antarctic ice sheet (the ice shelves) occurs in a dynamic environment with year‐to‐year changes in the distribution of icebergs and fast‐ice (the “icescape”). Dramatic events such as the collapse of glacier tongues are apparent in satellite images but oceanographic observations are insufficient to capture the synoptic impact of such events on the supply of oceanic heat to ice shelves. This study uses a 3D numerical model and semi‐ idealized experiments to examine whether the current high melting rates of ice shelves in the Amundsen Sea could be mitigated by certain icescape configurations. Specifically, the experiments quantify the impacts on oceanic heat supply of presence/absence of the Thwaites Glacier Tongue, Bear Ridge Iceberg Chain, tabular iceberg B22, and fast‐ice cover seaward of Pine Island Ice Shelf (PIS). The experiments reveal that future changes in the coastal icescape are unlikely to reverse the high ice shelf melting rates of the Amundsen Sea, and that icescape changes between 2011 and 2022 actually enhanced them slightly. Ice shelves such as Crosson and Thwaites are found to have multiple viable sources of oceanic heat whose relative importance may shift following icescape reconfigurations but the overall heat supply remains high. Similarly, the formation of a fast‐ice cover seaward of PIS slows down the cavity circulation (by 7%) but does not reduce its heat supply. [ABSTRACT FROM AUTHOR]

Jakobshavn Isbræ, a major outlet glacier in Greenland, lost its protective ice shelf in 2002 and has been speeding up and retreating since. We image its grounding line for the first time with a terrestrial radar interferometer deployed in 2016 and detect its migration at tidal frequencies. The southern half of the glacier develops a floating section (3 km × 3 km) that migrates in phase with the tidal difference, up to a distance of 2.8 km, far more than previously expected. We attribute the migration to kilometer‐scale seawater intrusions, 10–20 cm in height, with the tide. The intrusions reveal that the glacier bed may be up to 800 m deeper than expected on the south side, which illustrates that our knowledge of bed topography remains limited for this glacier. We expect seawater intrusions to cause rapid melt of basal ice and play a major role in the glacier evolution. Plain Language Summary: The transition boundary between grounded glacier ice and floating glacier ice, or grounding line, has never been mapped in much detail on the largest, fastest outlet glaciers of Greenland because available satellite radar imagery does not provide short enough repeat pass data. Here, we use a terrestrial radar interferometer which images the glacier every 2 min to map the grounding line repeatedly with differential interferometry. Surprisingly the glacier develops a small floating section on the south side where the grounding line migrates over considerable distances—0.5 to 2.8 km—during the tidal cycle, which is 10 times farther than previously expected from flotation. We attribute the migration to seawater intrusions over a bed 100–800 m deeper than previously known. Seawater intrusions will carry sufficient ocean heat to melt basal ice vigorously, a factor that has not been incorporated in modeling studies of this glacier. Key Points: We present the first mapping of the grounding line of Jakobshavn Isbr$\ae$, Greenland at ocean tidal frequenciesThe grounding line migrates over kilometers, far more than expected from flotation, which we attribute to kilometer‐size seawater intrusionsSeawater intrusions along a bed several hundred meters deeper than expected must be included in future glacier modeling studies [ABSTRACT FROM AUTHOR]

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

Glacial fjords often contain ice mélange, a frozen conglomeration of icebergs, brash ice, and sea ice, that have been postulated to influence both glacier dynamics and fjord circulation through coupled mechanical and thermodynamic processes. Ice mélange meltwater can alter stratification of the water column by releasing cool, fresh water across a range of depths in the upper layer of the fjord. This meltwater input can subsequently modify the depth at which the subglacial discharge plume reaches neutral buoyancy and therefore the underlying buoyancy-driven fjord circulation and heat exchange with warm ocean shelf waters. Despite a spate of recent modelling studies exploring these proposed feedbacks, we lack in situ observations quantifying changes to the water column induced by ice mélange meltwater. Here we use a novel dataset collected before and after the melt, breakup, and down-fjord transport of an ephemeral ice mélange in front of Kangilliup Sermia (Rink Isbræ) to directly investigate the extent to which ice mélange meltwater can modify glacier-adjacent water properties. We find that even a short-lived ice mélange (4 days) can cause substantial cooling (0.18 °C) and freshening (0.25 g kg-1) of the water column that leads to stratification change down to the depth of the outflowing discharge plume. We compare our observations to an adjacent fjord, Kangerlussuup Sermia, where ice mélange seldom forms in the summertime, and show that the presence or absence of ice mélange melt creates fundamental differences in their upper layer hydrography. These observations provide critical constraints for recent modelling studies that have suggested ice mélange meltwater needs to be included in ocean circulation models for glaciers with deep grounding lines and high ice fluxes, which are precisely the glaciers exhibiting the largest magnitude terminus retreats at present. [ABSTRACT FROM AUTHOR]

We present a novel and computationally efficient subglacial hydrological model that represents in a unified way both hard and soft bed rheologies as well as a dynamic switch between efficient and inefficient subglacial discharge. The subglacial model is dynamically linked to a regularized Coulomb friction law, allowing for a coupled evolution of the ice sheet on decadal to centennial time scales. The hydrological model is tested on an idealized marine ice sheet and subsequently applied to the drainage basin of Thwaites Glacier, West Antarctica, that is composed of a heterogeneous (hard/soft) bed. We find that subglacial hydrology embedded in the sliding law accelerates the grounding line retreat of Thwaites Glacier under present-day climatic conditions. Highest retreat rates are obtained for hard bed configurations and/or inefficient drainage systems.We show that the sensitivity is particularly driven by large gradients in effective pressure, more so than the value of effective pressure itself in the vicinity of the grounding line. Clearly, a better understanding of the subglacial system is needed with respect to both the spatial and temporal variability in effective pressure and bed rheological conditions. [ABSTRACT FROM AUTHOR]

Today, relatively warm Circumpolar Deep Water is melting Thwaites Glacier at the base of its ice shelf and at the grounding zone, contributing to significant ice retreat. Accelerating ice loss has been observed since the 1970s; however, it is unclear when this phase of significant melting initiated. We analyzed the marine sedimentary record to reconstruct Thwaites Glacier's history from the early Holocene to present. Marine geophysical surveys were carried out along the floating ice-shelf margin to identify core locations from various geomorphic settings. We use sedimentological data and physical properties to define sedimentary facies at seven core sites. Glaciomarine sediment deposits reveal that the grounded ice in the Amundsen Sea Embayment had already retreated to within ~45 km of the modern grounding zone prior to ca. 9,400 y ago. Sediments deposited within the past 100+ y record abrupt changes in environmental conditions. On seafloor highs, these shifts document ice-shelf thinning initiating at least as early as the 1940s. Sediments recovered from deep basins reflect a transition from ice proximal to slightly more distal conditions, suggesting ongoing grounding-zone retreat since the 1950s. The timing of ice-shelf unpinning from the seafloor for Thwaites Glacier coincides with similar records from neighboring Pine Island Glacier. Our work provides robust new evidence that glacier retreat in the Amundsen Sea was initiated in the mid-twentieth century, likely associated with climate variability. [ABSTRACT FROM AUTHOR]

The regular and robust mapping of grounding lines is essential for various applications related to the mass balance of marine ice sheets and glaciers, especially in Antarctica and Greenland. With Differential Interferometric Synthetic Aperture Radar (DInSAR) interferograms, it is possible to accurately capture the tide-induced bending of the ice shelf at a continent-wide scale and a temporal resolution of a few days. While current processing chains typically automatically generate differential interferograms, grounding lines are still primarily identified and delineated on the interferograms by a human operator. This method is time-consuming and inefficient, considering the volume of data from current and future SAR missions. We developed a pipeline that utilizes the Holistically-Nested Edge Detection (HED) neural network to delineate DInSAR interferograms automatically. We trained HED in a supervised manner using 421 manually annotated grounding lines for outlet glaciers and ice shelves on the Antarctic Ice Sheet. We also assessed the contribution of non-interferometric features like elevation, ice velocity and differential tide levels towards the delineation task. Our best-performing network generated grounding lines with a median distance of 186 m from the manual delineations. Additionally, we applied the network to generate grounding lines for undelineated interferograms, demonstrating the network's generalization capabilities and potential to generate high-resolution temporal and spatial mappings. [ABSTRACT FROM AUTHOR]

Ice rises or ice rumples act as ice-shelf pinning points that can have an important role in regulating the ice discharge of marine outlet glaciers. As an example, the observed recent gradual ungrounding of the ice shelf of West Antarctica's Thwaites Glacier from its last pinning points likely diminished the buttressing effect of the ice shelf and thus contributed to the destabilization of the outlet. Here we use an idealized experimental setting to simulate the response of an Antarctic-type, instability-prone marine outlet glacier to a successive ungrounding of its ice shelf from a topographic high and a subsequent re-grounding. We show that the glacier retreat down the landward down-sloping (retrograde) bed, induced by the loss in pinning-point buttressing, can be unstable and irreversible given a relatively deep subglacial bed depression. In this case, glacier retreat and re-advance show a hysteretic behavior and if the bed depression is sufficiently deep, the glacier does not recover from but remains locked in its collapsed state. Conversely, reversibility requires a sufficiently shallow bed depression. Based on a simple flux balance analysis, we argue that the combination of a deep bed depression and limited ice-shelf buttressing hampers grounding-line re-advance due to the dominant and highly non-linear influence of the bed depth on the ice discharge across the grounding line. We conclude that outlets that rest on a deep bed depression and are weakly buttressed, such as Thwaites Glacier, are more susceptible to abrupt and irreversible retreat than stronger buttressed glaciers on more moderate retrograde slope, such as Pine Island Glacier. In particular, our results suggest that the wide and deep marine bed depression in the interior of Thwaites Glacier's drainage basin might promote potential future unstable retreat and also represent a strong limitation for a possible re-advance of the glacier in case it would collapse. [ABSTRACT FROM AUTHOR]

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

We observe the evacuation of 11-year-old landfast sea ice in the Larsen B embayment on the East Antarctic Peninsula in January 2022, which was in part triggered by warm atmospheric conditions and strong offshore winds. This evacuation of sea ice was closely followed by major changes in the calving behaviour and dynamics of a subset of the ocean-terminating glaciers in the region. We show using satellite measurements that, following a decade of gradual slow-down, Hektoria, Green, and Crane glaciers sped up by approximately 20 %–50 % between February and the end of 2022, each increasing in speed by more than 100 ma-1. Circumstantially, this is attributable to their transition into tidewater glaciers following the loss of their ice shelves after the landfast sea ice evacuation. However, a question remains as to whether the landfast sea ice could have influenced the dynamics of these glaciers, or the stability of their ice shelves, through a buttressing effect akin to that of confined ice shelves on grounded ice streams. We show, with a series of diagnostic modelling experiments, that direct landfast sea ice buttressing had a negligible impact on the dynamics of the grounded ice streams. Furthermore, we suggest that the loss of landfast sea ice buttressing could have impacted the dynamics of the rheologically weak ice shelves, in turn diminishing their stability over time; however, the accompanying shifts in the distributions of resistive stress within the ice shelves would have been minor. This indicates that this loss of buttressing by landfast sea ice is likely to have been a secondary process in the ice shelf disaggregation compared to, for example, increased ocean swell or the drivers of the initial landfast sea ice disintegration. [ABSTRACT FROM AUTHOR]

The 79° North Glacier (Nioghalvfjerdsbrae, 79NG) is one of three remaining glaciers with a floating tongue in Greenland. Although the glacier has been considered exceptionally stable in the past, earlier studies have shown that the ice tongue has thinned in recent decades. By conducting high-resolution ground-based and airborne radar measurements in conjunction with satellite remote-sensing observations, we find significant changes in the geometry of 79NG. In the vicinity of the grounding line, a 500 m high subglacial channel has grown since ∼ 2010 and has caused surface lowering of up to 7.6 m a -1. Our results show extreme basal melt rates exceeding 150 m a -1 over a period of 17 d within a distance of 5 km from the grounding line, where the ice has thinned by 32 % since 1998. We find a heterogeneous distribution of melt rates, likely due to variability in water column thickness and channelization of the ice base. Time series of melt rates show a decrease in basal melting since 2018, indicating an inflow of colder water into the cavity below 79NG. We discuss the processes that have led to the changes in geometry and conclude that the inflow of warm ocean currents has led to the extensive thinning of 79NG's floating ice tongue near the grounding line over the last 2 decades. In contrast, we hypothesize that the growth of the channel results from increased subglacial discharge due to a considerably enlarged area of summer surface melt due to the warming of the atmosphere. [ABSTRACT FROM AUTHOR]

The new photon-counting laser altimetry satellite ICESat-2 was successfully launched on September 15, 2018 with an unprecedented ice surface elevation measurement accuracy of 2–4 cm. It is meaningful for accurately estimating volumetric changes in the Antarctic Ice Sheet. Cross-validation of different types of data, especially comparison with in-situ data, is important for the ice sheet mass balance results of a new satellite. This paper proposes an elevation-difference method of grid elevation change rate model, based on ICESat-2 ATL11 elevation time series data, to distinguish the linear change trend of ice sheet surface elevation from the elevation change resulting from periodic precipitation. In order to compare the estimated elevation-change rate with the flux computed from in-situ snow stake velocity measurements and GRACE-FO gravity survey data, we made corrections for firn air content, elastic, and glacier isostatic adjustment. Based on the ATL11 data from 2019 to 2023, our results show that the ice sheet change in Basin 11 along the CHINARE traverse is from ~0.019 m yr−1 to ~0.121 m yr−1, and the mass balance in the upstream of the traverse in Basin 11 is ~1.9±0.2 Gt yr−1. It is comparable to ~2.4±1.2 Gt yr−1 from GRACE-FO during the same time period. Furthermore, the flux across traverse-11.9±1.1 Gt yr−1 is comparable to that of ~-9.7±0.9 Gt yr−1 across the same flux gate during 1997–2009 which is calculated based on GNSS-derived ice velocity observations, considering the time period difference and uncertainties. [ABSTRACT FROM AUTHOR]

Adélie‐George V Land in East Antarctica, encompassing the vast Wilkes Subglacial Basin, has a configuration that could be prone to ice sheet instability: the basin's retrograde bed slope could make its marine terminating glaciers vulnerable to warm seawater intrusion and irreversible retreat under predicted climate forcing. However, future projections are uncertain, due in part to limited subglacial observations near the grounding zone. Here, we develop a novel statistical approach to characterize subglacial conditions from radar sounding observations. Our method reveals intermixed frozen and thawed bed within 100 km of the grounding‐zone near the Wilkes Subglacial Basin outflow, and enables comparisons to ice sheet model‐inferred thermal states. The signs of intermixed or near thawed conditions raises the possibility that changes in basal thermal state could impact the stability of Adélie‐George V Land, adding to the region's potentially vulnerable topographic configuration and sensitivity to ocean forcing driven grounding line retreat. Plain Language Summary: East Antarctica's Adélie‐George V Land has been relatively stable over the last few decades. However, this region contains the Wilkes Subglacial Basin, which has a downward‐sloping bed inland of the grounding zone. This could make irreversible retreat possible if warming seawater off the coast enters beneath the ice sheet. However, predicting the region's vulnerability is difficult, in part, because there is limited information about the conditions beneath the ice sheet. In this study, we develop a new statistical approach to synthesize radar sounding data and classify the conditions at the ice‐bed interface into frozen‐bed and thawed‐bed, which can then provide comparisons to ice sheet model output. We find that areas near the outflow of the Wilkes Subglacial Basin, critical in maintaining the stability of the region, might consist of mixed frozen‐bed and thawed‐bed or near‐thawed conditions on the scale of tens of kilometers across. This finding is important since the extent of basal thaw affects how easily ice can flow or slide over the bed. If parts of the bed are close to thawed, this could make Adélie‐George V Land more sensitive to climate forcing, possibly resulting in mass loss. Key Points: We develop an adaptable statistical framework using radar sounding data to classify the basal thermal state of ice sheetsApplied to the Adélie‐George V Coast, the framework reveals a mix of frozen and thawed‐bed, along with confidence in the classificationsAreas maintaining the region's stability have varied thermal states and we consider how this could increase sensitivity to climate forcing [ABSTRACT FROM AUTHOR]

Pine Island and Thwaites glaciers are the two largest contributors to sea level rise from Antarctica. Here we examine the influence of basal friction and melt in determining projected losses. We examine both Weertman and Coulomb friction laws with explicit weakening as the ice thins to flotation, which many friction laws include implicitly via the effective pressure. We find relatively small differences with the choice of friction law (Weertman or Coulomb) but find losses are highly sensitive to the rate at which the basal traction is reduced as the area above the grounding line thins. Consistent with earlier work on Pine Island Glacier, we find sea level contributions from both glaciers vary linearly with the melt volume averaged over time and space, with little influence from the spatial or temporal distribution of melt. Based on recent estimates of melt from other studies, our work simulations suggest that melt-driven combined sea-level rise contribution from both glaciers is unlikely to exceed 10 cm by 2200. We do not include other factors, such as ice shelf breakup that might increase loss, nor factors such as increased accumulation and isostatic uplift that may mitigate loss. [ABSTRACT FROM AUTHOR]

Geothermal heat flow (GHF) is the dominant factor affecting the basal thermal regime of ice sheet dynamics. But it is poorly defined for the Antarctic ice sheet. We compare the basal thermal state of the Totten Glacier catchment as simulated by eight different GHF datasets. We use a basal energy and water flow model coupled with a 3D full-Stokes ice dynamics model to estimate the basal temperature, basal friction heat and basal melting rate. In addition to the location of subglacial lakes, we use specularity content of the airborne radar returns as a two-sided constraint to discriminate between local wet or dry basal conditions and compare the returns with the basal state simulations with different GHFs. Two medium magnitude GHF distribution maps derived from seismic modelling rank well at simulating both cold- and warm-bed regions, the GHFs from Shen et al. (2020) and Shapiro and Ritzwoller (2004). The best-fit simulated result shows that most of the inland bed area is frozen. Only the central inland subglacial canyon, co-located with high specularity content, reaches the pressure melting point consistently in all the eight GHFs. Modelled basal melting rates in the slow-flowing region are generally 0–5 mm yr -1 but with local maxima of 10 mm yr -1 at the central inland subglacial canyon. The fast-flowing grounded glaciers close to the Totten ice shelf are lubricating their bases with meltwater at rates of 10–400 mm yr -1. [ABSTRACT FROM AUTHOR]

Satellite radar interferometry data reveals that the grounding line of Petermann Glacier, Greenland migrates by several kilometers during the tidal cycle, bringing pressurized, subsurface, warm ocean waters in regular contact with a large sector of grounded ice. We use the Massachusetts Institute of Technology general circulation model in two dimensions to calculate the ice melt rates as a function of grounding zone (GZ) length and ocean Thermal Forcing (TF). Ice melt rates are found to be higher in the GZ cavity than anywhere else in the ice shelf cavity. The melt rates increase sub‐linearly with the length of the GZ and ocean TF. The model results agree well with remote sensing estimates of ice melt. High basal ice melt rates in tidally flushed grounding zones imply that marine‐terminating glaciers are more sensitive to ocean TF than anticipated, which will increase their projected contribution to sea level rise. Plain Language Summary: The traditional view of ice melting in contact with ocean waters is that melt rates drop to zero at the grounding line, which is a semi‐fixed boundary at the junction between grounded ice and the ocean. In reality, the grounding line migrates by kilometers during the tidal cycle, more than 10 times beyond the range expected from hydrostatic equilibrium, which brings warm, saline water in rapid contact with broad sectors of grounded ice. We use an ocean model to calculate the melt rates caused by seawater intrusions. We find that the melt rates in the GZ are higher than anywhere else in the ice shelf cavity and increase as the GZ becomes wider and the ocean gets warmer. Ice melt in kilometer‐size grounding zones will reduce the basal resistance to flow and will increase the sensitivity of the glacier flow to ocean warming, hence projections of sea level rise from the glacier will go up. Key Points: First modeling of ice melt rates from seawater intrusions in the kilometer‐size grounding zone (GZ) of Petermann Glacier using an ocean modelModeled melt rates are highest in the GZ and increase linearly with GZ width and ocean thermal forcingHigh melt rates in kilometer‐size grounding zones imply a higher sensitivity of glaciers to ocean warming than anticipated [ABSTRACT FROM AUTHOR]

The Antarctic Ice Sheet represents the largest source of uncertainty in future sea level rise projections, with a contribution to sea level by 2100 ranging from -5 to 43 cm of sea level equivalent under high carbon emission scenarios estimated by the recent Ice Sheet Model Intercomparison for CMIP6 (ISMIP6). ISMIP6 highlighted the different behaviors of the East and West Antarctic ice sheets, as well as the possible role of increased surface mass balance in offsetting the dynamic ice loss in response to changing oceanic conditions in ice shelf cavities. However, the detailed contribution of individual glaciers, as well as the partitioning of uncertainty associated with this ensemble, have not yet been investigated. Here, we analyze the ISMIP6 results for high carbon emission scenarios, focusing on key glaciers around the Antarctic Ice Sheet, and we quantify their projected dynamic mass loss, defined here as mass loss through increased ice discharge into the ocean in response to changing oceanic conditions. We highlight glaciers contributing the most to sea level rise, as well as their vulnerability to changes in oceanic conditions. We then investigate the different sources of uncertainty and their relative role in projections, for the entire continent and for key individual glaciers. We show that, in addition to Thwaites and Pine Island glaciers in West Antarctica, Totten and Moscow University glaciers in East Antarctica present comparable future dynamic mass loss and high sensitivity to ice shelf basal melt. The overall uncertainty in additional dynamic mass loss in response to changing oceanic conditions, compared to a scenario with constant oceanic conditions, is dominated by the choice of ice sheet model, accounting for 52 % of the total uncertainty of the Antarctic dynamic mass loss in 2100. Its relative role for the most dynamic glaciers varies between 14 % for MacAyeal and Whillans ice streams and 56 % for Pine Island Glacier at the end of the century. The uncertainty associated with the choice of climate model increases over time and reaches 13 % of the uncertainty by 2100 for the Antarctic Ice Sheet but varies between 4 % for Thwaites Glacier and 53 % for Whillans Ice Stream. The uncertainty associated with the ice–climate interaction, which captures different treatments of oceanic forcings such as the choice of melt parameterization, its calibration, and simulated ice shelf geometries, accounts for 22 % of the uncertainty at the ice sheet scale but reaches 36 % and 39 % for Institute Ice Stream and Thwaites Glacier, respectively, by 2100. Overall, this study helps inform future research by highlighting the sectors of the ice sheet most vulnerable to oceanic warming over the 21st century and by quantifying the main sources of uncertainty. [ABSTRACT FROM AUTHOR]

One of the last remaining floating tongues of the Greenland ice sheet (GrIS), the Petermann Glacier ice shelf (PGIS), is seasonally shielded from warm Atlantic water (AW) by the formation of sea ice arches in the Nares Strait. However, continued decline of the Arctic sea ice extent and thickness suggests that arch formation is likely to become anomalous, necessitating an investigation into the response of PGIS to a year-round mobile and thin sea ice cover. We use a high-resolution unstructured grid 3-D ocean–sea ice–ice shelf setup, featuring an improved sub-ice-shelf bathymetry and a realistic PGIS geometry, to investigate in unprecedented detail the implications of transitions in the Nares Strait sea ice regime, that is, from a thick and landfast sea ice regime to a mobile, and further, a thin and mobile sea ice regime, with regard to PGIS basal melt. In all three sea ice regimes, basal melt near the grounding line (GL) presents a seasonal increase during summer, driven by a higher thermal driving. The stronger melt overturning increases the friction velocity slightly downstream, where enhanced friction-driven turbulent mixing further increases the thermal driving, substantially increasing the local melt. As the sea ice cover becomes mobile and thin, wind and (additionally in winter) convectively upwelled AW from the Nares Strait enter the PGIS cavity. While its effect on basal melting is largely limited to the shallower (<200 m) drafts during winter, in summer it extends to the GL (ca. 600 m) depth. In the absence of an increase in thermal driving, increased melting under the deeper (>200 m) drafts in winter is solely driven by the increased vertical shear of a more energetic boundary layer current. A similar behaviour is noted when transitioning from a mobile to a thin mobile sea ice cover in summer, when increases in thermal driving are negligible and increases in melt are congruent with increases in friction velocity. These results suggest that the projected continuation of the warming of the Arctic Ocean until the end of the 21st century and the accompanying decline in the Arctic sea ice extent and thickness will amplify the basal melt of PGIS, impacting the long-term stability of the Petermann Glacier and its contribution to the future GrIS mass loss and sea level rise. [ABSTRACT FROM AUTHOR]

Geomorphological records of past ice sheet change offer the opportunity to examine their centennial‐scale response to changing boundary conditions, which are not adequately captured in the satellite record. Here, we present the first reconstruction of ice surface lowering at Byrd Glacier, the largest outlet glacier of the Transantarctic Mountains. Using surface exposure ages from glacial erratic cobbles collected in two vertical transects along the Lonewolf Nunataks, we find the initial emergence of this set of nunataks occurred at ~15 ka, with a rapid pulse of thinning at ~8 ka. We compare our glacier thinning profiles with modelled ice sheet thickness and grounding line histories from two model ensembles to identify key processes responsible for ice sheet change. All model runs from the two ensembles predict grounding line retreat and inland thinning to occur in one rapid step from Last Glacial Maximum to present, in line with marine geology records, our exposure age data and derived glacier thinning rates. Experiments best matching the glacial thickness constraints, reconstructed from the surface exposure data, have faster basal sliding (i.e., promote greater sliding rates resulting in thinner ice). However, experiments best matching the timing and rapid rate of ice thinning derived from the same surface exposure data have higher basal friction. This apparent change in the modelled basal sliding regime, from when the ice surface is at maximum thickness, to the rapid thinning at ~8 ka, occurs as the grounding line retreats towards the Byrd Glacier and Ross Ice Shelf forms during the Holocene. This past context has implications for the stability of the modern grounding line of Byrd Glacier, which is characterised by high basal melt rates at the terminus—a process that has the potential to propagate glacier thinning far inland, impacting the overall (in)stability of the Byrd Glacier and Ross Ice Shelf. [ABSTRACT FROM AUTHOR]

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

Abstract. Grounding line discharge is a key component of the mass balance of the Antarctic Ice Sheet. Here we present an estimate of Antarctic Ice Sheet grounding line discharge from 1996 through to last month. We calculate ice flux at up to 100 m resolution through 16 algorithmically-generated flux gates, which are continuous around Antarctica. We draw on a range of ice velocity and thickness data to estimate grounding line discharge. For ice thickness, we use four bed topography datasets, two firn models and a temporally varying ice surface. For the ice velocity, we utilise a range of publicly-available ice velocity maps at resolutions ranging from 240 x 240 m to 1000 x 1000 m, as well as new, 100 x 100 m monthly velocity mosaics derived from intensity-tracking of Sentinel-1 image pairs, available since October 2014. The pixel-based ice fluxes and ice flux errors are integrated within all available Antarctic ice stream, ice shelf and glacier basins. Our dataset also includes the contributions to discharge from changes in ice thickness due to surface lowering, time-varying firn air content and surface mass change between the flux gates and grounding line. We find that Antarctic Ice Sheet grounding line discharge increased from 1,990 ± 23 Gt yr-1 to 2,205 ± 18 Gt yr-1 between 1996 and 2023, much of which was due to acceleration of ice streams in West Antarctica but with substantial contributions from ice streams in East Antarctica and glaciers on the Antarctic Peninsula. The uncertainties in our discharge dataset primarily result from uncertain bed elevation and flux gate location, which account for much of difference between our results and previous studies. It is our intention to update this discharge dataset each month, subject to continued Sentinel-1 acquisitions and funding availability. [ABSTRACT FROM AUTHOR]

Over the past decade, a wealth of research has been devoted to the detection of crevasses in glaciers and ice sheets via remote sensing and machine learning techniques. It is often argued that remotely sensed damage maps can function as early-warning signals for shifts in ice shelf conditions from intact to damaged states and can serve as an important tool for ice sheet modellers to improve future sea-level rise predictions. Here, we provide evidence for Filchner-Ronne and Pine Island ice shelves that remotely sensed damage maps are only weakly related to the ice rate factor field A derived by an ice-flow model when inverting for surface velocities. This technique is a common procedure in ice flow models, as it guarantees that any inferred changes in A relate to changes in ice flow measured through observations. The weak relationship found is improved when investigating heavily damaged shear margins, as observed on Pine Island Ice Shelf; yet, even in this setting, this association remains modest. Our findings suggest that many features identified as damage through remote sensing methods are not of direct relevance to present-day ice-shelf flow. While damage can clearly play an important role in ice-shelf processes and thus be relevant for ice-sheet behaviour and sea-level rise projections, our results imply that mapping ice damage directly from satellite observations may not directly help improve the representation of these processes in ice-flow models. [ABSTRACT FROM AUTHOR]

The glaciers of North Greenland are hosting enough ice to raise sea level by 2.1 m, and have long considered to be stable. This part of Greenland is buttressed by the last remaining ice shelves of the ice sheet. Here, we show that since 1978, ice shelves in North Greenland have lost more than 35% of their total volume, three of them collapsing completely. For the floating ice shelves that remain we observe a widespread increase in ice shelf mass losses, that are dominated by enhanced basal melting rates. Between 2000 and 2020, there was a widespread increase in basal melt rates that closely follows a rise in the ocean temperature. These glaciers are showing a direct dynamical response to ice shelf changes with retreating grounding lines and increased ice discharge. These results suggest that, under future projections of ocean thermal forcing, basal melting rates will continue to rise or remain at high level, which may have dramatic consequences for the stability of Greenlandic glaciers. North Greenland ice shelves have lost more than a third of their masses, thinning dramatically from below due to increased ocean temperatures. In response, grounding lines have retreated and the amount of ice discharged into the ocean has increased. [ABSTRACT FROM AUTHOR]

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

Soil and sediment mixing and associated permafrost processes are not widely studied or understood in the McMurdo Dry Valleys of Antarctica. In this study, we investigate the stability and depositional history of near-surface permafrost sediments to ∼ 3 m depth in the Pearse and lower Wright valleys using measured cosmogenic 10 Be and 26 Al depth profiles. In Pearse Valley, we estimate a minimum depositional age of ∼ 74 ka for the active layer and paleoactive-layer sediments (< 0.65 m). Combined depth profile modelling of 10 Be and 26 Al gives a depositional age for near-surface (< 1.65 m) permafrost in Pearse Valley of 180 +20/-40 ka, implying that the deposition of permafrost sediments predates MIS 5 advances of Taylor Glacier. Deeper permafrost sediments (> 2.09 m) in Pearse Valley are thus inferred to have a depositional age of > 180 ka. At a coastal, lower-elevation site in neighbouring lower Wright Valley, 10 Be and 26 Al depth profiles from a second permafrost core exhibit near-constant concentrations with depth and indicate the sediments are either vertically mixed after deposition or sufficiently young so that post-depositional nuclide production is negligible relative to inheritance. 26Al/10Be concentration ratios for both depth profiles range between 4.0 and 5.2 and are all lower than the nominal surface production rate ratio of 6.75, indicating that prior to deposition, these sediments experienced complex, yet similar, exposure–burial histories. Assuming a single-cycle exposure–burial scenario, the observed 26Al/10Be ratios are equivalent to a total minimum exposure–burial history of ∼ 1.2 Myr. In proximity to the depth profile core site, we measured cosmogenic 10 Be and 26 Al in three granite cobbles from thin, patchy drift (Taylor 2 Drift) in Pearse Valley to constrain the timing of retreat of Taylor Glacier. Assuming simple continuous exposure, our minimum, zero-erosion exposure ages suggest Taylor Glacier partially retreated from Pearse Valley no later than 65–74 ka. The timing of retreat after 65 ka and until the Last Glacial Maximum (LGM) when Taylor Glacier was at a minimum position remains unresolved. The surface cobble ages and permafrost processes reveal Taylor Glacier advances during MIS 5 were non-erosive or mildly erosive, preserving the underlying permafrost sediments and peppering boulders and cobbles upon an older, relict surface. Our results are consistent with U/Th ages from central Taylor Valley and suggest changes in moisture delivery over Taylor Dome during MIS 5e, 5c, and 5a appear to be associated with the extent of the Ross Ice Shelf and sea ice in the Ross Sea. These data provide further evidence of antiphase behaviour through retreat of a peripheral lobe of Taylor Glacier in Pearse Valley, a region that was glaciated during MIS 5. We suggest a causal relationship of cold-based glacier advance and retreat that is controlled by an increase in moisture availability during retreat of sea ice and perhaps the Ross Ice Shelf, as well as, conversely, a decrease during times of sea ice and Ross Ice Shelf expansion in the Ross Sea. [ABSTRACT FROM AUTHOR]

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

Some of the highest specific mass change rates in Antarctica are reported for the Antarctic Peninsula. However, the existing estimates for the northern Antarctic Peninsula (<70 ∘ S) are either spatially limited or are affected by considerable uncertainties. The complex topography, frequent cloud cover, limitations in ice thickness information, boundary effects, and uncertain glacial–isostatic adjustment estimates affect the ice sheet mass change estimates using altimetry, gravimetry, or the input-output method. Within this study, the first assessment of the geodetic mass balance throughout the ice sheet of the northern Antarctic Peninsula is carried out employing bi-static synthetic aperture radar (SAR) data from the TanDEM-X satellite mission. Repeat coverages from the austral winters of 2013 and 2017 are employed. Overall, coverage of 96.4 % of the study area by surface elevation change measurements and a total mass budget of -24.1±2.8 Gt a -1 are revealed. The spatial distribution of the surface elevation and mass changes points out that the former ice shelf tributary glaciers of the Prince Gustav Channel, Larsen A and B, and Wordie ice shelves are the hotspots of ice loss in the study area and highlights the long-lasting dynamic glacier adjustments after the ice shelf break-up events. The highest mass change rate is revealed for the Airy–Seller–Fleming glacier system at -4.9±0.6 Gt a -1 , and the highest average surface elevation change rate of -2.30±0.03 m a -1 is observed at Drygalski Glacier. The comparison of the ice mass budget with anomalies in the climatic mass balance indicates, that for wide parts of the southern section of the study area, the mass changes can be partly attributed to changes in the climatic mass balance. However, imbalanced high ice discharge drives the overall ice loss. The previously reported connection between mid-ocean warming along the southern section of the west coast and increased frontal glacier recession does not repeat in the pattern of the observed glacier mass losses, excluding in Wordie Bay. The obtained results provide information on ice surface elevation and mass changes for the entire northern Antarctic Peninsula on unprecedented spatially detailed scales and with high precision and will be beneficial for subsequent analysis and modeling. [ABSTRACT FROM AUTHOR]

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

Objectives: Reconstruction of ice flow velocity on the Antarctic Ice Sheet plays a critical role in estimating mass balance and contribution to global sea level rise. However, due to the low availability of optical satellite images before the 1990s, the time span of image pairs ranges from several years to more than ten years, leading to velocity overestimations. In this study, we correct the overestimated ice flow velocity based on the Lagrangian framework. After correction, we analyze the distribution of the overestimation values in the East Antarctic Ice Sheet (EAIS) according to the space and magnitude. Methods: We formulate a set of correction rules, including filtering points of ice flow velocity before correction and whether to retain results after correction. As a transition zone between grounded ice and ice shelf, ice flow velocity near the grounding line is one of the most important parameters of the ice sheet mass balance. In this paper, the overestimations of ice flow velocity are mainly studied from three aspects: First from different regions including inland region, grounding region and ice shelf region, then from different ocean sectors including Weddell Sea, West Indian Ocean, East Indian Ocean and Ross Sea, finally from different glaciers/ice shelves. Results: The average overestimation of ice flow velocity in EAIS is about 28 m/a. Both the grounding and ice shelf regions have larger overestimations than the inland region. The peak correction appears at ±25 km near the grounding line, and the maximum correction exceeds 300 m/a. As for ocean sectors, the overestimation of velocity is the largest in the East Indian Ocean. The average corrections of the Holmes Glacier, Ninnis Glacier, Publications Ice Shelf, and Totten Glacier in this sector exceed 50 m/a, much larger than that of EAIS. Conclusions: This paper shows that only a very small proportion of ice velocity points can meet the correction rules. The points show differences in regional distribution and they are concentrated mainly in the fast-flowing areas such as glaciers and ice shelves. Especially at the grounding zone and ice front areas of ice shelves, the ice flow velocity can be greatly overestimated. By analyzing the corrections and ice flow velocities in different ocean sectors, we found that the spatial acceleration trend of velocity in the East Indian Ocean is obvious and may be related to the interaction between the ocean and the ice shelves. The distribution of corrections varies among different ice shelves and glaciers. For example, Amery Ice Shelf has a larger correction at the confluence of three upstream glaciers, while Shirase Glacier has a larger correction in the inland region. The correction results in this paper have been applied to historical ice flow velocity products for East Antarctica from 1960s to 1980s. The overestimation correction of ice velocity is important for reducing the uncertainty in estimation and forecasting the contribution of Antarctic mass loss to the global sea level rise. [ABSTRACT FROM AUTHOR]

For better projections of sea level rise, two things are needed: an improved understanding of the contributing processes and their accurate representation in climate models. A major process is basal melting of ice shelves and glacier tongues by the ocean, which reduces ice sheet stability and increases ice discharge into the ocean. We study marine melting of Greenland's largest floating ice tongue, the 79° North Glacier, using a high‐resolution, 2D‐vertical ocean model. While our fjord model is idealized, the results agree with observations of melt rate and overturning strength. Our setup is the first application of adaptive vertical coordinates to an ice cavity. Their stratification‐zooming allows a vertical resolution finer than 1 m in the entrainment layer of the meltwater plume, which is important for the plume development. We find that the plume development is dominated by entrainment only initially. In the stratified upper part of the cavity, the subglacial plume shows continuous detrainment. It reaches neutral buoyancy near 100 m depth, detaches from the ice, and transports meltwater out of the fjord. Melting almost stops there. In a sensitivity study, we show that the detachment depth depends primarily on stratification. Our results contribute to the understanding of ice–ocean interactions in glacier cavities. Furthermore, we suggest that our modeling approach with stratification‐zooming coordinates will improve the representation of these interactions in global ocean models. Finally, our idealized model topography and forcing are close to a real fjord and completely defined analytically, making the setup an interesting reference case for future model developments. Plain Language Summary: The global increase of sea levels is a consequence of human‐induced climate change. It presents a threat to coastal regions and demands action to protect human life and infrastructure near the coast. Planning protective measures requires projections of sea level rise, computed with climate models. We present an approach to improve the simulation of an important contributor to sea level rise: melting of floating ice shelves by ocean circulation. Our modeling approach uses a vertical model grid that evolves over time. The temporal evolution depends on the density structure of the ocean. Large density differences appear just below an ice shelf, where fresh meltwater mixes with salty seawater. The adaptive grid of our model resolves this mixing process in great detail. This is important for an accurate computation of the melt rate and enables us to study in depth the ice shelf–ocean interactions. We study them at the glacier tongue of the 79° North Glacier, which is Greenland's largest ice shelf. The physical understanding gained from our simulations is also applicable to other floating glacier tongues and ice shelves. We suggest that using the presented model technique in global ocean models can improve projections of melting and sea level rise. Key Points: Melting of the 79° North Glacier ice tongue by turbulent ocean currents is studied with an idealized 2D‐vertical fjord modelThe subglacial plume behaves like an entraining plume close to the grounding line and like a detraining gravity current further downstreamA vertical resolution finer than 1 m is achieved in the subglacial plume by using adaptive vertical coordinates that zoom to stratification [ABSTRACT FROM AUTHOR]

Thwaites Glacier (TG) plays an important role in future sea-level rise (SLR) contribution from the West Antarctic Ice Sheet. Recent observations show that TG is losing mass, and its grounding zone is retreating. Previous modeling has produced a wide range of results concerning whether, when, and how rapidly further retreat will occur under continued warming. These differences arise at least in part from ill-constrained processes, including friction from the bed, and future atmosphere and ocean forcing affecting ice-shelf and grounding-zone buttressing. Here, we apply the Ice Sheet and Sea-level System Model (ISSM) with a range of specifications of basal sliding behavior in response to varying ocean forcing. We find that basin-wide bed character strongly affects TG's response to sub-shelf melt by modulating how changes in driving stress are balanced by the bed as the glacier responds to external forcing. Resulting differences in dynamic thinning patterns alter modeled grounding-line retreat across Thwaites' catchment, affecting both modeled rates and magnitudes of SLR contribution from this critical sector of the ice sheet. Bed character introduces large uncertainties in projections of TG under equal external forcing, pointing to this as a crucial constraint needed in predictive models of West Antarctica. [ABSTRACT FROM AUTHOR]