EGU General Assembly 2020
© Author(s) 2020. This work is distributed under
the Creative Commons Attribution 4.0 License.

Towards tackling ice-sheet ocean interaction with Finite Element Methods

Stefano Ottolenghi1 and Josefin Ahlkrona2
Stefano Ottolenghi and Josefin Ahlkrona
  • 1Stockholms universitet, Matematiska institutionen, Sweden (
  • 2Stockholms universitet, Matematiska institutionen, Sweden (

Ice sheet-ocean interaction is important to properly understand phenomena such as ice sheet melting and ocean circulation. While the long term goal of this project is to fully couple the ice and ocean in one single numerical framework, we here start by modelling the ocean. We use the full non-hydrostatic equations in order to accurately model the complex ocean dynamics near the ice sheets. As numerical method, we employ finite element methods due to their capability of representing a complex fjord geometry and locally refining the mesh in the areas which require more careful handling, and its strong mathematical foundation. This will allow to overcome classical problems such as representing a moving ice shelf in a discretized setting. We here present an example of modeled fjord circulation obtained simulating the model with the FEniCS computing platform.

How to cite: Ottolenghi, S. and Ahlkrona, J.: Towards tackling ice-sheet ocean interaction with Finite Element Methods, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-11279,, 2020


Display file