1. Create new computer models to simulate two-phase flow through a deforming thermo-mechanical visco-elasto-plastic lithosphere. Both mode-1 and mode-2 plasticity models will be taken into account. We plan to utilize PETSc for this purpose together with the DMSTAG framework.

2. Develop new approaches to couple state-of-the-art thermodynamic melting models with geodynamic models to simulate the evolving chemistry of crystallizing melts on its way from the mantle or lower crust to shallower levels in the lithosphere. Recently, hydrous mafic melting models have been developed that can, in principle, simulate the chemical differentiation of a hydrous mafic melt to more felsic compostions such as granite or granodiorite. Yet, coupling such thermodynamic models with geodynamic models remains a challenging task as the thermodynamic computations are too time-consuming to be performed on the fly. We will work on different approaches to overcome this.

3. Perform geodynamic inverse models of active magmatic systems, for which geophysics data is available. This allows interpreting the data in a mechanically consistent manner.

4. Perform systematic numerical simulations to better understand the physics of these models and compare it with observations.

As of now, work concentrates on the first three tasks. The software that will be developed in MAGMA will be released as open-source software, after we publish the methodology.