Ocean Thermal Forcing Research Articles

Seafloor roughness reduces melting of East Antarctic ice shelves

Heat delivered by the ocean circulation is melting the Antarctic ice sheet from below. This melt is largest where warm Circumpolar Deep Water accesses the continental shelf and reaches the ice shelf cavities. Future melt rate projections are based on ocean thermal forcing derived from climate models, which tend to be biased warm around Antarctica. The bias has been attributed to unresolved ocean processes that occur at scales poorly represented in models. Using a high-resolution model of the Denman Glacier region we show that seafloor roughness unresolved in climate models suppresses the impact of warm water on ice sheet melting. Seafloor roughness slows down the shelf circulation, reducing the presence of warm water over the shelf and the heat transport towards the ice cavities. As a result, the total meltwater discharge drops by 4 Gt year−1. Our results suggest a mechanism missing in global ocean and climate models that could reduce the spread in climate projections.

Heterogeneous impacts of ocean thermal forcing on ice discharge from Greenland's peripheral tidewater glaciers over 2000–2021

The Greenland Ice Sheet is losing mass at increasing rates. Substantial amounts of this mass loss occur by ice discharge which is influenced by ocean thermal forcing. The ice sheet is surrounded by thousands of peripheral, dynamically decoupled glaciers. The mass loss from these glaciers is disproportionately high considering their negligible share in Greenland’ overall ice mass. We study the relevance of ocean thermal forcing for ice discharge evolution in the context of this contrasting behaviour. Our estimate of ice discharge from the peripheral tidewater glaciers yields a rather stable Greenland-wide mean of 5.40 ± 3.54 Gt a−1 over 2000–2021. The evolutions of ice discharge and ocean thermal forcing are heterogeneous around Greenland. We observe a significant sector-wide increase of ice discharge in the East and a significant sector-wide decrease in the Northeast. Ocean thermal forcing shows significant increases along the northern/eastern coast, while otherwise unchanged conditions or decreases prevail. For East Greenland, this implies a clear influence of ocean thermal forcing on ice discharge. Similarly, we find clear influences at peripheral tidewater glaciers with thick termini that are similar to ice sheet outlet glaciers. At the peripheral glaciers in Northeast Greenland ice discharge evolution opposes ocean thermal forcing for unknown reasons.

Ocean cavity regime shift reversed West Antarctic grounding line retreat in the late Holocene

Recent geologic and modeled evidence suggests that the grounding line of the Siple Coast of the West Antarctic Ice Sheet (WAIS) retreated hundreds of kilometers beyond its present position in the middle to late Holocene and readvanced within the past 1.7 ka. This grounding line reversal has been attributed to both changing rates of isostatic rebound and regional climate change. Here, we test these two hypotheses using a proxy-informed ensemble of ice sheet model simulations with varying ocean thermal forcing, global glacioisostatic adjustment (GIA) model simulations, and coupled ice sheet-GIA simulations that consider the interactions between these processes. Our results indicate that a warm to cold ocean cavity regime shift is the most likely cause of this grounding line reversal, but that GIA influences the rate of ice sheet response to oceanic changes. This implies that the grounding line here is sensitive to future changes in sub-ice shelf ocean circulation.

Local forcing mechanisms challenge parameterizations of ocean thermal forcing for Greenland tidewater glaciers

Abstract. Frontal ablation has caused 32 %–66 % of Greenland Ice Sheet mass loss since 1972, and despite its importance in driving terminus change, ocean thermal forcing remains crudely incorporated into large-scale ice sheet models. In Greenland, local fjord-scale processes modify the magnitude of thermal forcing at the ice–ocean boundary but are too small scale to be resolved in current global climate models. For example, simulations used in the Ice Sheet Intercomparison Project for CMIP6 (ISMIP6) to predict future ice sheet change rely on the extrapolation of regional ocean water properties into fjords to drive terminus ablation. However, the accuracy of this approach has not previously been tested due to the scarcity of observations in Greenland fjords, as well as the inability of fjord-scale models to realistically incorporate icebergs. By employing the recently developed IceBerg package within the Massachusetts Institute of Technology general circulation model (MITgcm), we here evaluate the ability of ocean thermal forcing parameterizations to predict thermal forcing at tidewater glacier termini. This is accomplished through sensitivity experiments using a set of idealized Greenland fjords, each forced with equivalent ocean boundary conditions but with varying tidal amplitudes, subglacial discharge, iceberg coverage, and bathymetry. Our results indicate that the bathymetric obstruction of external water is the primary control on near-glacier thermal forcing, followed by iceberg submarine melting. Despite identical ocean boundary conditions, we find that the simulated fjord processes can modify grounding line thermal forcing by as much as 3 °C, the magnitude of which is largely controlled by the relative depth of bathymetric sills to the Polar Water–Atlantic Water thermocline. However, using a common adjustment for fjord bathymetry we can still predict grounding line thermal forcing within 0.2 °C in our simulations. Finally, we introduce new parameterizations that additionally account for iceberg-driven cooling that can accurately predict interior fjord thermal forcing profiles both in iceberg-laden simulations and in observations from Kangiata Sullua (Ilulissat Icefjord).

Modeling Ice Melt Rates From Seawater Intrusions in the Grounding Zone of Petermann Gletscher, Greenland

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

Rapid disintegration and weakening of ice shelves in North Greenland

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.

Local forcing mechanisms challenge parameterizations of ocean thermal forcing for Greenland tidewater glaciers

Local forcing mechanisms challenge parameterizations of ocean thermal forcing for Greenland tidewater glaciers

Exploring ice sheet model sensitivity to ocean thermal forcing and basal sliding using the Community Ice Sheet Model (CISM)

Abstract. Multi-meter sea level rise (SLR) is thought to be possible within the next few centuries, with most of the uncertainty originating from the Antarctic land ice contribution. One source of uncertainty relates to the ice sheet model initialization. Since ice sheets have a long response time (compared to other Earth system components such as the atmosphere), ice sheet model initialization methods can have significant impacts on how the ice sheet responds to future forcings. To assess this, we generated 25 different ice sheet spin-ups, using the Community Ice Sheet Model (CISM) at a 4 km resolution. During each spin-up, we varied two key parameters known to impact the sensitivity of the ice sheet to future forcing: one related to the sensitivity of the ice shelf melt rate to ocean thermal forcing (TF) and the other related to the basal friction. The spin-ups all nudge toward observed thickness and enforce a no-advance calving criterion, such that all final spin-up states resemble observations but differ in their melt and friction parameter settings. Each spin-up was then forced with future ocean thermal forcings from 13 different CMIP6 models under the Shared Socioeconomic Pathway (SSP)5-8.5 emissions scenario and modern climatological surface mass balance data. Our results show that the effects of the ice sheet and ocean parameter settings used during the spin-up are capable of impacting multi-century future SLR predictions by as much as 2 m. By the end of this century, the effects of these choices are more modest, but still significant, with differences of up to 0.2 m of SLR. We have identified a combined ocean and ice parameter space that leads to widespread mass loss within 500 years (low friction and high melt rate sensitivity). To explore temperature thresholds, we also ran a synthetically forced CISM ensemble that is focused on the Amundsen region only. Given certain ocean and ice parameter choices, Amundsen mass loss can be triggered with thermal forcing anomalies between 1.5 and 2 ∘C relative to the spin-up. Our results emphasize the critical importance of considering ice sheet and ocean parameter choices during spin-up for SLR predictions and suggest the importance of including glacial isostatic adjustment in ice sheet simulations.

Ocean‐Driven and Topography‐Controlled Nonlinear Glacier Retreat During the Holocene: Southwestern Ross Sea, Antarctica

AbstractRecent ice sheet mass loss in Antarctica has been attributed to an influx of warm ocean waters, which drove grounding‐line retreat and ice thinning. Episodic retreat and rapid thinning also occurred in the southwestern Ross Sea during the Holocene, which today accommodates cold ocean waters. We applied finite element ice‐flow modeling to investigate the roles of ocean temperature and bed topography in the deglaciation of this region. First, our experiments demonstrate that bed topography controlled the spatial pattern of grounding‐line retreat. Topographic pinning points limited the rate of ice loss until retreat progressed beyond a bathymetric threshold. Second, ocean thermal forcing determined the timing of this ice loss. Enhanced ocean‐driven melt is required during the Early‐to‐Mid Holocene to replicate geological records of deglaciation, possibly indicating that warm ocean waters were once present in this region. On multi‐centennial timescales, ocean temperature drove, while bed topography controlled, nonlinear rates of ice mass loss.

ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution using the Community Ice Sheet Model

Abstract. The future retreat rate for marine-based regions of the Antarctic Ice Sheet is one of the largest uncertainties in sea-level projections. The Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6) aims to improve projections and quantify uncertainties by running an ensemble of ice sheet models with atmosphere and ocean forcing derived from global climate models. Here, the Community Ice Sheet Model (CISM) is used to run ISMIP6-based projections of ocean-forced Antarctic Ice Sheet evolution. Using multiple combinations of sub-ice-shelf melt parameterizations and calibrations, CISM is spun up to steady state over many millennia. During the spin-up, basal friction parameters and basin-scale thermal forcing corrections are adjusted to optimize agreement with the observed ice thickness. The model is then run forward for 550 years, from 1950–2500, applying ocean thermal forcing anomalies from six climate models. In all simulations, the ocean forcing triggers long-term retreat of the West Antarctic Ice Sheet, especially in the Filchner–Ronne and Ross sectors. Mass loss accelerates late in the 21st century and then rises steadily for several centuries without leveling off. The resulting ocean-forced sea-level rise at year 2500 varies from about 150 to 1300 mm, depending on the melt scheme and ocean forcing. Further experiments show relatively high sensitivity to the basal friction law, moderate sensitivity to grid resolution and the prescribed collapse of small ice shelves, and low sensitivity to the stress-balance approximation. The Amundsen sector exhibits threshold behavior, with modest retreat under many parameter settings but complete collapse under some combinations of low basal friction and high thermal forcing anomalies. Large uncertainties remain, as a result of parameterized sub-shelf melt rates, simplified treatments of calving and basal friction, and the lack of ice–ocean coupling.

Ocean melting of the Zachariae Isstrøm and Nioghalvfjerdsfjorden glaciers, northeast Greenland

Zachariae Isstrøm (ZI) and Nioghalvfjerdsfjorden (79N) are marine-terminating glaciers in northeast Greenland that hold an ice volume equivalent to a 1.1-m global sea level rise. ZI lost its floating ice shelf, sped up, retreated at 650 m/y, and experienced a 5-gigaton/y mass loss. Glacier 79N has been more stable despite its exposure to the same climate forcing. We analyze the impact of ocean thermal forcing on the glaciers. A three-dimensional inversion of airborne gravity data reveals an 800-m-deep, broad channel that allows subsurface, warm, Atlantic Intermediate Water (AIW) (+1.[Formula: see text]C) to reach the front of ZI via two sills at 350-m depth. Subsurface ocean temperature in that channel has warmed by 1.3[Formula: see text]C since 1979. Using an ocean model, we calculate a rate of ice removal at the grounding line by the ocean that increased from 108 m/y to 185 m/y in 1979-2019. Observed ice thinning caused a retreat of its flotation line to increase from 105 m/y to 217 m/y, for a combined grounding line retreat of 13 km in 41 y that matches independent observations within 14%. In contrast, the limited access of AIW to 79N via a narrower passage yields lower grounded ice removal (53 m/y to 99 m/y) and thinning-induced retreat (27 m/y to 50 m/y) for a combined retreat of 4.4 km, also within 12% of observations. Ocean-induced removal of ice at the grounding line, modulated by bathymetric barriers, is therefore a main driver of ice sheet retreat, but it is not incorporated in most ice sheet models.

Constraining an Ocean Model Under Getz Ice Shelf, Antarctica, Using A Gravity‐Derived Bathymetry

AbstractGetz Ice Shelf, the largest producer of ice shelf meltwater in Antarctica, buttresses glaciers that hold enough ice to raise sea level by 22 cm. We present a new bathymetry of its sub‐ice shelf cavity using a three‐dimensional inversion of airborne gravity data constrained by multibeam bathymetry at sea and a reconstruction of the bedrock from mass conservation on land. The new bathymetry is deeper than previously estimated with differences exceeding 500 m in a number of regions. When incorporated into an ocean model, it yields a better description of the spatial distribution of ice shelf melt, specifically along glacier grounding lines. While the melt intensity is overestimated because of a positive bias in ocean thermal forcing, the study reveals the main pathways along which warm oceanic water enters the cavity and corroborates the observed rapid retreat of Berry Glacier along a deep channel with a retrograde bed slope.

Dynamics of Dalk Glacier in East Antarctica Derived from Multisource Satellite Observations Since 2000

Monitoring variability in outlet glaciers can improve the understanding of feedbacks associated with calving, ocean thermal forcing, and climate change. In this study, we present a remote-sensing investigation of Dalk Glacier in East Antarctica to analyze its dynamic changes. Terminus positions and surface ice velocities were estimated from Landsat and Sentinel-1 data, and the high-precision Worldview digital elevation model (DEM) was generated to determine the location of the potential ice rumple. We detected the cyclic behavior of glacier terminus changes and similar periodic increases in surface velocity since 2000. The terminus retreated in 2006, 2009, 2010, and 2016 and advanced in other years. The surface velocity of Dalk Glacier has a 5-year cycle with interannual speed-ups in 2007, 2012, and 2017. Our observations show the relationship between velocity changes and terminus variations, as well as the driving role of the ice rumple. Ice velocity often increases after calving events and continuous retreats. The loss of buttressing provided by an ice rumple may be a primary factor for increases in ice velocity. Given the restriction of the ice rumple, the surface velocity remains relatively stable when the glacier advances. The calving events may be linked to the unstable terminus caused by the ice rumple.

Neutral modes of surface temperature and the optimal ocean thermal forcing for global cooling

Inquiry into the climate response to external forcing perturbations has been the central interest of climate dynamics. But the understanding of two important aspects of climate change response—nonlinearity and regionality—is still lacking. Here a Green’s function approach is developed to estimate the linear response functions (LRFs) for both the linear and quadratic nonlinear response to ocean thermal forcing in a climate model, whereby the most excitable temperature modes, aka the neutral modes, can be identified for the current Earth climate. The resultant leading mode of the nonlinear response is characterized by a polar-amplified global cooling pattern, unveiling an intrinsic inclination of the modern climate towards cooling. Moreover, optimal forcing patterns are identified to most efficiently excite the corresponding neutral mode patterns. The forcing-response framework developed herein can be utilized to determine the optimal forcing patterns to inform solar geoengineering experiments and to interpret regional climate response and feedback in general.

Twenty-first century ocean forcing of the Greenland ice sheet for modelling of sea level contribution

Abstract. Changes in ocean temperature and salinity are expected to be an important determinant of the Greenland ice sheet's future sea level contribution. Yet, simulating the impact of these changes in continental-scale ice sheet models remains challenging due to the small scale of key physics, such as fjord circulation and plume dynamics, and poor understanding of critical processes, such as calving and submarine melting. Here we present the ocean forcing strategy for Greenland ice sheet models taking part in the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), the primary community effort to provide 21st century sea level projections for the Intergovernmental Panel on Climate Change Sixth Assessment Report. Beginning from global atmosphere–ocean general circulation models, we describe two complementary approaches to provide ocean boundary conditions for Greenland ice sheet models, termed the “retreat” and “submarine melt” implementations. The retreat implementation parameterises glacier retreat as a function of projected subglacial discharge and ocean thermal forcing, is designed to be implementable by all ice sheet models and results in retreat of around 1 and 15 km by 2100 in RCP2.6 and 8.5 scenarios, respectively. The submarine melt implementation provides estimated submarine melting only, leaving the ice sheet model to solve for the resulting calving and glacier retreat and suggests submarine melt rates will change little under RCP2.6 but will approximately triple by 2100 under RCP8.5. Both implementations have necessarily made use of simplifying assumptions and poorly constrained parameterisations and, as such, further research on submarine melting, calving and fjord–shelf exchange should remain a priority. Nevertheless, the presented framework will allow an ensemble of Greenland ice sheet models to be systematically and consistently forced by the ocean for the first time and should result in a significant improvement in projections of the Greenland ice sheet's contribution to future sea level change.

Impact of Iceberg Calving on the Retreat of Thwaites Glacier, West Antarctica Over the Next Century With Different Calving Laws and Ocean Thermal Forcing

AbstractThwaites Glacier, West Antarctica, has been a major contributor to global sea level rise over the past decades. Prior studies illustrated the critical role of ice shelf melt and iceberg calving based on cliff height in driving the retreat of Thwaites glacier. Here, we simulate its evolution with various calving laws and rates of frontal melt by the ocean in the absence of a buttressing ice shelf. Over the next century, we find that volume losses increase by 15–160% with a von Mises calving law compared to the case where the initial ice shelf is kept and the ice front is fixed at its current position, 10–20% with a buoyancy‐driven calving law, and 5–50% with frontal melt caused by ocean thermal forcing. Bed topography exerts the ultimate control on the evolution of Thwaites. In all simulations, once Thwaites Glacier retreats past the western subglacial ridge, the retreat becomes rapidly unstoppable.

Bathymetry of Southeast Greenland From Oceans Melting Greenland (OMG) Data

AbstractSoutheast Greenland has been a major participant in the ice sheet mass loss over the last several decades. Interpreting the evolution of glacier fronts requires information about their depth below sea level and ocean thermal forcing, which are incompletely known in the region. Here, we combine airborne gravity and multibeam echo sounding data from the National Aeronautics and Space Administration's Oceans Melting Greenland (OMG) mission with ocean probe and fishing boat depth data to reconstruct the bathymetry extending from the glacier margins to the edges of the continental shelf. We perform a three‐dimensional inversion of the gravity data over water and merge the solution with a mass conservation reconstruction of bed topography over land. In contrast with other parts of Greenland, we find few deep troughs connecting the glaciers to the sources of warm Atlantic Water, amidst a relatively uniform, shallow (350 m) continental shelf. The deep channels include the Kangerlugssuaq, Sermilik, Gyldenløve, and Tingmiarmiut Troughs.

Modeling the response of northwest Greenland to enhanced ocean thermal forcing and subglacial discharge

Abstract. Calving-front dynamics is an important control on Greenland's ice mass balance. Ice front retreat of marine-terminating glaciers may, for example, lead to a loss in resistive stress, which ultimately results in glacier acceleration and thinning. Over the past decade, it has been suggested that such retreats may be triggered by warm and salty Atlantic Water, which is typically found at a depth below 200–300 m. An increase in subglacial water discharge at glacier ice fronts due to enhanced surface runoff may also be responsible for an intensification of undercutting and calving. An increase in ocean thermal forcing or subglacial discharge therefore has the potential to destabilize marine-terminating glaciers along the coast of Greenland. It remains unclear which glaciers are currently stable but may retreat in the future and how far inland and how fast they will retreat. Here, we quantify the sensitivity and vulnerability of marine-terminating glaciers along the northwest coast of Greenland (from 72.5 to 76∘ N) to ocean forcing and subglacial discharge using the Ice Sheet System Model (ISSM). We rely on a parameterization of undercutting based on ocean thermal forcing and subglacial discharge and use ocean temperature and salinity from high-resolution ECCO2 (Estimating the Circulation and Climate of the Ocean, Phase II) simulations at the fjord mouth to constrain the ocean thermal forcing. The ice flow model includes a calving law based on a tensile von Mises criterion. We find that some glaciers, such as Dietrichson Gletscher or Alison Glacier, are sensitive to small increases in ocean thermal forcing, while others, such as Illullip Sermia or Cornell Gletscher, are remarkably stable, even in a +3 ∘C ocean warming scenario. Under the most intense experiment, we find that Hayes Gletscher retreats by more than 50 km inland by 2100 into a deep trough, and its velocity increases by a factor of 3 over only 23 years. The model confirms that ice–ocean interactions can trigger extensive and rapid glacier retreat, but the bed controls the rate and magnitude of the retreat. Under current oceanic and atmospheric conditions, we find that this sector of the Greenland ice sheet alone will contribute more than 1 cm to sea level rise and up to 3 cm by 2100 under the most extreme scenario.

Reconciling Drivers of Seasonal Terminus Advance and Retreat at 13 Central West Greenland Tidewater Glaciers

AbstractThe majority of Greenland tidewater glaciers undergo a seasonal cycle in terminus position, characterized by wintertime advance and summertime retreat. Understanding mechanisms that control seasonal cycles can help elucidate how tidewater glaciers regulate dynamic ice loss on longer timescales. However, controls on terminus position are numerous and complex. To address this, we compare time series of satellite‐derived terminus positions for tidewater glaciers in central west Greenland with observations of environmental forcings, including runoff at the grounding line, mélange presence, and, where available, ocean temperature in the proglacial fjord. We show that for most glaciers, seasonal terminus positions are more sensitive to glacial runoff than mélange or ocean thermal forcing. The strength of this relationship differs for two end‐member glacier types in the region, defined by their terminus geometry and dominant calving style. First, we find a strong relationship between magnitudes of runoff and terminus retreat at tidewater glaciers with shallow grounding lines (<400 m) that calve primarily through small‐magnitude serac failures. At these glaciers, subglacial plumes drive submarine melt and locally enhance retreat, causing heterogeneous position change across the terminus and local embayments where seasonal terminus changes are largest. In contrast, deep termini susceptible to buoyant flexure retreat sporadically through full ice thickness calving events less dependent on runoff. While less common, these glaciers deliver larger ice fluxes to the ocean. With predicted surface melt increases and diminished mélange coverage in a warming climate, our results reveal the impact of environmental forcings on diverse tidewater glacier systems in the region.

Effect of Topography on Subglacial Discharge and Submarine Melting During Tidewater Glacier Retreat

To first order, subglacial discharge depends on climate, which determines precipitation fluxes and glacier mass balance, and the rate of glacier volume change. For tidewater glaciers, large and rapid changes in glacier volume can occur independent of climate change due to strong glacier dynamic feedbacks. Using an idealized tidewater glacier model, we show that these feedbacks produce secular variations in subglacial discharge that are influenced by subglacial topography. Retreat along retrograde bed slopes (into deep water) results in rapid surface lowering and coincident increases in subglacial discharge. Consequently, submarine melting of glacier termini, which depends on subglacial discharge and ocean thermal forcing, also increases during retreat into deep water. Both subglacial discharge and submarine melting subsequently decrease as glacier termini retreat out of deep water and approach new steady state equilibria. In our simulations, subglacial discharge reached peaks that were 6–17% higher than preretreat values, with the highest values occurring during retreat from narrow sills, and submarine melting increased by 14% for unstratified fjords and 51% for highly stratified fjords. Our results therefore indicate that submarine melting acts in concert with iceberg calving to cause tidewater glacier termini to be unstable on retrograde beds. The full impact of submarine melting on tidewater glacier stability remains uncertain, however, due to poor understanding of the coupling between submarine melting and iceberg calving.