M2TRI_vep

M2TRI_vep is a 2-D geodynamic model applying the finite element method on unstructured triangular meshes to solve for the thermal and mechanical evolution of the lithosphere and the underlying mantle. Conservation equations of mass, momentum and energy are solved on a deforming Lagrangian mesh. M2TRI_vep features an incompressible, visco-elasto-plastic, non-Newtonian rheology and a free surface. The code includes modules for grain size reduction, serpentinization, parameterized hydrothermal cooling and the thermal effects of sill intrusions. It is written in MATLAB and through changes in geometry, initial and boundary conditions, and rheological properties the model can easily be adapted for other tectonic settings. A detailed model description can be found in: Glink, A., and Hasenclever, J. (2023). How Hydrothermal Cooling and Magmatic Sill Intrusions Control Flip-Flop Faulting at Ultraslow-Spreading Mid-Ocean Ridges. Geochemistry, Geophysics, Geosystems (accepted for publication