scholarly journals Coupling framework (1.0) for the ice sheet model PISM (1.1.1) and the ocean model MOM5 (5.1.0) via the ice-shelf cavity module PICO

2020 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Nicholas Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  
Keyword(s):  
Time Scales ◽  
General Circulation ◽  
Ice Sheet ◽  
Circulation Model ◽  
Ocean Model ◽  
Global Ocean ◽  
Energy Fluxes ◽  
Ice Shelf ◽  
Antarctic Ice Sheet ◽  
Wide Range

Abstract. The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degree, respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the entire Antarctic Ice Sheet and the Earth system over time spans on the order of centuries to millennia. In this study we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet model is calculated by PICO from modeled ocean temperatures and salinities at the depth of the continental shelf and, vice versa, the resulting mass and energy fluxes from the melting at the ice-ocean interface are transferred to the ocean model. Mass and energy fluxes are shown to be conserved to machine precision across the considered model domains. The implementation is computationally efficient as it introduces only minimal overhead. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions and time scales between the ice and ocean model in a generic way, and can thus be adopted to a wide range of model setups.

First results from coupling the Parallel Ice Sheet Model with the Modular Ocean Model via an Antarctic ice-shelf cavity module

2021 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  
Keyword(s):  
Time Scales ◽  
General Circulation ◽  
Ice Sheet ◽  
Circulation Model ◽  
Ocean Model ◽  
Global Ocean ◽  
Ice Shelf ◽  
Antarctic Ice Sheet ◽  
First Results ◽  
The Antarctic

<p>The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high resolution configurations, limiting these studies to individual glaciers or regions over short time scales of decades to a few centuries. To study global and long term interactions, we developed a framework to couple the dynamic ice sheet model PISM with the global ocean general circulation model MOM5 via the ice-shelf cavity module PICO. Since ice-shelf cavities are not resolved by MOM5, but parameterized with the box model PICO, the framework allows the ice sheet and ocean model to be run at resolution of 16 km and 3 degrees, respectively. We present first results from our coupled setup and discuss stability, feedbacks, and interactions of the Antarctic Ice Sheet and the global ocean system on millennial time scales.</p>

scholarly journals Coupling framework (1.0) for the PISM (1.1.4) ice sheet model and the MOM5 (5.1.0) ocean model via the PICO ice shelf cavity model in an Antarctic domain

2021 ◽  
Vol 14 (6) ◽  
pp. 3697-3714
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Nicholas Huiskamp ◽  
Stefan Petri ◽  
Torsten Albrecht ◽  
...  
Keyword(s):  
General Circulation ◽  
Ice Sheet ◽  
Ocean Model ◽  
Global Ocean ◽  
Cavity Model ◽  
Energy Fluxes ◽  
Ice Shelf ◽  
Antarctic Ice Sheet ◽  
Wide Range ◽  
The Antarctic

Abstract. The past and future evolution of the Antarctic Ice Sheet is largely controlled by interactions between the ocean and floating ice shelves. To investigate these interactions, coupled ocean and ice sheet model configurations are required. Previous modelling studies have mostly relied on high-resolution configurations, limiting these studies to individual glaciers or regions over short timescales of decades to a few centuries. We present a framework to couple the dynamic ice sheet model PISM (Parallel Ice Sheet Model) with the global ocean general circulation model MOM5 (Modular Ocean Model) via the ice shelf cavity model PICO (Potsdam Ice-shelf Cavity mOdel). As ice shelf cavities are not resolved by MOM5 but are parameterized with the PICO box model, the framework allows the ice sheet and ocean components to be run at resolutions of 16 km and 3∘ respectively. This approach makes the coupled configuration a useful tool for the analysis of interactions between the Antarctic Ice Sheet and the global ocean over time spans of the order of centuries to millennia. In this study, we describe the technical implementation of this coupling framework: sub-shelf melting in the ice sheet component is calculated by PICO from modelled ocean temperatures and salinities at the depth of the continental shelf, and, vice versa, the resulting mass and energy fluxes from melting at the ice–ocean interface are transferred to the ocean component. Mass and energy fluxes are shown to be conserved to machine precision across the considered component domains. The implementation is computationally efficient as it introduces only minimal overhead. Furthermore, the coupled model is evaluated in a 4000 year simulation under constant present-day climate forcing and is found to be stable with respect to the ocean and ice sheet spin-up states. The framework deals with heterogeneous spatial grid geometries, varying grid resolutions, and timescales between the ice and ocean component in a generic way; thus, it can be adopted to a wide range of model set-ups.

scholarly journals Competing influences of the ocean, atmosphere and solid earth on transient Miocene Antarctic ice sheet variability

2021 ◽  
Author(s):  
Lennert Bastiaan Stap ◽  
Constantijn J. Berends ◽  
Meike D. W. Scherrenberg ◽  
Roderik S. W. van de Wal ◽  
Edward G. W. Gasson
Keyword(s):  
General Circulation ◽  
Ice Sheet ◽  
Circulation Model ◽  
Climate Forcing ◽  
Ice Shelf ◽  
Antarctic Ice Sheet ◽  
Ice Volume ◽  
Bedrock Topography ◽  
Co2 Levels ◽  
Transient Simulations

Abstract. Benthic δ18O levels vary strongly during the warmer-than-modern early- and mid-Miocene (23 to 14 Myr ago), suggesting a dynamic Antarctic ice sheet (AIS). So far, however, realistic simulations of the Miocene AIS have been limited to equilibrium states under different CO2 levels and orbital settings. Earlier transient simulations lacked ice-sheet-atmosphere interactions, and used a present-day rather than Miocene Antarctic bedrock topography. Here, we quantify the effect of ice-sheet-atmosphere interactions, running IMAU-ICE using climate forcing from Miocene simulations by the general circulation model GENESIS. Utilising a recently developed matrix interpolation method enables us to interpolate the climate forcing based on CO2 levels (between 280 and 840 ppm) as well as ice sheet configurations (between no ice and a large ice sheet). We furthermore implement recent reconstructions of Miocene Antarctic bedrock topography. We find that the positive albedo-temperature feedback, partly compensated by the negative ice-volume-precipitation feedback, increases hysteresis in the relation between CO2 and ice volume (V). Together, these ice-sheet-atmosphere interactions decrease the amplitude of AIS variability caused by 40-kyr forcing CO2 cycles by 21 % in transient simulations. Thereby, they also diminish the contribution of AIS variability to benthic δ18O fluctuations. Furthermore, we show that under equal atmospheric and oceanic forcing, the amplitude of 40-kyr transient AIS variability becomes 10 % smaller during the early- and mid-Miocene, due to the evolving bedrock topography. Lastly, we quantify the influence of ice shelf formation around the Antarctic margins, by comparing simulations with Last Glacial Maximum (LGM) basal melt conditions, to ones in which ice shelf growth is prevented. Ice shelf formation increases hysteresis in the CO2-V relation, and amplifies 40-kyr AIS variability by 19 % using LGM basal melt rates, and by 5 % in our reference setting.

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

2017 ◽  
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.

scholarly journals Antarctic ice sheet fertilises the Southern Ocean

2014 ◽  
Vol 11 (10) ◽  
pp. 2635-2643 ◽  
Author(s):  
R. Death ◽  
J. L. Wadham ◽  
F. Monteiro ◽  
A. M. Le Brocq ◽  
M. Tranter ◽  
...  
Keyword(s):  
Southern Ocean ◽  
Ice Sheet ◽  
Ocean Model ◽  
Coastal Sediments ◽  
Global Ocean ◽  
Iron Cycling ◽  
Antarctic Ice Sheet ◽  
Significant Control ◽  
Subglacial Meltwater ◽  
The Impact

Abstract. Southern Ocean (SO) marine primary productivity (PP) is strongly influenced by the availability of iron in surface waters, which is thought to exert a significant control upon atmospheric CO2 concentrations on glacial/interglacial timescales. The zone bordering the Antarctic Ice Sheet exhibits high PP and seasonal plankton blooms in response to light and variations in iron availability. The sources of iron stimulating elevated SO PP are in debate. Established contributors include dust, coastal sediments/upwelling, icebergs and sea ice. Subglacial meltwater exported at the ice margin is a more recent suggestion, arising from intense iron cycling beneath the ice sheet. Icebergs and subglacial meltwater may supply a large amount of bioavailable iron to the SO, estimated in this study at 0.07–0.2 Tg yr−1. Here we apply the MIT global ocean model (Follows et al., 2007) to determine the potential impact of this level of iron export from the ice sheet upon SO PP. The export of iron from the ice sheet raises modelled SO PP by up to 40%, and provides one plausible explanation for seasonally very high in situ measurements of PP in the near-coastal zone. The impact on SO PP is greatest in coastal regions, which are also areas of high measured marine PP. These results suggest that the export of Antarctic runoff and icebergs may have an important impact on SO PP and should be included in future biogeochemical modelling.

scholarly journals Antarctic Ice Sheet fertilises the Southern Ocean

2013 ◽  
Vol 10 (7) ◽  
pp. 12551-12570 ◽  
Author(s):  
R. Death ◽  
J. L. Wadham ◽  
F. Monteiro ◽  
A. M. Le Brocq ◽  
M. Tranter ◽  
...  
Keyword(s):  
Southern Ocean ◽  
Ice Sheet ◽  
Ocean Model ◽  
Coastal Sediments ◽  
Global Ocean ◽  
Iron Cycling ◽  
Antarctic Ice Sheet ◽  
Significant Control ◽  
Subglacial Meltwater ◽  
The Impact

Abstract. Southern Ocean (SO) marine primary productivity (PP) is strongly influenced by the availability of iron in surface waters, which is thought to exert a significant control upon atmospheric CO2 concentrations on glacial/interglacial timescales. The zone bordering the Antarctic Ice Sheet exhibits high PP and seasonal plankton blooms in response to light and variations in iron availability. The sources of iron stimulating elevated SO PP are in debate. Established contributors include dust, coastal sediments/upwelling, icebergs and sea ice. Subglacial meltwater exported at the ice margin is a more recent suggestion, arising from intense iron cycling beneath the ice sheet. Icebergs and subglacial meltwater may supply a large amount of bioavailable iron to the SO, estimated in this study at 0.07–1.0 Tg yr−1. Here we apply the MIT global ocean model (Follows et al., 2007) to determine the potential impact of this level of iron export from the ice sheet upon SO PP. The export of iron from the ice sheet raises modelled SO PP by up to 40%, and provides one plausible explanation for very high seasonally observed PP in the near-coastal zone. The impact on SO PP is greatest in coastal regions, which are also areas of high observed marine PP. These results suggest that the export of Antarctic runoff and icebergs may have an important impact on SO PP and should be included in future biogeochemical modelling.

scholarly journals Response to Filchner–Ronne Ice Shelf cavity warming in a coupled ocean–ice sheet model – Part 1: The ocean perspective

Ocean Science ◽  
2017 ◽  
Vol 13 (5) ◽  
pp. 765-776 ◽  
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. The Regional Antarctic ice 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 in ice-shelf basal melt rates.

Coupling the Parallel Ice Sheet Model with the Modular Ocean Model via an Antarctic ice-shelf cavity module

2020 ◽  
Author(s):  
Moritz Kreuzer ◽  
Ronja Reese ◽  
Willem Huiskamp ◽  
Stefan Petri ◽  
Ricarda Winkelmann
Keyword(s):  
Time Scales ◽  
Ice Sheet ◽  
Mass Conservation ◽  
Ocean Model ◽  
Computationally Efficient ◽  
Ice Shelf ◽  
Current Mass ◽  
Modular Ocean Model ◽  
Interactive Modelling ◽  
The Antarctic

<p>Ocean-ice shelf interactions are the main drivers for the current mass loss from the Antarctic Ice Sheet. Recent studies have shown that increased continental meltwater input can enhance discharge through ice-ocean feedbacks. This raises the need for interactive modelling between ocean and ice-sheet systems to assess the consequences of additional freshwater input on the Antarctic region and beyond. While high-resolution simulations (1/4 degree or greater) can resolve detailed interactions between the ocean and ice shelf, the computational costs make them applicable only for regional studies or decadal to centennial time scales. In this study we present a framework for coupling a coarse resolution ocean model (MOM) to the Parallel Ice Sheet Model (PISM) via the Potsdam Ice-shelf Cavity mOdel (PICO). The intermediate model PICO approximates the overturning circulation in ice shelf cavities and includes ice-ocean boundary layer physics. We present this offline coupling approach and discuss the fluxes exchanged between the distinct model runs as well as energy and mass conservation. Using this flexible and computationally efficient framework, feedbacks between the ice and ocean can be analysed on a global spatial scale and paleoclimate time-scales.</p><p> </p>

scholarly journals The GPU version of LASG/IAP Climate System Ocean Model version 3 (LICOM3) under the heterogeneous-compute interface for portability (HIP) framework and its large-scale application

2021 ◽  
Vol 14 (5) ◽  
pp. 2781-2799
Author(s):  
Pengfei Wang ◽  
Jinrong Jiang ◽  
Pengfei Lin ◽  
Mengrong Ding ◽  
Junlin Wei ◽  
...  
Keyword(s):  
General Circulation ◽  
Large Scale ◽  
Vertical Direction ◽  
Circulation Model ◽  
Ocean General Circulation Model ◽  
Ocean Model ◽  
Climate System ◽  
Global Ocean ◽  
Processing Unit ◽  
Model Version

Abstract. A high-resolution (1/20∘) global ocean general circulation model with graphics processing unit (GPU) code implementations is developed based on the LASG/IAP Climate System Ocean Model version 3 (LICOM3) under a heterogeneous-compute interface for portability (HIP) framework. The dynamic core and physics package of LICOM3 are both ported to the GPU, and three-dimensional parallelization (also partitioned in the vertical direction) is applied. The HIP version of LICOM3 (LICOM3-HIP) is 42 times faster than the same number of CPU cores when 384 AMD GPUs and CPU cores are used. LICOM3-HIP has excellent scalability; it can still obtain a speedup of more than 4 on 9216 GPUs compared to 384 GPUs. In this phase, we successfully performed a test of 1/20∘ LICOM3-HIP using 6550 nodes and 26 200 GPUs, and on a large scale, the model's speed was increased to approximately 2.72 simulated years per day (SYPD). By putting almost all the computation processes inside GPUs, the time cost of data transfer between CPUs and GPUs was reduced, resulting in high performance. Simultaneously, a 14-year spin-up integration following phase 2 of the Ocean Model Intercomparison Project (OMIP-2) protocol of surface forcing was performed, and preliminary results were evaluated. We found that the model results had little difference from the CPU version. Further comparison with observations and lower-resolution LICOM3 results suggests that the 1/20∘ LICOM3-HIP can reproduce the observations and produce many smaller-scale activities, such as submesoscale eddies and frontal-scale structures.

Sign in / Sign up

Export Citation Format

Share Document