A quantitative study of the mechanisms governing dike propagation, dike arrest and sill formation

Abstract

Dikes and sills are the moving building blocks of the plumbing system of volcanoes and play a fundamental role in the accretionary processes of the crust. They nucleate, propagate, halt, resume propagation, and sometimes change trajectory with drastic implications for the outcome of eruptions ( Sigmundsson et al., 2010). Their dynamics is still poorly understood, in particular when different external influencing factors are interacting. Here we apply a boundary element model to study dike and sill formation, propagation and arrest in different scenarios. We model dikes as finite batches of compressible fluid magma, propagating quasi-statically in an elastic medium, and calculate their trajectories by maximising the energy release of the magma-rock system. We consider dike propagation in presence of density layering, of density plus rigidity layering, of a weakly welded interface between layers, under the action of an external stress field (of tectonic or topographic origin). Our simulations predict sill formation in several situations: i) when a horizontal weak interface is met by a propagating dike; ii) when a sufficiently high compressive tectonic environment is experienced by the ascending dike and iii) in case a dike, starting below a volcanic edifice, propagates away from the topographic load with a low dip angle. We find that dikes halt and stack when they become negatively buoyant and when they propagate with low overpressure at their upper tip toward a topographic load. Neutral buoyancy by itself cannot induce dikes to turn into sills, as previously suggested.