Coulomb failure of Earth's brittle crust controls growth, emplacement and shapes of igneous sills, saucer-shaped sills and laccoliths.

Abstract

Tabular intrusions are common features in the Earth's brittle crust. They exhibit a broad variety of shapes, ranging from thin sheet intrusions (sills, saucer-shaped sills, cone sheets), to more massive intrusions (domed and punched laccoliths, stocks). Such a diversity of intrusion shapes reflects different emplacement mechanisms caused by contrasting host rock and magma rheologies. Most current models of tabular intrusion emplacement assume that the host rock behaves purely elastically, whereas numerous observations show that shear failure plays a major role. In this study, we investigate the effects of the host rock's Coulomb properties on magma emplacement by integrating (1) laboratory models using dry Coulomb granular model hosts of variable strength (cohesion) and (2) limit analysis numerical models. Our results show that both sheet and massive tabular intrusions initiate as a sill, which triggers shear failure of its overburden along an inclined shear damage zone at a critical sill radius, which depends on the emplacement depth and the overburden's cohesion. Two scenarios are then possible: (1) if the cohesion of the overburden is significant, opening of a planar fracture along the precursory weakened shear damage zones to accommodate magma flow, leads to the formation of inclined sheets, or (2) if the cohesion of the overburden is negligible, the sill inflates and lifts up the overburden, which is dissected by several faults that control the growth of a massive intrusion. Finally, we derive a theoretical scaling that predicts the thickness-to-radius aspect ratios of the laboratory sheet intrusions. This theoretical prediction shows how sheet intrusion morphologies are controlled by a mechanical equilibrium between the flowing viscous magma and Coulomb shear failure of the overburden. Our study suggests that the emplacement of sheet and massive tabular intrusions are parts of the same mechanical regime, in which the Coulomb behavior of the Earth's brittle crust plays an essential role.