Ardnamurchan 3D cone-sheet architecture explained by a single elongate magma chamber

Magma transport in sheet intrusions of the Alnö carbonatite complex, central Sweden

Magma transport through the Earth’s crust occurs dominantly via sheet intrusions, such as dykes and cone-sheets, and is fundamental to crustal evolution, volcanic eruptions and geochemical element cycling. However, reliable methods to reconstruct flow direction in solidified sheet intrusions have proved elusive. Anisotropy of magnetic susceptibility (AMS) in magmatic sheets is often interpreted as primary magma flow, but magnetic fabrics can be modified by post-emplacement processes, making interpretation of AMS data ambiguous. Here we present AMS data from cone-sheets in the Alnö carbonatite complex, central Sweden. We discuss six scenarios of syn- and post-emplacement processes that can modify AMS fabrics and offer a conceptual framework for systematic interpretation of magma movements in sheet intrusions. The AMS fabrics in the Alnö cone-sheets are dominantly oblate with magnetic foliations parallel to sheet orientations. These fabrics may result from primary lateral flow or from sheet closure at the terminal stage of magma transport. As the cone-sheets are discontinuous along their strike direction, sheet closure is the most probable process to explain the observed AMS fabrics. We argue that these fabrics may be common to cone-sheets and an integrated geology, petrology and AMS approach can be used to distinguish them from primary flow fabrics.

Bimodal magmatism produced by progressively inhibited crustal assimilation

The origin of bimodal (mafic-felsic) rock suites is a fundamental question in volcanology. Here we use major and trace elements, high-resolution Sr, Nd and Pb isotope analyses, experimental petrology and thermodynamic modelling to investigate bimodal magmatism at the iconic Carlingford Igneous Centre, Ireland. We show that early microgranites are the result of extensive assimilation of trace element-enriched partial melts of local metasiltstones into mafic parent magmas. Melting experiments reveal the crust is very fusible, but thermodynamic modelling indicates repeated heating events rapidly lower its melt-production capacity. Granite generation ceased once enriched partial melts could no longer form and subsequent magmatism incorporated less fertile restite compositions only, producing mafic intrusions and a pronounced compositional gap. Considering the frequency of bimodal magma suites in the North Atlantic Igneous Province, and the ubiquity of suitable crustal compositions, we propose 'progressively inhibited crustal assimilation' (PICA) as a major cause of bimodality in continental volcanism.

Weibull-distributed dyke thickness reflects probabilistic character of host-rock strength

Magmatic sheet intrusions (dykes) constitute the main form of magma transport in the Earth’s crust. The size distribution of dykes is a crucial parameter that controls volcanic surface deformation and eruption rates and is required to realistically model volcano deformation for eruption forecasting. Here we present statistical analyses of 3,676 dyke thickness measurements from different tectonic settings and show that dyke thickness consistently follows the Weibull distribution. Known from materials science, power law-distributed flaws in brittle materials lead to Weibull-distributed failure stress. We therefore propose a dynamic model in which dyke thickness is determined by variable magma pressure that exploits differently sized host-rock weaknesses. The observed dyke thickness distributions are thus site-specific because rock strength, rather than magma viscosity and composition, exerts the dominant control on dyke emplacement. Fundamentally, the strength of geomaterials is scale-dependent and should be approximated by a probability distribution.

Understanding dyke thickness distributions is essential to quantify magma transport rates and improve eruption forecasting. Krumbholz et al. show that dyke thicknesses are Weibull-distributed and identify host-rock strength as the primary parameter that controls dyke emplacement.