# Paroxysmal eruptions tracked by variations of helium isotopes: inferences from Piton de la Fournaise (La Réunion island)

G. Boudoire, A. Rizzo, I. Arienzo, A. Di Muro

► To cite this version:

G. Boudoire, A. Rizzo, I. Arienzo, A. Di Muro. Paroxysmal eruptions tracked by variations of helium isotopes: inferences from Piton de la Fournaise (La Réunion island). Scientific Reports, 2020, 10 (1), 10.1038/s41598-020-66260-x. http://dx.doi.org/10.1038/s41598-020-66260-x

HAL Id: hal-02901340

https://uca.hal.science/hal-02901340

Submitted on 25 Feb 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Paroxysmal eruptions tracked by variations of helium isotopes: inferences from Piton de la Fournaise (La Réunion island)

G. Boudoire<fnref n="1" />,²✉, A. L. Rizzo²,³, I. Arienzo⁴ & A. Di Muro⁵,⁶

Helium (He) with its isotopes (³He, ⁴He) is a key tracer enabling the Earth's mantle and dynamics to be characterized. Enrichment in primordial helium (³He) has been detected in volcanic gases of numerous magmatic systems in different geodynamic settings. Despite past use to monitor volcanotectonic unrest, temporal ³He/⁴He variability in volcanic emissions is still poorly constrained. Here, we investigate noble gas chemistry of Piton de la Fournaise hotspot volcano, where temporal fluctuations of ³He/⁴He in response to the eruptive activity have never been studied. We compare the ³He/⁴He signature of volcanic gases and fluid inclusions and we highlight analogous evolution of the ³He/⁴He signature in both during the last decades of eruptive activity (1990–2017), even during the same eruption. We show that the maximum enrichment in ³He is found in magmatic fluids that fed the most voluminous eruptions which culminated in caldera collapse events. We argue that this enrichment in ³He mostly reflects a greater contribution of magmatic fluids from a primitive component of the mantle plume. These results emphasize that He isotopes may provide warnings of increases in deep magmatic contributions that potentially herald paroxysmal eruptions, as documented here at Piton de la Fournaise (2007) and also at Kilauea (2018).

Noble gases are known to provide important clues on Earth's genesis, reservoirs and current dynamics<fnref n="1" />⁻⁴. In particular, the isotopes of helium (He) are of particular interest due to the primordial origin of ³He with respect to ⁴He that is continuously produced by the radioactive decay of U-Th that have distinct proportions in Earth's reservoirs⁵,⁶. Coupled with other noble gases characterised by distinct physical-chemical properties, He is also useful to characterize mantle metasomatism, recycling of crustal material in the mantle, diffusive fractionation, melt evolution during ascent and aging-degassing⁷⁻¹¹. In volcanic contexts, long- to short-term temporal variability of He isotope signature has been documented in (i) gas trapped in fluid inclusions (FI) in mafic minerals and in (ii) volcanic gases. In the former case variations in ³He/⁴He were attributed to the variable influence of crustal and mantle components in magmatic melts¹²⁻¹⁸, whilst in the latter case it was considered to be related to time-changing contribution of magmatic fluids with respect to crustal sources or to the arrival of deep undegassed magma¹⁹⁻²⁵. In this respect, recent studies have stressed the relevance of ³He/⁴He monitoring in volcanic gases to detect a range of unrest conditions from magmatic to phreato-magmatic¹⁶,¹⁸,¹⁹,²¹⁻²³,²⁵,²⁶. With the exception of a few intra-plate volcanoes on Earth (e.g., Etna, El Hierro), most of the studies of ³He/⁴He monitoring volcanic gases have been carried out in subduction-related settings (e.g., Stromboli, Santorini, Turrialba, Ontake) where long quiescent periods and/or the emission of variably differentiated lavas, often free of gas-rich fluid inclusions in mafic minerals, are common¹⁸,²⁰⁻²³. This makes comparing ³He/⁴He in volcanic gases and FI challenging at active volcanoes, even if it could potentially shed light on the ongoing magmatic dynamics as well as on magma

<footnote n="1">¹Université Clermont Auvergne, CNRS, IRD, OPGC, Laboratoire Magmas et Volcans, 6 avenue Blaise Pascal, 63178, Aubière, France. ²Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Palermo, Via Ugo La Malfa 153, 90146, Palermo, Italy. ³Dipartimento di Fisica e Scienze della Terra, Università degli Studi di Ferrara, 44121, Ferrara, Italy. ⁴Istituto Nazionale Di Geofisica e Vulcanologia, Osservatorio Vesuviano, Via Diocleziano 328, 80124, Napoli, Italy. ⁵Université de Paris, Institut de Physique du Globe de Paris, CNRS, F-75005, Paris, France. ⁶Observatoire Volcanologique du Piton de la Fournaise, Institut de Physique du Globe de Paris, F-97418, La Plaine des Cafres, France. ✉e-mail: guillaume.boudoire@uca.fr</footnote>
residence and crystallization timing and allow the forecasting of unrest phases. Not surprisingly, detailed and innovative studies based on this approach have been totally missing until present.

Here, we aim to compare He systematics in FI in crystals from lavas and cumulates (Table S1) of Piton de la Fournaise (PdF) with (i) Sr isotopes of host rocks/crystals (Table S1) and (ii) He isotopes in volcanic gases from thermal springs (Table S2) from Piton des Neiges (PdN). Piton des Neiges represents the largest (and oldest) volcanic edifice of La Réunion island, with the only permanent CO₂-rich thermal springs<sup>10,27</sup>. The Piton de la Fournaise edifice has progressively built up on PdN flank since about 0.5 Ma<sup>28</sup>. Since the beginning of the 18th century, activity at PdF has been characterized by frequent and, on average, short-lived (hours-to-weeks, rarely months) eruptions<sup>29,30</sup>. The high frequency of PdF eruptions allows the investigation of fast <sup>3</sup>He/<sup>4</sup>He fluctuations in response to the eruptive activity over short time scales (e.g., days to years). La Réunion island is also set in a hotspot-related context<sup>31</sup> favouring the emission of basalts containing an abundance of mafic minerals that are enriched in gas entrapped in FI with respect to lavas emitted in subduction-related contexts<sup>10,21</sup>. This last point is crucial to investigate the composition of gases in FI and to compare it with volcanic gases. Geochemically, a distinction is classically made between eccentric and central lavas emitted at PdF<sup>32,33</sup>. The former are sporadically emitted along the North West rift zone connecting PdN and PdF volcanoes (NWRZ; Fig. 1) and ascend rapidly from the mantle and the mantle-crust underplating layer<sup>33</sup>. The latter are emitted at a high frequency (with an average of one eruption every 9 months) inside the Enclos Fouqué caldera and along the South-East and North-East rift zones (SERZ and NERZ; Fig. 1) and are classically related to the shallower, crustal plumbing system of PdF<sup>10,32–34</sup>. Note that the Enclos Fouqué caldera (8 × 13 km) results from the most violent explosive activity (“Bellecombe” events) of PdF of the last 5 kyr<sup>35</sup>. The endmembers of the three datasets used in this study (He isotopes in FI, Sr isotopes from host rocks/crystals and He isotopes in volcanic gases) are detailed in Table 1 (see Table S1 for full data set).

**Helium isotope variability in the plumbing system.** He concentrations (Fig. 2a) in FI in olivine range from $2.9 \times 10^{-11}$ to $1.0 \times 10^{-12} \text{mol} \cdot \text{g}^{-1}$ for eccentric products (NWRZ) and decreases to $6.9 \times 10^{-12}$ to $3.2 \times 10^{-13} \text{mol} \cdot \text{g}^{-1}$ for central cases (Enclos Fouqué, SERZ, NERZ). The correlation found here between the He and the total gas content (TGC) released by crystal crushing (Fig. S1) is consistent with the general view that central and eccentric products reflect entrapment pressures in the crust and in the upper mantle, respectively, and consequently record different extents of magmatic degassing<sup>10</sup> (Table 1). In addition, <sup>4</sup>He/<sup>40</sup>Ar*, which increases by fractional degassing at PdF during magma ascent<sup>10</sup>, displays higher values in central products ($3.7 \pm 1.7$, on average) than in eccentric ones ($1.5 \pm 0.5$, on average) (Table 1). These latter values are well above the ratio reported for primary degassed volatiles from the mantle beneath La Réunion island (<sup>4</sup>He/<sup>40</sup>Ar* = $0.3 \pm 0.2$)<sup>10</sup>. This increase further confirms the more evolved and degassed nature of central magmas as previously indicated by magma and mineral composition and fractional equilibrium degassing modeling<sup>10</sup> (Figs. 2a and S1). Helium isotope ratio (<sup>3</sup>He/<sup>4</sup>He) corrected for atmospheric contamination (i.e., Rc/Ra; see Methods) varies from 11.9 to 15.3 Ra in FI from crystals of cumulate rocks and from 12.5 to 14.4 Ra in those of lavas (Fig. 2b). Our new and extended sampling reveals a greater variability in Rc/Ra than previously reported on the island<sup>36,37</sup>. More importantly, we found on average, higher values in central crustal products (Rc/Ra = $13.8 \pm 0.7$, on average) than for eccentric ones (Rc/Ra = $12.8 \pm 0.5$, on average), which record a deeper origin<sup>10</sup>. We highlight that the whole range of <sup>3</sup>He/<sup>4</sup>He variability is not related to crystal size<sup>38</sup> (all crushed crystals being >0.5 cm). No noticeable difference exists between the <sup>3</sup>He/<sup>4</sup>He of the gas entrapped in FI from olivine or clinopyroxene crystals from wehrlites<sup>39</sup> (Table S1).

Thermal (28–36 °C) springs (Irénée, Véronique, Manès) at Cilaos (Fig. 1) are located on PdN volcano, 30 km from the PdF summit<sup>10,27,40,41</sup>. Volcanic gases from these thermal springs are CO₂-dominated<sup>27</sup> (>90,000 µmol/mol; Table S2). Ratios of <sup>4</sup>He/<sup>40</sup>Ar* vary from 1.3 to 1.4 (Irénée source only). These values are consistent with those from FI and are representative of the depth of the upper mantle, supporting the hypothesis of dominant exsolution and release from this depth<sup>10</sup> (Table 1). In these gases, the two sampling campaigns performed between 2016 and 2017 revealed Rc/Ra variations from 13.1 to 13.6 at Irénée, between 13.4 and 14.0 at Véronique, and between 13.4 and 13.9 at Manès. Irrespective of the slight variation in Rc/Ra values between the three sources, it is important to highlight that the <sup>3</sup>He/<sup>4</sup>He values measured in 2016–2017 are systematically higher (ΔRa = 0.7–1.5) than those reported in 1990 (i.e., Rc/Ra = 12.1 at Irénée, 12.7 at Véronique and 12.7 at Manès)<sup>27</sup>.

**Origin of helium isotope variability in fluid inclusions.** Previous studies, focused on noble gases from both eccentric and central products of PdF and PdN, have described a restricted range (ΔRa < 1) of <sup>3</sup>He/<sup>4</sup>He values at La Réunion for at least the last million years<sup>37</sup>. Conversely, our new data highlight a greater variability (ΔRa = 3.4; Fig. 2b). The low variability documented in previous studies could mostly reflect the effect of the former sampling strategy focused on olivine rich “oceanite” basalts. Those basalts are characterized by phenocrysts mixed with antecrysts<sup>28</sup>. When “bulk” crystals are crushed for the extraction of FI, this mixing could lead to the homogenization of the <sup>3</sup>He/<sup>4</sup>He signature of the different crystal populations, thus providing a narrower range of values. Conversely, we investigate products ranging from almost aphyric to olivine-rich basalts and cumulates. This sampling strategy may well explain why we observe a wider range of <sup>3</sup>He/<sup>4</sup>He ratios than previous studies. This variability could be explained by different processes:

(1) Shallow crustal contamination. Although this is a possibility, the highest Rc/Ra values reported for central-crustal products with respect to eccentric-mantelic cases rule out shallow crustal contamination as the main process responsible for our <sup>3</sup>He/<sup>4</sup>He variability. Similar conclusions were reported in recent studies based on Sr and Nd isotopes at PdF<sup>42,43</sup> (Fig. 3a).

(2) Kinetic effects due to diffusive fractionation during melting of the mantle source region<sup>9,44</sup>. Considering the partition coefficients of He isotopes<sup>39</sup>, such a fractionation process would only have played limited role in our case because it should have led to both a decrease in <sup>3</sup>He/<sup>4</sup>He and <sup>4</sup>He/<sup>40</sup>Ar* in time, which is not observed. In
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location map of sampled sites at La Réunion (PdF for Piton de la Fournaise and PdN for Piton des Neiges volcanoes). Crushed crystals from lavas (circles) and cumulates (squares) in the NW rift zone (NWRZ; in red) and in the Enclos Fouqué caldera (EF; in green) and related NE-SE rift zones (NERZ, SERZ; in green). The green arrow represents the chronology of the products from the August 2015 eruption described in this study. Volcanic gases sampled in Cilaos thermal springs at PdN (in blue). Analytical protocols are reported in the Methods.</figcaption></figure>

addition, part of the ³He/⁴He variability occurs during a single eruption over a very short (days) time span (e.g., August 2015; Figs. 4 and 7; Table S1).

(3) Post-eruptive radiogenic ingrowth of ⁴He or cosmogenic ³He addition due to exposure to cosmic rays⁴⁵. We can exclude such processes here due to the variability documented over the very short (eruptive) time spans. During the August 2015 eruption, the ³He/⁴He varied from 13.4 Ra (single measurement) on October 7, 2015 to 14.3 ± 0.2 Ra (on average) on October 26, 2015 (green arrow on Fig. 2b). Furthermore, the analytical method used in this work to release fluid inclusions (single-step crushing) minimizes the effect of post-eruptive ingrowth.
<table><caption>**Table 1.** Characteristics (sample type, gas contents and ratios, He-Sr isotopes) of the main end-members selected from the Table S1 and discussed in this study. Ages and inferred crystallization depths$^{10}$. In bold, minimum and maximum values for each parameter. 'TGC' for Total Gas Content.</caption><thead><tr><th colspan="5">Samples information</th><th colspan="4">Gas contents and ratios</th><th colspan="2">Isotopes ratios</th></tr><tr><th>Location</th><th>Type</th><th>Material</th><th>Eruption/Site</th><th>Emission date</th><th>Name</th><th>TGC (mol/g)</th><th><sup>4</sup>He (mol/g)</th><th><sup>4</sup>He/<sup>40</sup>Ar*</th><th>He/CO<sub>2</sub></th><th>Inferred crystallisation depth</th><th>Rc/Ra</th><th><sup>87</sup>Sr/<sup>86</sup>Sr</th></tr></thead><tbody><tr><td>Central</td><td>Lavas</td><td>Crystals</td><td>February 2005</td><td>February 2, 2005</td><td>REU0502-221</td><td><strong>9.6E-10</strong></td><td><strong>9.8E-13</strong></td><td><strong>4.9</strong></td><td><strong>1.0E-03</strong></td><td>Oceanic crust</td><td><strong>12.5</strong></td><td>0.704181</td></tr><tr><td>Central</td><td>Lavas</td><td>Crystals</td><td>August 2015</td><td>October 26, 2015</td><td>REU151026-1</td><td>1.1E-09</td><td>3.2E-13</td><td>3.5</td><td>3.9E-04</td><td></td><td><strong>14.3</strong></td><td>0.704165</td></tr><tr><td>Peripheral</td><td>Lavas</td><td>Crystals</td><td>Piton de Caille</td><td>&lt;5 kyr</td><td>PDC-1</td><td>1.0E-07</td><td>7.1E-12</td><td><strong>1.3</strong></td><td>7.0E-05</td><td>Upper mantle: underplating layer</td><td>13.8</td><td><strong>0.704104</strong></td></tr><tr><td>Peripheral</td><td>Lavas</td><td>Crystals</td><td>Trous Blancs</td><td>&lt;10 kry</td><td>TB</td><td>2.9E-08</td><td>4.9E-12</td><td>2.0</td><td>1.7E-04</td><td></td><td>13.4</td><td><strong>0.704261</strong></td></tr><tr><td>Peripheral</td><td>Cumulates</td><td>Crystals</td><td>Chisny</td><td>≈ 381 BP</td><td>ILR84-4</td><td><strong>7.8E-07</strong></td><td><strong>2.9E-11</strong></td><td>1.4</td><td>3.7E-05</td><td>Upper mantle: deep melt horizon</td><td>12.8</td><td></td></tr><tr><td>Peripheral</td><td>Spring waters</td><td>Free gas</td><td>Cilaos (Irénée)</td><td>April 12, 2017</td><td>IRE170412-1</td><td></td><td></td><td>1.4</td><td><strong>1.2E-05</strong></td><td>Upper mantle</td><td>13.6</td><td></td></tr><tr><td>Peripheral</td><td>Spring waters</td><td>Free gas</td><td>Cilaos (Véronique)</td><td>April 12, 2017</td><td>VER170412-1</td><td></td><td></td><td>1.8</td><td>5.7E-5</td><td></td><td>14.0</td><td></td></tr><tr><td>Peripheral</td><td>Spring waters</td><td>Free gas</td><td>Cilaos (Manès)</td><td>April 12, 2017</td><td>MAN170412-2</td><td></td><td></td><td>2.2</td><td>3.3E-5</td><td></td><td>13.9</td><td></td></tr></tbody></table>

(4) Diffusive loss of noble gases out of crystals. This process may fit part of our dataset leading to lower Rc/Ra and $^{4}$He/$^{40}$Ar* values$^{9,46}$ after a few years (less than 8 years for crystals smaller than 1 cm; Fig. 3a). However, it cannot explain the synchronous variability found in our volcanic gas measurements and fluid inclusions (Table S2).

(5) Magmatic degassing coupled to aging. This process could decrease the $^{3}$He/$^{4}$He values$^{6,12}$. However, the lowest Rc/Ra values documented for eccentric mantle-derived products, i.e., for the most undegassed products with respect to central crustal cases$^{10}$ are not consistent with the effect of magma degassing coupled to aging. We calculate that melts having degassed about 90% of their primary volatile content can decrease Rc/Ra from 15.3 to 11.9 only after 0.5 Myr of residence at the depth of the underplating layer$^{10}$ (Fig. 3a). This time span is unrealistic considering that magma residence time at PdF is on timescale of thousands of years at maximum$^{37,47}$. Conversely, residual material potentially ponding in the mantle melting region$^{48}$ for about 5 Myr, i.e., in the range of ages proposed for the early construction of La Réunion island$^{49}$, may decrease the Rc/Ra signature (Fig. 3a). In this case we would expect a progressive decrease of Rc/Ra values from the Older Series at Mauritius, through PdN to the youngest PdF products. However, this behaviour is not observed$^{50}$.

Thus, we propose that the $^{3}$He/$^{4}$He variability documented at PdF is related to the contribution of various extents of mixing of magmatic fluids coming from either a “high Rc/Ra” component that is identified in a few samples (e.g., “primitive” component in the following parts) or a “low Rc/Ra” component predominant in most samples from literature and our dataset (e.g., “steady-state” component in the following parts) (Fig. 3b). Similar cases of mixing between magmatic fluids coming from distinct components have been also reported at other hot-spot volcanoes$^{1,45,51,52}$. Our $^{4}$He/$^{40}$Ar* data provide a new insight into the nature of this mixing of magmatic fluids. We document a greater variability in $^{4}$He/$^{40}$Ar* for central crustal products with respect to eccentric cases and to previous studies performed at La Réunion$^{36,37,50,53}$. This is irreconcilable with the effect of fractional equilibrium degassing only$^{10}$ (Fig. 2a). Melt mixing and degassing$^{10}$ should be expected to maintain constant $^{4}$He/$^{40}$Ar* at constant pressure (Fig. 3b), unless primary melts originate from mantle portions that suffered different percentages of melt extraction. However, our results show that flushing by a gas phase exsolved at mantelic depths (e.g. (2) or (3) on Figs. 3b and 4a) is required to explain both the Rc/Ra and $^{4}$He/$^{40}$Ar* variability observed in our dataset. At PdF, gas flushing is sustained by the presence of extensive degassing at the mantle level beneath La Réunion island$^{10}$. Furthermore, the solubility of volatile elements in melts, together with the presence of multiple magma ponding zones and mixing steps at PdF$^{33}$, can produce a migration of the gas phase across the plumbing system. This can occur faster than the rate of melt ascent, making the gas more susceptible to generating secondary FI$^{10,33,54}$. This idea is supported by the $^{3}$He/$^{4}$He, which is similar in the gas phase released by from both thermal springs and secondary FI in crystals from lavas emitted in the same period (e.g., 2015–2017; Table 1; Fig. 4a).

**Paroxysmal eruptions linked to a predominant contribution of magmatic fluids from a “primitive” mantle component.** In historical times (post-1860 CE), discrete and short-lived effusive eruptions associated with low lava fountains, which quickly evolve into strombolian activity, represent the main eruptive dynamics at PdF$^{30,55}$. During the period 1998–2010, i.e., one of the most recent and best documented eruptive periods of PdF, this “classical” eruptive style produced eruptions with average erupted volumes and durations of 6.6 Mm$^3$ and 25 days, respectively$^{56}$. This activity contrasts with rarer “paroxysmal” eruptions linked to caldera collapses and to the emission of unusually large volumes of lava$^{30,35,55,57}$. Since 1930, two major “paroxysmal” eruptions have been observed$^{30,57}$: in 1931 (130 Mm$^3$)$^{57}$ and in 2007 (240 Mm$^3$)$^{57}$, both were associated with summit caldera collapses. Following the last “paroxysmal” eruption in 2007, a new phase of magmatic recharge of the central plumbing system began in 2014–2015$^{58}$.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Gas contents and isotopes ratios obtained by crystal crushing from lavas and cumulates. (a) $ ^{4}\text{He} $ concentration vs. $ ^{4}\text{He}/^{40}\text{Ar}^* $. FED for Fractional Equilibrium Degassing curves. See Methods for supplementary information on the modelled curves. (b) $ ^{4}\text{He} $ vs. Rc/Ra (see Methods). Same legend than in Fig. 1 for the samples (circles and squares for analyses from crystals in lavas and cumulates, respectively). Dark and bright grey fields for published data at La Réunion island for lavas at PdF and at PdN, respectively. White field for published data at Mauritius island. Helium contents as a function of the part of the plumbing system (crust, underplating layer ‘under. layer’, deep melt horizon): squares for the full ranges of contents calculated (Fig. S1) from TGC values reported at PdF$^{10}$. $ \sigma_{error} $ (bar in the corner) for the analytical uncertainty (on Rc/Ra). $ \sigma_{replicates} $ (individual bar for each dot) for the natural variability (standard deviation on the mean for each sample for which replicates were measured).</figcaption></figure>

Since 2016, frequent sampling of the Cilaos thermal springs at PdN has revealed that the $ ^{3}\text{He}/^{4}\text{He} $ values (Rc/Ra = 13.5 ± 0.3, on average) in gases have increased with respect to those measured in 1990 at the same sites$^{27}$ (Rc/Ra = 12.5 ± 0.3, on average; Fig. 4a). Importantly, our new sampling of eruptive products shows a similar increase in FI of products from syn- and post-April 2007 lavas (Rc/Ra = 14.0 ± 0.4, on average; Fig. 4b; Table 1) with respect to those from the 1986–2005 period (Rc/Ra = 12.8 ± 0.3, on average). The similarity of $ ^{3}\text{He}/^{4}\text{He} $ values in the gases related either to fluid inclusions or to thermal springs over time clearly indicates that ascent fluids are trapped in minerals in a short-time span that allows FI to preserve the same signature measured in volcanic gases. This result opens major prospects for future monitoring strategies. At PdF, the temporal evolution of Rc/Ra in gas released by crystal crushing from lavas emitted since 1931 shows two main events marked by Rc/Ra values greater than 13.7 (threshold value; see Methods): in 1931 and in April 2007 (and post-2007 products) (Fig. 4b). Together with products related to the Bellecombe's explosive events (5.5–3.0 kyr BP; Rc/Ra = 14.3 ± 0.6, on average), these high $ ^{3}\text{He}/^{4}\text{He} $ values are linked to paroxysmal eruptions$^{30,35,57}$ triggering caldera collapse and the emission of unusually high lava volumes (>100 Mm$^{3}$). The 2007 paroxysmal eruption was triggered by a major input of deep magma that provoked the extrusion of a large amount of magma from the crustal plumbing system of PdF$^{58}$. We thus infer that a predominant contribution (>95%) of the “primitive” mantle component leads to the production and ascent of large amount of melts and magmatic fluids with Rc/Ra >13.7 responsible of unusual paroxysmal eruptions. Conversely, the more frequent, small volume and short-lived eruptive activity at PdF would be related to a predominant contribution (>95%) of magmatic fluids from the “steady-state” component
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Noble gases isotopes ratios obtained by crustal crushing from lavas and cumulates. (a) $^4$He/$^{40}$Ar* vs. Rc/Ra modelling the effect of diffusion (green dashed lines; the star highlights the effect of the diffusion on a crystal of 0.5 or 1 cm of diameter after 2 or 8 years, respectively) and the combined effect of magma aging and degassing (black dashed lines). The extent of degassing is reported (e.g., percentage). (b) $^4$He/$^{40}$Ar* vs. Rc/Ra modelling the effect of degassing, gas flushing and melt mixing resulting from the melting/degassing of two distinct components of the mantle plume below La Réunion island. The part of the plumbing system where the volatile phase is exsolved (starting conditions) to generate gas flushing in the upper parts is characterised by its average $^4$He/$^{40}$Ar* value¹⁰: (0) oceanic crust, (1) mantle-crust underplating layer, (2) deep melt horizon and (3) mantle source. The extent of gas flushing is reported (e.g., percentages). See Methods for supplementary information on the modelled curves. Same legend than for Fig. 2. Dark grey field for previous data at PdF. The meaning of $\sigma_{error}$ and $\sigma_{replicates}$ is detailed in the Methods and the caption of Fig. 2.</figcaption></figure>

(Fig. 4b). This is the case of the February 2005 eruption (Rc/Ra from 12.5 to 13.1), i.e. the last eruption before the 2007 caldera collapse eruption from which we have get measurements of He isotopes.

**Changes in syn-eruptive magma dynamics deciphered using helium isotopes.**
$^3$He/$^4$He may also be used to track rapid changes in magma dynamics during a single eruption, as occurred in August 24 – October 31, 2015 (Fig. 5). This eruption, one of the longest (2 months) since 1930 (and involving $45 \pm 15$ Mm$^3$), evolved through three main phases each related to a distinct magma composition$^{58,59}$. Frequent sampling permitted the identification of a progressive change in magma composition over time towards a more mafic composition (see MgO content on Fig. 5a) and a progressive enrichment of lavas in olivine phenocrysts which began in early October 2015. In addition, three distinct behaviours in He isotopes data were observed during the eruption:

Lava sampled in October 7 had an intermediate $^3$He/$^4$He signature (Rc/Ra = 13.4; single measurement). Two interpretations may be considered. In the former case the Rc/Ra signature may be that of a high $^3$He/$^4$He melt (>13.7 Ra), like those emitted in April 2007, still ponding at shallow level ($^4$He/$^{40}$Ar* = 3.9) for many years and flushed by a gas phase related to steady-state degassing conditions (<12.8 Ra; Fig. 4a). This possibility is supported by petrological evidences that the 2008–2014 eruptions were fed by a 'April 2007'-derived melt that differentiated and cooled in the shallow part of the plumbing system$^{60}$. However, the arrival of a high $^3$He/$^4$He melt (>13.7 Ra) that progressively mixed with a steady-state melt already ponding at shallow level (<12.8 Ra) is also an arguable hypothesis considering the reverse zoning of a few host olivine phenocrysts$^{41}$. In both cases,
<figure><img src="image_6.png" /></figure>

**Figure 4.** Noble gases isotopes ratios obtained by crystal crushing from lavas and compared to volcanic gases from thermal springs at Cilaos. (a) Rc/Ra vs. $^4He/$$^{40}Ar$$*$ modelling the effect of magmatic fluids flushing and melts mixing starting from the highest Rc/Ra vs. $^4He/$$^{40}Ar$$*$ values. Thresholds of the >95%-contribution of either the “primitive” mantle component or the “steady-state” mantle component to Rc/Ra value determined using the maximum-likelihood method$^{59}$. See Methods for supplementary information on the modelled curves and the statistical method. Same legend than for Figs. 2 and 3. Blue field for the range of compositions from volcanic gases from thermal springs. Dark grey field for published data from lavas at PdF. The green arrow represents the chronology of the products from the August 2015 eruption. $^4He/$$^{40}Ar$$*$ variability as a function of the location of magma ponding in the plumbing system$^{10}$. (b) Time variations of Rc/Ra from both FI in olivine crystals from PdF lavas and from volcanic gases. During the last century$^{30,57}$, major eruptive events (caldera collapse and erupted volume >100 Mm$^3$) occured only in 1931 and 2007. The meaning of $\sigma_{error}$ and $\sigma_{replicates}$ is detailed in the Methods and the caption of Fig. 2.

this “mixing” process between magmatic fluids having distinct isotopic signature can be favoured by the presence of vertically extended magma storage zone, where magma ponding favours melt homogenisation and fluid re-equilibration during ascent$^{10}$ and fits most of our data (Fig. 4a). Such hybrid melts might be drained during the first phases of the eruption$^{58}$.

Since mid-October 2015, a progressive increase in Rc/Ra at quite constant $^4He/$$^{40}Ar$$*$ value was observed (October 16 and 23, 2015; Fig. 5). We interpret this behaviour as the effect of syn-eruptive refilling of the crustal part of the plumbing system resulting in the emission of more primitive melts that triggered two successive short-lived eruptive phases on October 22–24 and October 29–31 (Fig. 4a). The increase in soil CO$_2$ flux, plume SO$_2$ emissions and the emission of more mafic melts (MgO-rich) hosting zoned olivine crystals from mid-October 2015 support our assumption$^{58,59,61}$.

Lower $^3He/$$^4He$ values (Rc/Ra = 13.7; single measurement) but also highest $^4He/$$^{40}Ar$$*$ (= 6.4; single measurement) were obtained for crystals from lava erupted a few minutes before the end of the eruption. These values are consistent with the presence of cumulative olivine crystals entrained upon flushing of the upper part of the magma conduit$^{10,28,58}$, i.e., having higher $^4He/$$^{40}Ar$$*. We propose that such flushing of crystals may be favoured by an increase in the final magma output rate linked to closure of the dyke. This assumption is supported by a lava overflow that also preceeded the end of the eruption by a few minutes, and emptying of the lava lake inside the main eruptive vent (OVPF report). The lowest Rc/Ra values for this sample are associated with an increase
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Time evolution of noble gases isotopes ratios in gas-released by olivine crystals crushing from lavas during the August 2015 eruption. (a) Evolution of ³He/⁴He compared to hosting lava chemistry⁵⁸. (b) Evolution of ⁴He/⁴⁰Ar* and ⁴He/⁴⁰Ar* variability as a function of the location of magma ponding in the plumbing system¹⁰. Same legend as for Fig. 7. The meaning of σ<sub>error</sub> and σ<sub>replicates</sub> is detailed in the Methods and the caption of Fig. 2.</figcaption></figure>

of MgO contents and olivine crystals in the host lavas, as also documented for the mid-October products, but with an increase in Th content (Fig. 5a). This distinct behaviour of magma chemistry in this final sample of the August-October 2015 eruption mirrors the distinct behaviour observed for He isotope.

This eruption demonstrates that variations of ³He/⁴He may be documented on timescale of days and are closely linked to the dynamics of syn-eruptive magma transfer.

**Constraints on the “primitive” and “steady-state” mantle components.** Even if the ³He/⁴He variability that we document here may be related to various contributions of mantelic fluids from either a “primitive” component (data set for the “paroxysmal” eruptions) or a “steady-state” one (data set for the “classical” eruptions), the exact source of these two components is still unclear. This problem is even more intriguing considering that, if magmas at La Réunion are marked by the Dupal anomaly like most of the Ocean Island Basalts (OIBs) and Mid Ocean Ridge Basalts (MORBs) of the Indian Ocean⁶², the mantle source is considered extremely homogeneous with respect to other ocean basaltic islands²⁸,³⁷,⁵⁰,⁶³,⁶⁴. Only minor variations of Sr-Nd-Hf isotopes and incompatible trace elements ratios have been reported and mostly attributed to a primary feature of La Réunion mantle plume⁴²,⁴³,⁶⁴⁻⁶⁷. We argue here that the ³He/⁴He scale of variability (ΔRa greater than 3 units Ra) is also a primary feature of the mantle source beneath La Réunion island. Such an idea was previously developed for Hawaii, where ³He/⁴He fluctuations up to 8 Ra on a timescale of a few hundred years were attributed to the melting of primary heterogeneities of the mantle plume¹⁴.

In the Indian Ocean, lower ³He/⁴He values were measured at Rodrigues and at Mauritius (Younger Series) and attributed to ridge-hot spot interaction, i.e., related to a contribution of magmatic fluids with a lower ³He/⁴He sig-nature from the Central Indian Ridge (CIR) with respect to the mantle plume (Fig. 6a)⁴⁵,⁵⁰. At PdF, we do not find any correlation between ³He/⁴He and ⁸⁷Sr/⁸⁶Sr (Table 1; Fig. 6a) thus suggesting mixing between such distinct res-ervoirs¹,¹³,¹⁴,³⁷,⁵²,⁶⁷,⁶⁸. Instead, decade-scale fluctuations of Sr isotopes have been documented for recent PdF lavas and attributed to the fertility of the source and induced melt channelization⁴³ (Fig. 6b). In view of the absence of correlation between He and Sr isotopes, we propose that the ³He/⁴He decrease in some PdF products may be related to intrinsic features of the mantle plume beneath La Réunion island. This fits with other recent hypothesis showing that the composition of La Réunion mantle plume may derive from a mantle source that experienced an early-stage (Hadean) of differentiation generating distinct reservoirs or blobs within the mantle plume³⁷,⁴²,⁶⁴.
Similar episodic entrainment by thermal plumes of deep isolated mantle reservoirs of evolved versus primordial materials was previously reported in ocean island basalts⁶⁹. In this context, the suspected global ³He/⁴He increase from Older Series at Mauritius (not contaminated by the Central Indian Ridge) (Rc/Ra = 11.7 ± 0.1, on average; Fig. 2b) to present-day products at PdF (Rc/Ra = 13.3 ± 0.7, on average; Fig. 2b) may reflect such progressive melting of the more primitive component of the mantle plume⁵⁰.

On the other hand, the composition of some melts at PdF³³,⁶⁵ has suggested the presence of magmatic fluids enriched in U-Th, down to the upper mantle, which may be responsible for hydrothermal contamination in the PdF plumbing system¹⁰,³³,⁶⁵. The enrichment in U-Th in such magmatic fluids suggests a potential enrichment in ⁴He by radioactive decay in these fluids, which may lower the ³He/⁴He signature in some PdF products. For instance, we have reported a lower averaged Rc/Ra signature (Rc/Ra = 12.8 ± 0.5, on average) in eccentric products with respect to central crustal cases (Rc/Ra = 13.8 ± 0.7, on average), i.e., for products testifying to magma ponding at the depth of the mantle/crust underplating layer and even deeper (Fig. 2b). Such magma ponding, mixing and degassing are known to be important processes occurring in the upper mantle beneath La Réunion island³³,⁴¹. Thus, in addition to some mantle plume heterogeneities, we cannot exclude the role of contamination by ⁴He-U-Th-rich magmatic fluids in the lithospheric mantle able to reduce the ³He/⁴He signature of melts ponding in the plumbing system¹⁰,⁶⁵. Such decoupled ascent between ⁴He-U-Th-rich magmatic fluids and melts could be an additional argument supporting the apparent decoupling between He isotopes (gas in FI and volcanic gases) and Sr isotopes (melt and crystal phase) (Fig. 6).

## Conclusion

Through the study of the ³He/⁴He signature of volcanic gases and fluid inclusions at La Réunion island, we infer that a predominant contribution (>95%) of a primitive component of the mantle plume leads to the production and ascent of melts and magmatic fluids (Rc/Ra > 13.7) responsible for unusual “paroxysmal” eruptions, as observed for the most recent caldera collapse eruptions at Piton de la Fournaise (PdF) (Fig. 7). Conversely, the more frequent, small volume and short-lived eruptive activity at PdF would be related to (i) a predominant contribution (>95%) of magmatic fluids from a more degassed mantle component and/or (ii) greater residence times for melts that may be progressively contaminated by ⁴He-rich magmatic fluids coming from the lithosphere and producing steady-state magmatic conditions (Rc/Ra < 12.8). These steady-state conditions are well recorded in “classical” eruptions such as that of 2005 (Fig. 7).

Our results suggest that the eruptive activity at hotspot volcanoes may respond to fast changes in mantle dynamics and lead to “paroxysmal” eruptions with serious implications in terms of hazard⁷⁰. We emphasize that ³He/⁴He monitoring in volcanic gases may forecast the arrival of magmatic fluids triggering “paroxysmal” eruptions potentially related to caldera collapses, as observed at PdF (2007) and at Kilauea²⁵ (2018). Considering the synchronous variations of ³He/⁴He values observed in volcanic gases and in gas from FI, our study opens exciting perspectives for assessing magmatic unrest, as well as syn-eruptive changes in magma dynamics, beneath more hazardous volcanoes globally. Our study supports the idea that coupling petrology (fluid inclusions barometry and petrography) with gas geochemistry (He-Ar isotopes) is fundamental to understand the variability of helium isotopes and their use for volcano monitoring.

## Methods

**Sr isotopes measurements.** 8 analyses of Sr isotopes ratios (⁸⁷Sr/⁸⁶Sr) have been performed at INGV-OV (Istituto Nazionale di Geofisica e Vulcanologia – Osservatorio Vesuviano di Napoli, Italy). Isotopic compositions have been determined on ca. 0.1 grams of crystal powders, residues of the crushing for He isotopes analyses, after dissolution with ultrapure acids (HF-HNO₃-HCl mixtures). Sr was separated from the mineral matrix through conventional ion-exchange procedures (standard separation scheme⁷¹ for details) and measured statically by thermal ionization mass-spectrometry (TIMS, Thermo Finnigan Triton TI). During collection of isotopic data, replicate analyses of NIST-SRM 987 (SrCO₃) international reference standards were measured. Mean measured values for ⁸⁷Sr/⁸⁶Sr for the NIST-SRM 987 standard were 0.710219 ± 0.000018 (2σ, N = 101). Measured ⁸⁷Sr/⁸⁶Sr have been normalized for within-run isotopic fractionation to ⁸⁶Sr/⁸⁸Sr = 0.1194 and have been normalized to the recommended value of ⁸⁷Sr/⁸⁶Sr = 0.71025 for NIST-SRM 987 standard. On products sampled the same day during the August 2015 eruption (REU151031) we found a difference <0.000027 on ⁸⁷Sr/⁸⁶Sr signature (i.e. in the range of the 2σ uncertainty on measurements) between our analysis performed in olivine crystals and those previously performed from bulk rocks⁴³. Additionally, we have reported (Table S1) 8 analyses of Sr isotopes ratios (⁸⁷Sr/⁸⁶Sr) from literature performed on the same products that have been analysed in this study for noble gases⁴¹,⁴³.

**Noble gases measurements.** Samples analysed in our study from the central area spans a recent eruptive period (2002 to 2015). Eccentric ones were mostly erupted after the Bellecombe explosive events generating the Enclos Fouqué caldera³⁵ (≈5.5–3.0 kyr BP). Olivine crystals without impurities (as well as clinopyroxene crystals for two wehrlitic cumulates) were handpicked from 30 samples previously crushed and sieved in the fractions of 0.5 and 1 mm. Minerals were first cleaned ultrasonically in 6.5% HNO₃, second in deionized water and third in suprapure acetone. After being carefully weighed, the samples were loaded into a stainless steel crusher that allows to holding up to six samples simultaneously for noble gas extraction. FI were released by in-vacuo single-step crushing of minerals after that about 200 bar-pressure was applied by a hydraulic press. This procedure strongly reduces or avoids that cosmogenic ³He and radiogenic ⁴He that could possibly have grown or been trapped in the crystal lattice are released together with FI¹¹,⁷²⁻⁷⁵. However, since our samples are very recent, any cosmogenic effect can be thus excluded. We determined the Total Gas Content (TGC) during crushing and FI extraction by measuring the total gas pressure (CO₂ + N₂ + O₂ + noble gases). Then, a “cold finger” present in the crusher was immersed in liquid N₂ at –196 °C to remove CO₂ (and H₂O, if present). The residual pressure of N₂ + O₂ + noble gases was quantified and subtracted from the TGC. The gas mixture was then purified under
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Comparison between helium isotopes ($^{3}$He/$^{4}$He) and strontium isotopes ($^{87}$Sr/$^{86}$Sr) at PdF. (a) $^{87}$Sr/$^{86}$Sr vs. $^{3}$He/$^{4}$He. Dark grey field for published data from lavas at PdF$^{37}$. CIR for the Central Indian Ridge composition$^{72}$. Same legend than in Fig. 1 for our samples. (b) Time evolution of $^{3}$He/$^{4}$He and $^{87}$Sr/$^{86}$Sr since 1931 from lavas. $^{3}$He/$^{4}$He datasets before 2002 from previous studies$^{36,37}$. $^{3}$He/$^{4}$He after 2002 acquired in this study. $^{87}$Sr/$^{86}$Sr dataset completed$^{37}$ in this study for products emitted after 2002. New data are highlighted by filled circles. The meaning of $\sigma_{error}$ and $\sigma_{replicates}$ is detailed in the Methods and the caption of Fig. 2.</figcaption></figure>

getters in an ultra-high-vacuum ($10^{-9}$–$10^{-10}$ mbar) line in order to leave only noble gases. A cold finger with active charcoal immersed in liquid N$_{2}$ removed Ar, while He and Ne were separated by using a cold head preven- tively cooled at 10 K and then moved at 40 and 80 K by a temperature controller connected to a heater around the cold head in order to release He and Ne, respectively. A similar procedure was adopted in previous studies$^{11,75}$.

Helium ($^{3}$He and $^{4}$He) and neon ($^{20}$Ne) isotopes were measured separately by two different split-flight-tube mass spectrometers (Helix SFT-Thermo). Argon isotopes ($^{36}$Ar, $^{38}$Ar, and $^{40}$Ar) were analysed by a multicollector mass spectrometer (GVI Argus)$^{10}$. The analytical uncertainty in the measured concentrations of $^{4}$He, $^{20}$Ne and $^{40}$Ar was <0.1% (analytical uncertainty reported as $\sigma$). Typical blanks for $^{4}$He, $^{20}$Ne and $^{40}$Ar were <10$^{-16}$, <10$^{-16}$, and <10$^{-14}$ mol, respectively, which implies that blank incides on the sample peaks at the mass spectrometers of the above atomic mass units for <0.03%, <6%, and <0.4%, respectively. The protocol for the preparation, single-step crushing and analysis of FI was the same as that applied for other studies at the INGV, Sezione di Palermo (Italy)$^{10,11}$. We do not document relevant shift between noble gases measurements performed in this study and previous measurements performed in another laboratory on the same sample (i.e. ILR84-4)$^{36,37}$, which supports the accuracy of our measurements. Literature data at La Réunion island (PdF and PdN) and from Mauritius island are compared to our results$^{36,37,50,76}$. At Mauritius, only results from the Older Series are reported due to their direct link with the mantle plume with respect to Youngers Series being MORB-contaminated$^{45,50}$. In total, 72 measurements have been carried out for noble gases and TGC. A variable aliquot ranging from 0.08 g to 1.36 g of crystals was analysed in replicate measurements for the elemental and isotope composition of helium, neon, and argon in the gas trapped in FI (between 1 and 9 replicates per sample; the 9 replicates for PDC-1 repre- sent a total weight up to 9.26 g of crystals). Analytical results are reported in Table S1. The analytical uncertainty of single measurements is reported as “$\sigma_{error}$” for Rc/Ra and $^{40}$Ar/$^{36}$Ar in Table S1 and figures. The natural variability of replicates (e.g. reproducibility) is reported as “$\sigma_{replicates}$” (standard deviation on the mean of the measurements
<figure><img src="image_9.png" /><figcaption>**Figure 7.** Conceptual model of the magmatic plumbing system of PdF<sup>40,61</sup>. Model 1 for major fresh melt input leading to « paroxysmal » eruptions as observed in 1931 and 2007 caldera collapse eruptions (Rc/Ra >13.5 in both FI and volcanic gases). There, magmatic fluids are mainly related to a primitive component of the mantle plume (Rc/Ra >13.7). Model 2 for steady-state melt feeding leading to « classical » eruptions as observed in 2005 (Rc/Ra <12.8 in both FI and volcanic gases). There, magmatic fluids are mainly related to an old degassed buoyancy component of the mantle plume (Rc/Ra <12.8) and/or to the contribution of <sup>4</sup>He-rich magmatic fluids from the upper mantle. Pictures from the first author (CC BY 4.0).</figcaption></figure>

obtained for each replicates) for either the elemental concentrations or for the isotopic ratios. These analyses were
performed in the noble gas isotope laboratory of INGV, Sezione di Palermo (Italy).

Bubbling gaseous samples from Cilaos thermal springs were collected as follows: because thermal springs are captured by three different plexiglass cups, gas was directly pumped from the vent placed at the top of each cup into a 100 cm³ syringe and then injected into the sample containers (with volume about 20–30 cm³) by using three-way valve. Air was purged away from containers by means of several cycles of pumping and injection. Two-ways Pyrex bottles, with vacuum valves at both ends, were used as containers for measurements of major and reactive gases. Two-ways stainless steel tube, with swagelok valves at both ends, was instead employed to sample gas for noble gas analyses. In all samples, the concentrations of CO₂, CH₄, N₂, O₂, He, CO, and H₂ were routinely measured at INGV, Sezione di Palermo (Italy) by a Perkin Elmer Clarus 500 gas chromatograph equipped with a 3.5-m Carboxen 1000 column and double detector (hot-wire detector and flame ionization detector), with analytical errors of <3%<sup>18,26,77</sup>. Gas aliquots for noble gas analyses were introduced into three distinct stainless-steel ultra-high-vacuum lines for standard purification procedures to remove major and reactive species. In detail, each of the three preparation lines is dedicated to a single noble gas species, in our case He, Ne, and Ar, and is equipped with a pipette system connected to a manometer to introduce a known number of moles of gas sample. The purification procedure for helium, neon, and argon is in principle the same adopted for fluid inclusions above described and reported in previous studies for free gases<sup>11,21,77</sup>. The abundances and isotope compositions
of He were determined by a split flight tube mass spectrometer (Helix SFT-GVI). Neon abundance and isotope composition ($^{20}$Ne) was determined by a Helix MC Plus Thermo. The abundances and isotope compositions of Ar were measured in a multicollector mass spectrometer (Helix MC-GVI). The analytical errors of the He, Ne, and Ar-isotope analyses were less than 0.4%, 0.06% and 0.1%, respectively. Internal standard was purified from air, whose $^{3}$He/$^{4}$He reproducibility over 1 year of daily analysis was <3.5%, $^{20}$Ne reproducibility over 1 year of daily analysis was <4%, $^{40}$Ar/$^{36}$Ar reproducibility over 1 year of daily analysis was <3.5%. Typical blanks for He, Ne, and Ar were $<10^{-16}$, $<10^{-16}$ mol, and $<10^{-14}$ mol, respectively, being at least two orders of magnitude lower than sample signals at the mass spectrometers. Detailed measurements were performed in October 2016 and April 2017 for the 3 sources of Cilaos thermal springs (Irénée, Véronique, Manès). Historical (1990) data from the same sources are compared to our results$^{27}$.

**Noble gases corrections from the atmospheric component.** A correction from atmospheric contamination is necessary in order to obtain the so-called magmatic-derived $^{40}$Ar*. This has been calculated, accordingly with the following formula that assumes that all the measured $^{36}$Ar is of atmospheric origin:

$$40Ar^* = 40Ar_m - \left(\frac{40Ar}{36Ar}\right)_{air} \times 36Ar_m$$

where $^{40}$Ar* represents the corrected $^{40}$Ar ($^{40}$Ar/$^{36}$Ar = 295.5), and the “m” subscript indicates “measured”.

The $^{3}$He/$^{4}$He ratio is expressed as R/Ra (being Ra the He isotope ratio of air and equal to $1.39 \times 10^{-6}$) and was corrected from atmospheric contamination following the equation$^{78}$:

$$Rc/Ra = \frac{(R \times N - Ra \times Na)}{(N - Na)}$$

where Rc/Ra, R and Ra are $^{3}$He/$^{4}$He of the corrected, observed and atmospheric components, while N and Na are $^{4}$He/$^{20}$Ne measured and atmospheric components (Na = 0.318), respectively. The analytical uncertainty ($\sigma_{error}$ reported in Table S1) does not exceed 0.2 Ra on average.

**Noble gases modelling.** Various fractionation processes were used to model the noble gases variability described in our study. We here detail the equations and parameters used to perform this modelling.

### (1) Fractional equilibrium degassing

Fractional equilibrium degassing (FED) was modelled from the equation$^{79}$:

$$\left(\frac{^{4}He}{^{40}_{v}Ar^{*}}\right)_{FED} = \left(\frac{^{4}He}{^{40}_{l}Ar^{*}}\right)_{ini} \times \alpha \times F^{\alpha-1}$$

where $^{4}He/^{40}_{v}Ar^{*}$ is the ratio in the vapor phase following a FED trajectory, $^{4}He/^{40}_{l}Ar^{*}$ is the initial ratio in the melt, F is the fraction of residual gas (between 1 and 0) and $\alpha$ is the Ar/He relative solubility. In Fig. 2, $^{4}He/^{40}_{l}Ar^{*}$ is equal to 12.7 and 27 for eccentric (e.g. PDC-1; red) and central (e.g. ol; green) products and assuming a Ar/He solubility of 0.10 and 0.11, respectively$^{10}$. In Fig. S1, $^{4}He/^{40}_{l}Ar^{*}$ range from 1 to 5 for FED modeling in primary conditions$^{10}$ and assuming a Ar/He solubility of 0.10. The same equation is also used in Fig. S1 substituting $^{4}He/CO_{2}$ to $^{4}He/CO_{2}$ and assuming that $^{4}He/CO_{2}$ varies from $1.2 \times 10^{-5}$ to $2.2 \times 10^{-5}$ and the He/CO$_{2}$ solubility equal to 1.7.

### (2) Diffusion

To model noble gases diffusion across crystals we have made the assumption that crystals can be considered like perfect spheres using the equation$^{80}$:

$$\frac{C(t, r) - C_1}{C_0 - C_1} = 1 + \frac{2 \times a}{\pi \times r} \times \sum_{n=1}^{\infty} \frac{(-1)^n}{n} \times \sin\left(\frac{n \times \pi \times r}{a}\right) \times \exp\left(\frac{-D \times n^2 \times \pi^2 \times t}{a^2}\right)$$

where $C_1$ is the initial concentration of the sphere (here the studied gas content in the crystal, i.e. $^{3}$He, $^{4}$He or $^{40}$Ar*), $C_0$ is the surface concentration (here considered null for modelling noble gases diffusion out of the crystal)$^{9}$, a is the half-diameter of the sphere (here varying between 0.25 and 0.5 cm in accordance with the size of the crystals analysed in this study), t the time (in s) and r the distance from the center of the crystal (here r = a/2$^{1/3}$ in order to consider $C(t,r)$ as the average concentration of the sphere. The equation was resolved using a homemade script in Python 2.7 language using the SciPy, NumPy and math libraries. Diffusion coefficients D across olivine crystals were set at $1.61 \times 10^{-10}$, $1.85 \times 10^{-10}$ and $5.09 \times 10^{-11}$ cm$^2$/s for $^{4}$He, $^{3}$He and $^{40}$Ar, respectively$^{9,46}$. Initial conditions were fixed at Rc/Ra = 15.3 (maximum value in our dataset; e.g., CUM2) and $^{4}$He/$^{40}$Ar* varying from 4.5 (e.g., CUM2) to 10 inside the crystal and at 0 outside (i.e. negligible in the surrounding melt with respect to the crystal like proposed in previous modelings$^{7}$).

### (3) Aging

The effect of radiogenic $^{4}$He production from the decay of U and Th over time ($^{4}$He*) can be estimated according to the equation$^{81}$:

$$4He^* = 2.80 \times 10^{-8} \times [U] \times \left(4.35 + \frac{Th}{U}\right) \times t \left(cm^3STP/g\right)$$
where t is the time in Ma, [U] is the uranium concentration (0.5 ppm on average in PdF melts)⁴¹ and Th/U is the thorium-uranium atomic ratio (3.9 for La Réunion)⁴⁵,⁴⁷,⁸². We estimate the resulting ⁴He* production at 5.18 × 10⁻¹⁸ ol/g per year. Initial conditions were set at ⁴He/⁴⁰Ar* = 0.3 (primary conditions from the first exsolved phase)¹⁰ and Rc/Ra = 15.3 (maximum value in our dataset; e.g., CUM2). We have modelled (i) melt degassing (F varying from 0 to 1 to represent the extent of FED) without aging, (ii) melt degassing with aging (0.5 Myrs i.e. similar to the age of the beginning of the PdF subaerial volcanic activity)⁴⁹ and, (iii) melt without degassing but with the effect of aging (5 Myrs i.e. related to the end of the Older Series at Mauritius island and at the beginning of the volcanic activity at La Réunion island)⁴⁹,⁵⁰.

### (4) Mixing and statistical analysis

On Fig. 3, gas flushing was modelled as an arrival of a deeper magmatic gaseous phase mixing with the volatile phase already exsolved at the depth of magma ponding. We here modelled through isotopic and mass balances equations the effect of gas flushing of a gas phase similar to that of CUM2 (⁴He/⁴⁰Ar* = 4.5 i.e. ponding at crustal level), which have the highest ³He/⁴He ratio of our dataset (e.g., 15.3 Ra), with a volatile phase arriving from the deep melt horizon or the underplating layer with the lowest ³He/⁴He signature found in our dataset (e.g., 11.9 Ra) and ⁴He/⁴⁰Ar* = 0.3, 0.8 and 2.1¹⁰. The orange field shows the whole field of theoretical evolution of the gaseous phase considering the uncertainty relative to ⁴He/⁴⁰Ar* signatures¹⁰. The inverse process (e.g. gas flushing through the arrival of a volatile phase with the highest Rc/Ra signature in a gas phase marked by the lowest Rc/Ra signature) is also modelled with the red field and curves. On the contrary, melt mixing is expected to product a ³He/⁴He variability in the exsolved phase at constant ⁴He/⁴⁰Ar*. On Fig. 4, we have modelled the effect of gas flushing (with a ³He/⁴He signature of 12.5 Ra, i.e. similar to the minimum signature in our dataset for lava samples: REU0402-221 from the February 2005 eruption) on a volatile phase marked by the highest Rc/Ra signature (e.g., 14.4 Ra in REU0704-52 from the April 2007 eruption) and the highest ⁴He/⁴⁰Ar* value (e.g. 8.5 in REU0704-51 from the April 2007 eruption) in our dataset for lava samples.

Thresholds of the contribution (>95%) of either the “primitive” (³He/⁴He > 13.7 Ra) or “steady-state” (³He/⁴He < 12.8 Ra) mantle component were determined from the 47 analysis obtained from PdF lava samples (our study and literature; Table S1) with the maximum-likelihood method using a Gaussian Mixture Model (GMM) implementing an expectation-maximization (EM) algorithm developed in Python 2.7 language with the NumPy and sklearn libraries⁵⁹,⁸³.

Received: 21 October 2019; Accepted: 4 May 2020;

Published online: 17 June 2020

## References

1. Kurz, M. D., Jenkins, W. J., Schilling, J. G. & Hart, S. R. Helium isotopic variations in the mantle beneath the central North Atlantic Ocean. *Earth and Planetary Science Letters* 58(1), 1–14 (1982).

2. Ballentine, C. J. & Burnard, P. G. Production, release and transport of noble gases in the continental crust. *Reviews in mineralogy and geochemistry* 47(1), 481–538 (2002).

3. Moreira, M. & Allègre, C. J. Rare gas systematics on Mid Atlantic ridge (37–40 N). *Earth and Planetary Science Letters* 198(3-4), 401–416 (2002).

4. Moreira, M. Noble gas constraints on the origin and evolution of Earth's volatiles. *Geochemical Perspectives* 2(2), 229–230 (2013).

5. Jenkins, W. J., Edmond, J. M. & Corliss, J. B. Excess 3He and 4He in Galapagos submarine hydrothermal waters. *Nature* 272(5649), 156 (1978).

6. Brooker, R. A. et al. The 'zero charge' partitioning behaviour of noble gases during mantle melting. *Nature* 423(6941), 738 (2003).

7. Burnard, P., Harrison, D., Turner, G. & Nesbitt, R. Degassing and contamination of noble gases in Mid-Atlantic Ridge basalts. *Geochemistry, Geophysics, Geosystems* 4(1), 1–20 (2003).

8. Shaw, A. M., Hilton, D. R., Macpherson, C. G. & Sinton, J. M. The CO₂-He-Ar-H₂O systematics of the Manus back-arc basin: Resolving source composition from degassing and contamination effects. *Geochimica et Cosmochimica Acta* 68(8), 1837–1855 (2004).

9. Yamamoto, J. et al. Diffusive fractionation of noble gases in mantle with magma channels: origin of low He/Ar in mantle-derived rocks. *Earth and Planetary Science Letters* 280(1-4), 167–174 (2009).

10. Boudoire, G., Rizzo, A. L., Di Muro, A., Grassa, F. & Liuzzo, M. Extensive CO₂ degassing in the upper mantle beneath oceanic basaltic volcanoes: First insights from Piton de la Fournaise volcano (La Réunion Island). *Geochimica et Cosmochimica Acta* 235, 376–401 (2018).

11. Rizzo, A. L. et al. Geochemistry of noble gases and CO₂ in fluid inclusions from lithospheric mantle beneath Wilcza Góra (Lower Silesia, southwest Poland). *Frontiers in Earth Science*, https://doi.org/10.3389/feart.2018.00215 (2018).

12. Kurz, M. D., Kenna, T. C., Lassiter, J. C. & DePaolo, D. J. Helium isotopic evolution of Mauna Kea volcano: First results from the 1-km drill core. *Journal of Geophysical Research: Solid Earth* 101(B5), 11781–11791 (1996).

13. DePaolo, D. J., Bryce, J. G., Dodson, A., Shuster, D. L. & Kennedy, B. M. Isotopic evolution of Mauna Loa and the chemical structure of the Hawaiian plume. *Geochemistry, Geophysics, Geosystems* 2(7) (2001).

14. Mukhopadhyay, S., Lassiter, J. C., Farley, K. A. & Bogue, S. W. Geochemistry of Kauai shield-stage lavas: Implications for the chemical evolution of the Hawaiian plume. *Geochemistry, Geophysics, Geosystems* 4(1) (2003).

15. Nuccio, P. M., Paonita, A., Rizzo, A. & Rosciglione, A. Elemental and isotope covariation of noble gases in mineral phases from Etnean volcanics erupted during 2001–2005, and genetic relation with peripheral gas discharges. *Earth and Planetary Science Letters* 272(3), 683–690 (2008).

16. Martelli, M. et al. Noble-gas signature of magmas from a heterogeneous mantle wedge: The case of Stromboli volcano (Aeolian Islands, Italy). *Chemical Geology* 368, 39–53 (2014).

17. Correale, A. et al. A two-component mantle source feeding Mt. Etna magmatism: Insights from the geochemistry of primitive magmas. *Lithos* 184, 243–258 (2014).

18. Rizzo, A. L. et al. New mafic magma refilling a quiescent volcano: Evidence from He-Ne-Ar isotopes during the 2011–2012 unrest at Santorini, Greece. *Geochemistry, Geophysics, Geosystems* 16(3), 798–814 (2015a).

19. Caracausi, A. et al. Mount Etna: Geochemical signals of magma ascent and unusually extensive plumbing system. *Geophysical Research Letters* 30(2) (2003).

20. Rizzo, A. L. et al. The 2014 effusive eruption at Stromboli volcano (Italy): inferences from soil CO₂ flux and ³He/⁴He ratio in thermal waters. *Geophysical Research Letters* 42(7), 2235–2243 (2015b).
21. Rizzo, A. L. et al. Eruptive activity at Turrialba volcano (Costa Rica): Inferences from 3He/4He in fumarole gases and chemistry of the products ejected during 2014 and 2015. *Geochim. Geophys. Geosyst.* **17**, 4478–4494 (2016a).

22. Padrón, E. et al. Diffusive helium emissions as a precursory sign of volcanic unrest. *Geology* **41**(5), 539–542 (2013).

23. Sano, Y. et al. Ten-year helium anomaly prior to the 2014 Mt Ontake eruption. *Scientific reports* **5**, 13069 (2015).

24. Paonita, A., Caracausi, A., Martelli, M. & Rizzo, A. L. Temporal variations of helium isotopes in volcanic gases quantify pre-eruptive refill and pressurization in magma reservoirs: The Mount Etna case. *Geology* **44**(7), 499–502 (2016).

25. McMurtry, G. M., Dasilveira, L. A., Horn, E. L., DeLuze, J. R. & Blessing, J. E. High 3He/4He ratios in lower East Rift Zone steaming vents precede a new phase of Kilauea 2018 eruption by 8 months. *Scientific Reports* **9**, 11860 (2019).

26. Capasso, G., Carapezza, M. L., Federico, C., Inguaggiato, S. & Rizzo, A. Geochemical monitoring of the 2002–2003 eruption at Stromboli volcano (Italy): precursory changes in the carbon and helium isotopic composition of fumarole gases and thermal waters. *Bulletin of volcanology* **68**(2), 118–134 (2005).

27. Marty, B., Meynier, V., Nicolini, E., Griesshaber, E. & Toutain, J. P. Geochemistry of gas emanations: a case study of the Réunion Hot Spot, Indian Ocean. *Applied geochemistry* **8**(2), 141–152 (1993).

28. Albarède, F. et al. The geochemical regimes of Piton de la Fournaise volcano (Réunion) during the last 530 000 years. *Journal of Petrology* **38**(2), 171–201 (1997).

29. Roult, G. et al. A new comprehensive classification of the Piton de la Fournaise activity spanning the 1985–2010 period. *Search and analysis of short-term precursors from a broad-band seismological station. Journal of Volcanology and Geothermal Research* **241**, 78–104 (2012).

30. Michon, L. et al. Explosive activity of the summit cone of Piton de la Fournaise volcano (La Réunion island): a historical and geological review. *Journal of Volcanology and Geothermal Research* **264**, 117–133 (2013).

31. Kaneoka, I., Takaoka, N. & Upton, B. G. J. Noble gas systematics in basalts and a dunite nodule from Reunion and Grand Comore Islands, Indian Ocean. *Chemical Geology: Isotope Geoscience section* **59**, 35–42 (1986).

32. Valer, M., Schiano, P. & Bachèlery, P. Geochemical characteristics of the La Réunion mantle plume source inferred from olivine-hosted melt inclusions from the adventive cones of Piton de la Fournaise volcano (La Réunion Island). *Contributions to Mineralogy and Petrology* **172**(9), 74 (2017).

33. Boudoire, G. et al. Eruptive activity on the western flank of Piton de la Fournaise (La Réunion Island, Indian Ocean): insights on magma transfer, storage and evolution at an oceanic volcanic island. *Journal of Petrology*, https://doi.org/10.1093/petrology/egz04 (2019).

34. Di Muro, A. et al. Magma degassing at Piton de la Fournaise volcano. In *Active Volcanoes of the Southwest Indian Ocean* (Springer, Berlin, Heidelberg, 2016).

35. Ort, M. H., Di Muro, A., Michon, L. & Bachèlery, P. Explosive eruptions from the interaction of magmatic and hydrothermal systems during flank extension: the Bellecombe Tephra of Piton de La Fournaise (La Réunion Island). *Bulletin of Volcanology* **78**(1), 5 (2016).

36. Staudacher, T., Kurz, M. D. & Allègre, C. J. New noble-gas data on glass samples from Loihi Seamount and Hualalai and on dunite samples from Loihi and Reunion Island. *Chemical Geology* **56**(3–4), 193–205 (1986).

37. Graham, D., Lupton, J., Albarède, F. & Condomines, M. Extreme temporal homogeneity of helium isotopes at Piton de la Fournaise, Réunion Island. *Nature* **347**(6293), 545–548 (1990).

38. Protin, M., Blard, P. H., Marrocchi, Y. & Mathon, F. Irreversible adsorption of atmospheric helium on olivine: A lobster pot analogy. *Geochimica et Cosmochimica Acta* **179**, 76–88 (2016).

39. Heber, V. S., Brooker, R. A., Kelley, S. P. & Wood, B. J. Crystal-melt partitioning of noble gases (helium, neon, argon, krypton, and xenon) for olivine and clinopyroxene. *Geochimica et Cosmochimica Acta* **71**(4), 1041–1061 (2007).

40. Michon, L., Ferrazzini, V., Di Muro, A., Villeneuve, N. & Famin, V. Rift zones and magma plumbing system of Piton de la Fournaise volcano: How do they differ from Hawaii and Etna? *Journal of Volcanology and Geothermal Research* **303**, 112–129 (2015).

41. Boudoire, G. Architecture et dynamique des systèmes magmatiques associés aux volcans basaltiques: exemple du Piton de la Fournaise. Doctoral dissertation, La Réunion, 500pp. (2017).

42. Peters, B. J., Carlson, R. W., Day, J. M. & Horan, M. F. Hadean silicate differentiation preserved by anomalous 142 Nd/144 Nd ratios in the Réunion hotspot source. *Nature* **555**(7694), 89 (2018).

43. Vlastélic, I. et al. Control of source fertility on the eruptive activity of Piton de la Fournaise volcano, La Réunion. *Scientific reports* **8**(1), 14478 (2018).

44. Burnard, P., Graham, D. & Farley, K. Fractionation of noble gases (He, Ar) during MORB mantle melting: a case study on the Southeast Indian Ridge. *Earth and Planetary Science Letters* **227**(3), 457–472 (2004).

45. Füri, E. et al. Helium isotope variations between Réunion Island and the Central Indian Ridge (17–21 S): New evidence for ridge–hot spot interaction. *Journal of Geophysical Research: Solid Earth* **116**(B2) (2011).

46. Trull, T. W. & Kurz, M. D. Experimental measurements of 3He and 4He mobility in olivine and clinopyroxene at magmatic temperatures. *Geochimica et Cosmochimica Acta* **57**(6), 1313–1324 (1993).

47. Sigmarsson, O., Condomines, M. & Bachèlery, P. Magma residence time beneath the Piton de la Fournaise Volcano, Reunion Island, from U-series disequilibria. *Earth and Planetary Science Letters* **234**(1), 223–234 (2005).

48. Kerr, A. C., Pearson, D. G. & Nowell, G. M. Magma source evolution beneath the Caribbean oceanic plateau: New insights from elemental and Sr-Nd-Pb-Hf isotopic studies of ODP Leg 165 Site 1001 basalts. *Geological Society, London, Special Publications* **328**(1), 809–827 (2009).

49. Lénat, J. F. Construction of La Réunion. In *Active Volcanoes of the Southwest Indian Ocean* (Springer, Berlin, Heidelberg, 2016).

50. Hanyu, T. et al. Noble gas study of the Reunion hotspot: evidence for distinct less-degassed mantle sources. *Earth and Planetary Science Letters* **193**(1), 83–98 (2001).

51. Keller, R. A., Graham, D. W., Farley, K. A., Duncan, R. A. & Lupton, J. E. Cretaceous-to-recent record of elevated 3He/4He along the Hawaiian-Emperor volcanic chain. *Geochemistry, Geophysics, Geosystems* **5**(12) (2004).

52. Broadley, M. W., Sumino, H., Graham, D. W., Burgess, R. & Ballentine, C. J. Recycled components in mantle plumes deduced from variations in halogens (Cl, Br, and I), trace elements, and 3He/4He along the Hawaiian-Emperor seamount chain. *Geochemistry, Geophysics, Geosystems* **20**(1), 277–294 (2019).

53. Staudacher, T., Sarda, P. & Allègre, C. J. Noble gas systematics of Réunion island, Indian Ocean. *Chemical Geology* **89**(1–2), 1–17 (1990).

54. Liuzzo, M. et al. New evidences of CO₂ soil degassing anomalies on Piton de la Fournaise volcano and the link with volcano tectonics structures. *Geochemistry, Geophysics, Geosystems* **5**, 2002GC000502 (2015).

55. Michon, L., Villeneuve, N., Catry, T. & Merle, O. How summit calderas collapse on basaltic volcanoes: New insights from the April 2007 caldera collapse of Piton de la Fournaise volcano. *Journal of Volcanology and Geothermal Research* **184**(1–2), 138–151 (2009).

56. Staudacher, T., et al. Fifteen years of intense eruptive activity (1998–2013) at Piton de la Fournaise volcano: a review. In *Active Volcanoes of the Southwest Indian Ocean* (Springer, Berlin, Heidelberg, 2016).

57. Villeneuve, N. & Bachèlery, P. Revue de la typologie des éruptions au Piton de La Fournaise, processus et risques volcaniques associés. *Cybergeo: European Journal of Geography* (2006).

58. Coppola, D. et al. Shallow system rejuvenation and magma discharge trends at Piton de la Fournaise volcano (La Réunion Island). *Earth and Planetary Science Letters* **463**, 13–24 (2017).

59. Boudoire, G. et al. Small-scale spatial variability of soil CO₂ flux: Implication for monitoring strategy. *Journal of Volcanology and Geothermal Research* **366**, 13–26 (2018b).
60. Gurioli, L. et al. Integrating field, textural, and geochemical monitoring to track eruption triggers and dynamics: a case study from Piton de la Fournaise. *Solid Earth* 9(2), 431–455 (2018).

61. Boudoire, G. et al. New perspectives on volcano monitoring in a tropical environment: continuous measurements of soil CO₂ flux at Piton de la Fournaise (La Réunion Island, France). *Geophysical Research Letters* 44(16), 8244–8253 (2017).

62. Hart, S. R. A large-scale isotope anomaly in the Southern Hemisphere mantle. *Nature* 309(5971), 753 (1984).

63. Fisk, M. R., Upton, B. G. J., Ford, C. E. & White, W. M. Geochemical and experimental study of the genesis of magmas of Reunion Island, Indian Ocean. *Journal of Geophysical Research: Solid Earth* 93(B5), 4933–4950 (1988).

64. Bosch, D. et al. Pb, Hf and Nd isotope compositions of the two Réunion volcanoes (Indian Ocean): A tale of two small-scale mantle “blobs”? *Earth and Planetary Science Letters* 265(3), 748–765 (2008).

65. Pietruszka, A. J., Hauri, E. H. & Blichtert-Toft, J. Crustal contamination of mantle-derived magmas within Piton de la Fournaise Volcano, Réunion Island. *Journal of Petrology* 50(4), 661–684 (2009).

66. Di Muro, A. et al. The shallow plumbing system of Piton de la Fournaise Volcano (La Reunion Island, Indian Ocean) revealed by the major 2007 caldera-forming eruption. *Journal of Petrology* 55(7), 1287–1315 (2014).

67. Nauret, F., Famin, V., Vlastélic, I. & Gannoun, A. A trace of recycled continental crust in the Réunion hotspot. *Chemical Geology* 524, 67–76 (2019).

68. Li, M., McNamara, A. K. & Garnero, E. J. Chemical complexity of hotspots caused by cycling oceanic crust through mantle reservoirs. *Nature Geosciences* 7, 366–370 (2014).

69. Williams, C. D., Li, M., McNamara, A. K., Garnero, E. J. & van Soest, M. C. Episodic entraiment of deep primordial mantle material into ocean island basalts. *Nature Communications* 6, 8937, https://doi.org/10.1038/ncomms9937 (2015).

70. Harris A. & Villeneuve N. Newspaper reporting of the April 2007 eruption of Piton de la Fournaise, part 2: framing the hazard. *Journal of Applied Volcanology* 7, https://doi.org/10.1186/s13617-018-0072-2 (2018).

71. Arienzo, I. et al. Sr and Nd isotope analysis at the Radiogenic Isotope Laboratory of the Istituto Nazionale di Geofisica e Vulcanologia, Sezione di Napoli—Osservatorio Vesuviano. *Rapporti Tecnici INGV* 260, 1–18 (2013).

72. Kurz, M. K. Cosmogenic helium in a terrestrial igneous rock. *Nature* 320, 435–439 (1986).

73. Hilton, D. R., Fischer, T. P. & Marty, B. Noble Gases and Volatile Recycling at Subduction Zones. Rev. Mineral. *Geochemistry* 47, 319–37, https://doi.org/10.2138/rmg.2002.47.9 (2002).

74. Hilton, D. R., Hammerschmidt, K., Teufel, S. & Friedrichsen, H. Helium isotope characteristics of Andean geothermal fluids and lavas. *Earth Planet. Sci. Lett.* 120, 265–282, https://doi.org/10.1016/0012-821X(93)90244-4 (1993).

75. Faccini, B. et al. Subduction-related melt refertilisation and alkaline metasomatism in the Eastern Transylvanian Basin lithospheric mantle: Evidence from mineral chemistry and noble gases in fluid inclusions. Lithos, 464–465, https://doi.org/10.1016/j.lithos.2020.105516 (2020).

76. Mahoney, J. J. et al. Isotopic and geochemical provinces of the western Indian Ocean spreading centers. *Journal of Geophysical Research: Solid Earth* 94(B4), 4033–4052 (1989).

77. Rizzo, A. L. et al. Kolumbo submarine volcano (Greece): An active window into the Aegean subduction system. *Sci. Rep.* 6(28013) (2016b).

78. Sano, Y., Wakita, H., Ohsumi, T. & Kusakabe, M. Helium isotope evidence for magmatic gases in Lake Nyos, Cameroon. *Geophysical Research Letters* 14(10), 1039–1041 (1987).

79. Jambon, A., Weber, H. & Braun, O. Solubility of He, Ne, Ar, Kr and Xe in a basalt melt in the range 1250–1600 C. Geochemical implications. *Geochimica et Cosmochimica Acta* 50(3), 401–408 (1986).

80. Crank, J. The Mathematics of Diffusion: 2d Ed. (Clarendon Press., 1975).

81. Graham, D. W., Jenkins, W. J., Kurz, M. D. & Batiza, R. Helium isotope disequilibrium and geochronology of glassy submarine basalts. *Nature* 326(6111), 384 (1987).

82. Vlastélic, I., Lewin, E. & Staudacher, T. Th/U and other geochemical evidence for the Reunion plume sampling a less differentiated mantle domain. *Earth and Planetary Science Letters* 248(1–2), 379–393 (2006).

83. Benaglia, T., Chauveau, D., Hunter, D. & Young, D. mixtools: An R package for analyzing finite mixture models. *Journal of Statistical Software* 32(6), 1–29 (2009).

## Acknowledgements

M. Moreira and T. Staudacher are gratefully acknowledged for providing samples (ILR84-4) in order to test the accuracy of our analysis with respect to previous works in literature. We are also in debt with (i) Fausto Grassa and Marco Liuzzo for their help to Cilaos gases, (ii) Chiara Alessi and Mariagrazia Misseri for samples preparation as well as (iii) Mariano Tantillo and Mariagrazia Misseri for their assistance during noble gases isotope analysis. We thank the Editor and three anonymous Reviewers for their suggestions that greatly improved the clarity and quality of the paper. We are in debt with Andrew Harris for English suggestions. This research received funds from the Italian Ministero Istruzione Università e Ricerca (Project PRIN 2017LMNLAW) and the INGV of Palermo. We also acknowledge the French government IDEX-ISITE initiative 16-IDEX-0001 (CAP 20–25). This is Laboratory of Excellence ClerVolc contribution number 413.

## Author contributions

G.B. and A.L.R. conceptualized the project. G.B. and A.D.M. provided the samples. A.L.R. performed the noble gases analysis. I.A. performed the strontium isotopes analysis. G.B. and A.L.R. developed the interpretative model and G.B. wrote the paper with assistance of A.L.R.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-66260-x.

Correspondence and requests for materials should be addressed to G.B.

Reprints and permissions information is available at www.nature.com/reprints.

Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
<figure><img src="image_10.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020