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.
Computer models are necessary for understanding and predicting marine ice sheet behaviour. However, there is uncertainty over implementation of physical processes at the ice base, both for grounded and floating glacial ice. Here we implement several sliding relations in a marine ice sheet flow-line model accounting for all stress components and demonstrate that model resolution requirements are strongly dependent on both the choice of basal sliding relation and the spatial distribution of ice shelf basal melting.Sliding relations that reduce the magnitude of the step change in basal drag from grounded ice to floating ice (where basal drag is set to zero) show reduced dependence on resolution compared to a commonly used relation, in which basal drag is purely a power law function of basal ice velocity. Sliding relations in which basal drag goes smoothly to zero as the grounding line is approached from inland (due to a physically motivated incorporation of effective pressure at the bed) provide further reduction in resolution dependence.A similar issue is found with the imposition of basal melt under the floating part of the ice shelf: melt parameterisations that reduce the abruptness of change in basal melting from grounded ice (where basal melt is set to zero) to floating ice provide improved convergence with resolution compared to parameterisations in which high melt occurs adjacent to the grounding line.Thus physical processes, such as sub-glacial outflow (which could cause high melt near the grounding line), impact on capability to simulate marine ice sheets. If there exists an abrupt change across the grounding line in either basal drag or basal melting, then high resolution will be required to solve the problem. However, the plausible combination of a physical dependency of basal drag on effective pressure, and the possibility of low ice shelf basal melt rates next to the grounding line, may mean that some marine ice sheet systems can be reliably simulated at a coarser resolution than currently thought necessary.
Ekici, A.; Chadburn, S.; Chaudhary, N.; Hajdu, L. H.; Marmy, A.; Peng, S.; Boike, J.; Burke, E.; Friend, A. D.; Hauck, C.; +4 more
Ekici, A.; Chadburn, S.; Chaudhary, N.; Hajdu, L. H.; Marmy, A.; Peng, S.; Boike, J.; Burke, E.; Friend, A. D.; Hauck, C.; Krinner, G.; Langer, M.; Miller, P. A.; Beer, C.;
Project: EC | PAGE21 (282700), EC | GREENCYCLESII (238366), SNSF | The evolution of mountain... (136279)
Modeling soil thermal dynamics at high latitudes and altitudes requires representations of physical processes such as snow insulation, soil freezing and thawing and subsurface conditions like soil water/ice content and soil texture. We have compared six different land models: JSBACH, ORCHIDEE, JULES, COUP, HYBRID8 and LPJ-GUESS, at four different sites with distinct cold region landscape types, to identify the importance of physical processes in capturing observed temperature dynamics in soils. The sites include alpine, high Arctic, wet polygonal tundra and non-permafrost Arctic, thus showing how a range of models can represent distinct soil temperature regimes. For all sites, snow insulation is of major importance for estimating topsoil conditions. However, soil physics is essential for the subsoil temperature dynamics and thus the active layer thicknesses. This analysis shows that land models need more realistic surface processes, such as detailed snow dynamics and moss cover with changing thickness and wetness, along with better representations of subsoil thermal dynamics.