# What factors control superficial lava dome explosivity?

Georges Boudon, Hélène Balcone-Boissard, Benoît Villemant, Daniel J. Morgan

► To cite this version:

Georges Boudon, Hélène Balcone-Boissard, Benoît Villemant, Daniel J. Morgan. What factors control superficial lava dome explosivity?. Scientific Reports, 2015, 5, pp.14551. 10.1038/srep14551. hal-01217999

HAL Id: hal-01217999

https://hal.sorbonne-universite.fr/hal-01217999

Submitted on 20 Oct 2015

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
OPEN

# What factors control superficial lava dome explosivity?

Georges Boudon<fnref n="1" />, Hélène Balcone-Boissard<fnref n="2" />, Benoît Villemant<fnref n="2" /> & Daniel J. Morgan<fnref n="3" />

Received: 13 January 2015

Accepted: 10 August 2015

Published: 30 September 2015

Dome-forming eruption is a frequent eruptive style and a major hazard on numerous volcanoes worldwide. Lava domes are built by slow extrusion of degassed, viscous magma and may be destroyed by gravitational collapse or explosion. The triggering of lava dome explosions is poorly understood: here we propose a new model of superficial lava-dome explosivity based upon a textural and geochemical study (vesicularity, microcrystallinity, cristobalite distribution, residual water contents, crystal transit times) of clasts produced by key eruptions. Superficial explosion of a growing lava dome may be promoted through porosity reduction caused by both vesicle flattening due to gas escape and syn-eruptive cristobalite precipitation. Both processes generate an impermeable and rigid carapace allowing overpressurisation of the inner parts of the lava dome by the rapid input of vesiculated magma batches. The relative thickness of the cristobalite-rich carapace is an inverse function of the external lava dome surface area. Explosive activity is thus more likely to occur at the onset of lava dome extrusion, in agreement with observations, as the likelihood of superficial lava dome explosions depends inversely on lava dome volume. This new result is of interest for the whole volcanological community and for risk management.

Dome-forming eruption is a frequent eruptive style<fnref n="1" />; lava domes result from continuous or episodic slow extrusion of viscous lava. Most dome-forming eruptions produce highly microcrystallized and highly- to almost totally-degassed magmas which have a low explosive potential. During lava dome growth, recurrent collapses of unstable parts are the main destructive process of the lava dome, generating block- and ash-flows (concentrated pyroclastic flows) hereafter referred to as "concentrated pyroclastic density currents" (C-PDC)<fnref n="2" />. C-PDC's are channelized in valleys<fnref n="3" />⁻⁵ and have a high, but localized, damage potential that largely depends on the collapsed volume⁶. Sometimes, a dilute ash cloud surge develops at the top of the concentrated flow with an increased destructive effect, because it may overflow ridges and affect larger areas⁵. In addition, large lava dome collapses can induce a depressurization of volatile-bearing magma within the conduit, leading to vulcanian explosions<fnref n="3" />,⁵. By contrast, violent, laterally-directed explosions may occur at the base of a growing lava dome: these generate dilute and turbulent pyroclastic flows, or surges (péléean 'nuées ardentes'⁷) hereafter referred to as "dilute pyroclastic density currents" (D-PDC)<fnref n="2" />. They have a high velocity and their propagation is poorly dependent on the topography, leading to highly destructive effects⁷,⁸. Numerous studies on lava dome behaviors exist⁶,⁹⁻¹³, but the triggering of lava dome explosions—an important topic in hazard mitigation and risk assessment¹⁴—is poorly understood, which this article hopes to address.

A series of dome-forming eruptions are investigated through the petrology of their volcanic products. At Montagne Pelée (Martinique), during the first months of the 1902–1905 eruption, violent, superficial, laterally-directed explosions at the base of the growing lava dome generated D-PDC's which devastated the south-western flank of the volcano⁷,⁸. The May, 8th 1902 event, which occurred two days after the beginning of lava dome growth, was responsible for the death of 28 000 inhabitants. It was the first of a series of 7 eruptions, ending on August, 30th 1902. Following this explosive phase, the 1902–1905 activity

<footnote n="1">¹Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Univ. Paris Diderot, CNRS, F-75005, Paris, France.</footnote>

<footnote n="2">²Sorbonne Universités, UPMC Univ Paris 06, CNRS, Institut des Sciences de la Terre de Paris (iSTeP), 4 place Jussieu 75005 Paris, France.</footnote>

<footnote n="3">³Institute of Geophysics and Tectonics, School of Earth & Environment, University of Leeds, Leeds, LS2 9JT UK. Correspondence and requests for materials should be addressed to G.B. (email: boudon@ipgp.fr)</footnote>
moved toward a lava dome collapse style producing C-PDC's⁷. This eruptive style persisted throughout the 1929–1932 eruption¹⁵. The P1 eruption (650 y. BP), also started by the growth of a lava dome that was destroyed, as for the 1902 eruption, by two violent, superficial, laterally-directed explosions initiated at the base of the growing lava dome¹⁶. They generated D-PDC's which devastated the south-western flank of the volcano. This first phase was followed by a plinian phase, generating pumice fallout that covers the volcano flanks. The May, 8th 1902 D-PDC, one of the D-PDC from the P1 eruption and two C-PDC deposits from the 1929-1932 eruption as well as the two 1902–1905 and 1929–1932 lava domes were sampled.

During the on-going eruption of Soufrière Hills (Montserrat), a lava dome has been permanently present since 1995, alternating periods of quiescence with periods of lava dome growth and recurrent destruction⁵. Destruction occurs by collapses of variable volumes with attendant vulcanian explosions, together generating C-PDC⁵. This eruption destroyed the Castle Peak lava dome previously emplaced (1650 AD;¹⁷). The Castle Peak lava dome and a series of lava domes of the 1996–1997 period were sampled with their associated C-PDC deposits.

At Santa Maria-Santiaguito volcano (Guatemala), a plinian eruption occurred in 1902 and was followed 20 years later by a dome-forming eruption that is still on-going. A series of lava domes were built and produced more-or-less voluminous C-PDC channeled in the main valleys of the south-western flank of the volcano¹⁸,¹⁹. Different lava domes (Caliente: 1922–1925, El Monje: 1950's, El Brujo: 1970's), lava flows (1971–1972 and 1973–1975) and the C-PDC deposit from the large Caliente lava dome collapse of 1929 were sampled.

In La Chaîne des Puys (Massif Central, France), the Puy de Dôme eruption (~10,800 y. BP)²⁰ only produced C-PDCs whereas the Puy Chopine eruption (~9,700 y. BP)²⁰ generated a D-PDC during its first phase of activity, destroying the emerging lava dome. The lava dome and the C-PDC deposit from Puy de Dôme and the D-PDC deposit from Puy Chopine were sampled.

Finally, a sample of the lava dome from the 1991–1994 dome forming eruption of Unzen (Japan) is included in this study for comparison.

Petrological data from the plinian phase of the P1 eruption (Montagne Pelée), the 1902 plinian eruption of Santa Maria²¹,²² and the 1997 vulcanian explosions of Soufrière Hills²³ are also included for comparison of plinian or vulcanian eruptions with dome-forming eruptions at the same volcanoes during the active periods studied.

## Results

**Magma composition.** Magmas of Montagne Pelée and Montserrat are acid andesites (60–62 and 58–60 wt% SiO₂ respectively), those of Santiaguito are dacites (65–67 wt% SiO₂) and those of La Chaîne des Puys are trachytes (66–69 wt% SiO₂). All residual glasses are rhyolitic (70–79 wt% SiO₂).

**Textural characteristics: density, vesicularity, cristobalite content.** C-PDC clasts have a narrow and unimodal distribution of vesicularities (20–40%), whereas D-PDC clasts display a much larger range (10–75%) (Fig. 1). The most vesiculated clasts (vesicularity >50%) of D-PDC display characteristic pumiceous textures with subspherical and disconnected vesicles and only rare microlites (Fig. 2a,d). With decreasing vesicularity, vesicle number and size decrease, vesicle shapes become irregular with large vesicles concentrating in some areas, and both groundmass/vesicles and microlite/glass ratios increase (Fig. 2b,c,e). The less-vesiculated clasts show two types of textures: (i) in most clasts, vesicles are rare and crystalline silica precipitates abundant (Fig. 2c,f; Table 1). Crystalline silica exists as cristobalite, as identified by Raman spectrometry (Fig. 3a and supplementary material) and occurs either as cracked infillings in large vesicles (up to 30 µm in diameter) or as a pervasive form in small vesicles (down to 1 µm in diameter) (Figs 2c,f and 3b). The weight fraction of cristobalite in D-PDC clasts decreases with increasing vesicularity (Fig. 4a). No cristobalite is observed in the most vesiculated clasts. (ii) in a few clasts, a texture of sparse, small and irregular vesicles that are widely separated exists; these clasts are cristobalite-free.

Samples from the lava domes and from the C-PDC are poorly- to non-vesiculated (vesicularity < 50%), with vesicle patterns similar to those of the less-vesiculated clasts from the D-PDC. In some samples, coalescence phenomena create numerous irregular channels, sometimes several millimeters long and tens of micrometers wide. Cristobalite is present in all samples collected on lava domes in relatively high proportions (up to 28 area%; Fig. 4a; Table 1). Most C-PDC clasts contain cristobalite (up to 23 area%; Table 1) but its fraction is independent of the vesicularity (Fig. 4a).

**Residual water content.** Samples from the lava domes collected as clasts within the C-PDC deposits show low residual H₂O contents (H₂O_r < 0.6 wt%, bulk rock content corrected for phenocryst content; Table 1). In contrast, clasts collected from D-PDC deposits display larger H₂O_r ranges: 0–1.9 wt% at Montagne Pelée and 0–2.5 wt% at Puy Chopine (Fig. 4b; Table 1). Water content measured on whole-rock samples is considered as pure magmatic water (see methodology).

**Crystal transit time through diffusion profiles in magnetites.** Ti diffusion profiles in magnetite crystals from Montagne Pelée (May, 8th 1902 event) have been studied following the method developed on pyroxenes by Morgan and collaborators²⁴. This eruptive event has been chosen to estimate the transit
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Bulk rock vesicularity distribution (frequency %) of samples from Pyroclastic Density Current (PDC) deposits, as a function of vesicularity (in volume %). The distribution is based on at least 100 clasts for each deposit in the fraction 16–32 mm. The left column refers to D-PDCs (Montagne Pelée: May 8th, 1902, lava dome phase of P1 eruption; Puy Chopine). The right column refers to C-PDCs (Montagne Pelée, 2 C-PDCs from the 1929–1932 eruption; Puy de Dôme; Santiaguito, 1929). The C-PDC vesicularity distribution is unimodal for the four examples (centered around 27% and 22% for Montagne Pelée, 32%, for Puy de Dôme, and 17% for Santiaguito). The D-PDC vesicularity distribution is multimodal. In D-PDC deposits highly vesiculated clasts (>40%) are much less abundant in May 8th, 1902 of Montagne Pelée (7%) than in P1 eruption (15%) and Puy Chopine (30%).</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** BSE images of representative clasts from Montagne Pelée. Images (a-c): vesiculated to dense clasts from May, 8th, 1902 D-PDC deposit. (a) the most vesiculated clasts (bulk rock vesicularity >50%) of D-PDC display characteristic pumiceous textures with subspherical and disconnected vesicles, 10 to 100 µm in diameter, with thin bubble walls (<5 µm) and only rare microlites in glass. (b) intermediate clasts have decreasing bulk rock vesicularity, increasing microlite content and bubble wall thickness. (c) the densest clasts display heterogeneous textures: most clasts display rare, irregular and small (<20 µm) vesicles with 12 area % cristobalite precipitates. Cristobalite corresponds to the dark grey phase. (d) glassy vesiculated clast from May, 8th, 1902 D-PDC. (e) highly microcrystalline dense clast from 1929 C-PDC. (f) cristobalite precipitates - dark grey phase (dense clasts, May, 8th, 1902 D-PDC): these are found either as cracked infillings in some large vesicles (up to 30 µm) or as pervasive patches (few µm to 200 µm in size).</figcaption></figure>

time of magma in the conduit. The textural characteristics of the magnetite crystals depend on the vesic-
ularity of host clasts (Fig. 5a,d). In vesiculated clasts, magnetites display a normal zoning with Ti-rich
cores and Ti-poor rims that are in equilibrium with the residual melt (Fig. 5a,b). In dense clasts, two
types of textures are evidenced. In some dense, degassed clasts, only few magnetites exhibit the same
textural features as vesiculated clasts (~15%); but most crystals are exsolved (Fig. 5d) indicating that
they have re-equilibrated at temperatures below the solvus, and likely in oxidizing conditions²⁵. In other
dense clasts all magnetites are exsolved, indicating a sufficiently long storage time under the appropriate
conditions to allow complete exsolution. From over 100 crystals of magnetite separated in vesiculated
fragments, 13 diffusion profiles were usable, while over more than 150 separates in dense fragments only
three profiles were usable.

**Discussion**

Cristobalite is present in clasts from the PDC deposits in highly variable proportions and is heterogene-
ously distributed, which indicates that precipitation necessarily occurred before the collapse or explosion
of the growing lava dome and that cristobalite is not of post-eruptive origin, via hydrothermal circula-
tions or glass devitrification (Fig. 3 and Supplementary material). Cristobalite fills round-shaped voids of
any size, characteristic of former gas vesicles (Figs 2f and 4b). The clasts from D-PDC deposits display
a large range of textural features from highly degassed, poorly vesiculated and cristobalite-rich to highly
vesiculated, poorly degassed and cristobalite-free. Since the D-PDC's result from explosions at the base of
the growing lava dome<sup>7,8</sup>, clasts embody a representative sampling of the different zones of the destroyed
lava dome, showing a heterogeneous distribution of cristobalite. As cristobalite precipitation is mainly
controlled by pressure and temperature decrease<sup>26-28</sup>, it suggests that degassing at atmospheric pressure
and low temperature (lower at the lava-dome surface than the 400 °C below which cristobalite precipi-
tates)<sup>28,29</sup> may have produced an outer silicified carapace. The low vesicularity of the outer parts of a lava
dome may be due to a combination of three processes during lava dome growth: (i) connection, opening
and consequent extreme flattening of the vesicles due to the degassing in open-system during magma
ascent at shallow depth in the conduit<sup>9,30</sup>, (ii) rapid cooling of the magma at the surface, thus restricting
bubble growth<sup>31</sup>, (iii) cristobalite precipitation filling vesicles at low pressure<sup>27</sup>. In addition, slow degas-
sing at shallow depth induces melt crystallization, favoring the development of a rigid carapace delimited
by the brittle-ductile transition controlled by melt micro-crystallinity<sup>32</sup>. These processes progressively
<table><thead><tr><th></th><th>d measured kg.m<sup>-3</sup></th><th>Ves WR %</th><th>Ves glass %</th><th>H<sub>2</sub>O WR wt%</th><th>H<sub>2</sub>O r wt%</th><th>Silica area %</th><th></th><th>d measured kg.m<sup>-3</sup></th><th>Ves WR %</th><th>Ves glass %</th><th>H<sub>2</sub>O WR wt%</th><th>H<sub>2</sub>O r wt%</th><th>Silica area %</th></tr></thead><tbody><tr><td colspan="14">Montagne Pelée, Martinique</td></tr><tr><td colspan="14">P1, D-PDC</td></tr><tr><td>ME 1311fy 23</td><td>0,88</td><td>67</td><td>80</td><td>1,03</td><td>2,29</td><td>0</td><td>1902, Lava dome</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ME 1311fy 9</td><td>1,46</td><td>45</td><td>63</td><td>0,9</td><td>2,00</td><td>0</td><td>MH 804a</td><td>2,30</td><td>13</td><td>17</td><td>0,05</td><td>0,08</td><td>18</td></tr><tr><td>ME 1311fy 2</td><td>2,1</td><td>23</td><td>35</td><td></td><td></td><td>4</td><td>MG 804</td><td>2,37</td><td>11</td><td>13</td><td>0,05</td><td>0,08</td><td>23</td></tr><tr><td>ME 1311fy 29</td><td>2,3</td><td>18</td><td>28</td><td></td><td></td><td>0</td><td>MG 802</td><td>2,49</td><td>6</td><td>7</td><td>0,03</td><td>0,05</td><td>24</td></tr><tr><td>ME 1311fy 20</td><td>2,32</td><td>13</td><td>26</td><td>0,62</td><td>1,38</td><td>0</td><td>May 8<sup>th</sup> 1902, D-PDC</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ME 1311fy 30</td><td>2,3</td><td>15</td><td>24</td><td></td><td></td><td>0</td><td>ME 1314-1</td><td>2,31</td><td>13</td><td>21</td><td>0,21</td><td>0,39</td><td>21</td></tr><tr><td>ME 1311fy 14</td><td>2,4</td><td>13</td><td>21</td><td></td><td></td><td>0</td><td>ME 1314-2</td><td>1,86</td><td>30</td><td>42</td><td>0,21</td><td>0,39</td><td>0</td></tr><tr><td>ME 1311fy 6</td><td>2,4</td><td>12</td><td>19</td><td></td><td></td><td>8</td><td>ME 1314-3</td><td>2,28</td><td>14</td><td>22</td><td>0,03</td><td>0,05</td><td>12</td></tr><tr><td>ME 1311fy 18</td><td>2,5</td><td>11</td><td>17</td><td></td><td></td><td>0</td><td>ME 1314-4</td><td>1,64</td><td>38</td><td>52</td><td>0,49</td><td>0,91</td><td>0</td></tr><tr><td>ME 1311fy 32</td><td>2,5</td><td>10</td><td>16</td><td></td><td></td><td>10</td><td>ME 1314-5</td><td>1,05</td><td>60</td><td>72</td><td>0,68</td><td>1,27</td><td>0</td></tr><tr><td>ME 1311fy 10</td><td>2,6</td><td>7</td><td>11</td><td></td><td></td><td>10</td><td>ME 1314-6</td><td>1,30</td><td>51</td><td>64</td><td>0,62</td><td>1,16</td><td>0</td></tr><tr><td>ME 1311fx 6</td><td>2,4</td><td>12</td><td>20</td><td></td><td></td><td>0</td><td>ME 1314-7</td><td>2,27</td><td>14</td><td>23</td><td>0,28</td><td>0,52</td><td>0</td></tr><tr><td>ME 1311fx 25</td><td>2,4</td><td>14</td><td>22</td><td></td><td></td><td>0</td><td>ME 1314 y1 6</td><td>0,87</td><td>67</td><td>78</td><td>0,45</td><td>0,84</td><td>n.d.</td></tr><tr><td>ME 1311fx 30</td><td>0,97</td><td>64</td><td>78</td><td>1,11</td><td>2,47</td><td>0</td><td>ME 1314 y1 10</td><td>0,80</td><td>70</td><td>80</td><td>0,55</td><td>1,03</td><td>n.d.</td></tr><tr><td>ME 1311fx 10</td><td>1,79</td><td>32</td><td>50</td><td>0,81</td><td>1,80</td><td></td><td>ME 1314 y1 11</td><td>1,03</td><td>61</td><td>73</td><td>0,69</td><td>1,29</td><td>n.d.</td></tr><tr><td>ME 1311fx 8</td><td>1,93</td><td>27</td><td>45</td><td>0,81</td><td>1,80</td><td></td><td>ME 1314 y1 12</td><td>1,27</td><td>52</td><td>65</td><td>1,215</td><td>2,28</td><td>n.d.</td></tr><tr><td>ME 1311fx 12</td><td>2,35</td><td>15</td><td>24</td><td></td><td></td><td>3</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>ME 1311fx 21</td><td>2,44</td><td>12</td><td>19</td><td></td><td></td><td>8</td><td>1929, Lava dome</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1311E</td><td>1,75</td><td>34</td><td>52</td><td>0,39</td><td>0,87</td><td></td><td>MG 901</td><td>2,42</td><td>9</td><td>10</td><td>0,06</td><td>0,10</td><td>n.d.</td></tr><tr><td>1311F</td><td></td><td></td><td></td><td>0,28</td><td>0,62</td><td></td><td>MG 803</td><td>2,47</td><td>7</td><td>8</td><td>0,10</td><td>0,16</td><td>n.d.</td></tr><tr><td>1311F1</td><td>2,31</td><td>13</td><td>26</td><td>0,38</td><td>0,84</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1311F3</td><td>1,61</td><td>39</td><td>57</td><td>0,55</td><td>1,22</td><td></td><td>1929, C-PDC</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>901F1</td><td>2,36</td><td>11</td><td>23</td><td>0,4</td><td>0,89</td><td></td><td>ME 1317-1</td><td>2,53</td><td>4</td><td>8</td><td>0,06</td><td>0,11</td><td>n.d.</td></tr><tr><td>901F2</td><td>1,97</td><td>26</td><td>43</td><td>0,51</td><td>1,13</td><td></td><td>ME 1317-2</td><td>2,04</td><td>23</td><td>33</td><td>0,26</td><td>0,47</td><td>16</td></tr><tr><td>901F3</td><td>2,20</td><td>17</td><td>32</td><td></td><td></td><td></td><td>ME 1317-3</td><td>1,65</td><td>38</td><td>50</td><td>0,31</td><td>0,56</td><td>23</td></tr><tr><td>901F4</td><td>2,18</td><td>18</td><td>33</td><td>0,49</td><td>1,09</td><td>5</td><td>ME 1317-4</td><td>2,50</td><td>6</td><td>9</td><td>0,31</td><td>0,56</td><td>23</td></tr><tr><td colspan="14">French Massif Central</td></tr><tr><td colspan="14">Puy Chopine</td></tr><tr><td colspan="14">D-PDC</td></tr><tr><td>PC- 139</td><td>0,92</td><td>64</td><td>68</td><td>1,39</td><td>1,81</td><td>0,0</td><td>Puy de Dôme</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 133</td><td>1,05</td><td>59</td><td>64</td><td>1,47</td><td>1,91</td><td>n.d.</td><td>C PDC</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 91</td><td>1,14</td><td>56</td><td>60</td><td>0,64</td><td>0,83</td><td>n.d.</td><td>PdD 56</td><td>1,74</td><td>35</td><td>35</td><td>0,28</td><td>0,35</td><td>n.d.</td></tr><tr><td>PC- 130</td><td>1,22</td><td>52</td><td>57</td><td>1,27</td><td>1,65</td><td>n.d.</td><td>PdD 84</td><td>2,05</td><td>23</td><td>22</td><td>0,14</td><td>0,18</td><td>n.d.</td></tr><tr><td>PC- 19</td><td>1,35</td><td>47</td><td>52</td><td>1,17</td><td>1,52</td><td>0,0</td><td>PdD 41</td><td>1,94</td><td>27</td><td>27</td><td>0,1</td><td>0,13</td><td>7</td></tr><tr><td>PC- 11</td><td>1,33</td><td>48</td><td>53</td><td>0,8</td><td>1,04</td><td>n.d.</td><td>PdD 10</td><td>2,26</td><td>15</td><td>12</td><td>0,06</td><td>0,08</td><td>9</td></tr><tr><td>PC- 6</td><td>1,41</td><td>45</td><td>49</td><td>1,33</td><td>1,73</td><td>n.d.</td><td>PdD 44</td><td>1,65</td><td>38</td><td>39</td><td>0,38</td><td>0,48</td><td>10</td></tr><tr><td>PC- 99</td><td>1,46</td><td>43</td><td>48</td><td>0,89</td><td>1,16</td><td>0,0</td><td>Lava dome</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 58</td><td>1,65</td><td>36</td><td>40</td><td>0,82</td><td>1,06</td><td>0,0</td><td>PdD</td><td>2,01</td><td>24</td><td>24</td><td>0,26</td><td>0,33</td><td>26</td></tr><tr><td>PC- 29</td><td>1,62</td><td>37</td><td>41</td><td>0,85</td><td>1,10</td><td>n.d.</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 50</td><td>1,97</td><td>23</td><td>26</td><td>0,38</td><td>0,49</td><td>0,0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 70</td><td>1,93</td><td>25</td><td>28</td><td>0,27</td><td>0,35</td><td>2,0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 68</td><td>2,07</td><td>19</td><td>22</td><td>0,14</td><td>0,18</td><td>12,0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>PC- 56</td><td>2,07</td><td>19</td><td>22</td><td>0,21</td><td>0,27</td><td>6,0</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td colspan="14">lava dome</td></tr><tr><td>PC</td><td>2,1</td><td>18</td><td></td><td></td><td></td><td>n.d.</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

Continued
<table><caption>**Table 1.** Textural and geochemical data, from Montagne Pelée (Martinique), Chaîne des Puys, Santiaguito (Guatemala) and Soufrière Hills (Montserrat). **d measured (kg.m⁻³)**: the density is obtained by the “3 weights” methods (weighing uncoated clasts in air, then paraffin coated clast in air and in water) for clasts spanning the 16–32 mm sieve size fractions. **Yes WR (%)**: the whole rock vesicularity is calculated from density measurements using the Dense Rock Equivalent density. DRE density is specified for each system below the data; the DRE density is obtained by pycnometry measurement on fine-grained whole-rock powder on a minimum of 3 samples for each system. **Yes glass (%)**: The vesicularity is calculated from density measurements using the Dense Rock Equivalent density and corrected from phenocrysts content to refer to the initial melt. The crystallinity for each system is reported below the data. Vesglass is calculated from Vg/Vl by the formula: Vesglass = (Vg/Vl)/(1+Vg/Vl). The last two parameters are best suited for tracing the effect of degassing, because they only depend on initial melt content which is the only source of exsolved fluid (H₂O) and do not depend on the initial crystallinity of the magma which is highly variable in the different studied eruptions. **H₂O WR (wt%):** obtained by H₂ manometry⁴⁶. **H₂O r (wt%):** residual H₂O content calculated from H₂O WR corrected from phenocryst content. **Silica (area %):** obtained by chemical mapping (2 to 5 per sample) at magnification 100x of Si. The mean relative standard deviation is <5% for crystalline silica contents.</caption><thead><tr><th></th><th>d measured kg.m<sup>-3</sup></th><th>Ves WR %</th><th>Ves glass %</th><th>H<sub>2</sub>O WR wt%</th><th>H<sub>2</sub>O r wt%</th><th>Silica area %</th><th></th><th>d measured kg.m<sup>-3</sup></th><th>Ves WR %</th><th>Ves glass %</th><th>H<sub>2</sub>O WR wt%</th><th>H<sub>2</sub>O r wt%</th><th>Silica area %</th></tr></thead><tbody><tr><td>Santiaguito (Guatemala)</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Soufriere Hills (Montserrat)</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lava dome</td><td></td><td></td><td></td><td></td><td></td><td></td><td>Lava dome</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Castle Peak</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SM10a</td><td>2,44</td><td>8</td><td>8</td><td>0,1</td><td>0,13</td><td>0</td><td>MO1</td><td>2,44</td><td>11</td><td>16</td><td>0,2</td><td>0,40</td><td>n.d.</td></tr><tr><td>SM11</td><td>2,05</td><td>23</td><td>27</td><td>0,2</td><td>0,25</td><td>n.d.</td><td>MO2</td><td>2,44</td><td>11</td><td>16</td><td>0,2</td><td>0,40</td><td>19</td></tr><tr><td>SM12</td><td>2,38</td><td>10</td><td>11</td><td>0,2</td><td>0,22</td><td>26</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SM13</td><td>2,37</td><td>11</td><td>11</td><td>0,1</td><td>0,16</td><td>28</td><td>1995-1996</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SM15</td><td>2,07</td><td>22</td><td>26</td><td>0,3</td><td>0,46</td><td>n.d.</td><td>MVO24</td><td>2,42</td><td>12</td><td>17</td><td>0,03</td><td>0,06</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>MVO25</td><td>2,42</td><td>12</td><td>17</td><td>0,06</td><td>0,12</td><td>19</td></tr><tr><td>1929, C-PDC</td><td></td><td></td><td></td><td></td><td></td><td></td><td>MVO26</td><td>2,47</td><td>10</td><td>14</td><td>0,12</td><td>0,24</td><td>n.d.</td></tr><tr><td>SM9 3c1</td><td>1,92</td><td>28</td><td>33</td><td>0,4</td><td>0,51</td><td>n.d.</td><td>MVO27</td><td>2,14</td><td>22</td><td>32</td><td>0,07</td><td>0,14</td><td>0</td></tr><tr><td>SM9 3c3</td><td>2,42</td><td>9</td><td>9</td><td>0,1</td><td>0,17</td><td>n.d.</td><td>182a</td><td>1,97</td><td>28</td><td>40</td><td>0,2</td><td>0,40</td><td>n.d.</td></tr><tr><td>SM9 3c7</td><td>2,26</td><td>15</td><td>17</td><td>0,1</td><td>0,11</td><td>n.d.</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>SM8</td><td>1,89</td><td>29</td><td>34</td><td>0,2</td><td>0,34</td><td>8,0</td><td>1995-1996, C-PDC</td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Mo96-2</td><td>2,42</td><td>12</td><td>17</td><td>0,12</td><td>0,24</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Mo96-3C3</td><td>2,41</td><td>12</td><td>17</td><td>0,32</td><td>0,64</td><td>5</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Mo96-3C6</td><td>2,39</td><td>13</td><td>19</td><td>0,16</td><td>0,32</td><td>6</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>96-02A RT</td><td>2,37</td><td>14</td><td>20</td><td>0,1</td><td>0,20</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>96-01A</td><td>2,18</td><td>21</td><td>30</td><td>0,07</td><td>0,14</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>96-11B</td><td>2,36</td><td>14</td><td>20</td><td>0,18</td><td>0,36</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>96-11C</td><td>2,63</td><td>4</td><td>4</td><td>0,25</td><td>0,50</td><td>n.d.</td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>96-11C-Cb</td><td>2,63</td><td>4</td><td>4</td><td>0,25</td><td>0,50</td><td>n.d.</td></tr></tbody></table>

increase the viscosity and reduce the porosity up to annihilation, creating at the outer surface of the lava dome an impermeable, highly-silicified, mechanically resistant carapace.

To acquire a significant strength, the carapace must attain a thickness of few meters to several tens of meters in the absence of other consolidating effects²⁹. We suggest that silicification may reduce the threshold thickness to the lower range. This impermeable and resistant carapace mechanically isolates the core of the lava dome, thus preventing further volatile exsolution and induced crystallization of the melt confined in this core. The D-PDC's clasts span a wide range of vesicularity and H₂Oₜ content, from almost completely degassed, microcrystalline clasts to undegassed, microlite-free clasts (~2.5 wt% of H₂Oₜ for a glass vesicularity of ~80%; Figs 2d and 4b; Tables 1 and 2). The vesiculated, H₂O-rich clasts display textural characteristics similar to those of plinian clasts³³, contrary to the low vesiculated, H₂O-poor clasts, which are highly microcrystalline and cristobalite-rich. Thus volatile-rich and vesiculated clasts may represent the magma stored within the core of the lava dome and that evolves in a closed system degassing, leading to overpressurization of the upper volcanic edifice (Fig. 6). By contrast, the less-vesiculated, H₂O-poor clasts may represent magma that evolved in an open-system, degassing in the external part of the lava dome. Overpressures of 0.1–1 MPa are sufficient to surpass the tensile strength
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Cristobalite in clasts from Montagne Pelée and Puy de Dôme. (a) Raman spectra of crystalline silica patches and of interstitial crystalline silica in dense clasts from Montagne Pelée (May 8th, 1902) and Puy de Dôme: both consist of cristobalite (Micro-Raman spectroscopy, Zeiss Supra V55; UPMC, Paris). Y-axis corresponds to counts minus background for the silica precipitates measurements (not to scale). Minor peaks result from the laser excitation of surrounding material (bubble walls) and correspond to the groundmass phases (Feldspars, Pyroxene and glass). Qualitative details are given in the supplementary material. (b) BSE image of a dense clast from May, 8th, 1902 D-PDC of Montagne Pelée with patches of cristobalite (P) and interstitial cristobalite (I).</figcaption></figure>

of the carapace and to trigger the explosive destruction of the upper part of the lava dome<sup>6,34,35</sup>. The pressure distribution within the edifice may be significantly modified by the existence of shear stress at the conduit vent<sup>12</sup>: the largest gas overpressure may be located at the conduit wall (horizontal spreading at the vent) or at the center of the conduit (zero horizontal velocity as in the case of a pre-existing lava dome)<sup>36</sup>.

In addition, during magma ascent, as the rheology of the magma is controlled by both volatile behavior and crystal content, the supply system may have a non-linear behavior<sup>6,37</sup>. For silicic magmas, the largest part of the magma ascent occurs in a closed system: ascent rate increases with decreasing pressure in response to the decrease in water solubility and to the gas expansion. If, for any reason (variation in wall-rock permeability, wall rock-magma interactions, etc.) the system opens, gas escapes and the ascent rate is dramatically reduced because gas pressure re-equilibrates with the surroundings and melt crystallizes, increasing the bulk magma viscosity<sup>32</sup>. For such slow extrusion rates—which are the case in dome-forming eruptions—magma reaching the surface is generally highly degassed and highly crystalline (85 to 95% solid fraction). In addition lava dome acquire a more fragile behavior that in turn favors fracturing and gas escape. At high extrusion rates, the melt has less time to crystallize during ascent in the conduit, resulting in a more fluid-like behavior<sup>12</sup>. Local shear-induced fragmentation at the conduit walls decoupling the magma column from the wall rock may favor the segregation within the conduit of “small” batches of magma with higher ascent rates, which still evolve in closed-system degassing<sup>10,38</sup>. Such magma batches are less degassed and have a higher vesicularity and a lower density. When they intrude the shallow zones of a small, slowly-growing lava dome with a well-formed silicified carapace, they will induce significant pressure buildup under that carapace that may trigger superficial explosion. In addition, the intrusion of less-dense material reduces the load exerted on the conduit and may in turn increase the bulk magma ascent rate<sup>39</sup>. Explosions are laterally-directed because of the lowest strength in the lateral parts of the small lava dome, which has not developed a large talus pile around its base, and large horizontal variations of gas pressure and viscosity gradient that may act as a driving force for lateral gas escape<sup>9,11,30</sup>.

Diffusion modeling of Ti in magnetites allows estimation of transit time of magma in the conduit. Diffusion profiles are mostly acquired for magnetites from vesiculated clasts as they preserved the diffusive information compared to dense clasts that exhibit generally exsolved magnetites (Fig. 5a,b,d). Results evidence two most likely time intervals depending on the vesicularity of the clasts, a short and a
<figure><img src="image_5.png" /><figcaption>Figure 4. Bulk crystalline silica (cristobalite) and Residual H₂O content (H₂Oᵣ) of samples from lava domes and Pyroclastic Density Current (PDC) deposits, as a function of glass vesicularity (in volume %).</figcaption></figure>

Bulk rock vesicularity and H₂O content are corrected from phenocryst contents to refer to melt, entitled respectively glass vesicularity and residual H₂O content (H₂Oᵣ) (see methodology). (a) Cristobalite content vs vesicularity glass. Cristobalite contents are expressed in area % and measured by chemical mapping (see methodology). *Lava domes*. Samples are taken on the different lava domes: Montagne Pelée, Puy de Dôme, Soufrière Hills, Montserrat, Santiaguito. They display a low vesicularity and have the highest content of cristobalite (18–28%). C-PDCs. Clasts have a vesicularity from 10 to 50% and are cristobalite-rich with mean cristobalite contents up to ~23% for Montagne Pelée, and ~5–10% for Puy de Dôme, Montserrat and Santiaguito. The cristobalite content is independent of the glass vesicularity. D-PDCs. D-PDC samples display a large range of clast types. Beyond a glass vesicularity threshold value of ~40% no or low cristobalite is observed in clasts. Below this threshold, some clasts don't contain cristobalite whereas others show a negative correlation between vesicularity and cristobalite with a maximum value of 20% (Montagne Pelée). Clasts with glass vesicularities below the threshold, with or without cristobalite, may represent the silicified and rigid lava dome carapace whereas clasts with glass vesicularities above that threshold may represent the inner, less-degassed and vesicated magma. (b) Residual H₂O content (H₂Oᵣ) vs. glass vesicularity. Blue square domain: lava dome and C-PDC samples. Lava dome samples: H₂Oᵣ < 0.4 wt% and V < 40%. C-PDC clasts: H₂Oᵣ < 0.7 wt% and V < 50%. Red domain: D-PDC clasts. H₂Oᵣ ranges from 0.2 to 2.5 wt% and glass vesicularity from 10 to 80%. The pre-eruptive H₂O contents measured on melt inclusions for all these eruptions are significantly higher (>>5 wt%$^{16,21-23}$).

significantly longer interval of time. In vesicated clasts, Ti-zoning modeling suggests a bimodal distri-
bution of time intervals with peaks centered at ~8 days (n = 6) and ~95 days (n = 4). Three other crystals
give timescales below resolution limit (<4 days), suggesting very rapid transit times. In dense clasts, long
time intervals (~100 days, n = 2) were found and a single short time interval (~14 days). The long transit
times of similar durations found in both vesicated and dense clasts suggests a synchronous emplace-
ment and common history. By contrast, the short timescales (~8 days, Fig. 5c1, c2) are only evidenced
in vesicated clasts, suggesting a late phase of rapid ascent for these magmas. Most magnetites from
dense clasts have undergone exsolution, a reaction to slow cooling and oxidation²⁵, on the contrary to
<figure><img src="image_6.png" /><figcaption>**Figure 5. Magnetite textures at Montagne Pelée and Fe-Ti zoning pattern.** Samples are from May, 8th, 1902 D-PDC (Montagne Pelée). (a) BSE image of an unexsolved magnetite (green line: TiO₂ content profile in Fig. 5 B). The presence of a ring of melt inclusions helps to constrain the profile location. (b) Ti-diffusion profile modelling of (A), leading to a short timescale in vesiculated clast (~8 days). (c) Relative probability. Residence time data are from titanomagnetite grains recovered from dense and vesiculated pumices c1: Vesiculated clasts leading to “young” timescale (0–30 days) and c2: Dense clasts leading to “old” timescale (30+ days). The graphs show relative probability on an arbitrary vertical scale; populations are scaled by the number of crystals such that the area under the curve should equal unity and can be compared to individual crystal data in that graph. Uncertainties on individual diffusion data points (0.15 log units, 1-sigma) are Gaussian in log-time, meaning that they are asymmetric in linear time, as displayed. The height along the probability axis is scaled in order to conserve area, so that older crystals have a lower peak probability to match their wider uncertainty, and do not dominate signals when combining into population curves. Graph c1 shows the relative clustering of crystals from vesiculated pumices at short timescales. The single short-timescale crystal from dense pumice (c2) plots among the oldest of the crystals from the vesiculated pumice. We infer this to represent a relative increase in ascent velocity for the vesiculated pumice, which is also consistent with the pervasive alteration and exsolution of oxides in the dense pumice. Graphs c1 and c2 show that the older population of crystals in the vesiculated and dense products is effectively coincident and may reflect a similar origin and magma processing pathway. (d) BSE image of an exsolved magnetite.</figcaption></figure>

<table><caption>**Table 2. Dense Rock Equivalent density for studied magmas and the cristallinity (wt%).**</caption><thead><tr><th></th><th></th><th>DRE</th><th>Cristallinity (wt%)</th></tr></thead><tbody><tr><td rowspan="5">Montagne Pelée, Martinique</td><td>P1, D-PDC</td><td>2,65</td><td>55</td></tr><tr><td>1902 D-PDC</td><td>2,65</td><td>50</td></tr><tr><td>1902 lava dome</td><td>2,65</td><td>40</td></tr><tr><td>1929 C-PDC</td><td>2,65</td><td>40</td></tr><tr><td>1929 lava dome</td><td>2,65</td><td>40</td></tr><tr><td rowspan="2">Chaine des Puys (France)</td><td>Puy Chopine</td><td>2,56</td><td>23</td></tr><tr><td>Puy de Dôme</td><td>2,66</td><td>21</td></tr><tr><td rowspan="2">Santiaguito (Guatemala)</td><td>1929, C PDC</td><td>2,65</td><td>29</td></tr><tr><td>Lava dome</td><td>2,65</td><td>33</td></tr><tr><td rowspan="2">Montserrat</td><td>Lava dome</td><td>2,75</td><td>50</td></tr><tr><td>C PDC</td><td>2,75</td><td>48</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Residual H₂O content (H₂Oᵣ) vs Vesicularity glass of clasts from dome-forming eruptions; comparison with clasts from plinian and vulcanian eruptions. Magmatic H₂Oᵣ is calculated from H₂O content measured in bulk clast corrected from phenocrysts content (crystallinity). Lava domes, C-PDCs and D-PDCs: blue and red domains as in Fig. 4. Plinian and vulcanian clasts: **Orange circle**: 1902 plinian eruption of Santa Maria (Guatemala;²¹,²²); **Open circle**: plinian phase of the P1 eruption of Montagne Pelée (650 y. BP; Martinique;²¹,²²); **Green triangle**: vulcanian phase of Soufrière Hills (Montserrat;²²). Lines refer to closed system degassing models from initial rhyolitic melts containing 6.5% H₂O (**orange line**: plinian-type eruptions;²¹,²²) and 1% of H₂O (**green dotted line**: vulcanian explosions;²³) which fit the H₂Oᵣ - vesicularity evolution of clasts. All C-PDC fragments plot below these two closed system evolution lines. D-PDC clasts of every studied eruption span a very large domain between H₂O-poor, vesicle collapsed clasts and H₂O-rich, highly vesiculated clasts from plinian eruptions. This extreme heterogeneity indicates that before explosion, the different parts of the lava dome have suffered highly variable degassing conditions varying between pure closed system degassing, likely at the center of the lava dome or as new intruding magma batch, to highly degassed in open system (extreme H₂O loss and vesicle collapse), likely at the lava dome upper external carapace.</figcaption></figure>

magnetites from vesiculated clasts. The single short timescale found in dense clasts (~14 days, n = 1), is twice that of the average found in vesiculated clasts (~8 days, Fig. 5c1, c2). We believe this to be strongly suggestive, if not definitive, of faster ascent rates in the case of vesiculated melt. From these results, we propose a two-step magma ascent process. During dome-forming eruptions magma ascent rates are generally low, favoring open system degassing and leading to the emplacement of dense and highly degassed magmas, with exsolved magnetite crystals. In addition, oxidative exsolution of magnetites is favored by fluid migration related to large development of cristobalite precipitation characteristic of these magmas²⁵. In some cases magma batches may ascend more rapidly, at least during the later stages, leading to vesiculated, less-degassed melts. Since the specific textural characteristic of these melts (high vesicularity, low cristallinity, short diffusion times scales and lack of exsolution in magnetites) are preserved, it indicates that the arrival of such magma batches occurs shortly before the superficial explosion and quenching.

Observations show that the laterally-directed explosions of lava domes occur mainly at the onset of the eruption when the lava dome has still a small size (Montagne Pelée, 1902²⁷, ⁸, 650 BP¹⁶, Puy Chopine²⁰). On the basis of the volume of the D-PDC deposits, we can estimate that for these eruptions, the volume of the lava domes before explosion was probably not greater than 0.01 km³. Two factors may explain these observations.

i/the load of a small lava dome exerts a low pressure on its basement and the upper parts of the feeding conduit: the pressure buildup necessary to generate explosion is therefore low²¹,³⁷. Possibly, exsolved volatiles may migrate from ascending magma and concentrate within erupting magma beneath the impermeable carapace, leading to overpressurization³⁷. The low proportion of vesiculated clasts (vesicularity > 40%) in the D-PDC (Fig. 1) indicates that only a small batch of vesiculated magma is at the origin of the pressure buildup (Fig. 7A). The fraction of vesiculated clasts in lava dome explosion products is much lower in the May 8th 1902 eruption of Montagne Pelée (7%) than in the P1 eruption (15%) and Puy Chopine eruption (30%). These variations likely reflect variations in the depth of the explosion or the size of the lava dome. An explosion occurring within the lava dome requires less energy and likely a lower volume of vesiculated, less-degassed magma than an explosion occurring at the base of a lava dome or in the upper part of the feeding conduit. Lava dome growth is a self-defeating process: when the lava dome grows, the load progressively increases and the overpressure generated by exsolved gas becomes insufficient to overcome the load pressure and the lava dome strength to trigger an explosion³⁹.

ii/the size of the growing lava dome also controls cristobalite precipitation efficiency. To generate an impermeable, resistant carapace through cristobalite precipitation, the flux of silica-rich fluids has to be
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Evolution model for explosive (A) and “non-explosive” (B) dome-forming eruptions. (A) In the case of explosive activity, the lava dome is small (less than 0.01 km³), possessing a rigid and impermeable silicified carapace formed through H₄SiO₄ flux. A small batch of vesiculated magma is isolated by local shear-fragmentation along the conduit wall, and causes a pressure build-up when intruding the growing lava dome. (B) In the case of a large lava dome (>to 0.01 km³) crystalline silica precipitation occurs pervasively and deep-seated fractures allow magma extrusion without overpressurization of the lava dome. The texture of the lava dome is in that case similar to the dense and silicified part of the explosive lava dome (Fig. 7A, photos 2 and 3). See text for further discussion. The figure was drawn by G. Boudon.</figcaption></figure>

large relative to the lava dome size (Fig. 7B). These fluids are generated over the whole H₂O-saturated magma column and percolate through the lava dome at the surface within which they precipitate cristobalite²⁷. For a given flux of fluids corresponding to a given magma extrusion rate, the mean thickness of the silicified carapace decreases as the volume of the lava dome increases. In addition, a voluminous lava dome doesn't favor the permanence of a continuous impermeable carapace due to development of deep-seated fractures intersecting the outer parts of the lava dome. A low carapace strength less efficiently maintains overpressure inside the lava dome and degassing may occur passively through deep-seated fractures, reducing its explosive potential (Fig. 7b).

In recent decades, lava dome-forming eruptions have exhibited various type of explosivity. During the November 2010 eruptive phase of Merapi (Indonesia), explosive destruction of the emerging lava dome occurred one week after the beginning of its growth⁴⁰, likely due to the ascent of an undegassed batch of magma³⁸. Mount Lamington (Papua New Guinea) experienced similar explosive activity that destroyed the growing lava dome in the first days of the 1951 eruption⁴¹. In the case of cryptodome emplacement (Mt St. Helens, 1980 or Bezimianny, 1956), similar features have been described in clasts sampled in the blast deposits, with the coexistence of vesiculated clasts, with high H₂O₃ content (up to 1.5 wt%) and dense clasts which are highly crystalline, silica-rich and highly-degassed⁴². The presence of crystalline silica precipitates in the rocks surrounding the cryptodome likely favoured the isolation of a less-degassed batch of magma⁴³. It is thought, however, that the triggering factor of the explosion was the large decompression induced by the flank-collapse, because the overpressure reached inside the small cryptodome was likely insufficient to trigger it alone⁶. Similarly, large dome collapse leading to a sudden depressurization of the magma column may trigger fragmentation deeper in the conduit and produce vulcanian explosions, as at Soufrière Hills (Montserrat) or at Mount Unzen (Japan)¹,³,⁵. These vulcanian products have a higher vesicularity and a lower H₂O content (<0.6 or 0.7 wt% H₂O)²³, Fig. 6) than the vesiculated clasts from the D-PDCs and represent fragmentation conditions (pressure and water content) of magmas quenched a few hundred meters deep in the conduit.

Lava dome eruptions are highly controlled by shallow H₂O degassing processes and silica precipitation. Magma silicification (by cristobalite precipitation) and densification (by flattening of the vesicles) are systematic processes during dome-forming eruptions. But the conditions able to generate sufficiently large impermeable carapace at the periphery of a growing lava dome, allowing a possible explosion caused by the overpressure generated by the intrusion of a batch of undegassed magma, are rarely achieved. Laterally directed explosions at the base of a growing lava dome or in the superficial
part of the conduit may occur when two conditions are simultaneously fulfilled: a small size of the lava dome and a significant precipitation of vapor phase silica in voids and vesicles creating an impermeable carapace that prevents further gas loss and allows local overpressurization. Estimation of initial sizes of destroyed lava dome is difficult and can only be performed using well documented lava-dome eruptions such as the Montagne Pelée ones and through estimation of the volume of the D-PDC's. The proposed model for superficial explosivity of lava domes is based on the combination of textural, geochemical and petrological criteria that are obtained through a systematic study of a very detailed sampling of deposits resulting from both superficial explosions and gravitational dome collapses: Puy Chopine, Montagne Pelee, Montserrat and Santiaguito. In the latter two cases sampling was limited by access concerns due to current eruptive activities. These eruptions also cover a large range of magma compositions: from andesitic to trachytic. We believe that the resulting model is applicable to a variety of volcanic systems worldwide and rhyolitic melt compositions, although testing is necessary for other melt compositions.

Superficial, lava dome explosion-generating D-PDC's, demonstrably devastating events, embody one of the most hazardous aspects of dome-forming eruptions. Even if this type of event is less predictable because it dominantly occurs at the onset of such eruptions, and because no method exists to reliably detect potential overpressurization within a growing lava dome, the possibility of such activity must be taken into account during the management of future volcanic crises, especially during the early stage of dome-forming eruptions. This type of volcanic activity is hazardous regarding both pyroclastic activity and the effects of crystalline silica-rich volcanic ash dispersed into the atmosphere⁴⁴.

## Methodology

**Sampling.** For all the eruptions, the lava dome was sampled. For the dilute pyroclastic density currents (May 8th, 1902 and P1 eruption of Montagne Pelée, Puy Chopine in La Chaîne des Puys), and for the concentrated pyroclastic density currents (1929 eruption from Montagne Pelée (2 deposits sampled), Puy de Dome eruption, 1929 eruption of Santiaguito), we collected at least 100 clasts for each deposit in a fraction 16–32 mm, to have a statistical representation of the density distribution. In the concentrated pyroclastic density currents from the recent eruption of Soufrière Hills, Montserrat, only few clasts were sampled in the deposits.

**Density and Vesicularity of clasts.** Densities are measured on fragments taken on the lava domes and at least one hundred clasts from each associated PDC deposits by the “3 weighing method” (weighing uncoated clasts in air, then paraffin coated clast in air and in water). The bulk rock vesicularity is calculated from density measurements using the Dense Rock Equivalent density and corrected for phenocryst content to refer to the initial melt (glass vesicularity). The phenocryst weight fraction is estimated by point counting (Puy de Dôme and Puy Chopine samples⁴⁵) or by mass balance calculations on separate minerals and matrix using major and trace element data¹⁶.

**Residual H₂O content.** Water contents are measured by H₂ manometry (adapted from 46) on bulk samples spanning the whole density range of each eruption. Data are then corrected for phenocryst content to obtain the residual H₂O content of glass or matrix (H₂Oᵣ). Water measured is considered as pure magmatic water on selected samples, as ensured by step-heating during extraction⁴⁷. In addition, textural investigations never evidenced rehydration processes (Fig. 2). δD composition measurements of lava dome fragments and pumice clasts from the P1 eruption exclude significant contribution of meteoric water in this material¹⁶. Similar isotopic data on obsidian clasts from fallout deposits, lava domes and flows from the Western United States show that total water content is controlled by degassing processes in such samples of similar composition⁴⁸,⁴⁹.

**Textures (vesicularity and microcrystallinity) and crystalline silica (cristobalite) contents.** They were determined using SEM images and chemical mapping by Energy Dispersive X-ray Spectrometry (EDS) acquired with a Zeiss Supra V55 equipped with a Brucker Silicon drift Detector (UPMC, Paris)³³. Chemical mapping (2 to 5 maps per sample) of major elements (Si, Ti, Al, Fe, Mg, Ca, Na, K, P, Cl) at magnification 100× and 300× are used to estimate the contents in residual glass and crystals (phenocrysts + microlites) and cristobalite precipitates (Si maps). The mean relative standard deviation is ~15 area % for cristobalite contents.

**Diffusion modelling.** Diffusion chronometry estimates the time elapsed since the onset of a disequilibrium, for example when crystals are rapidly moved to an environment with different physical or chemical conditions (T, P, fO₂, composition), or when they grow rims that are not in equilibrium with their cores. The existence of partial, diffusion-mediated re-equilibration allows the timing of such particular phenomena to be determined by measuring the extent and lengthscale of re-equilibration. The rate at which diffusion occurs is strongly dependent on temperature, therefore when the crystal cools these compositional gradients can become “frozen in”. Modeling the measured composition (elements or isotopes) gradients allows estimation of the time between the last disequilibrium event and quenching. Such modeling using Fick's diffusion laws requires the knowledge of the diffusion coefficients and composition measurements at high spatial resolution and analytical precision. Changes in temperature, pressure and volatile content of the environment of crystals may significantly affect the growth and
composition of volcanic minerals such as magnetites, feldspars and pyroxenes²⁴,²⁵,⁵⁰ (Fig. 5a,b). Here the methodology developed for pyroxenes by Morgan on the examples of Oranaui eruption (New Zeland)²⁴ is adapted for magnetites. Ti zoning in magnetite is essentially temperature-dependent but may also be affected by sub-solidus exsolution (Fig. 5d). If magnetites are unexsolved, Ti content may be used to constrain T and/or fO₂ changes during the late stages of magma ascent in the conduit and magnetite crystallisation, shortly before quenching by explosion (timescales of the order of days to months,⁵⁰). At Montagne Pelée, vesiculated (V > 50%) and dense (V < 30%) clasts were crushed and the powder sieved. Magnetites were separated by hand-picking in the fraction 250–315 µm. Ti profiles of unexsolved magnetites were obtained using first SEM images to identify the best candidates for diffusion modeling. Ti zoning in crystals is difficult to recognize on the basis of SEM images or chemical mapping due to the low contrast in composition, but the occurrence of a ring of melt inclusions helps identification of suitable crystals (Fig. 5a). Then, quantitative Ti profiles were performed, together with major elements (Si, Ti, Al, Fe, Mn, Mg, Ca, Na, K, Co) by electron microprobe (Cameca, SX-Five, Camparis, UPMC, France). Analyses are performed at 15 keV and 20 nA with a focused beam of 3 µm. We used counting time of 10s for all elements except Ti(60s). Diffusion profiles are between 40 and 120 µm long, from core to rim with a 5 µm step. Then Ti diffusion profile may be fit by diffusion models at 875–900 °C and fO₂ of NNO+0.4–0.8⁵¹ using the diffusion data published giving a diffusion coefficient of 2.22854 10⁻¹⁷ m²/s⁵². Uncertainty is estimated by combining uncertainties in diffusion temperature and oxygen fugacity, and propagating this onto the timescale; the uncertainties are expressed as Gaussian in log-time, with a 1-sigma uncertainty value of ~0.15 log units, meaning that they display asymmetrically in linear time as shown in Fig. 5.

## References

1. Newhall, C. G. & Melson, W. G. Explosive activity associated with growth of volcanic domes. *J. Volcanol. Geotherm. Res.* **17**, 111–131 (1987).

2. Druitt, T. H. Pyroclastic density currents, in *The physics of explosive volcanic eruption*. Vol. 145 (eds Gilbert, J. S., Sparks, R. S. J.) 145–182 (Geological Society, London, Special Publications, 1998).

3. Nakada, S., Shimizu, H. & Ohta, K. Overview of the 1990–1995 eruption at Unzen Volcano. *J. Volcanol. Geotherm. Res.* **89**, 1–22 (1999).

4. Voight, B. et al. Historical eruptions of Merapi Volcano, Central Java, Indonesia, 1768–1998. *J. Volcanol. Geotherm. Res.* **100**, 69–138 (2010).

5. Sparks, R. S. J. & Young, S. R. The eruption of Soufrière Hills volcano, Montserrat (1995–1999): overview of scientific research, in *The Eruption of Soufrière Hills Volcano, Montserrat, from 1995 to 1999* (Druitt, T. H., Kokelaar, B. P., eds.), 115–152 (Geological Society, London, Memoirs 21, 2002).

6. Sparks, R. S. J. Causes and consequences of pressurization in lava domes. *Earth Planet. Sc. Lett.* **160**, 177–189 (1997).

7. Lacroix, A. *La Montagne Pelée et ses éruptions* (Masson, 1904).

8. Bourdier, J.-L., Boudon, G. & Bourgaud, A. Stratigraphy of the 1902 and 1929 nuée ardente deposits, Mount Pelée, Martinique. *J. Volcanol. Geotherm. Res.* **38**, 77–96 (1989).

9. Eichelberger, J. C., Carrigan, C. R., Westrich, H. R. & Price, R. H. Non-explosive silicic volcanism. *Nature* **323**, 598–602 (1986).

10. Fink, J. H. & Manley, C. R. Origin of pumiceous and glassy textures in rhyolite flows and domes, in *The emplacement of silicic domes and lava flows* (Fink, J. H., ed.), 77–88 (Geol. Soc. America Special Paper 212, 1987).

11. Gonnermann, H. M. & Manga, M. Explosive volcanism may not be an inevitable consequence of magma fragmentation. *Nature* **426**, 432–435 (2003).

12. Massol, H. & Jaupart, C. The generation of gas overpressure in volcanic eruptions. *Earth Planet. Sc. Lett.* **166**, 57–70 (1999).

13. Hale, A. J. & Mühlhaus, H. B. Modelling shearbands in a volcanic conduit: Implications for over-pressures and extrusion-rates. *Earth Planet. Sc. Lett.* **263**, 74–87 (2007).

14. Marzocchi, W., Newhall, C. & Woo, G. The scientific management of volcanic crises. *J. Volcanol. Geotherm. Res.* **247–248**, 181–189 (2012).

15. Perret, F. A. *The eruption of Mt. Pelée 1929–1932*. 458 (Publication of Carnegie Institution of Washington, 1937).

16. Villemant, B. & Boudon, G. Transition between dome-forming and plinian eruptive styles: H₂O and Cl degassing behavior. *Nature* **392**, 65–69 (1998).

17. Young, S. Summary of eruptive events and monitoring procedures at Soufrière Hills volcano Montserrat West Indies, 18 July 1995 to June 1996. *IAVCEI News.* **1/2** (1996).

18. Rose, W. I. Pattern and Mechanism of Volcanic Activity at the Santiaguito Volcanic Dome, Guatemala. *Bull. Volcanol.* **37**, 73–94 (1973).

19. Scott, J. A. J., Pyle, D. M., Mather, T. A. & Rose, W. I. Geochemistry and evolution of the Santiaguito volcanic dome complex, Guatemala. *J. Volcanol. Geotherm. Res.* **252**, 92–107 (2012).

20. Boivin, P. et al. *Volcanologie de la Chaine des Puys, 5ème édition* (Parc Naturel Régional de la Chaine des Puys, 2009).

21. Villemant, B., Boudon, G., Nougrigat, S., Poteaux, S. & Michel, A. H₂O and halogens in volcanic clasts: tracers of degassing processes during plinian and dome-forming eruptions. in *Volcanic Degassing*. Vol. 213 (Eds Oppenheimer, C. et al.) 63–79 (Geological Society of London, Special Publication, 2003).

22. Balcone-Boissard, H., Villemant, B. & Boudon, G. Behavior of halogens during the degassing of felsic magma. *Geochemistry, Geophysics and Geosystems* **11**, 477–485 (2010).

23. Villemant, B., Mouatt, J. & Michel, A. Andesitic Magma Degassing Investigated Through H₂O Vapour-Melt Partitioning of Halogens at Soufrière Hills Volcano, Montserrat (Lesser Antilles). *Earth Planet. Sci. Lett.* **269**, 212–229 (2008).

24. Allan, S. R., Morgan, D. J., Wilson, C. J. N. & Millet, M.-A. From mush to eruption in centuries: assembly of the super-sized Oruanui magma body. *Contrib. Mineral. Petrol.* **166**, 143–164 (2013).

25. Buddington, A. F. & Lindsley, D. H. Iron-Titanium Oxide Minerals and Synthetic Equivalents. *J. Petrol.* **5**, 310–357 (1964).

26. Fournier, R. O. & Rowe, J. J. The solubility of cristobalite along the three-phase curve, gas plus liquid plus cristobalite. *Amer. Mineral.* **47**, 897–902 (1962).

27. Martel, C. & Schmidt, B. C. Decompression experiments as an insight into ascent rates of silicic magmas. *Contrib. Mineral. Petrol.* **144**, 397–415 (2003).

28. De Hoog, J. C. M., van Bergen, M. J. & Jacobs, M. H. G. Vapour-phase crystallisation of silica from SiF₄-bearing volcanic gases: *Annals of geophysics* **48**, 775–785 (2005).
29. Bourgouin, L., Mühlhaus, H.-B., Hale, A. J. & Arsac, A. Studying the influence of a solid shell on lava dome growth and evolution using the level set method. *Geophys. J. Int.* **170**, 1431–1438 (2007).

30. Sparks, R. S. J. Dynamics of magma degassing. *in Volcanic Degassing*. Vol. 213 (Eds Oppenheimer, C. et al.) 5–22 (Geological Society of London, Special Publication, 2003).

31. Fink, J. H., Anderson, S. W. & Manley, C. R. Textural Constraints on Effusive Silicic Volcanism: Beyond the Permeable Foam Model. *J. Geophys. Res.* **97**, 9073–9083 (1992).

32. Lejeune, A.-M. & Richet, P. Rheology of crystal-bearing silicate melts: An experimental study at high viscosities. *J. Geophys. Res.* **100**, 4215–4229 (1995).

33. Balcone-Boissard, H., Boudon, G. & Villemant, B. Textural and geochemical constraints on eruptive style of the 79AD eruption at Vesuvius. *Bull. Volcanol.* **73**, 279–294 (2011).

34. Fink, J. H. & Kieffer, S. W. Estimate of pyoclastic flow velocities resulting from explosive decompression of lava domes. *Nature* **363**, 612–615 (1993).

35. Spieler, O. et al. The fragmentation threshold of pyroclastic rocks. *Earth Planet. Sc. Lett.* **226**, 139–148 (2004).

36. Holland, A. S. P., Watson, I. M., Phillips, J. C., Caricchi, L. & Dalton, M. P. Degassing processes during lava dome growth: Insights from Santiaguito lava dome, Guatemala. *J. Volcanol. Geotherm. Res.* **202**, 153–166 (2011).

37. Melnik, O. & Sparks, R. S. J. Controls on conduit magma flow dynamics during lava dome building eruptions. *J. Geophys. Res.* **110**, doi: 10.1029/2004JB003183 (2005).

38. Costa, F., Andreastuti, S., Bouvet de Maisonneuve, C. & Pallister, J. S. Petrological insights into storage conditions, and magmatic processes that yielded the centennial 2010 Merapi explosive eruption. *J. Volcanol. Geotherm. Res.* **261**, 209–235 (2013).

39. Taisne, B. & Jaupart, C. Magma degassing and intermittent lava dome growth. *Geophys. Res. Lett.* **35**, doi: 10.1029/2008GL035432 (2008).

40. Komorowski, J.-C. et al. Paroxysmal dome explosion during the Merapi 2010 eruption: Processes and facies relationships of associated high-energy pyroclastic density currents. *J. Volcanol. Geotherm. Res.* **261**, 260–294 (2013).

41. Taylor, G. A. The 1951 eruption of Mt. Lamington, Papua, Australia. *Bull. Min. Res., Geol. Geophys. Bull.* **38**, 1–112 (1958).

42. Neill, O. K. et al. Influence of pre-eruptive degassing and crystallization on the juvenile products of laterally directed volcanic explosions. *J. Volcanol. Geotherm. Res.* **198**, 264–274 (2010).

43. Komorowski, J.-C., Glicken, H. X. & Sheridan, M. F. Secondary electron imagery of microcracks and hackly fracture surfaces in sandstone-sized clasts from the 1980 Mount St. Helens debris avalanche deposits: Implication for particle-particle interactions. *Geology* **19**, 261–264 (1991).

44. Baxter, P. J. et al. Cristobalite in volcanic ash of the Soufriere Hills Volcano, Montserrat, British West Indies. *Science* **283**, 1142–1145 (1999).

45. Martel, C. et al. Trachyte phase relations and implications for magma storage conditions in the Chaîne des Puys (French Massif Central). *J. Petrology* **54**, 1071–1107 (2013).

46. Michel, A. & Villemant, B. Determination of halogens (F, Cl, Br, I), sulfur and water in seventeen geological reference materials. *Geostandard Newsletters* **27**, 163–171 (2003).

47. De Groat-Nelson, P. J., Cameron, B. I., Fink, J. H. & Holloway, J. R. Hydrogen isotope analysis of rehydrated silicic lavas: implications for eruption mechanisms. *Earth Planet. Sci. Lett.* **185**, 331–341 (2001).

48. Barnes, J. D. et al. Stable chlorine isotope behavior during volcanic degassing of H₂O and CO₂ at Mono Craters, CA. *Bull. Volcanol.* **76**, doi: 10.1007/s00445-014-0805-y (2014).

49. Taylor, B. E., Eichelberger, J. C. & Westrich, H. R. Hydrogen isotopic evidence of rhyolitic magma degassing during shallow intrusion and eruption. *Nature* **306**, 541–545 (1983).

50. Nakamura, M. Continuous mixing of crystal mush and replenished magma in the ongoing Unzen eruption. *Geology* **23**, 807–810, doi: 10.1130/0091-7613 (1995).

51. Martel, C. et al. Magma storage conditions and control of eruption regime in silica volcano: experimental evidence from Mt. Pelée. *Earth Planet. Sc. Lett.* **156**, 89–99 (1998).

52. Aragon, R., McCallister, R. H. & Harrison, H. R. Cation diffusion in titanomagnetites. *Contrib. Mineral. Petrol.* **85**, 174–185 (1984).

## Acknowledgments

We would like to thank A. Michel, O. Boudouma and S. Hidalgo for analytical support, P. Boivin providing the Puy de Dome lava sample, “le Volcan de Lemptégy” for access to the outcrop of Puy Chopine D-PDC deposit, Volcanological Observatory of Martinique (OVSM) for logistical support, O. Matias for his help on field on Santa Maria-Santiaguito, G. Ucciani for the development of a Matlab code for the treatment of the chemical maps. We thank also M. Manga for discussion on a preliminary version of the manuscript and J.F. Fink for his constructive review. IPGP contributions: N° 3668.

## Author Contributions

G.B. and H.B.B. defined and performed all data acquisition and interpretation, wrote the manuscript and prepared the figures. B.V. participated to the Raman data acquisition, all data interpretation and the writing of the manuscript text. D.J.M. was involved in the diffusion data acquisition, interpretation and the writing of the manuscript text. All authors reviewed the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Boudon, G. et al. What factors control superficial lava dome explosivity? *Sci. Rep.* **5**, 14551; doi: 10.1038/srep14551 (2015).

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

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/