May 2026 LIP of the Month

Large volcanic eruptions are mostly sourced above mobile basal mantle structures

Annalise Cucchiaro1, Nicolas Flament1, Maëlis Arnould2, Noel Cressie3

1 Environmental Futures, School of Earth, Atmospheric and Life Sciences, University of Wollongong, Wollongong, NSW 2522, Australia’. Email: A. Cucchiaro: ac576@uowmail.edu.au

2 Université Claude Bernard Lyon 1, École Normale Supérieure de Lyon, CNRS, Université Jean Monnet, LGL-TPE, UMR 5276, Villeurbanne, France

3 NIASRA, School of Mathematics and Physics, University of Wollongong, Wollongong, NSW 2522, Australia

 

Extracted and modified from the publication:

Cucchiaro, A., Flament, N., Arnould, M. and N. Cressie. Large volcanic eruptions are mostly sourced above mobile basal mantle structures. Nature Communications Earth and Environment6, 538 (2025). https://doi.org/10.1038/s43247-025-02482-z

 

Basal mantle structures, mantle plumes and LIPs

Large Igneous Provinces (LIPs) are the surface expression of mantle plumes (Bryan & Ernst, 2008; Ernst, 2014), which transport hot mantle material from the deep Earth to the surface. Seismic imaging of the lower mantle reveals two large continent-sized regions characterised by low velocity (Garnero & McNamara, 2008) called Large Low Velocity Provinces (LLVPs). Mantle plumes are thought to originate from these structures, but it is debated whether mantle plumes are generated preferentially in a plume generation zone along the edges of these basal mantle structures (Doubrovine et al., 2016; Torsvik et al., 2010), or from their interiors (Austermann et al., 2014; Davies, 2015), or from both. Some studies also propose that basal mantle structures have been stationary for hundreds of millions of years (Torsvik et al., 2010). However, mantle flow models that are either fully-dynamic or forced by plate tectonic histories both suggest that basal mantle structures have been mobile through time, deformed by sinking oceanic lithosphere over tens of millions of years (McNamara & Zhong, 2005). Here, we refer to the modelled equivalent of LLVPs as Big LOwer Mantle Basal Structures, or BLOBS.

Associations in space and time have been proposed to exist between 1/ LIPs, kimberlites and fixed LLVP edges (Torsvik et al., 2010), and 2/ LIPs, kimberlites and mobile BLOBS (Flament et al., 2022); however, previous studies have not explicitly considered mantle plumes; they are only implied as the mechanism in which LIPs are emplaced.

Here, we measure the distance between LIPs and plume conduits predicted by paleogeographically constrained mantle flow models at 1,040 km depth (Figure 1). We also measure the distance between LIPs and a suite of seismically-imaged LLVP (tomographic model cases T1–T4, from previous work (Auer et al., 2014; French & Romanowicz, 2014; Ritsema et al., 2011; Simmons et al., 2010), and between LIPs and BLOBS edges predicted by mantle flow model cases C1–C6 with varying basal layer density (between 0–2% greater than ambient mantle). We consider three LIP databases developed to assess distinct aspects of the volcanic record: D16 (Doubrovine et al., 2016), which consists of 26 LIPs associated with plume ‘heads’ representative of an updated selection proposed to be of deep origin; EY17 (Ernst & Youbi, 2017), which consists of 74 mafic and silicic LIPs associated with plume ‘heads’; and J18 (Johansson et al., 2018), which includes 167 volcanic products associated with both plume ‘heads’ and ‘tails’(Coffin et al., 2006).


Figure 1.Schematic of spatial relationships between basal mantle structures and volcanic eruptions with implicit mantle plumes as in previous work (Austermann et al., 2014; Davies et al., 2015; Doubrovine et al., 2016; Flament et al., 2022; Torsvik et al., 2010) (left) and with explicit mantle plumes as in this study (right), where blue arrows represent the distance between LIPs and respective features. The schematic is not drawn to scale.

 

Success metrics for LLSVPs and predicted BLOBS and plume conduits

We summarise angular distances (θ (°)) between large volcanic eruptions from 300 Ma and mantle structures (LLSVPs, BLOBS, and plumes) (Fig.2a) as cumulative distribution functions which we refer to as ‘sample distributions’ (Fig.2b). We repeat this for 1,000 random sets of LIPs with the same temporal resolution as the sample LIP database (Fig.2b). We then compute the p-value, denoted P, for a Monte Carlo significance test (Hope, 1968) based on the Kolmogorov-Smirnov statistic (Massey, 1951) (Fig.2b). This test allows us to compare the sample distance distribution with the random distance distribution. The null hypothesis of this test is that the distances between LIPs and mantle structures are the result of a sample of ‘uniform random locations’ where P < 0.05 indicates a departure from the null hypothesis (Casella & Berger, 2002). We report the mean of all angular distances within the sample distribution (θs) and the grand mean of all angular distances obtained from the 1,000 samples of uniform random locations (θr).


Figure 2.An example of the distance test between J18 LIPs, random J18 LIPs, and plume conduit locations at 60 Ma(a) angular distance to plume conduits (detected at 1,040km, shown as centroids - white disks) at 60 Ma for mantle flow model case C2, with LIPs in database J18 reconstructed to location of eruption (as black disks) and one set of random points with the same temporal resolution as sample points (as grey disks). (b) cumulative distribution function of sample LIPs (black curve) and 1,000 sets of random LIPs (grey curves) distance to plume conduits, with p-value (P) from a Monte Carlo significance test based on the Kolmogorov-Smirnov statistic.

 

Relationships between LIPs and model plumes

LIPs in the database J18 that contains both plume head and plume tail products are closer to plume conduits than random points for all mantle flow model cases (C1–C6; Fig. 3a,b). Through Monte Carlo significance testing we find a statistical-dependence relationship between modelled plume conduits and the J18 database for all model cases (Fig. 3c). We attribute this to the inclusion of plume ‘tail’ products in J18, which results in more numerous products of hotspot volcanism and clustering due to the motion of plates above a plume conduit. Coincidentally, plume tail products in database J18 are also more closely comparable to the plume conduits detected in our models. Conversely, databases EY17 and D16 containing silicic LIPs and selected plume ‘head’ products respectively, resulted in no statistical relationships for any model case (Fig.3c).


Figure 3. Summary of distances between volcanic eruptions and plume conduits at 1,040 km depth for all mantle flow model cases. (a), Average minimum angular distance of eruption locations to plume-conduit locations. (b), Grand mean of angular distances obtained from 1,000 samples of uniform random eruption locations to plume-conduit locations (J18, green discs; EY17 orange discs; D16 pink discs). (c), p-values for the Monte Carlo significance tests, with grey shading indicating departure from the null hypothesis (symmetrical logarithmic scale). Tick labels on the right of panel c indicate the BLOBS intrinsic density anomaly (δρ) in a given flow model.

 

Relationships between LIPs and: LLVPs and BLOBS

To assess the relationship between LIPs and basal structure edges, we consider LIPs above the interiors and exteriors of LLVPs and of BLOBS separately. For LIPs above the exterior of the LLVPs, we find a statistically significant relationship between LIPs in one database and LLVP edges in only one tomographic model (T1). There is no statistically significant relationship for any other permutation of LIPs in any database (D16, EY17, J18), and LLVP edges in any other tomographic model (T1–T4), either above the interior or the exterior of LLVPs. We attribute the success of tomographic model T1 to the overall smaller area of LLVPs, and to the presence of numerous small-scale basal mantle structures that are absent in other tomography models. For moving BLOBS, there is no statistical relationship between LIPs in any database and BLOBS interiors. However, there is a statistically significant relationship between J18 LIPs and BLOBS exteriors when BLOBS are between 1% and 1.6% denser than ambient mantle (Fig. 4c). This result supports a view in which plumes are not limited to a plume generation zone at the edge of fixed LLVPs (Doubrovine et al., 2016; Torsvik et al., 2010) but are instead anywhere above the interior of mobile BLOBS, or within ~5º to the outside of their edges due to plume tilting.


Figure 4. Angular distances and statistical significance for volcanic-eruption locations reconstructed above LLSVP and BLOBS exteriors for all model cases. (a), Mean minimum angular distance between large volcanic eruptions of three databases and fixed LLSVPs and mobile BLOBS edges. (b), Grand mean of angular distances obtained from 1,000 samples of uniform random eruption locations with the same temporal distribution as each respective large volcanic eruption database to fixed LLSVP edges or mobile BLOBS edges from 300 Ma. (c), p-values (P) for the Monte Carlo significance tests, with grey shading indicating departure from the null hypothesis (symmetrical logarithmic scale).

References:

Auer, L., Boschi, L., Becker, T. W., Nissen-Meyer, T., & Giardini, D. (2014). Savani: A variable resolution whole-mantle model of anisotropic shear velocity variations based on multiple data sets. Journal of Geophysical Research: Solid Earth, 119(4), 3006-3034. https://doi.org/10.1002/2013jb010773

Austermann, J., Kaye, B. T., Mitrovica, J. X., & Huybers, P. (2014). A statistical analysis of the correlation between large igneous provinces and lower mantle seismic structure. Geophysical Journal International, 197(1), 1-9. https://doi.org/10.1093/gji/ggt500

Bryan, S. E., & Ernst, R. E. (2008). Revised definition of Large Igneous Provinces (LIPs). Earth-Science Reviews, 86(1-4), 175-202. https://doi.org/10.1016/j.earscirev.2007.08.008

Casella, G., & Berger, R. L. (2002). Statistical inference. Duxbury Press: Belmont, CA.

Coffin, M. F., Duncan, R. A., Eldholm, O., Fitton, J. G., Frey, F. A., Larsen, H. C., Mahoney, J. J., Saunders, A. D., Schlich, R., & Wallace, P. J. (2006). Large igneous provinces and scientific ocean drilling: Status quo and a look ahead. Oceanography, 19(4), 150-160.

Davies, C. J. (2015). Cooling history of Earth’s core with high thermal conductivity. Physics of the Earth and Planetary Interiors, 247, 65-79.

Doubrovine, P. V., Steinberger, B., & Torsvik, T. H. (2016). A failure to reject: Testing the correlation between large igneous provinces and deep mantle structures with EDF statistics. Geochemistry, Geophysics, Geosystems, 17(3), 1130-1163. https://doi.org/10.1002/2015gc006044

Ernst, R. E. (2014). Large Igneous Provinces. Cambridge University Press, 653 p.

Ernst, R. E., & Youbi, N. (2017). How Large Igneous Provinces affect global climate, sometimes cause mass extinctions, and represent natural markers in the geological record. Palaeogeography, Palaeoclimatology, Palaeoecology, 478, 30-52. https://doi.org/10.1016/j.palaeo.2017.03.014

Flament, N., Bodur, O. F., Williams, S. E., & Merdith, A. S. (2022). Assembly of the basal mantle structure beneath Africa. Nature, 603(7903), 846-851. https://doi.org/10.1038/s41586-022-04538-y

French, S. W., & Romanowicz, B. A. (2014). Whole-mantle radially anisotropic shear velocity structure from spectral-element waveform tomography. Geophysical Journal International, 199(3), 1303-1327. https://doi.org/10.1093/gji/ggu334

Garnero, E. J., & McNamara, A. K. (2008). Structure and dynamics of Earth's lower mantle. Science, 320(5876), 626-628.

Hope, A. C. A. (1968). A simplified Monte Carlo significance test procedure. Journal of the Royal Statistical Society: Series B (Methodological), 30(3), 582-598.

Johansson, L., Zahirovic, S., & Müller, R. D. (2018). The interplay between the eruption and weathering of Large Igneous Provinces and the deep-time carbon cycle. Geophysical Research Letters, 45(11), 5380-5389. https://doi.org/10.1029/2017gl076691

Massey, F. J. (1951). The Kolmogorov-Smirnov test for goodness of fit. Journal of the American Statistical Association, 46(253), 68-78.

McNamara, A. K., & Zhong, S. (2005). Thermochemical structures beneath Africa and the Pacific Ocean. Nature, 437(7062), 1136-1139. https://doi.org/10.1038/nature04066

Ritsema, J., Deuss, A., van Heijst, H. J., & Woodhouse, J. H. (2011). S40RTS: a degree-40 shear-velocity model for the mantle from new Rayleigh wave dispersion, teleseismic traveltime and normal-mode splitting function measurements. Geophysical Journal International, 184(3), 1223-1236. https://doi.org/10.1111/j.1365-246X.2010.04884.x

Simmons, N. A., Forte, A. M., Boschi, L., & Grand, S. P. (2010). GyPSuM: A joint tomographic model of mantle density and seismic wave speeds. Journal of Geophysical Research: Solid Earth, 115(B12). https://doi.org/10.1029/2010jb007631

Torsvik, T. H., Burke, K., Steinberger, B., Webb, S. J., & Ashwal, L. D. (2010). Diamonds sampled by plumes from the core-mantle boundary. Nature, 466(7304), 352-355. https://doi.org/10.1038/nature09216