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: Earlier large-scale Greenland ice sheet sea-level projections (e.g. those run during the ice2sea and SeaRISE initiatives) have shown that ice sheet initial conditions have a large effect on the projections and give rise to important uncertainties. The goal of this initMIP-Greenland intercomparison exercise is to compare, evaluate, and improve the initialisation techniques used in the ice sheet modelling community and to estimate the associated uncertainties in modelled mass changes. initMIP-Greenland is the first in a series of ice sheet model intercomparison activities within ISMIP6 (the Ice Sheet Model Intercomparison Project for CMIP6), which is the primary activity within the Coupled Model Intercomparison Project Phase 6 (CMIP6) focusing on the ice sheets. Two experiments for the large-scale Greenland ice sheet have been designed to allow intercomparison between participating models of (1) the initial present-day state of the ice sheet and (2) the response in two idealised forward experiments. The forward experiments serve to evaluate the initialisation in terms of model drift (forward run without additional forcing) and in response to a large perturbation (prescribed surface mass balance anomaly); they should not be interpreted as sea-level projections. We present and discuss results that highlight the diversity of data sets, boundary conditions, and initialisation techniques used in the community to generate initial states of the Greenland ice sheet. We find good agreement across the ensemble for the dynamic response to surface mass balance changes in areas where the simulated ice sheets overlap but differences arising from the initial size of the ice sheet. The model drift in the control experiment is reduced for models that participated in earlier intercomparison exercises.

Description: Ice sheet numerical modeling is an important tool to estimate the dynamic contribution of the Antarctic ice sheet to sea level rise over the coming centuries. The influence of initial conditions on ice sheet model simulations, however, is still unclear. To better understand this influence, an initial state intercomparison exercise (initMIP) has been developed to compare, evaluate, and improve initialization procedures and estimate their impact on century-scale simulations. initMIP is the first set of experiments of the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6), which is the primary Coupled Model Intercomparison Project Phase 6 (CMIP6) activity focusing on the Greenland and Antarctic ice sheets. Following initMIP-Greenland, initMIP-Antarctica has been designed to explore uncertainties associated with model initialization and spin-up and to evaluate the impact of changes in external forcings. Starting from the state of the Antarctic ice sheet at the end of the initialization procedure, three forward experiments are each run for 100 years: a control run, a run with a surface mass balance anomaly, and a run with a basal melting anomaly beneath floating ice. This study presents the results of initMIP-Antarctica from 25 simulations performed by 16 international modeling groups. The submitted results use different initial conditions and initialization methods, as well as ice flow model parameters and reference external forcings. We find a good agreement among model responses to the surface mass balance anomaly but large variations in responses to the basal melting anomaly. These variations can be attributed to differences in the extent of ice shelves and their upstream tributaries, the numerical treatment of grounding line, and the initial ocean conditions applied, suggesting that ongoing efforts to better represent ice shelves in continental-scale models should continue.

Description: The mechanisms causing widespread flow acceleration of Jakobshavn Isbræ, West Greenland, remain unclear despite an abundance of observations and modeling studies. Here we simulate the glacier's evolution from 1985 to 2016 using a three-dimensional thermomechanical ice flow model. The model captures the timing and 90% of the observed changes by forcing the calving front. Basal drag in the trough is low, and lateral drag balances the ice stream's driving stress. The calving front position is the dominant control on changes of Jakobshavn Isbræ since the ice viscosity in the shear margins instantaneously drops in response to the stress perturbation caused by calving front retreat, which allows for widespread flow acceleration. Gradual shear margin warming contributes 5 to 10% to the total acceleration. Our simulations suggest that the glacier will contribute to eustatic sea level rise at a rate comparable to or higher than at present.

Description: Jakobshavn Isbræ, one of Greenland’s major outlet glaciers, displayed rapid changes since the mid-1990s. Its floating ice tongue broke up around 1997, followed by a rapid calving front retreat over tens of kilometres, possibly linked to a prior warming of the ocean waters adjacent to the fjord. Parallel to this major retreat, a quasi-simultaneous process of acceleration and thinning of the glacier has been observed, currently making it a major contributor to eustatic sea-level rise. However, the causal interplay between the various factors involved has not yet been fully understood. Numerical studies of Jakobshavn Isbræ so far are either 2-D plan view ice flow models with a fixed calving front or 2-D flowline models, which have to parameterize lateral stresses. Hence the interaction between changes in calving front position and ice dynamics could not be studied consistently. To overcome this limitation, we implemented an implicit boundary tracking method in the Ice Sheet System Model (ISSM). This tool allows us to freely evolve the calving front by prescribing a calving rate, the ice front velocity being therefore the sum of ice velocity and calving rate. A suite of sensitivity experiments perturbing an initial steady-state calving rate has been performed to study its impact on the dynamics of the glacier. Our numerical results suggest a high sensitivity of the glacier dynamics to the applied calving rate. Changes in calving rate quickly affect upstream areas of the ice stream through a combination of changes in calving front position, ice velocity, thickness, grounding and ungrounding, and surface gradient change. Consequently, acceleration triggered at the calving front quickly affects the entire drainage basin. Moreover, the model results suggest that the ice stream does not recover from a short duration (about a year) calving rate perturbation over timescales on the order of a century. We present selected results of the sensitivity experiments to support the discussion clarifying the causes of the current changes occurring at this dynamic ice stream.

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: Glacier-front dynamics is an important control on Greenland's ice mass balance. Warmer ocean waters trigger ice-front retreats of marine-terminating glaciers, and the corresponding loss in resistive stress leads to glacier acceleration and thinning. Here we present an approach to quantify the sensitivity and vulnerability of marine-terminating glaciers to ocean-induced melt. We develop a plan view model of Store Gletscher that includes a level set-based moving boundary capability, a parameterized ocean-induced melt, and a calving law with complete and precise land and fjord topographies to model the response of the glacier to increased melt. We find that the glacier is stabilized by a sill at its terminus. The glacier is dislodged from the sill when ocean-induced melt quadruples, at which point the glacier retreats irreversibly for 27 km into a reverse bed. The model suggests that ice-ocean interactions are the triggering mechanism of glacier retreat, but the bed controls its magnitude.

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: Calving is a major means of ice discharge of the Antarctic and Greenland Ice Sheets. The breaking off of icebergs changes the ice front configuration of marine terminating glaciers, which affects the stress regime of their upstream areas. Recent observations show the close correlation between the ice front position and the behaviour of many outlet glaciers. However, modelling of a glacier subject to calving poses various challenges. No universal calving rate parametrisation is known, and tracking of a moving ice front and the related boundary conditions in two or three spatial dimensions is non-trivial. Here, we present the theoretical and technical framework for a Level-Set Method, an implicit boundary tracking scheme, which we implemented into the Ice Sheet System Model (ISSM). The scheme allows us to study the dynamic response of a drainage basin to user-defined front ablation rates. We apply the method in a suite of experiments to Jakobshavn Isbræ, a major marine terminating outlet glacier of the western Greenland Ice Sheet. The model robustly reproduces the high sensitivity of the glacier to frontal ablation in form of calving. We find that enhanced calving is able to trigger significant acceleration of the ice stream. Upstream acceleration is sustained through a combination of various feedback mechanisms. However, lateral stress and ice influx into the trough are able to stabilise the ice stream. This study contributes to the present discussion on causes and effects of the continued changes occurring at Jakobshavn Isbræ, and emphasises that the incorporation of seasonal calving and dynamic lateral effects is key for realistic model projections of future global sea level rise on centennial time scales.