Description: Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Description: Recent observations highlight the influence of the calving front position on tidewater glaciers. Here, we present the theoretical and technical framework for a level-set method (LSM), which is now implemented into the Ice Sheet System Model. The LSM proves to be a robust tool to implicitly represent a dynamically evolving calving front in an Eulerian approach. We apply this method to a 3D thermodynamically coupled model of Jakobshavn Isbræ, West Greenland.

Description: Calving is a major mechanism of ice discharge of the Antarctic and Greenland ice sheets, and a change in calving front position affects the entire stress regime of marine terminating glaciers. The representation of calving front dynamics in a 2-D or 3-D ice sheet model remains non-trivial. Here, we present the theoretical and technical framework for a level-set method, an implicit boundary tracking scheme, which we implement into the Ice Sheet System Model (ISSM). This scheme allows us to study the dynamic response of a drainage basin to user-defined calving rates. We apply the method to Jakobshavn Isbræ, a major marine terminating outlet glacier of the West Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to calving, and we find that enhanced calving triggers significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of mechanisms. However, both lateral stress and ice influx stabilize the ice stream. This study provides new insights into the ongoing changes occurring at Jakobshavn Isbræ and emphasizes that the incorporation of moving boundaries and dynamic lateral effects, not captured in flow-line models, is key for realistic model projections of sea level rise on centennial timescales.

Description: As ice flows over a subglacial lake, the drop in bed resistance leads to an increase in ice velocities and a draw down of isochrones and cold ice. The ice surface flattens as it adjusts to the lack of resisting forces at the base. The rapid transition in velocity induces changes in ice viscosity and releases deformation energy that can raise the temperature locally. Recent studies of Antarctic subglacial lakes indicate that many lakes experience very fast and possibly episodic drainage, during which the lake size is rapidly reduced as water flows out. Questions that arise are what effect this would have on internal layers within the ice and whether such past drainage events could be inferred from isochrone structures downstream. Here, we study the effect of a subglacial lake on ice dynamics as well as the influence that such short timescale drainage would have on the internal layers of the ice. To this end, we use a full Stokes, polythermal ice flow model. An enthalpy-gradient method is used to account for the evolution of temperature and water content within the ice. We find that a rapid transition between slow-moving ice outside the lake, and full sliding over the lake, can release considerable amounts of deformational energy, with the potential to form a temperate layer at depth in the transition zone. In addition, we provide an explanation for a characteristic surface feature commonly seen at the edges of subglacial lakes, a hummocky surface depression in the transition zone between little to full sliding. We also conclude that rapid changes in the horizontal extent of subglacial lakes and slippery patches, compared to the average ice column velocity, can create a traveling wave at depth within the isochrone structure that transfers downstream with the advection of ice, thus indicating the possibility of detecting past drainage events with ice penetrating radar.