EGU2020-16485, updated on 09 Jan 2024
EGU General Assembly 2020
© Author(s) 2024. This work is distributed under
the Creative Commons Attribution 4.0 License.

Towards a staggered-grid finite difference code for modelling magmatic systems

Nicolas Berlie1, Boris Kaus1, Anton Popov1, and Patrick Sanan2
Nicolas Berlie et al.
  • 1Johannes Gutenberg University Mainz, Institute of Geosciences, Geodynamics, Mainz, Germany (
  • 2ETH Zürich, Department of Earth Sciences, Institute of Geophysics, Zürich, Switzerland

The in-depth behaviour of magmatic systems is still poorly constrained due to their lack of accessibility and the difficulty of finding good analogue representations. However, progresses in the field of geological numerical modelling can allow to better understand and interpret those constraints. The MAGMA project aims on developing tools and software for studying a range of magmatic processes in the lithosphere. On the way to building a general framework able to model the behaviour of a fluid-solid chemically coupled magmatic system, we present here the current development of a mechanical staggered-grid finite difference 2D code using robust analytical linear and non-linear solvers via the PETSc infrastructure, able to run in parallel on highly performant computers. The mesh is assembled using the recently developed DMStag framework, which is part of PETSc. This code solves the Stokes equations for elasto-visco-plastic rheologies, including tensile plasticity essential to the development of dyke structures. Here, we will present the equations and implementations used, and show initial results and benchmarks.

How to cite: Berlie, N., Kaus, B., Popov, A., and Sanan, P.: Towards a staggered-grid finite difference code for modelling magmatic systems, EGU General Assembly 2020, Online, 4–8 May 2020, EGU2020-16485,, 2020.

This abstract will not be presented.