SAM: Simplified Analytical Model of Dyke Pathways in Three Dimensions

Cite as:

Mantiloni, Lorenzo; Rivalta, Eleonora; Davis, Timothy (2023): SAM: Simplified Analytical Model of Dyke Pathways in Three Dimensions. GFZ Data Services. https://doi.org/10.5880/GFZ.2.1.2023.001

Status

I N R E V I E W : Mantiloni, Lorenzo; Rivalta, Eleonora; Davis, Timothy (2023): SAM: Simplified Analytical Model of Dyke Pathways in Three Dimensions. GFZ Data Services. https://doi.org/10.5880/GFZ.2.1.2023.001

Abstract

SAM ("Simplified Analytical Model") is a MatLab-based software that allows for fast and flexible simulations of three-dimensional dyke pathways in an elastic medium. The model was first introduced in "Mechanical modeling of pre-eruptive magma propagation scenarios at calderas" (Mantiloni, L. et al. 2023). In SAM, dykes are modelled as penny-shaped cracks of fixed radius, opening against the local direction of the least-compressive principal stress. The direction of propagation is determined by the gradient of the external stress normal to the crack's plane and the buoyancy force of the magma filling the dyke, calculated at a set of observation points along the crack's tipline. The model can also include a uniform internal pressure within the dyke and compute the stress intensity factor along the crack's tipline, comparing it to the fracture toughness of the host rock to determine if the dyke will advance. SAM needs a model for the stress field of the host rock as input, as well as magma and rock densities, rock elastic properties, the dyke's radius and the number of observation points. The model may be applied to simulate dyke pathways in realistic volcanic settings with different stress sources, and can perform large numbers of simulations in little time. The model does not, however, account for any viscous flow of magma within the dyke, nor the velocity of dyke propagation. Dykes cannot change shape or area during the propagation, and are always bound to be oriented normally to the local least-compressive principal stress axis. This repository also includes data and parameters of the synthetic scenarios discussed in "Mechanical modeling of pre-eruptive magma propagation scenarios at calderas".

Authors

Mantiloni, Lorenzo;GFZ German Research Centre for Geosciences, Potsdam, Germany;Potsdam University, Institute for Geosciences, Potsdam, Germany
Rivalta, Eleonora;GFZ German Research Centre for Geosciences, Potsdam, Germany;Università di Bologna, Dipartimento di Fisica e Astronomia, Bologna, Italy
Davis, Timothy;University of Oxford, Department of Earth Sciences, Oxford, UK

Contact

Mantiloni, Lorenzo (PhD Candidate) ; GFZ German Research Centre for Geosciences, Potsdam, Germany; Potsdam University, Institute for Geosciences, Potsdam, Germany;