Outgassing behaviour during highly explosive basaltic eruptions

Explosivity of basaltic eruptions is related to the ef ﬁ ciency in which exsolved gas can separate from the melt during ascent, which is controlled by magma permeability. However, basaltic pyroclasts from eruptions of varying explosivity can show similar permeability, indicating a possible complex relationship between permeability, outgassing and eruptive style. Here, we provide 3D measurements of basaltic pyroclasts using X-ray microtomography. We investigate the role of permeability and outgassing on magma ascent dynamics by using a numerical conduit model. Among the permeable parameters, bubble number density and friction coef ﬁ cient largely affect explosivity. However, for fast ascending basaltic magmas, gas-melt coupling is maintained independent of magma permeability. In this case, magma storage conditions may determine eruptive style, driving rapid magma ascent, crystallisation and bubble nucleation, producing a highly explosive eruption. Monitoring parameters which reveal pre-eruptive conditions may assist hazard mitigation, particularly for basaltic systems which exhibit a wide range in eruptive style.


B
asaltic volcanoes exhibit a considerable range in eruptive style, where the intensity of an eruption directly influences the impact on proximal communities and the environment 1 .Eruptive style is governed by interdependent conduit processes during magma ascent, such as crystallisation, gas exsolution and expansion 1 .If outgassing, i.e., the process by which gas can decouple from a low viscosity melt, is facilitated during ascent, an effusive eruption is likely 2,3 .However, if coupling between the gas and melt phases is maintained, gas escape is restricted, increasing the potential for magma fragmentation and an explosive eruption 2 .Explosive-effusive transitions in volcanic activity have been linked to the development of permeability and outgassing efficiency during magma ascent 1,3 .
During magma ascent, outgassing can occur through different regimes, such as bubbly flow, slug flow, annular flow and permeable flow regimes [4][5][6][7][8][9][10] .The slug flow regime is commonly assumed as the main mechanism for transient Strombolian eruptions, whereas bubbly flow and annular flow, are usually associated with more continuous lava fountaining activity 7,10,11 .The permeable flow regime, instead, is not associated with a specific style of activity, as effusive, fountaining or highly explosive (sub-Plinian or Plinian) eruptions can occur 3 depending on the efficiency of gas flow through permeable pathways within the ascending magma.The efficiency of outgassing through the permeable flow regime during magma ascent is controlled by the formation and maintenance of permeable pore networks within the magma [12][13][14] .
Permeable outgassing can be described using Forchheimer's Law 2,6 : where P is the pressure, Z is the flow direction, U the volume flux of gas with respect to the melt (i.e., the gas-slip velocity), and µ g and ρ g the viscosity and density of the gas phase, respectively.The viscous (k v ) and inertial (k i ) permeabilities account for energy loss during gas flow through viscous and inertial dissipation mechanisms 15 .
Together they account for the control of pore network geometry on permeability, where k v describes the bulk pore space (porosity and pore aperture) and k i the pore distribution and tortuosity of connections 15 .These properties are linked to k v and k i through the Kozeny-Carman relationship 2,16 .Overall, permeability is affected by the geometry of the pore network and shows a correlation with porosity 17,18 .Constraints on the three-dimensional (3D) geometry of the pore space are required to effectively describe the relationship between porosity and permeability.X-ray computed microtomography (μCT) is an important technique for the 3D textural characterisation of volcanic rocks 16,[19][20][21] .For explosive eruptions, pyroclasts preserve information on degassing at fragmentation, assuming that they have not been subjected to deformation, and cooling of the pyroclast was fast enough to prevent further vesiculation and degassing following ejection from the vent 2,16,22 .
Textural studies of pyroclasts examining vesicle (pore) sizes, shapes, and their distributions can provide insight into bubble nucleation and growth during magma ascent 18,22,23 .However, magmas can show a considerable variation in permeability for the same value of porosity, encompassing several orders of magnitude [16][17][18]24,25 , as permeability also depends on the anisotropy, tortuosity, and the morphology of connections between pores 2,14,[25][26][27] . By usng X-ray μCT, we can directly visualise the volumetric pore network within the clast and provide a characterisation of the 3D geometric properties of permeable pathways such as connectivity, tortuosity, and the throat-pore size ratio 28 .
Several studies have investigated magma permeability and established complex porosity-permeability relationships which depend on pore network properties 6,16,17,25,27,29 .However, recent studies show that permeability data for pyroclasts from basaltic eruptions of varying explosivity overlap 10,29 , indicating that there may be a more complex relationship between magma permeability, outgassing and eruptive style and intensity.
Basaltic magmas, which commonly ascend slowly within the conduit, usually have sufficient time to allow decoupling of the gas and melt phases, producing lava flows and mild lava fountaining activities 10,[29][30][31] .However, the fast ascent rates estimated for basaltic Plinian and sub-Plinian eruptions (with average rates ranging between 4 and 75 m s −1 ) 32,33 likely restricts the time available for gas-melt decoupling during magma ascent, promoting, depending on the magma viscosity, fragmentation and an explosive eruption or high-intensity fountain activity 3,23,29,[34][35][36][37][38] .This suggests that in order to better understand the influence of permeability and porosity on eruptive style and intensity, it is fundamental to also consider other processes which occur during magma ascent, using a holistic and comprehensive approach.
We report 3D quantification of the pore network for pyroclasts of the basaltic Plinian eruptions of Las Sierras-Masaya volcano (Nicaragua) and pyroclasts of the 122 BC Plinian eruption and recent lava fountain activity at Etna volcano (Italy), acquired using phase-contrast synchrotron X-ray μCT 28,39 .Our sample suite includes vesicular scoria lapilli from three basaltic Plinian eruptions: the Fontana Lapilli (FL), the Masaya Triple Layer (MTL) and the Etna 122 BC eruptions (Supplementary Table 1) [40][41][42] .To compare different styles of activity, one pyroclast from the 23 November 2013 lava fountain of the New South-East Crater of Etna 43,44 and one pyroclast from the 24-31 July fountaining phase (lower vent) of the 2001 flank eruption of Etna 45,46 were also examined.Andronico et al. 47 classified the eruptive style of activity at Etna based on volcanological observations, where the 2001 eruption is a sustained lava fountain (LF) and the 23 November 2013 event a large-scale lava fountain (LSLF).LSLFs are a more intense eruptive style compared to LFs, where the erupted tephra volume is one order of magnitude larger and ash dispersal reaches up to 100s of km from the vent.Instead, the smaller eruption plume produced by LFs leads to ash and lapilli fallout in a more restricted area, reaching up to 10s of km from the vent.By examining the pyroclasts of these two eruptions, we can also compare our results from fountain activity of different intensities at Etna.By using X-ray μCT, we can enhance the contrast between microstructural features of interest within the sample volume and directly quantify, with high accuracy, the key parameters required for describing magma permeability, as outlined by the Kozeny-Carman relationship, such as the connected porosity, tortuosity, and throat-pore size ratio.We incorporate our data into a numerical conduit model 48 , to simulate outgassing during magma ascent under several different ascent rates.With this model, we can examine the complex, interdependent conduit processes that govern the explosive-effusive transition at basaltic systems and the impact of magma permeability on conduit flow dynamics and the style and intensity of basaltic eruptions.

Results and discussion
Vesicle morphology, distributions and number densities.FL pyroclasts show a high number of spherical vesicles with regular shaped margins, whilst Etna 122 BC pyroclasts show more irregular, polylobate vesicles (Fig. 1).MTL pyroclasts show both textures.For the Etna pyroclasts, vesicle size decreases and their shapes become more irregular as eruption intensity increases, from the Etna 2001 and 2013 fountains to the Etna 122 BC Plinian eruption.
Vesicle size distributions (VSDs) of Plinian clasts (Fig. 2) can be described by a power-law trend with an exponent >1 when bubble volumes exceed 10 6 μm 3 .Overall, the exponent of the power-law fit of the VSD generally increases with explosivity, with the Etna 2001 fountain clast showing a lower exponent compared to the Plinian clasts (Fig. 2).A few clasts (FL_LG_2; MTL2_a0_2) show mixed power-law-exponential distributions, with the exponential trend observed between 10 2 and 10 6 μm 3 size classes (Fig. 2 and Supplementary Data 1).
FL pyroclasts have the highest bubble number densities (N d ), ranging between 2635 and 11,804 mm −3 (Table   Bubble nucleation and growth processes.Variations in vesicle morphology (Fig. 1) may relate to the variable crystallinity of the clasts across the sample suite.Previous studies find that clasts of the FL eruption show considerable textural heterogeneity, where crystal fraction varies between 0.03 and 0.19 in crystal-poor domains and between 0.45 and 0.64 in crystal-rich domains 36,49 .Similarly, clasts of the MTL eruption have been found to show crystal-poor (where crystal fraction varies between 0.2 and 0.25) and crystal-rich (0.35 and 0.5) domains 50 .Crystalline clasts of unit C of the Etna 122 BC eruption have crystal fractions reaching 0.69-0.72 34 , greater than the crystal fractions of lava fountain clasts such as the Etna 2001 eruption (0.15-0.23) 46 and Etna 2013 eruption (0.37).The irregular, polylobate vesicles observed in samples of the Etna 122 BC eruption and MTL clasts may relate to the highly crystallised groundmass (Fig. 1 and Supplementary Fig. 1) 34,36,50 , as bubbles deform around the rigid, crystalline network.Instead, the spherical, regular morphology of vesicles in FL clasts may indicate that these clasts represent the crystal-poor endmember.Compared to the Etna 122 BC clasts, the more crystal-poor Etna fountain clasts show spherical-subspherical vesicles with more regular morphology, also as a result of their low crystal content.
The observed VSD likely reflects bubble nucleation and growth processes.For crystal-poor clasts such as those of the FL eruption,

(l) Etna_2001
N(>V) = 1 x 10 7 V -0.8 r 2 = 0.9934   1.The solid black lines show where vesicle volume distributions can be described by a power-law trend.The dashed black curves show where vesicle size distributions can be described by an exponential trend.Most vesicle size distributions can be fit with a powerlaw trend, although samples FL_LG_2 and MTL2_2 have mixed power-law-exponential distributions.
the observed VSD was not likely modulated by crystallisation and bubble coalescence.However, we recognise that interbubble film drainage times may be short for lower viscosity, crystal-poor basaltic melts and that consequently, some bubbles may have been able to coalesce 51 and form larger bubbles.Instead, for highly crystalline pyroclasts such as those of the Etna 122 BC eruption, coalescence, likely facilitated by crystallisation, is dominant, producing irregular, polylobate vesicles with high connectivity (Fig. 3).Both the FL and Etna 122 BC pyroclasts produce VSDs which can be described by power-law trends at large vesicle volumes, despite their considerable differences in vesicle morphology.
A VSD power-law fit indicates continuous, multiple bubble nucleation events due to the maintenance of volatile supersaturation, coalescence, and non-equilibrium degassing behaviour 18,20,[52][53][54][55] .Continuous nucleation can occur in H 2 O-CO 2 -bearing basaltic melts under disequilibrium conditions where CO 2 supersaturation is maintained 55 .The lower exponent of the power-law fit for the 2001 Etna clast, and from other studies of Strombolian activity 19,20 suggests a lower number of nucleation events compared to the Plinian case.
The transition towards equilibrium degassing, where vesicles are of similar size and separated by comparable distance, requires timescales of minutes to hours for basaltic melts, and is reflected by a shift from power-law to exponential VSDs 20,56 .As most FL, MTL and Etna 122 BC pyroclasts show power-law VSDs, rapid, continuous bubble nucleation and growth under disequilibrium conditions likely occurred.The mixed power-law-exponential VSD observed in some FL and MTL clasts suggests that also some bubble growth and coalescence occurred alongside continuous  bubble nucleation 56 .Mixed power-law-exponential VSDs have been observed in erupted products of the 2019 Stromboli paroxysms, suggesting minor bubble growth and coalescence, continuous, multiple bubble nucleation events and disequilibrium degassing promoted by rapid magma ascent, resulting in the enhanced explosivity of these events 57 .
The high N d of Plinian clasts may arise from the high decompression rates estimated for basaltic Plinian eruptions of 0.1-2 MPa s − 1 3,32,33,38,49 .Our calculated N d are consistent with decompression experiments on H 2 O-CO 2 bearing basaltic melts at decompression rates of 0.078 MPa s −1 , which produce N d between 10 3 and 10 4 mm −3 55, 56 .Instead, experiments performed with lower decompression rates between 0.017 and 0.039 MPa s −1 produce N d lower than 10 3 mm − 3 56,58 .In general, the decompression rate affects the volatile supersaturation of the melt and, in turn, controls N d when bubbles nucleate via a homogeneous mechanism 8 , by nucleating directly from the melt.In this case, volatiles exsolve without facilitation from pre-existing crystals and other phases 59 .High decompression rates increase volatile supersaturation, promoting nucleation and restricting the time available for bubble growth and coalescence 8,56 .Our Plinian clasts have VSDs characterised by power-law trends, suggesting rapid and continuous bubble nucleation.Spherical vesicles in crystal-poor FL pyroclasts suggest bubble nucleation was the dominant process; however, vesicle morphology in crystal-rich MTL and Etna 122 BC pyroclasts was likely modified by coalescence, facilitated by the more crystalline matrix, following bubble nucleation.
It has also been shown that N d increases with eruption intensity 35,60 .A higher N d implies fewer connections between bubbles, favouring an increase in overpressure and the likelihood of magma fragmentation 9 .The N d of Plinian clasts share more similarities with large-scale lava fountain and paroxysm 61,62 clasts than clasts from Strombolian activity (Fig. 4), suggesting that high nucleation rates may favour a Plinian (or sub-Plinian) eruption or a large-scale lava fountain.Calculated N d for clasts of Plinian and large-scale lava fountain activity are also comparable to the N d estimated for clasts from the 2019 Stromboli paroxysms 57 .The higher N d of large-scale lava fountain and Strombolian paroxysm pyroclasts, when compared with clasts from more typical lava fountain and Strombolian activity (Fig. 4a), further suggests the connection between eruption intensity and N d .
Connectivity, tortuosity and the throat-pore size ratio.Magma permeability depends on the 3D geometry of the pore network, including connectivity, the throat-pore size ratio and tortuosity of connections 14 .FL clasts show the largest variation in porosity (0.21 and 0.79), compared to the MTL clasts (0.43-0.69) and the Etna 122 BC clasts (0.62-0.68).The Etna 2013 and Etna 2001 lava fountain clasts tend to show lower porosities (0.58-0.59) (Fig. 4 and Table 1).
Connectivity (C) is calculated following Colombier et al. 14 : where Φ is the porosity and Φ c the connected porosity.MTL and Etna 122 BC clasts show high connectivity (0.98-1), whereas FL clasts have the widest range in connectivity (0.05 and 1).The throat-pore size ratio (f tb ) represents the ratio between the radius of the throat and vesicle (average of the radii of the two connected vesicles) (Supplementary Methods).Being a nondimensional number, this value is independent from the size of the pores and provides general information on the pore network within the sample.In terms of outgassing and permeability, the throat size is the relevant parameter to consider, which can be calculated as the mean throat-pore size ratio multiplied by the average bubble radius (r b ).There is a narrow range (0.16-0.27) in f tb across clasts of Plinian and fountain eruptions (Table 1).The mean f tb for the Etna 122 BC clasts overlaps with clasts of the 2001 and 2013 fountain eruptions.The mean f tb calculated for the basaltic Plinian FL, MTL and Etna 122 BC clasts are comparable to frothy pumices of the silicic Kos Plateau Tuff eruption, which range between 0.1 and 0.2 16 .The mean f tb calculated for crystalrich clasts is higher than crystal-poor clasts at similar values of Φ, although these values overlap when considering the standard deviation.
The average tortuosity (Supplementary Data 2) and tortuosity factor (m) show greater variation across the clasts compared to f tb and N d .Tortuosity is the ratio between the actual length of the flow path and the length of a path following a straight line in the direction of flow (represented by a value of 1) (Supplementary Methods) and therefore increases as the flow path becomes more convoluted 31 .Calculated tortuosity varies across the sample suite, although crystal-poor clasts generally show values of tortuosity ~1.5.The tortuosity was used to calculate the tortuosity factor (m), following Degruyter et al. 16 : The above expression shows the relationship between τ and Φ c through Archie's law, where τ behaves as a fractal with respect to Φ c 63,64 .The mean m calculated for FL clasts varies between 1.16 and 4.09, whilst MTL clasts show a greater range in mean m, between 1.79 and 6.The Etna 122 BC clasts show a higher mean m (5.79-6.46)than the 2013 and 2001 clasts (2.14-2.28).
Controls on magma permeability.Pyroclasts with high porosity generally show high connectivity (Figs. 3 and 4b).However, vesicles are more connected, complex and irregular in shape in more crystalline clasts (Fig. 1).High crystallinities may promote coalescence of isolated, small-medium sized bubbles, favouring the formation of connected pathways [65][66][67] .Indeed, our crystalrich MTL and Etna 122 BC clasts show high connectivity.Generally, lava fountain clasts show a wide range in connectivity for the same value of porosity (Fig. 4b), which may reflect the variety of clast types erupted during fountain eruptions with variable quench and cooling rates, possibly enabling post-fragmentation vesiculation and an increase in connectivity 10,14 .Instead, shrinkage and densification of large bubbles in Strombolian clasts may decrease porosity 10,68 .
Although there is only a small difference in connectivity between our clasts from Plinian and less explosive activity, basaltic Plinian clasts generally have lower porosity and higher connectivity than fountain and Strombolian activity (Fig. 4).In general, large-scale fountain clasts show higher porosities than Plinian clasts (Fig. 4b), suggesting a greater influence of bubble expansion and/or coalescence.The spherical-subspherical vesicles observed in the 2001 and 2013 Etna lava fountain clasts suggest bubble expansion, rather than coalescence, promoted the formation of larger bubbles.
Smaller values of f tb decrease magma permeability, increasing the likelihood of an explosive eruption 2,9 .The low values of f tb calculated for the Plinian and large-scale lava fountain clasts suggest that magma permeability was low and gas-melt coupling was maintained during these high-intensity eruptions.However, the low f tb also calculated for the low-intensity Etna 2001 lava fountain clast suggests that f tb may have less of a control on explosivity than other parameters such as N d .
In general, crystal-rich pyroclasts tend to show a higher mean tortuosity factor (m) (Table 1 and Supplementary Data 2).Although crystallisation can promote connectivity, crystallisation can also produce more tortuous outgassing pathways as bubbles deform to accommodate the rigid crystalline network 31 , reducing magma permeability as the greater path length and curvature increases flow resistance from viscous and inertial effects 6 .Indeed, the crystal-rich Etna 122 BC and MTL clasts show the highest values of m.Rapid syn-eruptive crystallisation has been suggested for the MTL and Etna 122 BC eruptions 33,50 , where the crystallisation timescale may be sufficiently fast to restrict coalescence and the development of permeable outgassing pathways during magma ascent.Basaltic Plinian clasts show similar values of m to those calculated for silicic frothy pumices (3-10.7) 16.Overall, our results suggest that m increases with explosivity and crystal content.
Our dataset confirms not only the overlap in porosity and pore connectivity observed for basaltic Plinian and less explosive basaltic pyroclasts 10 , but also for f tb .This result is consistent with works which find similar values of magma permeability when considering different eruptive styles at basaltic systems 29 .
Sensitivity analyses.Our results suggest that the characteristics of the pore network (i.e., N d , f tb and m) extrapolated from basaltic pyroclasts are not sufficient to distinguish the different styles of activity at basaltic volcanoes.To investigate the role of permeable outgassing at basaltic volcanoes using a more holistic and comprehensive approach, we performed numerical simulations of magma ascent dynamics for different test case scenarios: Etna 122 BC, MTL crystal-poor and MTL crystal-rich.We performed different sensitivity analyses using the 1D steady-state magma ascent model for a cylindrical conduit of La Spina et al. 48, varying the 4 characteristic parameters reported in Eqs.(5-7): N d , m, f tb and f 0 .We use a range of 10 9 −10 14 m −3 for N d , 1.01-10 for m, 0.03-1 for f tb and 10 −5 to 1 for f 0 .These ranges are constrained by the data obtained from this study on natural pyroclasts (Supplementary Table 2 and Table 1) and Polacci et al. 20 for basaltic eruptions of lower explosivity.For the Etna and MTL (crystal-rich and crystal-poor) eruptions we performed sensitivity analyses using a conduit radius (r) of 35 and 25 m, respectively, to produce mass eruption rates consistent with estimates from the corresponding eruption deposits 33,40 .Furthermore, for the MTL crystal-rich case only, we performed 3 different sensitivity analyses with three different conduit radii: 5, 10 and 25 m.Since the mass eruption rate (and thus magma ascent rate) is highly affected by the conduit radius, investigating different radii will allow us to understand the role of the permeable network in outgassing for different magma ascent rates.
The description of the numerical model and the initial conditions for each test case study are provided in the Supplementary Methods section of the Supplementary Information.In our numerical simulations, we distinguish between two types of solutions: explosive and non-explosive solutions.An explosive solution is obtained when fragmentation is achieved during magma ascent within the conduit, and this can be seen as an example of sub-Plinian or Plinian eruptions, depending on the mass eruption rate.Instead, a non-explosive solution is obtained when fragmentation is not achieved within the conduit, and these solutions may be interpreted as representing different styles of   62 .Compared with the classification of Andronico et al. 47 for fountain eruptions at Etna, the erupted volume and eruption rate of high and intermediate fountains at Kīlauea 83 are more representative of large-scale lava fountains (LSLF) and low fountains more representative of lava fountains (LF).Literature data were obtained from the following references: Polacci et al. 19,20 , Stovall et al. 61 and Parcheta et al. 62 .Data from Stovall et al. 61 and Parcheta et al. 62 represent textural measurements performed on 2D images.b Porosity vs connectivity.Grey circles represent literature data on fountain eruptions but are not coloured according to fountain height.Literature data were obtained from the following references: Colombier et al. 14 , Kawabata et al. 84 and Polacci et al. 66 .
activity: from lava flows (where gas is able to decouple efficiently from the magma) to lava fountaining of varying intensity (where there is a greater degree of coupling between gas and magma) 3 .For the MTL crystal-poor case (Supplementary Fig. 2), <1% of the solutions are explosive.Non-explosive solutions are found for all input values, which have a mass eruption rate of 10 8 kg s −1 and exit velocity on the order of 100 m s −1 .These solutions may represent large-scale lava fountains, which typically have high mass eruption rates and ascent velocities 3 .
For the crystal-rich MTL sensitivity analysis (r = 25 m), 77% of the solutions are explosive, with a higher percentage at high N d , high m, low f tb and high f 0 (Fig. 5).Similarly, clasts of the FL and MTL basaltic Plinian eruptions show high N d ≥ 10 12 m −3 , high m ≤ 6 and low f tb ≤ 0.27 (Fig. 5a-c and Table 1).Similar trends are observed for the MTL sensitivity analyses performed with lower conduit radii (Supplementary Fig. 3); however, there is a lower percentage of explosive solutions for a conduit radius of 10 m (73%) and 5 m (63%).Explosive solutions (r = 25 m), where magma fragments within the conduit, have mass eruption rates of 10 6 −10 8 kg s −1 and exit velocities on the order of 100 m s −1 (Fig. 6), consistent with an explosive eruption of Plinian intensity 40,41 .Non-explosive solutions (r = 25 m) are found at low N d , low m, high f tb and low f 0 (Fig. 5), and have mass eruption rates of ~10 6 kg s −1 , but with exit velocities on the order of 1 m s −1 or less (Fig. 6).These are representative of effusive eruptions.
At lower r, the mass eruption rate of explosive solutions decreases to 10 5 −10 7 kg s −1 (r = 10 m) and 10 4 −10 6 kg s −1 (r = 5 m), whilst there is an increase in the number of explosive solutions with an exit pressure of 10 5 Pa.In this case, the gas-ash cloud ejected from the vent of the conduit will not experience further acceleration, resulting in a lower ash column compared to the overpressurized gas-ash cloud.We also observe a larger range of exit velocities for explosive solutions (ranging from ~50 to ~130 m s −1 ) compared to that obtained for a conduit radius of 25 m (Fig. 6).The mass eruption rate of non-explosive solutions decreases from 10 6 to 10 3 kg s −1 as r decreases from 25 to 5 m.
The Etna 122 BC sensitivity analysis shows a lower percentage of explosive solutions (59%).However, a comparable trend is observed, with explosive solutions found at high N d , m and f 0 and low f tb (Fig. 5e-h).Clasts of the Etna 122 BC Plinian eruption are comparable and have high N d ≥ 10 12 m −3 , high m ~6 and low f tb ≤ 0.24 (Fig. 5e-g and Table 1).N d and f 0 show the strongest correlation with explosivity (Fig. 5e, h); however, ~20% of the solutions at N d > 10 12 m −3 are non-explosive.
Numerical solutions show comparable mass eruption rates for explosive and non-explosive solutions (between ~10 7 and ~10 8 kg s −1 ), even though non-explosive solutions tend to The results are expressed as a percentage as the sum of the frequencies is equal to 100.Results from the MTL crystal-rich sensitivity analysis using a conduit radius of 25 m are highlighted in blue.Results from the Etna 122 BC sensitivity analysis are highlighted in green.Explosive solutions are represented by filled squares, whilst non-explosive solutions are represented by white diamonds with a coloured outline.The grey circles represent the total number of solutions, accounting for both the explosive and non-explosive solutions.The sum of the frequencies is 100; therefore, where there is a decrease in the total number of solutions, there is an increase in the number of simulations where the conduit model could not find a solution which satisfied the boundary conditions.The vertical grey bars represent the range in N d , m, f tb and f 0 obtained from X-ray μCT data from this study (Table 1) for clasts of the FL, MTL and Etna eruptions.
have lower values than explosive solutions.Exit velocities for explosive solutions range between ~70 and ~140 m s −1 .Instead, non-explosive solutions form two distinct groups characterised by different exit velocities, where one group of solutions has exit velocities <~40 m s −1 and the other group has higher exit velocities ranging between ~90 and ~130 m s −1 (Fig. 6).These two groups may represent lower intensity and large-scale lava fountain activity respectively 3 .
The effect of the Forchheimer parameters on conduit dynamics.MTL, Etna 122 BC and FL Plinian clasts show high N d ≥ 10 12 m −3 , f tb < 0.27 and m ≤ 6.46 (Table 1).The sensitivity analysis confirms that these ranges may promote an explosive eruption (Fig. 5).Indeed, high N d , high m and low f tb are factors which can reduce magma permeability and outgassing efficiency 2,6,9 .Our numerical results indicate that, within the ranges investigated, N d and f 0 have the greatest impact on magma permeability and outgassing, and thus, on the style of activity.For the MTL crystal-rich sensitivity analyses, non-explosive solutions are not found at N d > 10 13 m −3 , and non-explosive solutions are rarely found at f 0 > 10 −1 , occurring only at the lowest N d (Supplementary Fig. 4).Instead, for the Etna 122 BC case, two groups of non-explosive solutions can be observed (Supplementary Fig. 4).At N d < ~10 12 m −3 and f 0 < 10 −1 , non-explosive solutions are found with variable f tb and m.A second group of non-explosive solutions are found at high N d and f 0 < 10 −1 , instead with a low f tb , indicating a higher degree of gas-melt coupling.Similarly, clasts of the Etna 2001 lava fountain and 2013 large-scale lava fountain show low values of m and f tb , but a large difference in N d , where N d of the 2013 fountain reach 10 12 m −3 (Table 1).These two groups of solutions are likely those characterised by different exit velocities and exit pressure, representing large-scale and lower-intensity fountain activity, respectively.A high N d , low f 0 and low f tb may also produce a large-scale lava fountain (Supplementary Fig. 4), where magma ascends at high rates within the conduit, maintaining gas-melt coupling 3 .The Forchheimer parameters also influence output parameters such as exit velocity and mass flow rate.For the Etna 122 BC sensitivity analysis (Fig. 7), f 0 is the most important control on exit velocity and exit gas volume fraction of explosive solutions, whilst f tb influences the mass flow rate and fragmentation depth.For non-explosive solutions (Supplementary Fig. 5), N d has the greatest control on the exit velocity and exit gas volume fraction of non-explosive solutions, further supporting the influence of N d on the exit velocity of lava fountain activity.Instead, m has a more important control on mass flow rate, whilst f 0 has the least influence on the outcome of the style of non-explosive solutions.However, for explosive solutions, the parameters with the greatest control on the flow regime differ for explosive and non-explosive solutions.
For the MTL sensitivity analysis (r = 25 m) (Fig. 7), N d , f 0 and f tb are the most important factors which influence the exit gas volume fraction, mass flow rate, exit velocity and fragmentation depth of explosive solutions.For non-explosive solutions, m has the greatest control, succeeded by N d .A similar trend is observed for the MTL sensitivity analysis with a conduit radius of 5 m (Supplementary Fig. 6).
The role of magma ascent rates in outgassing.If low permeability and gas-melt coupling are maintained during fast ascent, such as in Plinian and sub-Plinian eruptions, magma failure can occur with increasing bubble growth and expansion 1,54 .Instead, a combination of low ascent velocities and permeability development allows gas bubbles to decouple from the melt, move, and coalesce, increasing the likelihood of a Strombolian or effusive eruption 19,69 .Estimated syn-eruptive magma ascent rates show a clear distinction between explosive and effusive eruptions 1 , indicating that the outgassing timescale must be considered alongside the development of permeability during ascent.
The overlap in porosity, connectivity and permeability observed for clasts of basaltic Plinian and Strombolian activity 10,29,31 could relate to the high decompression and ascent rates which occur during a Plinian eruption, which may restrict outgassing even for magmas with high permeability.Numerical simulations of the 2010 Merapi eruption 27 also find that varying the percolation threshold and magma permeability has less influence on eruptive style than decreasing the conduit radius, which reduces the ascent velocity and facilitates the transition from an explosive to effusive eruption.
Comparing the numerical results obtained from the MTL crystal-rich sensitivity analyses, those performed at lower r show greater sensitivity to the parameters describing outgassing, with an expansion of the parameter space where non-explosive solutions would occur as r decreases (Supplementary Fig. 4).However, we also observe an increase in the number of simulations with decreasing r, for which the boundary conditions at the vent of the conduit are not satisfied (either the atmospheric pressure or choked flow condition) for the inlet boundary conditions fixed.For the MTL crystal-rich sensitivity analysis, assuming a 25 m conduit radius, this occurs for ~3% of the simulations, increasing to ~6% and ~18% of the simulations as r is decreased to 10 and 5 m, respectively.Analysing these simulations, we found that, for the inlet conditions assumed and for some combinations of the Forchheimer parameters, strong decoupling between the gas and magma allowed gas to escape completely from the system, leaving the degassed magma behind.In these cases, magma was not able to reach the surface, and thus the numerical model cannot find a solution which satisfies the right boundary conditions at the vent of the conduit.These simulations occur mostly at low N d , low m and high f tb , conditions which enhance decoupling (Supplementary Fig. 3).The increase in the number of these simulations with decreasing r indicates that magma ascent rate plays an important role in outgassing.If magma ascent rate is low, efficient outgassing may occur, resulting in an effusive eruption or even no eruption at all.Instead, for the same combinations of the Forchheimer parameters (N d , m, f tb and f 0 ), but at high magma ascent rates, a highly explosive eruption (Plinian or sub-Plinian) or large-scale fountain activity may occur, even for a magma with high permeability.
The results of the sensitivity analyses demonstrate that magma ascent rate has a strong control on the style of activity (Fig. 8).Generally, explosive solutions are found at higher Reynolds number (Re), or at low values of u rel /u mix , where u rel is the relative velocity between the gas and melt phases and u mix the mixture velocity during magma ascent.The Reynolds number provides an indication of the magma ascent rate (Supplementary Information), whereas u rel /u mix provides an estimation of the outgassing efficiency (where lower values of u rel /u mix correspond to a lower outgassing efficiency).We remark that the outgassing efficiency is strongly correlated with the Reynolds number, since both are a function of the ascent rate.For the Etna 122 BC case (Fig. 8e-h) the numerical results can be divided into three main groups, depending on the viscous permeability (k v ).For k v < 10 −11 m 2 , explosive simulations are always obtained.These simulations are characterised by high ascent rate (high Reynolds number) and low outgassing efficiency (low u rel /u mix ).For 10 −11 <k v < 10 −10 m 2 most of the simulations are explosive, with a cluster of non-explosive simulations obtained at higher outgassing efficiency (high u rel / u mix ).Given the high ascent rates of the non-explosive simulations in this group (Reynolds number ~10 −1 -10 0 ), these simulations are most likely representative of large-scale lava fountaining activity 3 .For k v > 10 −10 m 2 , instead, both explosive and non-explosive simulations occur, which can be grouped according to the Reynolds number (or the outgassing efficiency).Specifically, explosive simulations can be obtained for Reynolds numbers >~10 -2 , whereas for Reynolds numbers <~10 -1 non-explosive simulations occur.
A similar behaviour can be observed for the MTL sensitivity analyses.Specifically, using a conduit radius of 25 m (Fig. 8a-d), large-scale lava fountaining activity is never obtained, as we observe non-explosive solutions only at low ascent rates (low Reynolds number) and high outgassing efficiency (high u rel / u mix ).For k v < 10 −10 m 2 explosive solutions always occur.However, for k v > 10 -10 m 2 explosive or non-explosive solutions are found, where the outcome of the simulation depends on the ascent rate but not on the viscous or inertial permeability.Thus our sensitivity analyses show (Fig. 8 and Supplementary Fig. 7) that, at high magma permeability (k v > 10 -10 m 2 ), the style of activity appears to be not directly affected by the viscous and inertial permeability, but rather by the magma ascent rate, as explosive solutions are still observed at high values of magma permeability.
Several studies have shown that strong gas-melt coupling and restricted outgassing may promote magma acceleration and largescale fountain activity 10 or a highly explosive eruption 37 .This is supported by our observations of clasts of basaltic Plinian eruptions and large-scale lava fountain activity, which have similar N d , C and f tb (Fig. 4 and Table 1).Our numerical results confirm that a high degree of gas-melt decoupling is unable to occur even at high magma permeabilities under fast magma ascent rates.Fast magma ascent rates, therefore, may inhibit outgassing independent of magma permeability, producing eventually highly explosive eruptions.
The effect of initial conditions on eruptive style.As <1% of the MTL crystal-poor simulations are explosive, and 23% of the MTL crystal-rich (r = 25 m) simulations are non-explosive, initial storage conditions may also be important in promoting Plinian activity at basaltic systems.Initial temperature has been proposed as an important control on the eruptive style at basaltic volcanoes 3,33,[48][49][50] , driving rapid syn-eruptive crystallisation under conditions of high undercooling.The difference between the initial conditions for the crystal-poor and crystal-rich MTL sensitivity analysis are only a 20 °C decrease in temperature, and a 7 vol.%increase in crystal content 50 .As the input ranges for N d , m, f tb , f 0 are the same for both sensitivity analyses, the difference in the outcome of the simulations for the crystal-poor and crystal-rich case may also be attributed to the initial temperature and crystal content in the magma reservoir.Magma storage conditions can influence the explosivity of silicic eruptions, with explosive eruptions promoted within a specific range of initial H 2 O content and crystallinity 70 .Common pre-eruptive conditions have also been found for basaltic Plinian eruptions, such as moderate storage temperatures and low volatile contents 33,41,48 .The results of our sensitivity analysis highlight that eruptive style is more sensitive to changes in magma storage conditions, such as pre-eruptive temperature and crystallinity, rather than the properties which control magma permeability.Figure 9 summarises the results of the sensitivity analysis, supported by our observations and measurements of natural pyroclasts.At high magma ascent rates, with decreasing temperature and increasing crystal content, outgassing is inhibited, and magma viscosity increases 33,49,50 , leading to brittle fragmentation.A highly explosive eruption may transition to lava fountain activity with an increase in temperature or decrease in crystal content under high magma ascent rates.For basaltic systems such as Etna which show a wide range in eruptive style, small changes in outgassing efficiency can lead to an increase in explosivity 3,48 .
The driving mechanisms which may promote basaltic Plinian activity may, therefore, relate to rheological changes induced by rapid syn-eruptive crystallisation, which, when combined with rapid magma ascent, maintains gas-melt coupling whilst bubbles nucleate, promoting magma fragmentation 10,23,29,31,33,34,49,50,71,72 .Even if a permeable pore network is able to develop, our results show that fast magma ascent inhibits outgassing independent of magma permeability, preventing gas escape and accumulating overpressure.

Conclusions
By using X-ray microCT, we quantify the 3D pore network of clasts from basaltic Plinian eruptions, finding that clasts are characterised by bubble number densities exceeding 10 12 m -3 , low throat-pore size ratios ≤ 0.27 and tortuosities reaching 6.5.We demonstrate that parameters such as the bubble number density and friction coefficient have considerable influence on outgassing and the eruptive style at basaltic systems.However, by combining our data on vesicle micro-textures with a numerical conduit model, we demonstrate that bubble nucleation and growth processes alone cannot explain the occurrence of basaltic Plinian eruptions.Instead, we find that fast magma ascent rates inhibit outgassing independent of magma permeability during highly explosive basaltic eruptions, as gas-melt coupling is maintained.Small changes in magma storage conditions, such as pre-eruptive temperature and crystallinity, can have considerable influence on eruptive style, driving rapid crystallisation and bubble nucleation during magma ascent, or enabling the establishment of connected, permeable pore networks that facilitate outgassing.Close monitoring of pre-eruptive conditions may, therefore, provide crucial information for hazard assessment and risk mitigation, particularly for basaltic systems that exhibit a wide range of eruptive styles.

Methods
Phase-contrast synchrotron X-ray μCT.Scoria lapilli from the FL and MTL eruptions 41,50 are from the localities described by Bamber et al. 49,50 respectively.Scoria lapilli from the Etna 122 BC eruption were taken from a proximal site (37°40'19.4°N,15°02'09.36°E)and belong to unit C of Coltelli et al. 42 .Further details on the samples used in this study are provided in the Supplementary Information.
Fig. 8 Results of the sensitivity analysis showing how the Reynolds number and degree of outgassing vary with the viscous and inertial permeability.Panels (a-d) present the results of the sensitivity analysis using the MTL initial conditions and a conduit radius of 25 m, showing how the Reynolds number (Re) varies with a the viscous permeability (k v ) and b the inertial permeability (k i ).Also shown is how u rel /u mix varies with c k v and d k i , where u rel represents the relative velocity between the gas and melt phases and u mix the velocity of the mixture.Results of the sensitivity analysis using the Etna 122 BC initial conditions are presented in (e-h).Each symbol represents the outcome of a single simulation, where explosive solutions are indicated by the coloured squares and non-explosive solutions by the white diamonds.All values are obtained from the fragmentation level for explosive solutions and at the vent for non-explosive solutions.For the Etna 122 BC sensitivity analysis, two groups of non-explosive solutions are observed, with different ranges of the Reynolds number and u rel /u mix .
The synchrotron X-ray μCT experiments were carried out at the SYRMEP beamline of the Elettra synchrotron facility in Basovizza (Trieste, Italy).Experiments were performed in propagation-based phase-contrast mode 73 using a filtered polychromatic (filters: 1.5 mm Si plus 1.0 mm Al), near-parallel X-ray beam with a mean beam energy of 27 keV.The sample was mounted on an air-bearing rotation stage and rotated continuously over a total rotation angle of 180°.The exposure time for each radiograph (projection) was 2 s and 1800 projections were acquired for each scan.The total scan duration was ~1 h.The sample-to-detector distance was set to 250 mm for the MTL clasts and 160 mm for the FL and Etna clasts.A water-cooled, 16-bit, 2048 × 2048 pixels, sCMOS macroscope camera (Hamamatsu, Orca Flash 4.0) coupled via a high-numerical aperture optics to a 17-μm-thick GGG:Eu scintillator screen, was used to acquire the scans.The effective pixel size was set to 0.9 μm, corresponding to a field of view of approximately 1.8 mm × 1.8 mm.
The SYRMEP Tomo Project (STP v.1.2.5) software 74 was used to reconstruct the 2D axial slices from the projections and remove eventual ring artefacts 75 .A single-distance phase-retrieval algorithm 76 was applied to projections prior to reconstruction to enhance the contrast between vesicles and glass, whilst minimising the effect of phase-contrast artefacts at the edge of vesicles 77 .The δ/ β parameter (the ratio between the real and imaginary parts of the refraction index of the investigated material) was set to 80 for the MTL clasts, and 100 for the FL and Etna clasts (ad hoc adjustment).The axial slices, reconstructed via a filtered back-projection algorithm, were imported into ImageJ 78 , converted into the 8-bit raw format and stacked to produce complete 3D data (virtual volumes).
Vesicles were segmented by applying the 3D Otsu automatic thresholding algorithm 79 with manual refinement following visual inspection of the sample using Avizo, to ensure the selection of thin pore-melt films.Porosity (Φ), connected porosity (Φ c ), bubble number density (N d ), tortuosity (τ) and the throat-pore size ratio (f tb ) were calculated using modules available in the Avizo software (v.2019.1;Thermo Fisher Scientific, USA) (Supplementary Information).Tortuosity was used to calculate the tortuosity factor (m), following Degruyter et al. 16 .The 3D visualisation of raw (after tomographic reconstruction) and processed images was obtained through a volume rendering procedure, using Avizo and the commercial software package VGStudio MAX 2.0 (Volume Graphics, Germany).
Sensitivity analyses using the 1D steady-state conduit model.We used the 1D steady-state numerical conduit model of La Spina et al. 3,9,48 to perform sensitivity analyses investigating the effect of the different outgassing parameters, as outlined by the Forchheimer equation and the Kozeny-Carman relationship, on the outcome of numerical solutions of magma ascent.The model simulates conduit processes such as crystallisation, outgassing, conduit wall friction and magma fragmentation.A complete description of the model and its governing equations is provided in the Supplementary Methods section of the Supplementary Information.The results of the numerical simulations are provided in Supplementary Data 3 and 4.
During ascent, magma can fragment within the conduit, transitioning from a mixture of melt, crystals, and gas (dissolved and exsolved) to a continuous mixture of gas and solid particles.In this study, magma fragmentation is simulated using the strain-rate criterion of Papale 80 , where fragmentation occurs when the Deborah number (De) (the ratio between the Maxwell relaxation timescale of the fluid and the timescale of deformation) exceeds 0.01 3 .
Following La Spina et al. 9 , above the fragmentation level, interphase friction is described using the model of Yoshida and Koyaguchi 12 .Below the fragmentation level, outgassing is described by a formulation of Forchheimer's Law, following Degruyter et al. 2 .The conservation of momentum for the liquid and gas phases is used to write the equation for the mixture and relative velocity 3 .In the conduit model, mechanical decoupling between the gas and the melt is controlled by the relaxation parameter τ (f) , and is defined as follows: where C D is a drag coefficient and r a the average size of the fragmented particles.Above the fragmentation level, ϕ f is equal to 1, where ϕ f is the fragmentation efficiency.Below the fragmentation level, ϕ f is equal to 0. The subscript l and g denote the liquid and gas phase, respectively, whilst u and ρ represent the phase velocity and phase mass density.The viscous (k v ) and inertial (k i ) permeabilities are calculated following Degruyter et al. 2 : where f tb is the throat-pore size ratio, m is the tortuosity factor, f 0 is the friction coefficient, and α g is the exsolved gas volume fraction.The average bubble size r b is calculated as follows: where N d is the bubble number density.The effect of the different parameters on outgassing are evaluated by varying the values of N d , m, f tb and f 0 .The opensource DAKOTA toolkit (Design Analysis Kit for Optimization and Terascale Applications) 81 , developed at Sandia National Laboratories, was used to perform the sensitivity analyses.For each analysis, 10,000 simulations were performed, varying the parameters within defined ranges informed by the data obtained by this study and Polacci et al. 20 (Supplementary Table 2).There are no available constraints on f 0 for basaltic magmas.Degruyter et al. 2 considered a range between 10 and 100, estimated for rhyolite compositions.However, f 0 may be much lower for lowviscosity basaltic magmas 9 .To consider this possibility, we assumed a wide range for the friction coefficient.
Sensitivity analyses were performed for both the Etna 122 BC 82 and MTL 50 melt compositions.As there are greater constraints on the vent location for the MTL eruption 49 , we chose to simulate the MTL eruption for our sensitivity analysis.The vent location for the FL eruption is less constrained and thought to lie outside of the present Masaya caldera 41 .As we are investigating the impact of outgassing on highly explosive basaltic eruptions, we chose to simulate the Plinian 122 BC eruption for the Etna case.As our study focuses on the effects of the Forchheimer parameters on explosivity, the remaining input parameters for the sensitivity analyses are constrained using literature data (Supplementary Table 3).For the MTL sensitivity analyses, both the crystal-poor and crystal-rich endmembers are considered 49 .
range, ranging between 2329 and 5111 mm −3 and 3316 and 3434 mm −3 , respectively, whilst the large-scale Etna 2013 fountain clast has a comparable N d of 4028 mm −3 .Instead, the low fountain Etna 2001 clast has a lower N d of 716 mm −3 compared to pyroclasts of more intense activity.

Fig. 1
Fig.1Representative orthoslices and 3D renderings of pyroclasts.Pyroclasts (a-c) FL_D_1, (d-f) MTL2_a0_1, (g-i) Etna122_C3_1, (j-l) Etna2013_1 and (m-o) Etna2001_1 are shown.Panels (b), (e), (h), (k) and (n) show renderings (voxel size = 0.9 µm), which focus on the matrix glass and crystals, whilst panels (c), (f), (i), (l) and (o) show renderings which focus on the segmented bubbles.All sample volumes are 3 mm 3 .The height of each cylinder is 1.35 mm and the diameter is 1.69 mm.There is a change in vesicle morphology from the spherical, smaller bubbles of crystal-poor clast FL_D_1, to the irregular, polylobate bubbles of crystal-rich clasts MTL2_a0_1 and Etna122_C3_1.The clast Etna2013_1 shows smaller vesicles with more irregular margins than clast FL_D_1, whilst clast Etna2001_1 shows the largest vesicles.Volume renderings were produced using VGStudio MAX 2.0.The renderings presented here are selected as representative examples of the pyroclasts analysed for these eruptions.

Fig. 2
Fig. 2 Cumulative vesicle size distributions.Vesicle size distributions of clasts of the a-d Fontana Lapilli, e-h Masaya Triple Layer and i-l Etna eruptions are presented, where N d represents the bubble number density and is calculated as the number of vesicles per mm 3 of the total melt volume.Sample names correspond to those presented in Table1.The solid black lines show where vesicle volume distributions can be described by a power-law trend.The dashed black curves show where vesicle size distributions can be described by an exponential trend.Most vesicle size distributions can be fit with a powerlaw trend, although samples FL_LG_2 and MTL2_2 have mixed power-law-exponential distributions.

aFig. 3 3D
Fig. 3 3D renderings showing the variation in porosity and connectivity for clasts of the Fontana Lapilli and Etna 122 BC eruptions.Both porosity (Φ) and connectivity (C) increase from the crystal-poor clasts a FL_E_1 and b FL_D_1 to the crystal-rich clast c Etna122_C3_1.Connected porosity is highlighted in purple.Connectivity increases with crystallinity, to a maximum value of 1 for the most crystalline clast.Vesicles change from a spherical to more irregular, polylobate morphology as crystallinity and connectivity increase.All sample volumes are 3 mm 3 .The height of each cylinder is 1.35 mm and the diameter is 1.69 mm.Volume renderings were produced using the Avizo software.

Fig. 5
Fig.5Results of the sensitivity analyses showing the percentage of solutions with the varying Forchheimer parameters.The number of explosive/nonexplosive solutions with varying input parameters a, e bubble number density (N d ), b, f tortuosity (m), c, g throat-pore size ratio (f tb ) and d, h friction coefficient (f 0 ) is shown.The results are expressed as a percentage as the sum of the frequencies is equal to 100.Results from the MTL crystal-rich sensitivity analysis using a conduit radius of 25 m are highlighted in blue.Results from the Etna 122 BC sensitivity analysis are highlighted in green.Explosive solutions are represented by filled squares, whilst non-explosive solutions are represented by white diamonds with a coloured outline.The grey circles represent the total number of solutions, accounting for both the explosive and non-explosive solutions.The sum of the frequencies is 100; therefore, where there is a decrease in the total number of solutions, there is an increase in the number of simulations where the conduit model could not find a solution which satisfied the boundary conditions.The vertical grey bars represent the range in N d , m, f tb and f 0 obtained from X-ray μCT data from this study (Table1) for clasts of the FL, MTL and Etna eruptions.

Fig. 6
Fig.6Results of the sensitivity analyses showing how the output parameters melt exit velocity, mass eruption rate and exit pressure vary with explosivity.The MTL crystal-rich sensitivity analyses are highlighted in blue and are shown for an input conduit radius (r) of a-c 25 m, d-f 10 m and g-i 5 m.Results from the j-l Etna 122 BC sensitivity analysis are also shown in green.Explosive solutions are represented by filled squares, whilst nonexplosive solutions are represented by white diamonds with a coloured outline.The range in melt exit velocity, mass eruption rate and exit pressure for the non-explosive solutions is also highlighted by the vertical grey bar.For the Etna 122 BC sensitivity analysis, two groups of non-explosive solutions can be observed, which are characterised by different ranges of melt exit velocity and exit pressure.

Fig. 7
Fig. 7 Sobol index showing the influence of the input parameters N d , m, f tb and f 0 on the output parameters for explosive solutions.The results of the a-d MTL crystal-rich and e-h Etna 122 BC sensitivity analyses are shown, illustrating the influence of the input parameters bubble number density (N d ), tortuosity factor (m), throat-pore size ratio (f tb ) and the friction coefficient (f 0 ) on exit gas volume fraction, mass flow rate, exit velocity and the depth of fragmentation.

Fig. 9 A
Fig. 9 A diagram illustrating how the initial conditions and resulting conduit processes may influence eruptive style at basaltic systems.Different preeruptive conditions and the resulting variation in crystallisation and outgassing within the conduit may promote a fountain activity, b large-scale fountain activity or a c Plinian eruption at basaltic systems.As the conduit radius increases, large-scale fountain activity or a Plinian eruption may occur.Increasing initial crystal content and decreasing initial temperature may promote a transition from a large-scale lava fountain to an explosive Plinian eruption.The degree of coupling between the gas and melt phases increases from the fountain to the large-scale fountain and Plinian examples.For the large-scale fountain and Plinian eruption examples, high magma ascent rates maintain gas-melt coupling during ascent and restrict outgassing.Insets show representative vesicle textures from pyroclasts of this study of the three types of activity.Pyroclast textures show increasing bubble number density (N d ) as eruption intensity increases to the large-scale fountain and Plinian cases.However, highly crystalline clasts of Plinian eruptions show a vesicle network with a higher tortuosity factor (m) compared to fountain activity.

Table 1
Results of the 3D quantitative textural analysis on digital volumes reconstructed from synchrotron-based X-ray microtomographic images.
Porosity, bubble number density and connectivity, for basaltic clasts representing a range of explosivity.a Porosity vs N d , calculated for clasts of the FL, MTL and Etna eruptions.An increase in porosity would likely result from the expansion and/or coalescence of bubbles, whilst N d would increase with increasing nucleation, as indicated by the interpretive arrows on the plot.Coloured symbols represent data from this study, with filled symbols representing crystal-rich clasts, whilst half-filled symbols correspond to crystal-poor clasts.Squares represent data from clasts of Plinian eruptions, whilst circles represent data from clasts of fountain activity.Open triangles represent literature data on Strombolian eruptions at Stromboli, whilst open circles (black, dark blue and light blue) represent literature data from clasts of the 1959 Kīlauea Iki and 1969 Mauna Ulu fountaining eruptions of Kīlauea, coloured according to fountain height.High fountains refer to fountains reaching >400 m in height, intermediate fountains between 100 and 400 m in height and low fountains as <100 m in height