Simulating interactive ice sheets in the multi-resolution AWI-ESM 1.2: A case study using SCOPE 1.0

Abstract. Interactions between the climate and the cryosphere have the potential to induce strong non-linear transitions in the Earth's climate. These interactions influence both the atmospheric circulation, by changing the ice sheet's geometry, as well as the oceanic circulation, by modification of the water mass properties. Furthermore, the waxing and waning of large continental ice sheets influences the global albedo, altering the energy balance of the Earth System and inducing climate-ice sheet feedbacks on a global scale as evident in Pleistocene glacial-interglacial cycles. To date, few fully comprehensive models exist, that do not only contain a coupled atmosphere/land/ocean component, but also consider interactive cryosphere physics. Yet, on glacial-interglacial and tectonic time scales, as well as in the Anthropocene, ice sheets are not in equilibrium with the climate, and prescribed fixed ice sheet representations in the model can principally be only an approximation to reality. Only climate models, that contain interactive ice sheets, can produce simulations of the Earth's climate which include all feedbacks and processes related to atmosphere-land-ocean-ice interactions. Previous fully coupled models were limited either by low spatial resolution or an incomplete representation of ice sheet processes, such as iceberg calving, surface ablation processes, and ocean/ice-shelf interactions. Here, we present the newly developed AWI-Earth System Model (AWI-ESM), which tackles some of these problems. Our modelling toolbox is based on the AWI-climate model, including atmosphere and vegetation components suitable for paleoclimate studies, a multi-resolution global ocean component which can be refined to simulate regions of interest at high resolution, and an ice sheet component suitable for simulating both ice sheet and ice shelf dynamics and thermodynamics. We describe the currently implemented coupling between these components, present first results for the Mid-Holocene and Last Interglacial, and introduce further ideas for scientific applications for both future and past climate states with a focus on the Northern Hemisphere. Finally, we provide an outlook on the potential of such fully coupled Earth System models in improving representation of climate-ice sheet feedbacks in future paleoclimate studies with this model.

Download Full-text

Simulating Interactive Ice Sheets in the Multi-Resolution AWI-ESM: A case study using the SCOPE Coupler

10.5194/egusphere-egu21-15682 ◽

2021 ◽

Author(s):

Paul Gierz ◽

Lars Ackermann ◽

Christian Rodehacke ◽

Uta Krebs-Kanzow ◽

Christian Stepanek ◽

...

Keyword(s):

Climate Model ◽

Ice Sheets ◽

Ice Sheet ◽

Global Ocean ◽

Earth System ◽

Last Interglacial ◽

Ice Shelf ◽

Coupled Models ◽

Fully Coupled ◽

Paleoclimate Studies

Interactions between the climate and the cryosphere have the&#160;potential to induce strong non-linear transitions in the Earth's&#160;climate. These interactions influence both the atmospheric circulation, by changing the ice sheet's geometry, as well as the oceanic circulation, by modification of the water mass properties. Furthermore, the waxing and waning of large continental ice sheets influences the global albedo, altering the energy balance of the Earth System and inducing climate-ice sheet feedbacks on a global scale as evident in Pleistocene glacial-interglacial cycles. To date, few fully comprehensive models exist, that do not only contain a coupled atmosphere/land/ocean component, but also consider interactive cryosphere physics. Yet, on glacial-interglacial and tectonic time scales, as well as in the Anthropocene, ice sheets are not in equilibrium with the climate, and prescribed fixed ice sheet representations in the model can principally be only an approximation to reality. Only climate models, that contain interactive ice sheets, can produce simulations of the Earth's climate which include all feedbacks and processes related to atmosphere-land-ocean-ice interactions. Previous fully coupled models were limited either by low spatial resolution or an incomplete representation of ice sheet processes, such as iceberg calving, surface ablation processes, and ocean/ice-shelf interactions. Here, we present the newly developed AWI-Earth System Model (AWI-ESM), which tackles some of these problems. Our modelling toolbox is based on the AWI-climate model, including&#160;atmosphere and vegetation components suitable for paleoclimate studies, a multi-resolution global ocean component which can be refined to simulate regions of interest at high resolution, and an ice sheet component suitable for simulating both ice sheet and ice shelf dynamics and thermodynamics. We describe the currently implemented coupling between these components, present first results for the Mid-Holocene and Last Interglacial, and introduce further ideas for scientific applications for both future and past climate states with a focus on the Northern Hemisphere. Finally, we provide an outlook on the potential of such fully coupled Earth System models in improving representation of climate-ice sheet feedbacks in future paleoclimate studies with this model.

Download Full-text

Asynchronous Antarctic and Greenland ice-volume contributions to the last interglacial sea-level highstand

10.5194/egusphere-egu2020-1513 ◽

2020 ◽

Author(s):

Eelco Rohling ◽

Fiona Hibbert

Keyword(s):

Sea Level Rise ◽

Sea Level ◽

Global Climate ◽

Ice Sheets ◽

Ice Sheet ◽

Greenland Ice Sheet ◽

Climate Forcing ◽

Potential Contribution ◽

Last Interglacial ◽

Ice Shelf

Sea-level rise is among the greatest risks that arise from anthropogenic global climate change. It is receiving a lot of attention, among others in the IPCC reports, but major questions remain as to the potential contribution from the great continental ice sheets. In recent years, some modelling work has suggested that the ice-component of sea-level rise may be much faster than previously thought, but the rapidity of rise seen in these results depends on inclusion of scientifically debated mechanisms of ice-shelf decay and associated ice-sheet instability. The processes have not been active during historical times, so data are needed from previous warm periods to evaluate whether the suggested rates of sea-level rise are supported by observations or not. Also, we then need to assess which of the ice sheets was most sensitive, and why. The last interglacial (LIG; ~130,000 to ~118,000 years ago, ka) was the last time global sea level rose well above its present level, reaching a highstand of +6 to +9 m or more. Because Greenland Ice Sheet (GrIS) contributions were smaller than that, this implies substantial Antarctic Ice Sheet (AIS) contributions. However, this still leaves the timings, magnitudes, and drivers of GrIS and AIS reductions open to debate. I will discuss recently published sea-level reconstructions for the LIG highstand, which reveal that AIS and GrIS contributions were distinctly asynchronous, and that rates of rise to values above 0 m (present-day sea level) reached up to 3.5 m per century. Such high pre-anthropogenic rates of sea-level rise lend credibility to high rates inferred by ice modelling under certain ice-shelf instability parameterisations, for both the past and future. Climate forcing was distinctly asynchronous between the southern and northern hemispheres as well during the LIG, explaining the asynchronous sea-level contributions from AIS and GrIS. Today, climate forcing is synchronous between the two hemispheres, and also faster and greater than during the LIG. Therefore, LIG rates of sea-level rise should likely be considered minimum estimates for the future.

Download Full-text

The Framework For Ice Sheet–Ocean Coupling (FISOC) V1.1

Geoscientific Model Development ◽

10.5194/gmd-14-889-2021 ◽

2021 ◽

Vol 14 (2) ◽

pp. 889-905

Author(s):

Rupert Gladstone ◽

Benjamin Galton-Fenzi ◽

David Gwyther ◽

Qin Zhou ◽

Tore Hattermann ◽

...

Keyword(s):

Dynamic Model ◽

Ice Sheets ◽

Ice Sheet ◽

Coupled Processes ◽

Earth System ◽

Ice Shelf ◽

Time Step ◽

Ice Dynamics ◽

Ocean Coupling ◽

The Impact

Abstract. A number of important questions concern processes at the margins of ice sheets where multiple components of the Earth system, most crucially ice sheets and oceans, interact. Such processes include thermodynamic interaction at the ice–ocean interface, the impact of meltwater on ice shelf cavity circulation, the impact of basal melting of ice shelves on grounded ice dynamics and ocean controls on iceberg calving. These include fundamentally coupled processes in which feedback mechanisms between ice and ocean play an important role. Some of these mechanisms have major implications for humanity, most notably the impact of retreating marine ice sheets on the global sea level. In order to better quantify these mechanisms using computer models, feedbacks need to be incorporated into the modelling system. To achieve this, ocean and ice dynamic models must be coupled, allowing runtime information sharing between components. We have developed a flexible coupling framework based on existing Earth system coupling technologies. The open-source Framework for Ice Sheet–Ocean Coupling (FISOC) provides a modular approach to coupling, facilitating switching between different ice dynamic and ocean components. FISOC allows fully synchronous coupling, in which both ice and ocean run on the same time step, or semi-synchronous coupling in which the ice dynamic model uses a longer time step. Multiple regridding options are available, and there are multiple methods for coupling the sub-ice-shelf cavity geometry. Thermodynamic coupling may also be activated. We present idealized simulations using FISOC with a Stokes flow ice dynamic model coupled to a regional ocean model. We demonstrate the modularity of FISOC by switching between two different regional ocean models and presenting outputs for both. We demonstrate conservation of mass and other verification steps during evolution of an idealized coupled ice–ocean system, both with and without grounding line movement.

Download Full-text

A Framework for Ice Sheet – Ocean Coupling (FISOC) V1.1

10.5194/gmd-2020-206 ◽

2020 ◽

Author(s):

Rupert Gladstone ◽

Benjamin Galton-Fenzi ◽

David Gwyther ◽

Qin Zhou ◽

Tore Hattermann ◽

...

Keyword(s):

Dynamic Model ◽

Ice Sheets ◽

Ice Sheet ◽

Coupled Processes ◽

Earth System ◽

Ice Shelf ◽

Time Step ◽

Ice Dynamics ◽

Ocean Coupling ◽

The Impact

Abstract. A number of important questions concern processes at the margins of ice sheets where multiple components of the Earth System, most crucially ice sheets and oceans, interact. Such processes include thermodynamic interaction at the ice-ocean interface, the impact of melt water on ice shelf cavity circulation, the impact of basal melting of ice shelves on grounded ice dynamics, and ocean controls on iceberg calving. These include fundamentally coupled processes in which feedback mechanisms between ice and ocean play an important role. Some of these mechanisms have major implications for humanity, most notably the impact of retreating marine ice sheets on global sea level. In order to better quantify these mechanisms using computer models, feedbacks need to be incorporated into the modelling system. To achieve this ocean and ice dynamic models must be coupled, allowing run time information sharing between components. We have developed a flexible coupling framework based on existing Earth System coupling technologies. The open-source Framework for Ice Sheet – Ocean Coupling (FISOC) provides a modular approach to online coupling, facilitating switching between different ice dynamic and ocean components. FISOC allows fully synchronous coupling, in which both ice and ocean run on the same time-step, or semi-synchronous coupling in which the ice dynamic model uses a longer time step. Multiple regridding options are available, and multiple methods for coupling the sub ice shelf cavity geometry. Thermodynamic coupling may also be activated. We present idealised simulations using FISOC with a Stokes flow ice dynamic model coupled to a regional ocean model. We demonstrate the modularity of FISOC by switching between two different regional ocean models and presenting outputs for both. We demonstrate conservation of mass and other verification steps during evolution of an idealised coupled ice – ocean system, both with and without grounding line movement.

Download Full-text

Dynamical processes involved in the retreat of marine ice sheets

Journal of Glaciology ◽

10.3189/172756501781832269 ◽

2001 ◽

Vol 47 (157) ◽

pp. 271-282 ◽

Cited By ~ 95

Author(s):

Richard C.A. Hindmarsh ◽

E. Le Meur

Keyword(s):

Ice Sheets ◽

Ice Sheet ◽

Ice Shelf ◽

Ice Streams ◽

Antarctic Ice Sheet ◽

West Antarctic Ice Sheet ◽

West Antarctic ◽

Grounding Line ◽

Neutral Equilibrium ◽

Generally True

AbstractMarine ice sheets with mechanics described by the shallow-ice approximation by definition do not couple mechanically with the shelf. Such ice sheets are known to have neutral equilibria. We consider the implications of this for their dynamics and in particular for mechanisms which promote marine ice-sheet retreat. The removal of ice-shelf buttressing leading to enhanced flow in grounded ice is discounted as a significant influence on mechanical grounds. Sea-level rise leading to reduced effective pressures under ice streams is shown to be a feasible mechanism for producing postglacial West Antarctic ice-sheet retreat but is inconsistent with borehole evidence. Warming thins the ice sheet by reducing the average viscosity but does not lead to grounding-line retreat. Internal oscillations either specified or generated via a MacAyeal–Payne thermal mechanism promote migration. This is a noise-induced drift phenomenon stemming from the neutral equilibrium property of marine ice sheets. This migration occurs at quite slow rates, but these are sufficiently large to have possibly played a role in the dynamics of the West Antarctic ice sheet after the glacial maximum. Numerical experiments suggest that it is generally true that while significant changes in thickness can be caused by spatially uniform changes, spatial variability coupled with dynamical variability is needed to cause margin movement.

Download Full-text

Last Interglacial climate and sea-level evolution from a coupled ice sheet–climate model

Climate of the Past ◽

10.5194/cp-12-2195-2016 ◽

2016 ◽

Vol 12 (12) ◽

pp. 2195-2213 ◽

Cited By ~ 22

Author(s):

Heiko Goelzer ◽

Philippe Huybrechts ◽

Marie-France Loutre ◽

Thierry Fichefet

Keyword(s):

Surface Temperature ◽

Sea Level ◽

Ice Sheets ◽

Ice Sheet ◽

Greenland Ice Sheet ◽

Last Interglacial ◽

A Minor ◽

The Antarctic ◽

The Impact ◽

And Sea Level

Abstract. As the most recent warm period in Earth's history with a sea-level stand higher than present, the Last Interglacial (LIG, ∼  130 to 115 kyr BP) is often considered a prime example to study the impact of a warmer climate on the two polar ice sheets remaining today. Here we simulate the Last Interglacial climate, ice sheet, and sea-level evolution with the Earth system model of intermediate complexity LOVECLIM v.1.3, which includes dynamic and fully coupled components representing the atmosphere, the ocean and sea ice, the terrestrial biosphere, and the Greenland and Antarctic ice sheets. In this setup, sea-level evolution and climate–ice sheet interactions are modelled in a consistent framework.Surface mass balance change governed by changes in surface meltwater runoff is the dominant forcing for the Greenland ice sheet, which shows a peak sea-level contribution of 1.4 m at 123 kyr BP in the reference experiment. Our results indicate that ice sheet–climate feedbacks play an important role to amplify climate and sea-level changes in the Northern Hemisphere. The sensitivity of the Greenland ice sheet to surface temperature changes considerably increases when interactive albedo changes are considered. Southern Hemisphere polar and sub-polar ocean warming is limited throughout the Last Interglacial, and surface and sub-shelf melting exerts only a minor control on the Antarctic sea-level contribution with a peak of 4.4 m at 125 kyr BP. Retreat of the Antarctic ice sheet at the onset of the LIG is mainly forced by rising sea level and to a lesser extent by reduced ice shelf viscosity as the surface temperature increases. Global sea level shows a peak of 5.3 m at 124.5 kyr BP, which includes a minor contribution of 0.35 m from oceanic thermal expansion. Neither the individual contributions nor the total modelled sea-level stand show fast multi-millennial timescale variations as indicated by some reconstructions.

Download Full-text

Response to Filchner-Ronne Ice Shelf cavity warming in a coupled ocean–ice sheet model. Part I: The ocean perspective

10.5194/os-2017-41 ◽

2017 ◽

Cited By ~ 1

Author(s):

Ralph Timmermann ◽

Sebastian Goeller

Keyword(s):

World Ocean ◽

Ice Sheet ◽

Coupled Model ◽

Freezing Temperature ◽

Ocean Model ◽

Climate Forcing ◽

Global Ocean ◽

Ice Shelf ◽

Marginal Seas ◽

Grounding Line

Abstract. A Regional Antarctic and Global Ocean (RAnGO) model has been developed to study the interaction between the world ocean and the Antarctic ice sheet. The coupled model is based on a global implementation of the Finite Element Sea-ice Ocean Model (FESOM) with a mesh refinement in the Southern Ocean, particularly in its marginal seas and in the sub-ice shelf cavities. The cryosphere is represented by a regional setup of the ice flow model RIMBAY comprising the Filchner-Ronne Ice Shelf and the grounded ice in its catchment area up to the ice divides. At the base of the RIMBAY ice shelf, melt rates from FESOM's ice-shelf component are supplied. RIMBAY returns ice thickness and the position of the grounding line. The ocean model uses a pre-computed mesh to allow for an easy adjustment of the model domain to a varying cavity geometry. RAnGO simulations with a 20th-century climate forcing yield realistic basal melt rates and a quasi-stable grounding line position close to the presently observed state. In a centennial-scale warm-water-inflow scenario, the model suggests a substantial thinning of the ice shelf and a local retreat of the grounding line. The potentially negative feedback from ice-shelf thinning through a rising in-situ freezing temperature is more than outweighed by the increasing water column thickness in the deepest parts of the cavity. Compared to a control simulation with fixed ice-shelf geometry, the coupled model thus yields a slightly stronger increase of ice-shelf basal melt rates.

Download Full-text

Viscosity and elasticity: a model intercomparison of ice-shelf bending in an Antarctic grounding zone

Journal of Glaciology ◽

10.1017/jog.2017.15 ◽

2017 ◽

Vol 63 (240) ◽

pp. 573-580 ◽

Cited By ~ 7

Author(s):

CHRISTIAN T. WILD ◽

OLIVER J. MARSH ◽

WOLFGANG RACK

Keyword(s):

Young’S Modulus ◽

Ocean Circulation ◽

Young's Modulus ◽

Viscoelastic Model ◽

Ice Sheet ◽

Global Ocean ◽

Elastic Model ◽

Ice Shelf ◽

Bending Stresses ◽

Global Ocean Circulation

ABSTRACTGrounding zones are vital to ice-sheet mass balance and its coupling to the global ocean circulation. Processes here determine the mass discharge from the grounded ice sheet, to the floating ice shelves. The response of this transition zone to tidal forcing has been described by both elastic and viscoelastic models. Here we examine the validity of these models for grounding zone flexure over tidal timescales using field data from the Southern McMurdo Ice Shelf (78° 15′S, 167° 7′E). Observations of tidal movement were carried out by simultaneous tiltmeter and GPS measurements along a profile across the grounding zone. Finite-element simulations covering a 64 d period reveal that the viscoelastic model fits best the observations using a Young's modulus of 1.6 GPa and a viscosity of 1013.7 Pa s (≈ 50.1 TPa s). We conclude that the elastic model is only well-constrained for tidal displacements >35% of the spring-tidal amplitude using a Young's modulus of 1.62 ± 0.69 GPa, but that a viscoelastic model is necessary to adequately capture tidal bending at amplitudes below this threshold. In grounding zones where bending stresses are greater than at the Southern McMurdo Ice Shelf or ice viscosity is lower, the threshold would be even higher.

Download Full-text

Suppression of marine ice sheet instability

Journal of Fluid Mechanics ◽

10.1017/jfm.2018.742 ◽

2018 ◽

Vol 857 ◽

pp. 648-680 ◽

Cited By ~ 3

Author(s):

Samuel S. Pegler

Keyword(s):

Steady State ◽

Rapid Assessment ◽

Ice Sheets ◽

Ice Sheet ◽

Stable Steady State ◽

Ice Shelf ◽

West Antarctic Ice Sheet ◽

West Antarctic ◽

Grounding Line ◽

The Stability

A long-standing open question in glaciology concerns the propensity for ice sheets that lie predominantly submerged in the ocean (marine ice sheets) to destabilise under buoyancy. This paper addresses the processes by which a buoyancy-driven mechanism for the retreat and ultimate collapse of such ice sheets – the marine ice sheet instability – is suppressed by lateral stresses acting on its floating component (the ice shelf). The key results are to demonstrate the transition between a mode of stable (easily reversible) retreat along a stable steady-state branch created by ice-shelf buttressing to tipped (almost irreversible) retreat across a critical parametric threshold. The conditions for triggering tipped retreat can be controlled by the calving position and other properties of the ice-shelf profile and can be largely independent of basal stress, in contrast to principles established from studies of unbuttressed grounding-line dynamics. The stability and recovery conditions introduced by lateral stresses are analysed by developing a method of constructing grounding-line stability (bifurcation) diagrams, which provide a rapid assessment of the steady-state positions, their natures and the conditions for secondary grounding, giving clear visualisations of global stabilisation conditions. A further result is to reveal the possibility of a third structural component of a marine ice sheet that lies intermediate to the fully grounded and floating components. The region forms an extended grounding area in which the ice sheet lies very close to flotation, and there is no clearly distinguished grounding line. The formation of this region generates an upsurge in buttressing that provides the most feasible mechanism for reversal of a tipped grounding line. The results of this paper provide conceptual insight into the phenomena controlling the stability of the West Antarctic Ice Sheet, the collapse of which has the potential to dominate future contributions to global sea-level rise.

Download Full-text

Marine Ice Sheet Instability and Ice Shelf Buttressing Influenced Deglaciation of the Minch Ice Stream, Northwest Scotland

10.5194/tc-2018-116 ◽

2018 ◽

Cited By ~ 1

Author(s):

Niall Gandy ◽

Lauren J. Gregoire ◽

Jeremy C. Ely ◽

Christopher D. Clark ◽

David M. Hodgson ◽

...

Keyword(s):

Ice Sheets ◽

Ice Sheet ◽

Last Deglaciation ◽

Ice Shelf ◽

Dynamical Processes ◽

Ice Dynamics ◽

Surface Mass ◽

Ice Shelves ◽

Ice Stream ◽

The Last Deglaciation

Abstract. Uncertainties in future sea level projections are dominated by our limited understanding of the dynamical processes that control instabilities of marine ice sheets. A valuable case to examine these processes is the last deglaciation of the British-Irish Ice Sheet. The Minch Ice Stream, which drained a large proportion of ice from the northwest sector of the British-Irish Ice Sheet during the last deglaciation, is well constrained, with abundant empirical data which could be used to inform, validate and analyse numerical ice sheet simulations. We use BISICLES, a higher-order ice sheet model, to examine the dynamical processes that controlled the retreat of the Minch Ice Stream. We simulate retreat from the shelf edge under constant "warm" surface mass balance and subshelf melt, to isolate the role of internal ice dynamics from external forcings. The model simulates a slowdown of retreat as the ice stream becomes laterally confined at a "pinning-point" between mainland Scotland and the Isle of Lewis. At this stage, the presence of ice shelves became a major control on deglaciation, providing buttressing to upstream ice. Subsequently, the presence of a reverse slope inside the Minch Strait produces an acceleration in retreat, leading to a "collapsed" state, even when the climate returns to the initial "cold" conditions. Our simulations demonstrate the importance of the Marine Ice Sheet Instability and ice shelf buttressing during the deglaciation of parts of the British-Irish Ice Sheet. Thus, geological data could be used to constrain these processes in ice sheet models used for projecting the future of our contemporary ice sheets.

Download Full-text