Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

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

## Abstract

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.

## Introduction

Low viscosity magmas are the most common erupted on Earth1,2. 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 challenging3,4,5,6. The style and efficiency of magma fragmentation directly controls the type of products erupted and, in turn, the hazards they pose7,8,9. 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 hot10,11,12,13,14 for a longer time. Such lavas pose an additional risk to communities and infrastructure down flank3. 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 airspace15. Despite the diverse array of hazards that such mafic eruptions pose, too little is known about pyroclast transport and the associated eruption dynamics7,16,17 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 expression18,19,20. Eruption intensities and styles of activity can also fluctuate over short timescales and distances21,22—this presents challenges to hazard management23 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 locations21. Furthermore, the physical processes governing explosive eruption of low viscosity magmas are fundamentally different to those operating in the now relatively well-studied silicic systems24,25,26,27. Most commonly their fragmentation processes are fluid dynamic in nature such that breakage does not occur simply by decompression upon crossing the glass transition7,27,28,29. 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 relaxation30 transforms irregular pyroclast morphologies into more spherical shapes. The consequence is that pyroclast properties are susceptible to modification after the initial fragmentation12,29,30,31,32,33. 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;34,35,36 (2) early products are frequently buried or rafted by subsequent eruption episodes37,38,39 and (3) rheomorphic processes and in-flight secondary fragmentation processes destroy primary fragmentation features1,11,17,29,33,40,41. 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, Canada42,43,44,45, which is reported to have caused up to 2,000 fatalities making it the deadliest eruption in Canada’s history46,47. 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 Province48 and is notable for a 32-km-long basanite-trachybasalt lava flow (~0.5 km3 covering ~36 km2; Fig. 1). The tephra cone and valley-filling lavas overlay intercalated sandstone to siltstone and mudstone from the late Jurassic Bowser Lake group49,50. Detailed volcanological mapping shows that Tseax volcano consists of a ~65 m high, 350-400 m diameter tephra cone (2.8 ± 0.4 × 106 m3) situated within a 550-600 m diameter oxidised horseshoe-shaped spatter rampart44. Another smaller (~20 m high, ~55 m diameter, 1.7 ± 0.1 ×104 m3) 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 mounds44 (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 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. 2bd) 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 this dull appearance to be a consequence of microlitic51 and possibly further nanolitic52 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}_{{{{{\rm{g}}}}}}$$. Additionally, on some clasts, larger ($$\gtrsim$$ 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 titanomagnetite53. Plagioclase crystals show albite and polysynthetic 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 composition53. 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−3 with a mean bulk density of 967 kg m−3. Inflated pyroclast bulk density varied between 311 and 1235 kg m−3 with a mean bulk density of 755 kg m−3. Using Helium pycnometry, the dense rock equivalent of the pyroclasts was calculated to be 2957 kg m−3 (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 $$\gtrsim$$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 100.3–101.2 Pa s over the explored temperature range of 1267–1531 °C. We also measured the 1 atmosphere subliquidus viscosity where it rose to between 101.3–103.4 Pa s over the 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 (GRD54 and Raman55,56) 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 (Tg). The predicted values of Tg corresponding to the temperature where $$\eta$$ ~ 1012 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 >Tg processes can occur57,58,59. 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 trajectory30,32,40,58. The large pyroclasts (e.g., spatter) may exit the fountain hot (>Tg) allowing dynamic processes to continue modifying their interiors (i.e., coalescence, vesiculation)12,13. Smaller pyroclasts tend to be cooled to below Tg 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 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 sphere30,60 (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 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.

Thus, for a standard thermally insulated lava fountain30,61,62, only pyroclasts that reside in the hot ($$\gtrsim 700^\circ$$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−1. 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 Tg for short periods of time, pyroclasts could exit the fountain and reach the deposit above Tg, potentially forming inflated pyroclasts. Such high gas temperatures ($$\gtrsim 700^\circ$$C) surrounding lava fountains are unusual, although they have been documented33 during the 2018 eruption of Kīlauea 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 substantially18. These intra-episode height variations are often due to changes in the mass eruption rate40,58 but other factors such as proximal lava ponding around the vent can also play a role34,58. Inflated pyroclasts may form instead of scoria if the lava fountain temporarily pauses or stops and the eruptive activity transitions to low-energy bubble bursts63 (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 Kīlauea, Hawaii21. 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 Tg 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 Tg—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, Tg (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 Tg increases ΔT = TeTg thereby extending the time available for pyroclast modification after primary magmatic fragmentation. Thirdly, the iron-rich pyroclasts were subject to rapid cooling—conditions that favour microlite crystallisation52. For these melt compositions it has been shown53 that plagioclase and oxides will rapidly crystalize at temperatures just below the Tseax eruption temperature, Te. 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 droplets13.

• 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 > Tg.

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

• The eruption/episode volume was small and the deposit accumulation rate low, below the threshold that would produce rheomorphic lava flows14.

• 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.

## 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 microprobe 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, Fe2O3 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 Accupyc 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 & Wilson66. 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 (100.30−103.42 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 samples54. 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 basalts67.

We used the composition of the remelted glasses of TS-S71 to calculate the temperature dependence of the liquid viscosity using the GRD model54. The corresponding VFT function68,69,70 is:

$${{\log }}\,\eta =\,-4.55+\,\frac{5904.6}{T\left(K\right)-591.8}$$
(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 Recktenwald60 as used by Porritt et al.30. Here the pyroclasts are modelled as spherical and lose heat from their surfaces by forced convection, $${F}_{{{{{\rm{c}}}}}}$$:

$${F}_{{{{{\rm{c}}}}}}={h}_{{{{{\rm{c}}}}}}({T}_{{{{{\rm{s}}}}}}-{T}_{{{{{\rm{a}}}}}})$$
(2)

where $${h}_{{{{{\rm{c}}}}}}$$ is the heat transfer coefficient, $${T}_{{{{{\rm{s}}}}}}$$ is the surface temperature of the pyroclast and $${T}_{{{{{\rm{a}}}}}}$$ is the ambient temperature. Given that small microphenocrysts are observed in the inflated pyroclasts, we set $${T}_{{{{{\rm{s}}}}}}$$ to be 1244 °C, just below the liquidus temperature. Various values for $${T}_{{{{{\rm{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}_{{{{{\rm{c}}}}}}$$ is encapsulated within the Nusselt number, Nu, a dimensionless heat loss parameter defined as:

$${{{{{\rm{Nu}}}}}}=\frac{{h}_{{{{{\rm{c}}}}}}d}{{k}_{{{{{\rm{a}}}}}}}$$
(3)

where $$d$$ is the pyroclast diameter and $${k}_{{{{{\rm{a}}}}}}$$ is the thermal conductivity of the air (0.0257 Wm−1 K−1). We calculated Nu based on previous experimental work that used spheres71 and natural volcanic pyroclasts72. 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:

$${{{{{\rm{Re}}}}}}=\frac{\triangle {Vd}{\rho }_{{{{{\rm{a}}}}}}}{{\mu }_{{{{{\rm{a}}}}}}}$$
(4)

where $$\triangle V$$ is the differential velocity between the ejected pyroclast and the surrounding air, $${\rho }_{a}$$ is the air density and $${\mu }_{{{{{\rm{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−3 and 1.75 × 10−5 Pa s, respectively. We assume the average pyroclast differential velocity, $$\triangle V$$ to be ~100 m s−1 for a lava fountain30,73 and ~0.1 m s−1 for a bubble burst event63 considering drag and the entire ballistic motion. Second, the Prandtl number is given by:

$${{{{{\rm{Pr }}}}}}=\frac{{\mu }_{{{{{\rm{a}}}}}}{C}_{{{{{\rm{p}}}}}}}{{k}_{{{{{\rm{a}}}}}}}$$
(5)

where $${C}_{{{{{\rm{p}}}}}}$$ is the air heat capacity, here taken as a constant value of 1005 J kg−1 K−1.

Previous experimental work71,72 found that for low Reynolds numbers, $$<$$ 2.0 × 105, the Nusselt number is given by the following relationship:

$${{{{{\rm{Nu}}}}}}=2+x{{{{{{\rm{Re}}}}}}}^{1/2}{{{\Pr }}}^{1/3}$$
(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−3; Fig. 5). The largest particle Reynolds number used in Fig. 8 was 1.4 × 104, 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−1 K−1 and 1 × 10−5 m2 s−1, respectively74,75,76,77,78.

## 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 Recktenwald60. Our adapted version can be accessed via Github/Zenodo (https://doi.org/10.5281/zenodo.6381629).

## 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 Kīlauea 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. Kilauea’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 Kīlauea 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. & Hollendonner, L. Melt stripping and agglutination of pyroclasts during the explosive eruption of low viscosity magmas. Nat. Commun. 13, 1–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 Kilauea 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 Kilauea 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 (1989).

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’asanskwhl 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 Cordilleran 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 22, 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 Kilauea Volcano, Hawaii: Espisodes 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, Kilauea, 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, Kilauea, 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 Kilauea Tephra Units, Kilauea 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).

77. Dobran, F. Volcanic Processes: Mechanisms in Material Transport (Springer Science & Business Media, 2012).

78. Bartlett, R. W. Magma convection, temperature distribution, and differentiation. Am. J. Sci. 267, 1067–1082 (1969).

## 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 information

Authors

### 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.

### Corresponding author

Correspondence to Thomas J. Jones.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Peer review

### Peer review information

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

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Jones, T.J., Le Moigne, Y., Russell, J.K. et al. Inflated pyroclasts in proximal fallout deposits reveal abrupt transitions in eruption behaviour. Nat Commun 13, 2832 (2022). https://doi.org/10.1038/s41467-022-30501-6

• Accepted:

• Published:

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

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.