Evolution of ocean-induced ice melt beneath Zachariæ Isstrøm, Northeast Greenland combining observations and an ocean general circulation model from 1978 to present

Zachariæ Isstrøm, a major ice stream in northeast Greenland, lost a large fraction of its ice shelf during the last decade. We study the evolution of subaqueous melting of its floating section from 1978 to present. The ice shelf melt rate depends on thermal forcing from warm, salty, subsurface ocean waters of Atlantic origin (AW), the mixing of AW with fresh, buoyant subglacial discharge at the calving margin, and the shape of the sub-ice-shelf cavity. Subglacial discharge doubled as a result of enhanced ice sheet runoff caused by warmer air temperatures. Ocean thermal forcing has increased due to enhanced advection of AW. Using an Eulerian method, MEaSUREs ice velocity, Operation IceBridge (OIB) ice thickness, and RACMO2.3 surface balance data, we evaluate the ice shelf melt rate in 1978, 1999 and 2010. The melt rate doubled from 1999 to 2010. Using a Lagrangian method with World View imagery, we map the melt rate in detail from 2011 to 2016. We compare the results with 2D simulations from the Massachusetts Institute of Technology general circulation model (MITgcm), at a high spatial resolution (20-m horizontal and 40-m vertical grid spacing), using OIB ice thickness and sub-ice-shelf cavity for years 1978, 1996, 2010 and 2011, combined with in-situ ocean temperature/salinity data from Ocean Melting Greenland (OMG) 2017. We find that winter melt rates are 2 3 times smaller than summer rates and melt rates increase by one order magnitude during the transition from ice shelf termination to near-vertical calving wall termination. As the last remaining bits of floating ice shelf disappear, ice-ocean interaction will therefore play an increasing role in driving the glacier retreat into its marine-based basin. This work was performed under a contract with NASA Cryosphere Program at UC Irvine and Caltech's Jet Propulsion Laboratory.