Search

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]

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]

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]

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]

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]

The seafloor depths under the Cook Ice Shelf and Ninnis Glacier Tongue have not been directly measured, despite their importance for understanding ocean circulation and ice shelf change. We model the bathymetry underneath the floating ice and surrounding ocean using airborne gravity data. Our model is constrained by few ship‐based seafloor measurements near the ice front and by ice‐base measurements over areas of grounded ice from radar data. Localized basins (∼1,400 m deep) are found beneath both ice shelves. The shallowest modeled bathymetry (∼200 m) represents the offshore extension of Cape Freshfield. Near the grounding line, seafloor depths are found to be deeper than the observed depth of the modified Circumpolar Deep water in the region (<350 m), key factor for basal melt analyses. From transit flight gravity anomalies, we suggest the relocation of the mapped edge of the continental shelf and a narrowing of the Cook Shelf Depression. Plain Language Summary: The knowledge of how deep the ocean floor is under the floating ice shelves that connect to grounded ice sheets, is crucial for understanding how ocean water circulates and interacts with the overlying ice. We present a new bathymetric model of the seafloor beneath two ice shelves located in East Antarctica: Cook Ice Shelf and Ninnis Glacier Tongue. Both ice shelves are inaccessible to ships due to heavy sea ice conditions, so the data used in our model were collected from airborne surveys. Our bathymetry model shows new information on the depth and shape of the seafloor that will help understanding the ocean circulation in the area and how this might impact ice thickness changes. Key Points: High resolution bathymetry model of Cook Ice Shelf, Ninnis Glacier Tongue, and surrounding open ocean from airborne gravity inversionNew bathymetry model improves the understanding of water pathways between the ice shelves and the continental shelf edgeTransit flight gravity anomalies suggest relocation of the edge of the continental shelf northwards of the currently attributed position [ABSTRACT FROM AUTHOR]

Thwaites Ice Shelf (TWIS), the floating extension of Thwaites Glacier, West Antarctica, is changing rapidly and may completely disintegrate in the near future. Any buttressing that the ice shelf provides to the upstream grounded Thwaites glacier will then be lost. Previously, it has been argued that this could lead to onset of dynamical instability and the rapid demise of the entire glacier. Here we provide the first systematic quantitative assessment of how strongly the upstream ice is buttressed by TWIS and how its collapse affects future projections. By modeling the stresses acting along the current grounding line, we show that they deviate insignificantly from the stresses after ice shelf collapse. Using three ice‐flow models, we furthermore model the transient evolution of Thwaites Glacier and find that a complete disintegration of the ice shelf will not substantially impact future mass loss over the next 50 years. Plain Language Summary: As ice is lost, global sea levels rise. Currently glaciers and ice sheets contribute about half of observed sea level rise (SLR). Some ice is lost through surface melting, but sea levels can also be impacted by ice that flows directly into the ocean. This loss of ice, from the grounded interior into the surrounding ocean, is particularly important for the Antarctic Ice Sheet. There the grounded sections of the ice sheet are sometimes surrounded by floating seaward extensions referred to as ice shelves. These ice shelves can restrict, that is, buttress, upstream flow. Here we estimate the buttressing provided to one of the largest and most dynamical glaciers of the Antarctic Ice Sheet, the Thwaites Glacier. The glacier is currently abutted by an ice shelf, which has become progressively smaller and weaker over the last few decades. It has been suggested that a disintegration of the ice shelf might lead to greatly enhanced ice flow. However, we find that the ice shelf does not significantly buttress upstream flow, and removing the ice shelf has little effect on SLR over the next 50 years. Key Points: Thwaites Ice Shelf provides limited buttressing to the grounded upstream Thwaites GlacierRemoval of the ice shelf changes the contribution of the West Antarctic Ice Sheet to sea level rise by at most a few millimeters over 50 yearsObservations and modeling suggest that Thwaites Glacier is currently unsteady and stable [ABSTRACT FROM AUTHOR]

The fastest projected rates of sea level rise appear in models which include "the marine ice cliff instability (MICI)," a hypothesized but mostly unobserved process defined by rapid fracture and wastage of terminal ice cliffs that outpaces viscous relaxation and ice‐shelf formation. Crane Glacier's response to the Larsen B Ice Shelf collapse has been invoked as evidence of MICI in the observational record. Using available remote sensing data, we analyze Crane's retreat, arrest, and regrowth over the last two decades. Much of Crane's terminus retreat occurred in floating, not grounded ice. Retreat accelerated by at least 54% during the 2 years following ice shelf collapse. This is inconsistent with contemporaneous regional forcing that promoted ice shelf growth during this period, but consistent with a geometrically controlled positive feedback. We infer a maximum possible cliff height of 111 m, which according to process models, could enable cliff calving assuming damaged ice. Plain Language Summary: The behavior of Antarctic glaciers will largely determine the pace and magnitude of future sea level rise. But the projections made by ice sheet models are uncertain, in part due to the uncertain response of Antarctica to the future loss of its floating ice shelves. It has been hypothesized that ice shelf breakup could trigger a self‐sustaining mechanism of ice loss whereby ice cliffs collapse under their own weight. This idea is controversial because it has not been unambiguously observed in modern glacier systems. We show that after the loss of its ice shelf, Crane Glacier experienced a 2‐year period of accelerating ice loss. While not an unambiguous validation of the instability hypothesis, what we do observe is consistent with a positive feedback and unforced retreat. Models of ice cliff failure that assume glacier ice has pre‐existing weaknesses are consistent with the behavior we observe at Crane, suggesting that if ice was sufficiently damaged at the Crane terminus, retreat by cliff failure was theoretically possible. Key Points: Retreat of Crane Glacier's terminus after the collapse of the Larsen B Ice Shelf accelerated between November 2002 and 20046.74 km of the 10.30 km retreat occurred in floating ice, with ice cliff failure possible during retreat assuming damaged iceCrane's retreat into a narrow fjord and sea‐ice growth re‐established buttressing stresses, slowed calving, and reversed terminus retreat [ABSTRACT FROM AUTHOR]

We mapped the terminus position for every marine‐terminating glacier in the Northern Hemisphere for 2000, 2010, and 2020, including the Greenland Ice Sheet, to provide the first complete measure of their variability. In total, these 1,704 glaciers lost an average of 389.7 ± 1.6 km2 a−1 (total 7,527 ± 31 km2) from 2000 to 2020 with 123 glaciers becoming no longer marine‐terminating over this period. Overall, 85.3% of glaciers retreated, 2.5% advanced, and the remaining 12.3% did not change outside of uncertainty limits. Outlet glaciers of the Greenland Ice Sheet are responsible for 61.9% of total area loss, although their rate of retreat was 34% less in 2010–2020 than 2000–2010. Glaciers with the largest area loss terminate in ice shelves or ice tongues, are surge‐type, have an unstable basal geometry, or have an unusually wide calving margin. Plain Language Summary: North of the equator, 1,704 glaciers touched the ocean in 2000. Here, we present the first analysis to document the frontal position of every one of these glaciers in 2000, 2010, and 2020. We found that 85.3% retreated and are now reduced in area. Only 2.5% of glaciers advanced or increased in area. The remaining 12.3% did not change within uncertainty limits. Total area losses were 389.7 ± 1.6 km2 per year (total 7,527 ± 31 km2) over the 20‐year period. Glaciers flowing from the Greenland Ice Sheet accounted for over 60% of total area losses. We found wide variations in the response of glaciers to similar changes in air and ocean temperature and sea ice concentrations, showing that environmental conditions alone cannot explain why some glaciers retreated more than others. Instead, unique glacier characteristics are the most important factor in controlling the variability of terminus retreat. Glaciers with floating ice at their front (ice shelves or ice tongues), those that undergo periodic changes in their flow velocity (surges), those which have a weak connection to their beds, and glaciers that are unusually wide, experienced the largest area loss from 2000 to 2020. Key Points: There were 1704 marine‐terminating glaciers in 2000 in the Northern Hemisphere, of which 85.3% retreated and 2.5% advanced from 2000 to 2020Tidewater glaciers lost a total area of 7,527 ± 31 km2 from 2000 to 2020, with the Greenland Ice Sheet responsible for 61.9% of total lossesVariations in retreat are best explained by glacier characteristics: ice shelves/tongues, surging, basal geometry, and calving width [ABSTRACT FROM AUTHOR]

Tidal modulation of ice streams and their adjacent ice shelves is a real‐world experiment to understand ice‐dynamic processes. We observe the dynamics of Priestley Glacier, Antarctica, using Terrestrial Radar Interferometry (TRI) and GNSS. Ocean tides are predominantly diurnal but horizontal GNSS displacements also oscillate semi‐diurnally. The oscillations are strongest in the ice shelf and tidal signatures decay near‐linearly in the TRI data over >10 km upstream of the grounding line. Tidal flexing is observed >6 km upstream of the grounding line including cm‐scale uplift. Tidal grounding line migration is small and <40% of the ice thickness. The frequency doubling of horizontal displacements relative to the ocean tides is consistent with variable ice‐shelf buttressing demonstrated with a visco‐elastic Maxwell model. Taken together, this supports previously hypothesized flexural ice softening in the grounding‐zone through tides and offers new observational constraints for the role of ice rheology in ice‐shelf buttressing. Plain Language Summary: Temperatures in Antarctica's interior are below zero, so that ice continuously accumulates through snowfall. This mass gain is balanced by ice transport from the interior toward the coast. There, floating ice shelves form on the ocean and ice is eventually lost through iceberg calving and ocean induced melting. Small changes in any of these processes impact global ocean circulation and mean sea level. The speed of the ice varies with ocean tides. We use a specialized instrument that detects this variability and tidal flexing area‐wide and every few hours. We find that although low‐ and high‐tide only occur once a day in this area, ice flow is largest twice a day. This can be explained because ice softens when it is flexed by the ocean tides so that it can flow episodically faster. This is important because the deformation of ice shelves determines to an extent the stability of the entire ice sheet. Key Points: Terrestrial radar interferometry and GNSS identify dynamic changes at Priestley Glacier, Antarctica, forced by diurnal ocean tidesHorizontal displacements vary semi‐diurnally, are strongest in the shelf and decay near‐linearly >10 km upstream of the grounding lineGrounding line migration is small and flexural softening in ice‐shelf shear margins can explain the semi‐diurnal ice flow component [ABSTRACT FROM AUTHOR]