Modelling magma-induced surface uplift and dynamic fracturing around laccoliths on the Moon, Mars, and Earth

Dome-shaped uplifted and fractured terrain observed at the surface of the Moon and Mars includes floor-fractured impact craters. Such deformation features are inferred to form by the emplacement and inflation of sill- and laccolith-shaped magma bodies in the shallowest 1-2 km of a planetary body’s rocky crust. Only the final surface deformation features can be observed from space, modelling helps to understand the emplacement dynamics and the deformation of the overlying rock. A mismatch exists, however, between the complex mechanical response of host rocks to magma-induced stresses observed on Earth in exposed volcanic plumbing systems and the linearly elastic deformation assumed by most of the often-used numerical models. We have implemented simulations of the inflation of a laccolith intrusion in a particle-based host medium in the two-dimensional (2D) Discrete Element Method (DEM). Our approach allows us to investigate magma-induced, highly discontinuous, deformation and dynamic fracturing and visualizes the localization of subsurface strain. We systematically varied a range of numerical model parameters that govern host rock strength (bond cohesion, bond tensile strength, bond elastic modulus), and specific gravity known for the Moon, Mars and Earth. For equal rock stiffness and amounts of intruded magma, our model results show that we can expect more vertical surface displacement on the Moon due to the lower gravity there compared to Mars, and Earth. Rock toughness and rock stiffness control the amount of fracturing more than gravity does. We also tested how host rock strengths in our 2D DEM model could be upscaled from intact strengths of rock samples collected at Earth analogue sites, or by implementing a digital crack network that simulates the highly fractured conditions of the intensively impacted Lunar and Martian crusts. Our results show that laccolith inflation in pre-cracked host rocks results in higher surface displacements and a higher amount of magma-induced cracking in broader fractured zones. We expect that our model results will induce a better understanding of the emplacement and architecture of shallow magmatic intrusions below magma-induced uplifted terrain and floor-fractured craters on the Moon and Mars.