• European Marine Science
• Other research productsclose
• DK close
• US close
• The Cryosphere (TC) close

Predicting the climate for the future and how it will impact ice sheet evolution requires coupling ice sheet models with climate models. However, before we attempt to develop a realistic coupled setup, we propose, in this study, to first analyse the impact of a model simulated climate on an ice sheet. We undertake this exercise for a set of regional and global climate models. Modelled near surface air temperature and precipitation are provided as upper boundary conditions to the GRISLI (GRenoble Ice Shelf and Land Ice model) hybrid ice sheet model (ISM) in its Greenland configuration. After 20 kyrs of simulation, the resulting ice sheets highlight the differences between the climate models. While modelled ice sheet sizes are generally comparable to the observed one, there are considerable deviations among the ice sheets on regional scales. These deviations can be explained by biases in temperature and precipitation near the coast. This is especially true in the case of global models. But the deviations between the climate models are also due to the differences in the atmospheric general circulation. To account for these differences in the context of coupling ice sheet models with climate models, we conclude that appropriate downscaling methods will be needed. In some cases, systematic corrections of the climatic variables at the interface may be required to obtain realistic results for the Greenland ice sheet (GIS).

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.

The Northeast Greenland Ice Stream (NEGIS) currently drains more than 10 % of the Greenland Ice Sheet area and has recently undergone significant dynamic changes. It is therefore critical to accurately represent this feature when assessing the future contribution of Greenland to sea level rise. At present, NEGIS is reproduced in ice sheet models by inferring basal conditions using observed surface velocities. This approach helps estimate conditions at the base of the ice sheet but cannot be used to estimate the evolution of basal drag in time, so it is not a good representation of the evolution of the ice sheet in future climate warming scenarios. NEGIS is suggested to be initiated by a geothermal heat flux anomaly close to the ice divide, left behind by the movement of Greenland over the Icelandic plume. However, the heat flux underneath the ice sheet is largely unknown, except for a few direct measurements from deep ice core drill sites. Using the Ice Sheet System Model (ISSM), with ice dynamics coupled to a subglacial hydrology model, we investigate the possibility of initiating NEGIS by inserting heat flux anomalies with various locations and intensities. In our model experiment, a minimum heat flux value of 970 mW m−2 located close to the East Greenland Ice-core Project (EGRIP) is required locally to reproduce the observed NEGIS velocities, giving basal melt rates consistent with previous estimates. The value cannot be attributed to geothermal heat flux alone and we suggest hydrothermal circulation as a potential explanation for the high local heat flux. By including high heat flux and the effect of water on sliding, we successfully reproduce the main characteristics of NEGIS in an ice sheet model without using data assimilation.