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

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.