European Ice Sheet Modelling Initiative (EISMINT) model intercomparison experiments with first-order mechanics

Author(s):  
Fuyuki Saito ◽  
Ayako Abe-Ouchi ◽  
Heinz Blatter
2014 ◽  
Vol 26 (6) ◽  
pp. 674-686 ◽  
Author(s):  
C.J. Fogwill ◽  
C.S.M. Turney ◽  
N.R. Golledge ◽  
D.H. Rood ◽  
K. Hippe ◽  
...  

AbstractDetermining the millennial-scale behaviour of marine-based sectors of the West Antarctic Ice Sheet (WAIS) is critical to improve predictions of the future contribution of Antarctica to sea level rise. Here high-resolution ice sheet modelling was combined with new terrestrial geological constraints (in situ14C and 10Be analysis) to reconstruct the evolution of two major ice streams entering the Weddell Sea over 20 000 years. The results demonstrate how marked differences in ice flux at the marine margin of the expanded Antarctic ice sheet led to a major reorganization of ice streams in the Weddell Sea during the last deglaciation, resulting in the eastward migration of the Institute Ice Stream, triggering a significant regional change in ice sheet mass balance during the early to mid Holocene. The findings highlight how spatial variability in ice flow can cause marked changes in the pattern, flux and flow direction of ice streams on millennial timescales in this marine ice sheet setting. Given that this sector of the WAIS is assumed to be sensitive to ocean-forced instability and may be influenced by predicted twenty-first century ocean warming, our ability to model and predict abrupt and extensive ice stream diversions is key to a realistic assessment of future ice sheet sensitivity.


2020 ◽  
Author(s):  
Sophie Nowicki ◽  
Antony J. Payne ◽  
Heiko Goelzer ◽  
Helene Seroussi ◽  
William H. Lipscomb ◽  
...  

Abstract. Projection of the contribution of ice sheets to sea-level change as part of the Coupled Model Intercomparison Project – phase 6 (CMIP6) takes the form of simulations from coupled ice-sheet-climate models and standalone ice sheet models, overseen by the Ice Sheet Model Intercomparison Project for CMIP6 (ISMIP6). This paper describes the experimental setup for process-based sea-level change projections to be performed with standalone Greenland and Antarctic ice sheet models in the context of ISMIP6. The ISMIP6 protocol relies on a suite of polar atmospheric and oceanic CMIP-based forcing for ice sheet models, in order to explore the uncertainty in projected sea-level change due to future emissions scenarios, CMIP models, ice sheet models, and parameterizations for ice-ocean interactions. We describe here the approach taken for defining the suite of ISMIP6 standalone ice sheet simulations, document the experimental framework and implementation, as well as present an overview of the ISMIP6 forcing to be used by participating ice sheet modeling groups.


2019 ◽  
Vol 13 (5) ◽  
pp. 1441-1471 ◽  
Author(s):  
Hélène Seroussi ◽  
Sophie Nowicki ◽  
Erika Simon ◽  
Ayako Abe-Ouchi ◽  
Torsten Albrecht ◽  
...  

Abstract. 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.


2012 ◽  
Vol 5 (4) ◽  
pp. 963-974 ◽  
Author(s):  
A. M. Dolan ◽  
S. J. Koenig ◽  
D. J. Hill ◽  
A. M. Haywood ◽  
R. M. DeConto

Abstract. During the mid-Pliocene warm period (3.264 to 3.025 million years ago), global mean temperature was similar to that predicted for the next century and atmospheric carbon dioxide concentrations were slightly higher than today. Sea level was also higher than today, implying a reduction in the extent of the ice sheets. Thus, the mid-Pliocene warm period (mPWP) provides a unique testing ground to investigate the stability of the Earth's ice sheets and their contribution to sea level in a warmer-than-modern world. Climate models and ice sheet models can be used to enhance our understanding of ice sheet stability; however, uncertainties associated with different ice-sheet modelling frameworks mean that a rigorous comparison of numerical ice sheet model simulations for the Pliocene is essential. As an extension to the Pliocene Model Intercomparison Project (PlioMIP; Haywood et al., 2010, 2011a), the Pliocene Ice Sheet Modelling Intercomparison Project (PLISMIP) will provide the first assessment as to the ice sheet model dependency of ice sheet predictions for the mPWP. Here we outline the PLISMIP experimental design and initialisation conditions that have been adopted to simulate the Greenland and Antarctic ice sheets under present-day and warm mid-Pliocene conditions. Not only will this project provide a new benchmark in the simulation of ice sheets in a past warm period, but the analysis of model sensitivity to various uncertainties could directly inform future predictions of ice sheet and sea level change.


2018 ◽  
Author(s):  
Matthew J. Hoffman ◽  
Mauro Perego ◽  
Stephen F. Price ◽  
William H. Lipscomb ◽  
Douglas Jacobsen ◽  
...  

Abstract. We introduce MPAS-Albany Land Ice (MALI), a new, variable resolution land ice model that uses unstructured Voronoi grids on a plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework for developing variable resolution Earth System Model components and the Albany multi-physics code base for solution of coupled systems of partial-differential equations, which itself makes use of Trilinos solver libraries. MALI includes a three-dimensional, first-order momentum balance solver ("Blatter-Pattyn") by linking to the Albany-LI ice sheet velocity solver, as well as an explicit shallow ice velocity solver. Evolution of ice geometry and tracers is handled through an explicit first-order horizontal advection scheme with vertical remapping. Evolution of ice temperature is treated using operator splitting of vertical diffusion and horizontal advection and can be configured to use either a temperature or enthalpy formulation. MALI includes a mass-conserving subglacial hydrology model that supports distributed and/or channelized drainage and can optionally be coupled to ice dynamics. Options for calving include "eigencalving", which assumes calving rate is proportional to extensional strain rates. MALI is evaluated against commonly used exact solutions and community benchmark experiments and shows the expected accuracy. We report first results for the MISMIP3d benchmark experiments for a Blatter-Pattyn type model and show that results fall in-between those of models using Stokes flow and L1L2 approximations. We use the model to simulate a semi-realistic Antarctic Ice Sheet problem for 1100 years at 20 km resolution. MALI is the glacier component of the Energy Exascale Earth System Model (E3SM) version 1, and we describe current and planned coupling to other components.


1990 ◽  
Vol 14 ◽  
pp. 347-347
Author(s):  
V.L. Mazo

Tidewater glaciers and large ice sheets, e.g. the Antarctic ice sheet and a late-Würm Arctic ice sheet, are complex but single dynamic systems composed of terrestrial, marine and floating parts. Morphology and dynamics of the different parts are different. The terrestrial parts are convex and their dynamics are controlled by shear stress only (the longitudinal stress is zero); the floating parts are concave and their dynamics are controlled by longitudinal stress only (the shear stress is zero). To connect the different parts we should consider transitional zones where shear and longitudinal stresses are comparable.To describe glacier and ice-sheet dynamics, longwave approximation of the first order is used. In this approximation it is impossible to connect terrestrial and floating parts dynamically, only morphologically and kinematically. It means that the first-order longwave approximation is not sufficient.If the transitional zone between the terrestrial and floating parts is long in comparison to ice thickness (in hydrodynamics the term “weak” is used) we can do the next step in the longwave approximation to describe the single dynamical system consisting of the terrestrial and floating parts and the weak transitional zones (ice streams). It is a purely hydrodynamical approach to the problem without ad hoc hypothesis.The presented model is a non-stationary three-dimensional hydrodynamic model of glaciers and ice sheets interacted with ocean, involving the conditions of ice continuity and dynamic equilibrium, ice rheology, and boundary conditions on the free surface (dynamic and kinematic) and on the bed (ice freezing or sliding). Longwave approximation is used to reduce the three-dimensional model to a two-dimensional one. The latter consists of (1) evolution equations for grounded and floating parts and weak transitional zones; (2) boundary conditions on the fronts (e.g. the conditions of calving); (3) equations governing the junctions of the parts (the most important junction is the grounded line) with the conditions connecting the evolution equations.


2018 ◽  
Vol 11 (9) ◽  
pp. 3747-3780 ◽  
Author(s):  
Matthew J. Hoffman ◽  
Mauro Perego ◽  
Stephen F. Price ◽  
William H. Lipscomb ◽  
Tong Zhang ◽  
...  

Abstract. We introduce MPAS-Albany Land Ice (MALI) v6.0, a new variable-resolution land ice model that uses unstructured Voronoi grids on a plane or sphere. MALI is built using the Model for Prediction Across Scales (MPAS) framework for developing variable-resolution Earth system model components and the Albany multi-physics code base for the solution of coupled systems of partial differential equations, which itself makes use of Trilinos solver libraries. MALI includes a three-dimensional first-order momentum balance solver (Blatter–Pattyn) by linking to the Albany-LI ice sheet velocity solver and an explicit shallow ice velocity solver. The evolution of ice geometry and tracers is handled through an explicit first-order horizontal advection scheme with vertical remapping. The evolution of ice temperature is treated using operator splitting of vertical diffusion and horizontal advection and can be configured to use either a temperature or enthalpy formulation. MALI includes a mass-conserving subglacial hydrology model that supports distributed and/or channelized drainage and can optionally be coupled to ice dynamics. Options for calving include “eigencalving”, which assumes that the calving rate is proportional to extensional strain rates. MALI is evaluated against commonly used exact solutions and community benchmark experiments and shows the expected accuracy. Results for the MISMIP3d benchmark experiments with MALI's Blatter–Pattyn solver fall between published results from Stokes and L1L2 models as expected. We use the model to simulate a semi-realistic Antarctic ice sheet problem following the initMIP protocol and using 2 km resolution in marine ice sheet regions. MALI is the glacier component of the Energy Exascale Earth System Model (E3SM) version 1, and we describe current and planned coupling to other E3SM components.


2020 ◽  
Author(s):  
Baptiste Vandecrux ◽  
Ruth Mottram ◽  
Peter L. Langen ◽  
Robert S. Fausto ◽  
Martin Olesen ◽  
...  

Abstract. Perennial snow, or firn, covers 80 % of the Greenland ice sheet and has the capacity to retain part of the surface meltwater, buffering the ice sheet’s contribution to sea level. Multi-layer firn models are traditionally used to simulate the firn processes and estimate meltwater retention. We present the output from nine firn models, forced by weather-station-derived mass and energy fluxes at four sites representative of the dry snow, percolation, ice slab and firn aquifer areas. We compare the model outputs and evaluate them against in situ observations. Models that explicitly account for deep meltwater percolation overestimate percolation depth and consequently firn temperature at the percolation and ice slab sites although they accurately simulate the recharge of the firn aquifer. Models using Darcy's law and a bucket scheme compare favourably to observations at the percolation site but only the Darcy models accurately simulate firn temperature and thus meltwater percolation at the ice slab site. We find that Eulerian models, that transfer firn through fixed layers, smooth sharp gradients in firn temperature and density over time. From the model spread, we find that simulated densities (respectively temperature) have an uncertainty envelope of ±60 kg m−3 (resp. ±14 °C) in the dry snow area and up to ±280 kg m−3 (resp. ±15–18 °C) at warmer sites.


Sign in / Sign up

Export Citation Format

Share Document