# Inflated pyroclasts in proximal fallout deposits reveal abrupt transitions in eruption behaviour

Thomas J. Jones, Yannick Le Moigne, James K. Russell, Glyn Williams-Jones, Daniele Giordano, Donald B. Dingwell

► To cite this version:

Thomas J. Jones, Yannick Le Moigne, James K. Russell, Glyn Williams-Jones, Daniele Giordano, et al. Inflated pyroclasts in proximal fallout deposits reveal abrupt transitions in eruption behaviour. Nature Communications, 2022, 13, 327, p. 53-78. 10.1038/s41467-022-30501-6. insu-03825086

HAL Id: insu-03825086

https://insu.hal.science/insu-03825086

Submitted on 21 Oct 2022

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>

ARTICLE

https://doi.org/10.1038/s41467-022-30501-6

OPEN

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

Check for updates

# Inflated pyroclasts in proximal fallout deposits reveal abrupt transitions in eruption behaviour

Thomas J. Jones ID<fnref n="1" />✉, Yannick Le Moigne²,³, James K. Russell ID<fnref n="4" />, Glyn Williams-Jones ID², Daniele Giordano⁵ & Donald B. Dingwell ID⁶

During explosive eruption of low viscosity magmas, pyroclasts are cooled predominantly by forced convection. Depending on the cooling efficiency relative to other timescales, a spectrum of deposits can be formed. Deposition of hot clasts, above their glass transition temperature, can form spatter mounds, ramparts and clastogenic lava flows. Clasts may also be deposited cold, producing tephra cones and blankets. Thus, the deposit and pyroclast type can provide information about eruption dynamics and magma properties. Here we examine pyroclasts from Tseax volcano, British Columbia, Canada. These newly identified inflated pyroclasts, are fluidal in form, have undergone post-depositional expansion, and are found juxtaposed with scoria. Detailed field, chemical and textural observations, coupled with high temperature rheometry and thermal modelling, reveal that abrupt transitions in eruptive behaviour — from lava fountaining to low-energy bubble bursts — created these pyroclastic deposits. These findings should help identify transitions in eruptive behaviour at other mafic volcanoes worldwide.

<footnote n="1">¹Department of Earth, Ocean and Ecological Sciences, University of Liverpool, Liverpool L69 3GP, UK. ²Centre for Natural Hazards Research, Department of Earth Sciences, Simon Fraser University, Burnaby, BC, Canada. ³Laboratoire Magmas et Volcans, Université Clermont-Auvergne, Clermont-Ferrand, France.</footnote>

<footnote n="4">⁴Department of Earth, Ocean and Atmospheric Sciences, University of British Columbia, Vancouver, BC, Canada. ⁵Dipartimento di Scienze della Terra, Università degli studi di Torino, 10125 Torino, Italy. ⁶Department of Earth and Environmental Sciences, Ludwig-Maximilians-Universität, 80333 Munich, Germany. ✉email: thomas.jones@liverpool.ac.uk</footnote>
Low viscosity magmas are the most common erupted on Earth¹,². The associated eruptions can be long-lived, causing large economic losses and can be highly variable both in magnitude and style, thus, making hazard management challenging³⁻⁶. The style and efficiency of magma fragmentation directly controls the type of products erupted and, in turn, the hazards they pose⁷⁻⁹. For example, Hawaiian style fountaining of fluidal clasts frequently produces vent proximal spatter ramparts and mounds but can also yield coeval lava flows where accumulation rates are high enough and clasts remain hot¹⁰⁻¹⁴ for a longer time. Such lavas pose an additional risk to communities and infrastructure down flank³. If fragmentation is highly efficient and clasts are cooled to a greater extent, the resulting widespread tephra blankets pose a risk to a larger geographic area and in some cases may even cause closure of airspace¹⁵. Despite the diverse array of hazards that such mafic eruptions pose, too little is known about pyroclast transport and the associated eruption dynamics⁷,¹⁶,¹⁷ in such eruptions.

Explosive activity at mafic volcanoes, especially in its milder forms (i.e., lava fountaining, spattering events and discrete bubble bursts) can be highly variable in its surface expression¹⁸⁻²⁰. Eruption intensities and styles of activity can also fluctuate over short timescales and distances²¹,²²—this presents challenges to hazard management²³ and hinders systematic observation of eruptive phenomena. Fluctuations in intensity and style can occur between eruptive episodes, within a single episode, and simultaneously at different spatial locations²¹. Furthermore, the physical processes governing explosive eruption of low viscosity magmas are fundamentally different to those operating in the now relatively well-studied silicic systems²⁴⁻²⁷. Most commonly their fragmentation processes are fluid dynamic in nature such that breakage does not occur simply by decompression upon crossing the glass transition⁷,²⁷⁻²⁹. Instead, their low melt viscosities allow surface tension-driven reshaping and bubble nucleation, growth, and coalescence to operate on syn- and post-eruptive timescales. A common example includes achneliths such as Pele's tears where in-flight surface tension-driven relaxation³⁰ transforms irregular pyroclast morphologies into more spherical shapes. The consequence is that pyroclast properties are susceptible to modification after the initial fragmentation¹²,²⁹⁻³³. As a result, the eruptive deposits derived from low viscosity magmas are challenging to interpret as: (1) vent structures evolve and migrate during an eruption;³⁴⁻³⁶ (2) early products are frequently buried or rafted by subsequent eruption episodes³⁷⁻³⁹ and (3) rheomorphic processes and in-flight secondary fragmentation processes destroy primary fragmentation features¹,¹¹,¹⁷,²⁹,³³,⁴⁰,⁴¹. These multiple factors hinder our understanding of mafic eruption dynamics solely by viewing their associated products.

Here we examine pyroclasts from the ~1700s CE eruption of Tseax in northwest British Columbia, Canada⁴²⁻⁴⁵, which is reported to have caused up to 2,000 fatalities making it the deadliest eruption in Canada's history⁴⁶,⁴⁷. The pyroclasts are fluidal in form and have undergone post-depositional expansion and are therefore termed inflated pyroclasts herein. We performed a set of detailed field, chemical, and textural observations coupled with high temperature rheometry and thermal modelling to reveal the unique and abrupt transitions in eruptive behaviour that created these pyroclasts.

## Results and discussion

**Tseax volcano.** Tseax volcano, located in the Nisga'a Memorial Lava Bed Provincial Park in northwest British Columbia, Canada, is the southernmost volcanic centre of the Northern Cordillera Volcanic Province⁴⁸ and is notable for a 32-km-long basaltic trachybasalt lava flow (~0.5 km³ covering ~36 km²; Fig. 1).

<figure><img src="image_4.png" /><figcaption>Fig. 1 Location and description of the area surrounding Tseax volcano.</figcaption></figure>

a Tseax (red hexagon) in the context of Neogene-Quaternary volcanic centres and complexes of the Northern Cordilleran Volcanic Province (NCVP). b Map of Tseax and the 32-km-long valley-filling lava flows (yellow). c Volcanological map of Tseax cone, the Satellite cones and spatter rampart, all shaded in red. Yellow stars indicate locations of inflated pyroclast samples. Red contours denote tephra deposit thickness. Modified after Le Moigne et al.⁴⁴. d Oblique 3D view of Tseax and surrounding area generated by digital photogrammetry. Modified after Le Moigne et al.⁴⁴.

The tephra cone and valley-filling lavas overlay intercalated sandstone to siltstone and mudstone from the late Jurassic Bowser Lake group⁴⁹,⁵⁰. Detailed volcanological mapping shows that Tseax volcano consists of a ~65 m high, 350-400 m diameter tephra cone ($2.8 \pm 0.4 \times 10^6$ m³) situated within a 550-600 m diameter oxidised horseshoe-shaped spatter rampart⁴⁴. Another smaller (~20 m high, ~55 m diameter, $1.7 \pm 0.1 \times 10^4$ m³) highly oxidised tephra cone, unnamed but referred to here as Satellite cone, is located 470 m to the north of Tseax and in close proximity to a group of small tephra mounds⁴⁴ (Fig. 1).

**Deposit description.** The highly vesicular, fluidal inflated pyroclasts were found at two separate vent proximal locations at Tseax volcano (yellow stars in Fig. 1c). The southern site is a
<figure><img src="image_5.png" /><figcaption>**Fig. 2** Deposit scale observations of the Tseax Satellite cone. **a** Satellite scoria cone (height ~20 m) with inflated pyroclasts only found within ~5 m of the summit. The black arrow points to the location of the other figure panels. **b–d** Variably oxidised scoria cone deposits at the summit in the surface deposits. Some inflated pyroclasts are identified by the white arrows. These clasts have a smooth outer surface, a fluidal morphology and are commonly rounded in shape. **e** Close up of an inflated pyroclast (~5 cm) surrounded by scoria. **f** Close up of a dissected block (30 cm in length) found at the summit of the Satellite cone which contains several inflated pyroclasts that are fluidal in nature and infill the space between scoria clasts. Further deposit photographs can be found in Fig. S1.</figcaption></figure>

spatter rampart wall, containing rare inflated pyroclasts and slightly rheomorphic spatter clasts that have been variably oxidised. At this location, the inflated pyroclasts make up <2% of the deposit by volume and are not concentrated in any stratigraphic horizon. The northern site features much better exposure and is the primary location sampled and investigated in this study. The inflated pyroclasts occur within a ~5 m radius of the top of the Satellite cone (Fig. 2a). The near-summit region of the Satellite cone has no pit or crater. The deposit is variably oxidised and dominated by moderately to well-sorted scoria clasts, with the inflated pyroclasts contributing <5% of the deposit by volume. No clear deposit stratigraphy was identified. In places the deposit is agglutinated, however individual pyroclasts can be identified (Fig. 2b–d) and in most cases the inflated pyroclasts can be separated from the deposit by hand whilst remaining intact.

Although partially agglutinated, there is no evidence of compactional welding at the Satellite cone. The inflated pyroclasts are only observed at vent proximal locations and fill the available void space in the deposit (Fig. 2e, f). In some cases, they are observed to have drained by gravity within the deposit or very occasionally are deformed by the pyroclast(s) directly above, forming squeeze-out textures. Further deposit photographs can be found in Fig. S1.

**Inflated pyroclast morphology and surfaces.** The individual inflated pyroclasts are commonly spherical in shape (Fig. 3a), but more elongated morphologies are also found. Their surfaces have a dull, matt appearance, rather than the shiny, glassy surfaces commonly observed on Pele's tears, for example. We interpret
<figure><img src="image_6.png" /><figcaption>**Fig. 3** Inflated pyroclast textures and detailed features. **a** Close up photograph of an inflated pyroclast within the scoria cone deposit. The clast surface contains a few small (≤0.5 mm) holes. **b–d** Scanning electron micrographs of inflated pyroclast surfaces showing: **b** circular non-brittle holes on the pyroclasts surface; **c** smaller pieces of vesicular scoria stuck to the inflated pyroclast and **d** multiple layers or skins. Further photographs of individual pyroclasts can be found in Fig. S2.</figcaption></figure>

this dull appearance to be a consequence of microlitic⁵¹ and possibly further nanolitic⁵² crystallisation (see next sub-section). The surfaces of the inflated pyroclasts are pierced by a few randomly distributed (sub)circular holes that are on the order of 100 µm in diameter (Fig. 3b, c). Images of these holes, obtained by scanning electron microscopy (SEM), reveal that they have smooth edges and do not show clear evidence of brittle rupture, suggesting that they formed when the pyroclast was above the glass transition temperature, T_g. Additionally, on some clasts, larger (≥ 1 mm) surface ruptures are observed revealing a second skin underneath (Fig. 3d). Further photographs of individual pyroclasts can be found in Fig. S2.

**Inflated pyroclast petrography.** The Tseax volcanic deposits are weakly porphyritic containing phenocrysts, glomerocrysts (>0.5 mm) and microphenocrysts (0.2–0.5 mm) of plagioclase, olivine and titanomagnetite. The dominant phenocryst phase is plagioclase (70%) followed by olivine (25%) and oxide (>5%). Clinopyroxene is never observed as a phenocryst phase. For the inflated pyroclasts, the focus of this study, we used thin section observations to record the internal texture and petrography of the pyroclast (Fig. 4). The oxides have previously been identified as titanomagnetite⁵³. Plagioclase crystals show albite and poly-synthetic twinning and a composition of ~An35 is obtained using the Michel–Levy Method. One large and fractured (~ 1 mm long) plagioclase is present (Fig. 4f) as well as a fragment of a plagioclase xenocryst (Fig. 4h). The olivine crystals are euhedral to subhedral, have sizes up to 0.45 mm in length, frequently exhibit hopper textures (Fig. 4e) and intracrystalline fractures are rarely filled with iddingsite. Titanomagnetite occurs as rare subhedral microphenocrysts. The sample is hyper-vesiculated (~80% vesicles), with vesicles that range in size from 0.1 mm to as large as the thin section slide (~40 mm).

The glomerocrysts consist of small aggregates of 5–12 phenocrysts of plagioclase (dominant) + olivine + oxide (Fig. 4d, i). Glomerocrysts exhibit poikilitic olivine with oxide oikocrysts and the plagioclase is interstitial to olivine and plagioclase. This suggests the following crystallisation sequence: (1) oxide, (2) olivine and (3) plagioclase. The matrix is highly crystallised and mainly composed

of plagioclase microlites, olivine and oxides. Note that there is a large plagioclase microlite population even on the thinnest part of the sample (Fig. 4g). Interstitial glass is rare, and the plagioclase laths are randomly oriented suggesting that bubble growth occurred in the absence of these microlites. Only phenocrysts occasionally distort bubble shapes (Fig. 4i) and thus were present in the melt during the dominant period of bubble growth and coalescence.

**Pyroclast chemical and physical properties.** All Tseax eruption products are alkali, Fe-Ti-rich basanite-to-trachybasalt in composition⁵³. Bulk major element chemical compositions determined by X-Ray fluorescence measurements of representative samples are shown in Table S1.

Due to agglutination within the proximal fallout deposit, it was not possible to measure a full pyroclast size distribution. The axis dimensions of the five largest and five smallest inflated pyroclasts can be found in Table S2. Bulk density (i.e., combined density of rock and vesicles comprising the entire pyroclast) was measured using the Archimedes method (see Methods) on the proximal scoria and the inflated pyroclasts that could be easily removed from the deposit (Fig. 5; Supplementary Data 1). Scoria bulk density varied between 621 and 1429 kg m⁻³ with a mean bulk density of 967 kg m⁻³. Inflated pyroclast bulk density varied between 311 and 1235 kg m⁻³ with a mean bulk density of 755 kg m⁻³. Using Helium pycnometry, the dense rock equivalent of the pyroclasts was calculated to be 2957 kg m⁻³ (Fig. S3). This DRE density was used to calculate the bulk vesicularity distributions shown in Fig. 5. The scoria has a bulk vesicularity ranging from 50% to 78% with a mean of 66% and the measured inflated pyroclasts range in vesicularity from 58% to 89% with a mean of 74%.

To document the internal textures of these inflated pyroclasts and to assess how they differ from the adjacent regular scoria, a sample containing both pyroclast types was scanned using micro-X-Ray computed tomography (XRCT). Representative image slices are shown in Fig. 6a and b for the scoria and inflated pyroclast, respectively. The proximal scoria clast (Fig. 6a) has a high vesicle number density, with most vesicles sub-mm in size. The vesicles (i.e., empty pore space), shown in black, range from spherical to highly irregular in shape. Dense, crystalline phases are shown in white. (Micro)phenocrysts in the scoria forming the satellite cone contribute a total modal abundance of <16% with typical mineral abundances of 3% olivine, 3–5% plagioclase and 1–3% oxides. In the sample shown in Fig. 6a, the small ~100 µm equant crystals are interpreted to be iron oxides and the larger ~1 mm elongate laths, plagioclase. The inflated pyroclast (Fig. 6b) has a much lower bubble number density and a coarser bubble size distribution. The inflated pyroclast is dominated volumetrically by large round bubbles at its centre. These larger bubbles are enclosed by a thin ≥50 µm crystal-bearing quenched melt film. The bubbles within this film are often rounded in shape and sometimes intersect the pyroclast exterior. In rare instances, the bubbles in the film are observed to intersect both the large internal bubble and the exterior surface of the pyroclast to form a narrowly interconnected gas-escape pathway (cf. Fig. 3b, c).

**Melt viscosity.** We measured the high-temperature, anhydrous viscosity of the Tseax melt (Table 1; Fig. 7). As detailed in Methods, the starting materials were obtained by direct fusion of the whole rock sample TS-S71 (Table S1) and were measured using a concentric cylinder apparatus in equilibrium with air at 1 atmosphere. The superliquidus anhydrous liquid viscosities ranged from 10<sup>0.3</sup>–10<sup>1.2</sup> Pa s over the explored temperature range of 1267–1531 °C. We also measured the 1 atmosphere subliquidus viscosity where it rose to between 10<sup>1.3</sup>–10<sup>3.4</sup> Pa s over the
<figure><img src="image_7.png" /><figcaption>**Fig. 4** Microphotographs of an inflated pyroclast (Sample: IP_211015_A; see Fig. S2h). **a–c** PPL, phenocrysts of plagioclase, olivine and oxide in a microcrystalline matrix (plagioclase+olivine+oxide+ minor interstitial glass). **d, e** XPL, same as **a–c**. **f** XPL, a 1 mm long plagioclase phenocryst surrounded by three large vesicles and several smaller ones (0.1–0.2 mm in diameter). **g** PPL, close up of a thin boundary between two large vesicles. Note the dense population of plagioclase microlites. **h** XPL, a plagioclase microphenocryst and microcrystalline matrix between two vesicles. A broken fragment of a plagioclase xenocryst is in the upper left corner of the microphotograph. **i** A glomerocryst surrounded by three large vesicles. Glomerocrysts of plagioclase +olivine+oxide are outlined in a dashed yellow line in (d, i). Phenocryst phases are labelled. Abbreviations: pl: plagioclase; ol: olivine; ox: oxide (previously identified as titanomagnetite⁵³).</figcaption></figure>

temperature range 1123–1243 °C. Micro-Raman spectroscopy was performed on the experimental glass to check for the presence of nanolites that may have formed during the quench; the absence of nanolites confirms that superliquidus viscosity measurements represent melt viscosities.

Superliquidus measurements of viscosity (white circles; Fig. 7) define a smooth trend against reciprocal temperature, whilst subliquidus values of viscosity (grey circles; Fig. 7) show a sharp rise in viscosity due to crystallisation and show greater variance. The experimental dataset can be compared to model temperature-dependent curves predicted from the compositions of the remelted Tseax samples. The two sets of viscosity model parameters (GRD⁵⁴ and Raman⁵⁵,⁵⁶) for the remelted glass are reported in Table 1 and reproduce the experimental measurements of melt viscosity well (± 0.25 log units; dashed lines Fig. 7). Furthermore, we compared these datasets to the GRD model for natural inflated pyroclast (green line; Fig. 7) compositions. Again, the near uniform chemical composition between the products yields no significant deviations in viscosity. These model curves allow for extrapolation to lower temperatures approaching the glass transition temperature ($T_g$). The predicted values of $T_g$ corresponding to the temperature where $\eta \sim 10^{12} \text{ Pa s vary}$

between 660 °C and 676 °C below which crystallisation cannot occur.

**Origins of inflated pyroclasts.** The proximal deposits at Tseax feature an intriguing juxtaposition of scoria and inflated pyroclasts that is not easily ascribed to a standard model for a thermally insulated lava fountain. In this standard model, temperature variations are predominantly controlled by the amount of entrained ambient air. A cooler fountain exterior, with a large amount of entrained ambient air, envelopes a hotter, thermally insulated, incandescent fountain interior where $>T_g$ processes can occur⁵⁷⁻⁵⁹. In this conventional scenario, pyroclasts exiting the fountain cool rapidly as they pass through the cooler exterior; the extent of cooling being dictated by their size and trajectory³⁰,³²,⁴⁰,⁵⁸. The large pyroclasts (e.g., spatter) may exit the fountain hot ($>T_g$) allowing dynamic processes to continue modifying their interiors (i.e., coalescence, vesication)¹²,¹³. Smaller pyroclasts tend to be cooled to below $T_g$ except for those with very short trajectories as found in very proximal deposits.

In the highly proximal tephra deposits, the inflated pyroclasts represent a substantially smaller volume fraction of the deposit than does scoria. However, the inflated pyroclasts are found over
<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /><figcaption>**Fig. 5 Bulk vesicularity histograms.** Data shown for **a** the proximal scoria and **b** the inflated pyroclasts. The number of samples (n) measured is reported for each dataset. The box and whisker plots display the data range (whiskers) and the first quartile, median and third quartile (box). The corresponding bulk vesicularity can be found in the Supplementary Data 1 file.</figcaption></figure>

a wide range of sizes (Fig. 2). They are not restricted to the largest sized pyroclasts but can be found as small as ~0.2 cm (post inflation). Furthermore, there is clear macro- and microscale evidence for the bubble growth, which drives the inflation, to occur within the deposit, after deposition. Thus, the inflated pyroclasts can only occur under one of two distinct eruption scenarios: (1) extreme thermal heterogeneity within the lava fountain or (2) abrupt transitions from lava fountaining to discrete bubble bursts. These two hypotheses will now be evaluated for the Tseax eruption.

Extreme thermal heterogeneities within the parent lava fountain could account for the contrasting thermal histories of the scoria and inflated pyroclasts, where a small volume of thermally insulated lava proceeds to form the inflated pyroclasts. To assess the feasibility of this scenario at Tseax, we used a 1D transient heat conduction model for a sphere³⁰,⁶⁰ (see Methods). The observed field deposits require that all inflated pyroclasts, irrespective of their size, landed and inflated in the deposit above Tg. Furthermore, no breadcrust textured surfaces were observed

<figure><img src="image_10.png" /><figcaption>**Fig. 6 X-ray computed tomography derived images.** **a** Image slice through regular scoria found adjacent to the inflated pyroclasts within the same deposit. This scoria clast has a bulk vesicularity of 76.5%. **b** Image slice through an inflated pyroclast. This inflated pyroclast has a bulk vesicularity of 89.7%. Note the difference in scale bar length between the two images.</figcaption></figure>

indicating that the entire inflated pyroclast, including its surface was kept well above Tg during transport and emplacement. Given that pyroclasts loose heat as function of their size, we modelled the cooling of a pyroclast 2 mm in diameter exposed to different ambient temperatures (Fig. 8a). This small particle diameter is consistent with the smallest inflated pyroclasts observed in the deposit (Table S2) and therefore provides the most robust hypothesis test. As shown in Fig. 8a, pyroclasts are only able to stay hot enough to form inflated pyroclasts when the ambient temperature is above the glass transition temperature. When the ambient temperature drops below the glass transition temperature (cooler lines; Fig. 8a), the pyroclast becomes effectively solid in <0.2 s and therefore cannot inflate in the deposit.
<table><caption>Table 1 Experimental measurements of viscosity of Tseax lava based on remelt of sample TS-S71 (see Table S1) compared to values predicted by VFT models (i.e., A<sub>VFT</sub>, B<sub>VFT</sub>, C<sub>VFT</sub>) for temperature dependence of melt viscosity.</caption><thead><tr><th rowspan="2">State</th><th colspan="3">Experimental measurements</th><th colspan="2">Predicted log η (Pa s)</th></tr><tr><th>T (°C)</th><th>10,000/T(K)</th><th>log η (Pa s)</th><th>GRD<sup>a</sup></th><th>R<sub>Global</sub><sup>b</sup></th></tr></thead><tbody><tr><td rowspan="15">Superliquidus melt</td><td>1531</td><td>5.54</td><td>0.302</td><td>0.32</td><td>0.21</td></tr><tr><td>1508</td><td>5.62</td><td>0.361</td><td>0.42</td><td>0.31</td></tr><tr><td>1484</td><td>5.69</td><td>0.423</td><td>0.52</td><td>0.41</td></tr><tr><td>1460</td><td>5.77</td><td>0.488</td><td>0.62</td><td>0.51</td></tr><tr><td>1435</td><td>5.85</td><td>0.561</td><td>0.74</td><td>0.62</td></tr><tr><td>1412</td><td>5.94</td><td>0.640</td><td>0.85</td><td>0.74</td></tr><tr><td>1388</td><td>6.02</td><td>0.722</td><td>0.97</td><td>0.85</td></tr><tr><td>1363</td><td>6.11</td><td>0.810</td><td>1.10</td><td>0.98</td></tr><tr><td>1340</td><td>6.20</td><td>0.932</td><td>1.23</td><td>1.11</td></tr><tr><td>1315</td><td>6.29</td><td>1.03</td><td>1.37</td><td>1.24</td></tr><tr><td>1292</td><td>6.39</td><td>1.13</td><td>1.52</td><td>1.39</td></tr><tr><td>1268</td><td>6.49</td><td>1.24</td><td>1.67</td><td>1.54</td></tr><tr><td rowspan="15">Subliquidus melt (with crystallisation)</td><td>1244</td><td>6.59</td><td>1.35</td><td>A<sub>VFT</sub></td><td>-4.55</td><td>-4.55</td></tr><tr><td>1244</td><td>6.59</td><td>1.43</td><td>B<sub>VFT</sub></td><td>5905</td><td>5785</td></tr><tr><td>1220</td><td>6.70</td><td>1.30</td><td>C<sub>VFT</sub></td><td>592</td><td>590</td></tr><tr><td>1220</td><td>6.70</td><td>1.56</td><td>T<sub>g</sub> (°C)</td><td>676</td><td>667</td></tr><tr><td>1220</td><td>6.70</td><td>1.78</td><td>Fragility</td><td>44.0</td><td>44.5</td></tr><tr><td>1220</td><td>6.70</td><td>2.01</td><td></td><td></td><td></td></tr><tr><td>1196</td><td>6.81</td><td>2.19</td><td></td><td></td><td></td></tr><tr><td>1196</td><td>6.81</td><td>2.56</td><td></td><td></td><td></td></tr><tr><td>1172</td><td>6.92</td><td>2.88</td><td></td><td></td><td></td></tr><tr><td>1172</td><td>6.92</td><td>2.60</td><td></td><td></td><td></td></tr><tr><td>1148</td><td>7.04</td><td>2.63</td><td></td><td></td><td></td></tr><tr><td>1124</td><td>7.16</td><td>2.88</td><td></td><td></td><td></td></tr><tr><td>1124</td><td>7.16</td><td>3.42</td><td></td><td></td><td></td></tr></tbody></table>

Raman spectrum of remelted Tseax glass (TS-S71) has a normalised peak ratio value of 0.1595.

<sup>a</sup>Model values predicted by GRD<sup>54</sup> model.

<sup>b</sup>Values predicted by Raman-based model<sup>55,56</sup> for melt viscosity.

Thus, for a standard thermally insulated lava fountain<sup>30,61,62</sup>, only pyroclasts that reside in the hot (≥700°C) interior for the entirety of the transport path can form inflated pyroclasts. For the modelled conditions at Tseax this is highly unlikely; firstly, a 2 mm pyroclast would not be able to settle vertically within an opposing gas fountain velocity of ~100 m s<sup>-1</sup>. Secondly, pyroclasts are likely to enter the cooler exterior either at the fountain top or edges, from this point transport and in-deposit inflation must occur on timescales < 0.2 s, this is infeasible. However, if ambient gas temperatures surrounding the fountain are locally raised to temperatures above T<sub>g</sub> for short periods of time, pyroclasts could exit the fountain and reach the deposit above T<sub>g</sub>, potentially forming inflated pyroclasts. Such high gas temperatures (≥700°C) surrounding lava fountains are unusual, although they have been documented<sup>33</sup> during the 2018 eruption of Kilauea and are most likely when mass eruption rates are high, the vent is wide, and hosts ponded lava. We therefore cannot rule out the production of inflated pyroclasts via extreme thermal anomalies at other localities, but suggest that it is highly unlikely at Tseax.

Hawaiian style fountains are inherently unsteady and during a single eruption episode the fountain height can vary substantially<sup>18</sup>. These intra-episode height variations are often due to changes in the mass eruption rate<sup>40,58</sup> but other factors such as proximal lava ponding around the vent can also play a role<sup>34,58</sup>. Inflated pyroclasts may form instead of scoria if the lava fountain temporarily pauses or stops and the eruptive activity transitions to low-energy bubble bursts<sup>63</sup> (similar to bubble bursting at a lava lake or very weak Strombolian activity). Such eruptive behaviour transitions are not unprecedented and similar abrupt (but regular) changes in eruptive activity have been witnessed during the 2021 eruption of Fagradalsfjall, Iceland and

during the 2018 eruption of Kilauea, Hawaii<sup>21</sup>. The cause of such transitions remains uncertain, but it is likely due to reorganisation of gas in the shallow volcanic plumbing system. This forms our second hypothesis and is again tested through thermal modelling (Fig. 8b).

During low-energy bubble bursts, the differential velocity between the ejected clast and the surrounding air is low compared to lava fountain events. Lower differential velocities correspond to lower heat transfer coefficients meaning that cooling is less effective and pyroclasts can remain hot, and above the glass transition temperature, for longer. These conditions are shown in Fig. 8b, where a 2 mm diameter pyroclast (including the exterior surface) remains above T<sub>g</sub> for ~1 s: an order of magnitude longer than the lava fountain scenario (cf. Fig. 8a). Furthermore, if the full range of inflated pyroclast sizes are considered (sequential coloured lines in Fig. 8b), we show that the majority of inflated pyroclasts required 10's or 100's of seconds to reach T<sub>g</sub>—ample time for the pyroclasts to be ejected, reach the deposit, and for the bubbles to grow and coalesce to produce the textures observed. We therefore contend that repeated rapid transitions from lava fountaining (producing the regular scoria) to discrete low-energy bubble bursts provided the in-deposit conditions required to form the inflated pyroclasts at Tseax (Fig. 9).

The melt properties also played a role in the formation of inflated pyroclasts. The Tseax magma has a low viscosity and a predicted abnormally low glass transition temperature, T<sub>g</sub> (660–676 °C; Fig. 7). These physical properties are key to the formation of inflated pyroclasts. Firstly, the low melt viscosity allows for bubble growth and coalescence to occur on timescales faster than pyroclast cooling. Secondly, a low T<sub>g</sub> increases ΔT = T<sub>e</sub> – T<sub>g</sub> thereby extending the time available for pyroclast modification after primary magmatic fragmentation. Thirdly, the
<figure><img src="image_11.png" /><figcaption>**Fig. 7 Arrhenius plot describing the viscosity, η (Pa s)-temperature, T (K) relationship for natural and remelted samples of Tseax lava.**</figcaption></figure>

Experimentally measured values of melt viscosity (white symbols) above liquidus temperatures ($T_L$ -1245 °C) and the viscosity of melt-crystal mixtures (grey symbols). The GRD model⁵⁴ values of melt viscosity for the composition of the natural inflated pyroclast are shown by the green labelled solid line (Table S1). Also shown are model viscosity curves for the remelted glass sample (TS-S71) predicted by the GRD model (thick dashed line) and the Raman-based model⁵⁵,⁵⁶ (fine dashed line). Model values of $T_g$ taken as the temperature corresponding to ~ $10^{12}$ Pa s vary between 660 and 676 °C (grey shaded field). The corresponding rheometry data can be found in Table 1.

iron-rich pyroclasts were subject to rapid cooling—conditions that favour microlite crystallisation⁵². For these melt compositions it has been shown⁵³ that plagioclase and oxides will rapidly crystalize at temperatures just below the Tseax eruption temperature, $T_e$. Supported by our thin section observations, we suggest that extensive plagioclase microlites formed towards the end of pyroclast inflation to generate the non-glassy surface texture observed (Fig. 3a). Furthermore, these microlites provided a supporting framework, and a rise in viscosity (Fig. 7), that allowed the inflated pyroclasts to retain their shape, rather than collapsing, as interior bubbles breached the exterior surface of the pyroclast allowing egress of pressurised gasses (Fig. 3b–d).

**The critical conditions required to produce inflated pyroclasts.**

Here, we have documented the conditions that lead to inflated pyroclasts at Tseax but contend that inflated pyroclasts could occur at a range of volcanoes worldwide, provided a series of critical conditions are met. Thus, when recognised in deposits, inflated pyroclasts indicate the following eruption conditions:

* Pyroclast velocities were small, and certainly lower than the critical impact velocity causing pyroclasts to highly deform on impact and further fragment into a series of droplets¹³.

* The eruption temperatures were hot enough and/or the flight trajectories were short enough and/or in-flight cooling was limited such that clasts land in the deposit with temperatures > $T_g$.

* Melt viscosities were low to allow bubbles to grow, relax and coalesce on short timescales.

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

**Fig. 8 Heat conduction model for pyroclasts.** The coloured lines track the surface temperature of the pyroclast. The black dashed lines represent the corresponding ambient air temperature, and the grey shaded region denotes the temperatures below the glass transition temperature, $T_g$. Both $T_g$ values of 660 °C and 676 °C are marked to represent the range determined by high-temperature rheometry. a Cooling of a 2 mm diameter pyroclast within a lava fountain, the ambient temperatures ranging from 30°C (light yellow) to 900 °C (dark red) represent the thermal zonation in the fountain, from the cold edge to the hotter, insulated interior. b Cooling of 2 mm (dark green), 10 mm (light green) and 40 mm (blue) diameter pyroclasts generated from a low-velocity bubble burst. Note, the change of the x-axis to a log scale.

* The eruption/episode volume was small and the deposit accumulation rate low, below the threshold that would produce rheomorphic lava flows¹⁴.

* If vesicles are observed intersecting the pyroclast exterior, microlites and/or nanolites may have formed preventing interior bubble collapse.

* If scoria and inflated pyroclasts are intermixed within the deposit (rather than a single horizon of inflated pyroclasts), then transitions between fountaining and bubble bursts occurred multiple times.

We have shown that when juxtaposed with normal scoria, inflated pyroclasts are indicators of abrupt transitions in eruptive style—repeated rapid transitions between lava fountaining activity and low-energy bubble bursts— transitions that could be considered when hazard planning and managing access to these (basaltic) eruptions that are popular with tourists.
<figure><img src="image_13.png" /><figcaption>**Fig. 9** Summary cartoon of the eruptive processes that operated at the Satellite cone, Tseax volcano. The diagrams are schematically drawn and are not to scale. **a** Standard lava fountain producing regular scoria. **b** Repeated transitions to a low-energy bubble burst ejecting low viscosity lava fragments that land in close proximity to the vent. **c** Close up of the proximal deposits showing the interior bubble growth and in-deposit inflation of the pyroclasts. **d** Sketch of a cross-section through an inflated pyroclast with a crystal framework preventing interior bubble collapse.</figcaption></figure>

## Methods

**Deposit sampling.** The samples were collected during three field seasons (Summers of 2016, 2017 and 2019) at Tseax volcano (55.11085 °N, 128.89944 °W), in northwest British Columbia, Canada. The inflated pyroclasts sampled for this study were collected within a 5 m radius of the summit of the Satellite cone (Fig. 1). No deposit stratigraphy was identified so the inflated pyroclasts were sampled from the present-day surface, or from the upper 0.5 m of deposit when an erosional surface exposed a small section. If inflated pyroclasts were too agglutinated to be successfully removed from the deposit a larger sample was taken for later, careful separation in the lab. During transport, the samples were individually wrapped in bubble wrap to prevent breakage, comminution, and surface modifications.

**Scanning electron microscopy analysis.** Six representative inflated pyroclasts were cut into ~2 cm pieces and mounted on scanning electron microscope (SEM) stubs using adhesive carbon tape to reduce sample charging. All sample stubs were then carbon coated using a sputter coater. Analysis and image acquisition were then performed using a Philips XL30 SEM in scanning electron mode with a 15 kV accelerating voltage, a 35 µA beam current, and an average working distance of 11 mm.

**Electron microscopy analysis.** The major element compositions of glasses derived from remelting of whole rock sample (TS-S71; see viscosity) were measured using a

Cameca SX100 electron probe microanalyser (EPMA) at the Department for Earth and Environmental Sciences, LMU Munich, Germany. Operating conditions were an acceleration voltage of 15 kV and beam current of 4 nA, using a defocused 10 µm beam to minimise alkali loss. The peak time of our analyses was 10 s and a background time of 5 s (on each side of the peak) for all elements. Albite, periclase, apatite, wollastonite, bustamite, $Fe_2O_3$ and ilmenite were used as standards. Standard deviations based on replicate measurements are < 2.5% for all analysed elements. Precision and accuracy were checked by analysing the reference glasses VG-2 (basalt) and VG-568 (rhyolite)$^{64,65}$ at the start of each analytical session. No significant alkali loss was observed relative to tabulated values of standards.

**X-Ray fluorescence.** X-Ray fluorescence measurements of the bulk major element chemical compositions of representative samples were measured by ALS Minerals in North Vancouver, British Columbia (Table S1). Ferrous iron (FeO wt.%) content was measured by volumetric analysis. The detection limit is 0.01 wt.% for all the major oxides. The sample suite includes: (i) samples of highly oxidised tephra from the satellite cone (SC) where the inflated pyroclasts are found, (ii) tephra from the main cone (MC) which overlies SC tephra, (iii) a sample of lava (TS-S71), (iv) an inflated pyroclast (IP211015A) and (v) the glass produced by remelting the lava in support of the viscosity measurements.
**Dense rock equivalent, or skeletal density measurements.** Using a mortar and pestle, several inflated pyroclasts were crushed to a fine powder, finer than the smallest vesicles. The dense rock equivalent, or skeletal density of the inflated pyroclasts was calculated using an analytical balance to measure mass and a Micrometrics Acupyc II 1340 helium pycnometer to measure volume. Five aliquots of the fine pyroclast powder were measured ten times for mass and volume, then plotted as mass against volume. A linear regression was then fitted through the five datasets and the origin, with the slope determining the density. These data are shown in Fig. S3.

**Bulk density measurements.** A total of 40 inflated pyroclasts and 50 scoria clasts were measured for bulk density using the Archimedes techniques outlined in Houghton & Wilson⁶⁶. The bulk vesicularity was calculated from the bulk density using the dense rock equivalent density measured independently by helium pycnometry (as detailed above). Three repeat measurements on the same inflated pyroclast revealed a maximum absolute uncertainty of ± 1% vesicularity.

**Viscosity measurements.** We measured the anhydrous, 1 atmosphere viscosity of the remelted TS-S71 sample (Tables 1 and S1). Viscosity of the Tseax melt was measured using standard rheological measurements. High-temperature viscosity measurements were made, after normal calibration procedures, at the Department for Earth and Environmental Sciences, LMU Munich, Germany. A concentric cylinder apparatus was used for sample homogenisation and determination of anhydrous liquid viscosities (10³⁰.³⁰–10³⁴.⁴² Pa s) at both superliquidus (1531–1267 °C) and subliquidus (1244–1124 °C) temperatures.

The measured liquid viscosity of the Tseax melt is one of the lowest values amongst those measured for natural samples⁵⁴. The subliquidus measurements show that the onset of an apparent crystallisation (<1244 °C) has a measurable rheological effect manifest by a drift in viscosity values with lower temperatures. The high temperature marking the onset of crystallisation is not typical of other basaltic melts including MORB, OIB, or more alkaline, iron-poorer basalts⁶⁷.

We used the composition of the remelted glasses of TS-S71 to calculate the temperature dependence of the liquid viscosity using the GRD model⁵⁴. The corresponding VFT function⁶⁸⁻⁷⁰ is:

$$ \log \eta = -4.55 + \frac{5904.6}{T(K) - 591.8} \qquad (1) $$

which predicts the measured values well (cf. Fig. 7).

**X-ray computed tomography.** One ~4.5 cm sample containing both scoria and an adhering inflated pyroclast was scanned using a Scanco Medical µCT100 device housed at the Centre for High-Throughput Phenogenomics, Faculty of Dentistry, University of British Columbia. Acquisition conditions were a 55 kV accelerating voltage, a 200 µA current and a 0.5 mm aluminium filter. The total exposure time was 9.7 h. The voxel dimension (i.e., the resolution) is 11.4 µm.

**Thermal models.** We modelled the cooling of inflated pyroclasts using a 1D transient heat conduction model based on the MATLAB code by Recktenwald⁶⁰ as used by Porritt et al.³⁰. Here the pyroclasts are modelled as spherical and lose heat from their surfaces by forced convection, $F_c$:

$$ F_c = h_c(T_s - T_a) \qquad (2) $$

where $h_c$ is the heat transfer coefficient, $T_s$ is the surface temperature of the pyroclast and $T_a$ is the ambient temperature. Given that small microphenocrysts are observed in the inflated pyroclasts, we set $T_s$ to be 1244 °C, just below the liquidus temperature. Various values for $T_a$ were used ranging from 900 °C and 30 °C which represent the core of a vigorous lava fountain and background air temperature, respectively. The heat transfer coefficient, $h_c$ is encapsulated within the Nusselt number, Nu, a dimensionless heat loss parameter defined as:

$$ \text{Nu} = \frac{h_c \cdot d}{k_a} \qquad (3) $$

where d is the pyroclast diameter and $k_a$ is the thermal conductivity of the air (0.0257 Wm⁻¹ K⁻¹). We calculated Nu based on previous experimental work that used spheres⁷¹ and natural volcanic pyroclasts⁷². The equations are based on two further dimensionless groups, the particle Reynolds number, Re, and the Prandtl number, Pr. First, the particle Reynolds number is given by:

$$ \text{Re} = \frac{\Delta V d \rho_a}{\mu_a} \qquad (4) $$

where ΔV is the differential velocity between the ejected pyroclast and the surrounding air, $ρ_a$ is the air density and $μ_a$ is the air viscosity. The temperature dependence of air density and viscosity is assumed to be negligible, and we use constant values of 1.225 kg m⁻³ and 1.75 × 10⁻⁵ Pa s, respectively. We assume the average pyroclast differential velocity, ΔV to be ~100 m s⁻¹ for a lava fountain³⁰,⁷³ and ~0.1 m s⁻¹ for a bubble burst event⁶³ considering drag and the entire ballistic

motion. Second, the Prandtl number is given by:

$$ \text{Pr} = \frac{\mu_a C_p}{k_a} \qquad (5) $$

where $C_p$ is the air heat capacity, here taken as a constant value of 1005 J kg⁻¹ K⁻¹.

Previous experimental work⁷¹,⁷² found that for low Reynolds numbers, < 2.0 × 10⁵, the Nusselt number is given by the following relationship:

$$ \text{Nu} = 2 + x \text{Re}^{1/2} \text{Pr}^{1/3} \qquad (6) $$

where $x = 2.2 \times 10^{-4} \rho + 0.31$ and $\rho$ is the density of the pyroclast, taken here to be the mean bulk density of the proximal scoria (967 kg m⁻³; Fig. 5). The largest particle Reynolds number used in Fig. 8 was 1.4 × 10⁴, thus Eq. 6 is valid. The cooling model also uses the melt thermal conductivity and diffusivity; here we assume these to be constant and use typical values of 2 W m⁻¹ K⁻¹ and 1 × 10⁻⁵ m² s⁻¹, respectively⁷⁴⁻⁷⁸.

## Data availability

The rheometry, bulk density, geochemical and field data generated or analysed in this study are provided in the main article, the Supplementary Information, and Supplementary Data 1. The original X-Ray Computed Tomography data are available from the authors on request.

## Code availability

The 1D transient heat conduction code is written in MATLAB and adapted from Recktenwald⁶⁰. Our adapted version can be accessed via Github/Zenodo (https://doi.org/10.5281/zenodo.6381629).

Received: 19 August 2021; Accepted: 28 April 2022;

Published online: 20 May 2022

## References

1. Valentine, G. A. & Gregg, T. K. P. Continental basaltic volcanoes- processes and problems. *J. Volcanol. Geotherm. Res.* 177, 857–873 (2008).

2. Sigurdsson, H. Volcanic episodes and rates of volcanism. In *Encyclopedia of Volcanoes* 271–279 (Academic Press San Diego, CA, 2000).

3. Neal, C. A. et al. The 2018 rift eruption and summit collapse of Kilauea Volcano. *Science* 363, 367–374 (2019).

4. Heggie, T. W. Geotourism and volcanoes: health hazards facing tourists at volcanic and geothermal destinations. *Travel Med. Infect. Dis.* 7, 257–261 (2009).

5. Andronico, D., Scollo, S. & Cristaldi, A. Unexpected hazards from tephra fallouts at Mt Etna: the 23 November 2013 lava fountain. *J. Volcanol. Geotherm. Res.* 304, 118–125 (2015).

6. Giordano, D. et al. Thermo-rheological magma control on the impact of highly fluid lava flows at Mt. Nyiragongo. *Geophys. Res. Lett.* 34, L06301 (2007).

7. Jones, T. J., Reynolds, C. D. & Boothroyd, S. C. Fluid dynamic induced break-up during volcanic eruptions. *Nat. Commun.* 10, 3828 (2019).

8. Giordano, D. & Dingwell, D. Viscosity of hydrous Etna basalt: implications for Plinian-style basaltic eruptions. *Bull. Volcanol.* 65, 8–14 (2003).

9. Coltelli, M., Del Carlo, P. & Vezzoli, L. Discovery of a Plinian basaltic eruption of Roman age at Etna volcano, Italy. *Geology* 26, 1095–1098 (1998).

10. Turbeville, B. N. Tephra fountaining, rheomorphism, and spatter flow during emplacement of the Pitigliano Tuffs, Latera caldera, Italy. *J. Volcanol. Geotherm. Res.* 53, 309–327 (1992).

11. Parcheta, C. E., Houghton, B. F. & Swanson, D. A. Hawaiian fissure fountains 1: decoding deposits—episode 1 of the 1969–1974 Mauna Ulu eruption. *Bull. Volcanol.* 74, 1729–1743 (2012).

12. Jones, T. J., Houghton, B. F., Llewellin, E. W., Parcheta, C. E. & Höltgen, L. Spatter matters – distinguishing primary (eruptive) and secondary (non-eruptive) spatter deposits. *Sci. Rep.* 8, 9179 (2018).

13. Sumner, J. M., Blake, S., Matela, R. J. & Wolff, J. A. Spatter. *J. Volcanol. Geotherm. Res.* 142, 49–65 (2005).

14. Rader, E. & Geist, D. Eruption conditions of spatter deposits. *J. Volcanol. Geotherm. Res.* 304, 287–293 (2015).

15. Tesche, M. et al. Volcanic ash over Scandinavia originating from the Grimsvötn eruptions in May 2011. *J. Geophys. Res. Atmos.* 117, D09201 (2012).

16. Giordano, D., Polacci, M., Papale, P. & Caricchi, L. Rheological control on the dynamics of explosive activity in the 2000 summit eruption of Mt. Etna. *Solid Earth* 1, 61–69 (2010).
17. Taddeucci, J. et al. Fracturing and healing of basaltic magmas during explosive volcanic eruptions. *Nat. Geosci.* **14**, 248–254 (2021).

18. Houghton, B. F. et al. Stronger or longer: discriminating between Hawaiian and Strombolian eruption styles. *Geology* **44**, 163–166 (2016).

19. Gaudin, D. et al. Integrating puffing and explosions in a general scheme for Strombolian-style activity. *J. Geophys. Res. Solid Earth* **122**, 1860–1875 (2017).

20. Spina, L., Cannata, A., Morgavi, D. & Perugini, D. Degassing behaviour at basaltic volcanoes: New insights from experimental investigations of different conduit geometry and magma viscosity. *Earth Sci. Rev.* **192**, 317–336 (2019).

21. Houghton, B. F. et al. The birth of a Hawaiian fissure eruption. *J. Geophys. Res. Solid Earth* **126**, e2020JB020903 (2020).

22. Orr, T. R. et al. Kilaua'e's 5–9 March 2011 Kamoamoa fissure eruption and its relation to 30+ years of activity from Pu'u'O'o. In *Hawaiian Volcanoes: From Source to Surface* **208**, 393 (John Wiley & Sons, 2015).

23. Patrick, M. R. et al. The cascading origin of the 2018 Kilauea eruption and implications for future forecasting. *Nat. Commun.* **11**, 1–13 (2020).

24. Dingwell, D. B. & Webb, S. L. Structural relaxation in silicate melts and non-Newtonian melt rheology in geologic processes. *Phys. Chem. Miner.* **16**, 508–516 (1989).

25. Dingwell, D. B. Volcanic dilemma: flow or blow? *Science* **273**, 1054 (1996).

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

27. Jones, T. J., Cashman, K. V., Liu, E. J., Rust, A. C. & Scheu, B. Magma fragmentation: a perspective on emerging topics and future directions. *Bull. Volcanol.* **84**, 1–8 (2022).

28. Namiki, A. & Manga, M. Transition between fragmentation and permeable outgassing of low viscosity magmas. *J. Volcanol. Geotherm. Res.* **169**, 48–60 (2008).

29. Jones, T. J., Russell, J. K., Brown, R. J. & Hollenderdon, L. Melt stripping and agglutination of pyroclasts during the explosive eruption of low viscosity magmas. *Nat. Commun.* **13**, 11–11 (2022).

30. Porritt, L. A., Russell, J. K. & Quane, S. L. Pele's tears and spheres: examples from Kilauea Iki. *Earth Planet. Sci. Lett.* **333**, 171–180 (2012).

31. Stovall, W. K., Houghton, B. F., Gonnermann, H., Fagents, S. A. & Swanson, D. A. Eruption dynamics of Hawaiian-style fountains: the case study of episode 1 of the Kilaua Iki 1959 eruption. *Bull. Volcanol.* **73**, 511–529 (2011).

32. Stovall, W. K., Houghton, B. F., Hammer, J. E., Fagents, S. A. & Swanson, D. A. Vesiculation of high fountaining Hawaiian eruptions: episodes 15 and 16 of 1959 Kilaua Iki. *Bull. Volcanol.* **74**, 441–455 (2012).

33. Namiki, A., Patrick, M. R., Manga, M. & Houghton, B. F. Brittle fragmentation by rapid gas separation in a Hawaiian fountain. *Nat. Geosci.* **14**, 1–6 (2021).

34. Jones, T. J., Llewellin, E. W., Houghton, B. F., Brown, R. J. & Vye-Brown, C. Proximal lava drainage controls on basaltic fissure eruption dynamics. *Bull. Volcanol.* **79**, 81 (2017).

35. Bruce, P. M. & Huppert, H. E. Thermal control of basaltic fissure eruptions. *Nature* **342**, 665–667 (1979).

36. Thorarinsson, S., Steinthorsson, S., Einarsson, T., Kristmannsdottir, H. & Oskarsson, N. The eruption on Heimaey, Iceland. *Nature* **241**, 372–375 (1973).

37. Brown, R. J., Thordarson, T., Self, S. & Blake, S. Disruption of tephra fall deposits caused by lava flows during basaltic eruptions. *Bull. Volcanol.* **77**, 90 (2015).

38. Reynolds, P., Brown, R. J., Thordarson, T. & Llewellin, E. W. The architecture and shallow conduits of Laki-type pyroclastic cones: insights into a basaltic fissure eruption. *Bull. Volcanol.* **78**, 1–18 (2016).

39. Thordarson, T. & Self, S. The Laki (Skaftár Fires) and Grimsvötn eruptions in 1783–1785. *Bull. Volcanol.* **55**, 233–263 (1993).

40. Head, J. W. & Wilson, L. Basaltic pyroclastic eruptions: influence of gas-release patterns and volume fluxes on fountain structure, and the formation of cinder cones, spatter cones, rootless flows, lava ponds and lava flows. *J. Volcanol. Geotherm. Res.* **37**, 261–271 (1989).

41. Edwards, M. J., Pioli, L., Harris, A. J. L., Gurioli, L. & Thivet, S. Magma fragmentation and particle size distributions in low intensity mafic explosions: the July/August 2015 Piton de la Fournaise eruption. *Sci. Rep.* **10**, 1–14 (2020).

42. Sutherland Brown, A. Aiyansh lava flow, British Columbia. *Can. J. Earth Sci.* **6**, 1460–1468 (1969).

43. Hanson, G. Reconnaissance between Skeena River and Stewart, British Columbia. In *Geological Survey of Canada, Summary Report* 29–45, https://doi.org/10.4095/103192 (1923).

44. Le Moigne, Y., Williams-Jones, G., Russell, J. K. & Quane, S. Physical volcanology of Tseax Volcano, British Columbia, Canada. *J. Maps* **16**, 363–375 (2020).

45. Williams-Jones, G. et al. The age of the Tseax volcanic eruption, British Columbia, Canada. *Can. J. Earth Sci.* **57**, 1238–1253 (2020).

46. Hickson, C. J. & Edwards, B. R. Volcanoes and volcanic hazards in Canada. In *A Synthesis Of Geological Hazards in Canada* **548**, 1e248 (Natural Resources Canada, 2001).

47. Nisga'a Tribal Council and B.C. Parks Committee. *Anhluut'ukwsim Laxmihl Angwinga'asankswhl Nisga'a: Nisga'a Memorial Lava Bed Park Master Plan*. (Grant Copeland & Associates, 1997).

48. Edwards, B. R. & Russell, J. K. Distribution, nature, and origin of Neogene—Quaternary magmatism in the northern Cordillera volcanic province, Canada. *Geol. Soc. Am. Bull.* **112**, 1280–1295 (2000).

49. Van Der Heyden, P., Woodsworth, G. J. & Snyder, L. D. *Reconnaissance Geological Mapping in Southwest Nass River Map Area, British Columbia* (Natural Resources Canada, Geological Survey of Canada, 2000).

50. Evenchick, C. A. et al. Geology, Nass River, British Columbia. *Geol. Surv. Canada, Open File* 5705. https://doi.org/10.4095/225792 (2008).

51. Klug, C. & Cashman, K. V. Vesiculation of May 18, 1980, Mount St. Helens magma. *Geology* **29**, 468–472 (1994).

52. Di Genova, D. et al. In situ observation of nanolite growth in volcanic melt: a driving force for explosive eruptions. *Sci. Adv.* **6**, eabb0413 (2020).

53. Le Moigne, Y., Vigouroux, N., Russell, J. K. & Williams-Jones, G. Magmatic origins and storage conditions for the historic eruption of Tseax Volcano, British Columbia, Canada. *Chem. Geol.* **588**, 120648 (2022).

54. Giordano, D., Russell, J. K. & Dingwell, D. B. Viscosity of magmatic liquids: a model. *Earth Planet. Sci. Lett.* **271**, 123–134 (2008).

55. Giordano, D. & Russell, J. K. Towards a structural model for the viscosity of geological melts. *Earth Planet. Sci. Lett.* **501**, 202–212 (2018).

56. Giordano, D. et al. A calibrated database of Raman spectra for natural silicate glasses: implications for modelling melt physical properties. *J. Raman Spectrosc.* **51**, 1822–1838 (2020).

57. Parfitt, E. A., Wilson, L. & Neal, C. A. Factors influencing the height of Hawaiian lava fountains: implications for the use of fountain height as an indicator of magma gas content. *Bull. Volcanol.* **57**, 440–450 (1995).

58. Wilson, L., Parfitt, E. A. & Head, J. W. Explosive volcanic eruptions—VIII. The role of magma recycling in controlling the behaviour of Hawaiian-style lava fountains. *Geophys. J. Int.* **121**, 215–225 (1995).

59. Wolfe, E. W., Neal, C., Bamks, N. & Duggan, T. Puu Oo Eruption of Kilaua'e Volcano, Hawaii: Episodes 1 Through 20, January 3, 1983, Through June 8, 1984. (US. Geol. Surv. Prof. Pap. 1463, 1988).

60. Recktenwald, G. Transient, one-dimensional heat conduction in a convectively cooled sphere. *MATLAB Code* https://www.webcitation.org/60nDyv3Yy (2006).

61. Head, J. W. & Wilson, L. Lava fountain heights at Pu'u'O'o, Kilaua'e, Hawaii: indicators of amount and variations of exsolved magma volatiles. *J. Geophys. Res. Solid Earth* **92**, 13715–13719 (1987).

62. Wolff, J. A. & Sumner, J. M. Lava fountains and their products. In *Encyclopedia of Volcanoes* 321–329 (eds Sigurdsson, H., Houghton, B. F., McNutt, S. R., Rymer, H. & Stix, J.) (Encyclopedia of volcanoes. Academic, San Diego, 2000).

63. Mintz, B. G. et al. *Patterns of Bubble Bursting and Weak Explosive Activity in an Active Lava Lake*—Halema 'uma'u, Kilaua'e, 2015. (United States Geological Survey, 2021).

64. Helz, R. T., Clague, D. A., Mastin, L. G. & Rose, T. R. *Electron Microprobe Analyses of Glasses from Kilaua'e Tephra Units, Kilaua'e Volcano, Hawaii* (US Department of the Interior, US Geological Survey, 2014).

65. Jarosewich, E., Nelen, J. A. & Norberg, J. A. Reference samples for electron microprobe analysis. *Geostand. Newsl.* **4**, 43–47 (1980).

66. Houghton, B. F. & Wilson, C. J. N. A vesicularity index for pyroclastic deposits. *Bull. Volcanol.* **51**, 451–462 (1989).

67. Kolzenburg, S., Giordano, D., Hess, K. U. & Dingwell, D. B. Shear rate-dependent disequilibrium rheology and dynamics of basalt solidification. *Geophys. Res. Lett.* **45**, 6466–6475 (2018).

68. Fulcher, G. S. Analysis of recent measurements of the viscosity of glasses. *J. Am. Ceram. Soc.* **8**, 339–355 (1925).

69. Vogel, H. Das temperaturabhängigkeitsgesetz der viskosität von flüssigkeiten. *Phys. Z.* **22**, 645–646 (1921).

70. Tammann, G. & Hesse, W. Die Abhängigkeit der Viscosität von der Temperatur bie unterkühlten Flüssigkeiten. *Z. Anorg. Allg. Chem.* **156**, 245–257 (1926).

71. Achenbach, E. Heat transfer from spheres up to Re=6x10^6. In *International Heat Transfer Conference Digital Library* p. 341–346. (Hemisphere Publishing Corp., 1978).

72. Stroberg, T. W., Manga, M. & Dufek, J. Heat transfer coefficients of natural volcanic clasts. *J. Volcanol. Geotherm. Res.* **194**, 214–219 (2010).

73. Taddeucci, J., Edmonds, M., Houghton, B., James, M. R. & Vergniolle, S. Hawaiian and Strombolian eruptions. In *The Encyclopedia of Volcanoes* (Second Edition) 485–503 (Elsevier, 2015).

74. Huppert, H. E. & Sparks, R. S. J. The generation of granitic magmas by intrusion of basalt into continental crust. *J. Petrol.* **29**, 599–624 (1988).

75. Jellinek, A. M. & Kerr, R. C. Mixing and compositional stratification produced by natural convection: 2. Applications to the differentiation of basaltic and silicic magma chambers and komatiite lava flows. *J. Geophys. Res. Solid Earth* **104**, 7203–7218 (1999).

76. Jones, T. J., Llewellin, E. W. & Mader, H. M. The use of a shear thinning polymer as a bubbly magma analogue for scaled laboratory experiments. *J. Volcanol. Geotherm. Res.* 106768. https://doi.org/10.1016/j.jvolgeores.2020.106768 (2020).
## Acknowledgements

The authors thank Dr. Deana Nyce of the Wilp Wilxoʻoskwhl Nisga'a Institute (WWNI) for her continued and enthusiastic support. This work was supported by the WWNI, the Nisga'a Lisims Government, a Royal Society research grant (RGS\R1\211015) awarded to T.J.J. and NSERC Discovery grants RGPIN-2018-03841 and RGPIN-2016-03953 awarded to J.K.R. and G.W.J. respectively. This research was carried out as part of B.C. Parks/ NLG permit 108312.

## Author contributions

The project was conceived by T.J.J., Y.L., J.K.R. and G.W.J. The sampling, methodology and analysis was completed by T.J.J., Y.L., J.K.R., D.G. and D.B.D. All authors contributed to the writing of the paper after an initial draft was prepared by T.J.J.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-30501-6.

**Correspondence** and requests for materials should be addressed to Thomas J. Jones.

**Peer review information** *Nature Communications* thanks Wendy Stovall and the anonymous reviewer(s) for their contribution to the peer review of this work.

**Reprints and permission information** is available at http://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_14.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 Creative 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 permitted 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) 2022

<footnote n="77">77. Dobran, F. *Volcanic Processes: Mechanisms in Material Transport* (Springer Science & Business Media, 2012).</footnote>

<footnote n="78">78. Bartlett, R. W. Magma convection, temperature distribution, and differentiation. *Am. J. Sci.* **267**, 1067–1082 (1969).</footnote>