
==== Front
Sci Adv
Sci Adv
sciadv
advances
Science Advances
2375-2548
American Association for the Advancement of Science

39150999
ado2585
10.1126/sciadv.ado2585
Research Article
Earth, Environmental, Ecological, and Space Sciences
SciAdv r-articles
Geology
Physical Sciences
Geology
Direct observation of degassing during decompression of basaltic magma
Direct observation of magma degassing
https://orcid.org/0000-0002-0242-4320
Bonechi Barbara Conceptualization Data curation Formal analysis Investigation Validation Visualization Writing - original draft Writing - review & editing 1 *
https://orcid.org/0000-0003-3318-8700
Polacci Margherita Conceptualization Data curation Formal analysis Funding acquisition Investigation Methodology Project administration Resources Supervision Validation Visualization Writing - original draft Writing - review & editing 1
https://orcid.org/0000-0002-5293-7300
Arzilli Fabio Conceptualization Investigation Methodology Resources Validation 2
https://orcid.org/0000-0001-6736-7884
La Spina Giuseppe Conceptualization Formal analysis Investigation Methodology Writing - original draft Writing - review & editing 3
https://orcid.org/0000-0002-3717-2151
Hazemann Jean-Louis Investigation Methodology Resources Writing - review & editing 4
https://orcid.org/0000-0003-4931-9912
Brooker Richard A. Conceptualization Funding acquisition Investigation Methodology Project administration Resources Supervision Validation Writing - review & editing 5
https://orcid.org/0000-0001-6708-1085
Atwood Robert Data curation Formal analysis Investigation Methodology Resources Software 6
https://orcid.org/0000-0003-4537-210X
Marussi Sebastian Investigation Methodology 7
https://orcid.org/0000-0002-3898-8881
Lee Peter D. Conceptualization Funding acquisition Methodology Writing - review & editing 7 8
https://orcid.org/0000-0002-5781-2152
Wogelius Roy A. Formal analysis Investigation Resources Writing - review & editing 1
https://orcid.org/0009-0008-7224-537X
Fellowes Jonathan Formal analysis Investigation Methodology Resources Validation 1
https://orcid.org/0000-0003-3779-4812
Burton Mike R. Conceptualization Funding acquisition Investigation Methodology Project administration Supervision Visualization Writing - review & editing 1
1 Department of Earth and Environmental Sciences, University of Manchester, Manchester, UK.
2 School of Science and Technology, Geology Division, University of Camerino, Camerino, Italy.
3 Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo, Catania, Italy.
4 Université Grenoble Alpes, CNRS, Grenoble INP, Institut Néel, Grenoble, France.
5 School of Earth Sciences, University of Bristol, Bristol, UK.
6 Diamond Light Source, Harwell Science and Innovation Campus, Harwell, Oxfordshire, UK.
7 Department of Mechanical Engineering, University College London, London, UK.
8 Research Complex at Harwell, Rutherford Appleton Laboratory, Harwell, Oxfordshire, UK.
* Corresponding author. Email: barbara.bonechi@manchester.ac.uk
16 8 2024
16 8 2024
10 33 eado258524 1 2024
11 7 2024
Copyright © 2024 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution License 4.0 (CC BY).
2024
The Authors
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution license, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Transitions in eruptive style during volcanic eruptions strongly depend on how easily gas and magma decouple during ascent. Stronger gas-melt coupling favors highly explosive eruptions, whereas weaker coupling promotes lava fountaining and lava flows. The mechanisms producing these transitions are still poorly understood because of a lack of direct observations of bubble dynamics under natural magmatic conditions. Here, we combine x-ray radiography with a novel high-pressure/high-temperature apparatus to observe and quantify in real-time bubble growth and coalescence in basaltic magmas from 100 megapascals to surface. For low-viscosity magmas, bubbles coalesce and recover a spherical shape within 3 seconds, implying that, for lava fountaining activity, gas and melt remain coupled during the ascent up to the last hundred meters of the conduit. For higher-viscosity magmas, recovery times become longer, promoting connected bubble pathways. This apparatus opens frontiers in unraveling magmatic/volcanic processes, leading to improved hazard assessment and risk mitigation.

Direct observation of degassing and bubble coalescence in magmas reveals processes that control the style of volcanic eruptions.

http://dx.doi.org/10.13039/100010661 Horizon 2020 Framework Programme 101025887 http://dx.doi.org/10.13039/100011889 Diamond Light Source MG28538-1 http://dx.doi.org/10.13039/100014013 UK Research and Innovation MR/V023985/1 NextGenerationEU DM737/2021 RCUK NERC NE/N018575/1 French “grand emprunt” Equipex ANR-11-EQPX-36 Research Complex at Harwell CopyeditorJudith Urtula
==== Body
pmcINTRODUCTION

Volatile exsolution, expansion, and outgassing during magma ascent play a key control on the intensity and style of eruptions (1). The vesiculation process is a consequence of the supersaturation of dissolved volatiles within magma (caused, for example, by decreasing pressure during magma ascent or by magma crystallization), producing volatile exsolution (2–4). Bubbles of initially supercritical fluid, mostly containing water and lower amounts of carbon dioxide, halogens, sulfur, and other volatiles, nucleate and grow during magma ascent (5). During ascent, multiple processes occur: decompression-induced nucleation and growth of bubbles, coalescence of these bubbles with one another, and potentially outgassing, through channels formed by a network of coalescing bubbles (6–8). Ultimately, all of these processes play major roles in determining the eruptive style by influencing whether magma explosively fragments, produces lava fountains, or erupts effusively in lava flows (7, 9, 10). Understanding bubble coalescence is therefore crucial, given its role in controlling permeability (11, 12), the amount and rate of outgassing (13), and hence eruption dynamics (7).

As bubbles nucleate and grow, whether during eruptive magma ascent or experiments in the laboratory, interbubble melt films are thinned to the point of rupture, and bubble coalescence starts to occur (7, 14–20). Bubble coalescence in volcanology has been investigated by several studies, mainly focusing on the mechanisms that lead toward coalescence (7, 15, 21–28). In general, coalescence is a function of the viscosity of the surrounding melt, and it is an important process affecting bubble dynamics and, hence, magma ascent and transport to the surface (13). Bubble coalescence rates (i.e., how many times coalescence events occur per unit time) depend on the timescale of approach between two bubbles (29), thinning of the melt film separating two bubbles to a critical value, film rupture, and relaxation (21, 22, 27, 28, 30). Several studies focused on estimating the value of the critical film thickness in silicate melts (21, 31), as well as the timescales of film rupture once the critical thickness is reached (21, 22, 25–28, 30–32). The estimated rupture timescales range from 1 to 104 s, as a function of melt viscosity, bubble size, and the critical film thickness (21, 30). Bubble deformation plays also a key role on coalescence rate and on bubble population (23, 24, 33), and as a result, they can affect magma ascent dynamics. Other mathematical and numerical studies investigated more closely the dynamics of two merging bubbles in a viscous fluid, showing that the Ohnesorge number (a dimensionless number that relates the viscous forces to inertial and surface tension forces) plays a critical role on how quickly two merged bubbles recover the spherical shape (34–37). These studies, however, are not calibrated for magmatic fluids; thus, more studies are necessary to constrain these dynamics for volcanic systems.

Following the pioneering study of Sparks (1), degassing and bubble dynamics in silicic melts have been widely investigated by means of experiments (38–42) and numerical models (20, 43–47). Many studies have focused on the effect of decompression rate on bubble nucleation in silicic melts (48–54), while others on the parameters controlling bubble growth (1, 16, 38, 43, 45, 55, 56). Accurate insights into the exsolution process were obtained with experiments performed to study vesiculation in silicic melts at both 1 atm (42, 55, 57, 58) and high pressure (39, 41, 59). Recently, further developments concerned the use of in situ four-dimensional (4D) x-ray tomographic microscopy (where sequences of 3D tomographic scans are collected rapidly and continuously, creating a time series of 3D scans) to study vesiculation of silicate melts at high temperature, but at atmospheric pressure (60, 61). However, none of these studies were able to investigate the vesiculation process in basaltic melts in real time at pressures and temperatures comparable to those of an ascending basaltic magma from depth. Here, we combined x-ray synchrotron radiography with a novel x-ray transparent Internally Heated Pressure Vessel (IHPV) apparatus to simulate magma storage and ascent within the crust at pressures ≤100 MPa and temperatures ≤1180°C. With this apparatus, we performed in situ vesiculation experiments to study bubble growth and coalescence in a basaltic magma in real time at magmatic pressures and temperatures under water-saturated conditions.

Our experiments provide visualization and quantification of timescales of bubble formation (i.e., bubble growth, expansion, and coalescence) in real time, confirming and empirically validating theoretical and modeling results for bubble growth and expansion. The experimental results offer an improved understanding of coupling and decoupling between magma and volatiles during ascent in the conduit, providing insights into processes leading to eruptive style transitions and, ultimately, having fundamental implications for hazard assessment and risk mitigation in areas of active basaltic volcanism. Although basaltic volcanoes are usually characterized by effusive and mildly explosive Strombolian and lava fountaining activity (8), some basaltic volcanoes also produce highly explosive Plinian eruptions (62–68), with a much higher risk for population safety and critical infrastructure, as well as larger environmental impacts.

RESULTS AND DISCUSSION

High-pressure, high-temperature x-ray radiography experiments

Experiments were performed in situ at beamline I12-JEEP, Diamond Light Source, Harwell, UK, combining a novel x-ray transparent IHPV apparatus (fig. S1) with fast synchrotron x-ray radiography. We used a hydrous basaltic glass (table S1) from the 2001 Mt. Etna eruption as the starting material (see Materials and Methods).

We performed decompression experiments at superliquidus and subliquidus conditions to study the vesiculation process at different viscosities and crystallinities. To investigate vesiculation kinetics in pure basaltic melts with 1 wt % of H2O dissolved, we performed three decompression experiments at superliquidus (Superliq_Dec) conditions (i.e., no crystals) with decompression rates of 0.05 and 0.08 MPa s−1 representative of basaltic magma ascent rates during fountaining activity (69). The temperature was kept constant at 1180°C during decompression, above the bulk liquidus that is represented by clinopyroxene (table S2). To investigate vesiculation kinetics in crystal-bearing basaltic melts with 0.5, 1, and 2 wt % dissolved H2O, we performed five decompression experiments at subliquidus conditions (Subliq_Dec). This set of experiments is characterized by an initial cooling, with a cooling rate of 0.75°C s−1, at isobaric conditions (50 and 75 MPa) from 1180°C to different target temperatures (1050° to 1080°C) to promote crystallization before decompression. After subliquidus conditions were reached, the system was isothermally decompressed to 0.1 MPa with a decompression rate of 0.08 MPa s−1 to simulate different magma ascent rates during basaltic explosive or fountaining activities (69).

An additional decompression experiment at superliquidus condition was performed to investigate vesiculation kinetics in a hydrous rhyolitic melt (with 0.2 wt % of H2O); this experiment gives us the opportunity to compare vesiculation kinetics of basaltic and rhyolitic melts and to extrapolate the role of viscosity and crystals on bubble kinetics and dynamics. In particular, this allows us to compare experiments with a high-viscosity crystal-free magma (rhyolite at superliquidus temperature) with those with a less viscous melt but with similar bulk viscosity, due to the presence of crystals (Subliq_Dec experiments). Once all the experiments performed reached 0.1 MPa, the temperature was dropped to the atmospheric one with a continuous cooling of 0.75°C s−1 (fig. S2).

Radiographic images show that the vesiculation process of hydrous basaltic melts is quite different in the Superliq_Dec and the Subliq_Dec experiments (table S2). In the first group, we observe the growth of single bubbles that during the decompression paths tend to coalesce forming an individual bubble moving upwards to the top of the system setup (movies S1A and S2). A different behavior is observed in the rhyolite sample (Rhyo) in which bubbles remain confined in the melt, reflecting the higher viscosity compared with basalt. In this case, bubbles expand but do not separate from the melt as individual bubbles because of the higher viscosity (movie S1B). The Subliq_Dec experiments behave like the Rhyo as a consequence of the initial cooling (0.75°C s−1) that favors crystallization of microlites and a consequent increase in viscosity. In these runs, bubble expansion remains confined in the melt (movies S1C and S3), and outgassing is clearly observed through channels of interconnected bubbles, with a “breathing” pattern where groups of bubbles expand and then release gas through a pathway. While in the Superliq_Dec experiments it is possible to observe the nucleation of spherical bubbles that can easily grow and coalesce, in the Rhyo and Subliq_Dec samples, instead, bubbles are always deformed (e.g., bubbles repulsing each other to accommodate their volume increase), and coalescence occurs but to a lesser extent. In the Subliq_Dec samples, the presence of microlites seems to have a relevant impact on bubble behavior, because bubble expansion occurs following a less regular pattern than that observed in the Rhyo sample (movie S1), in which crystals are absent (fig. S6) because of its superliquidus conditions. The irregular bubble expansion seems dictated by the presence of microlites that physically impede bubble growth confining expansion in the less crystalline portion of the melt. The differences observed in the radiographic images are also visible in the textures of the recovered samples (see Supplementary Text and table S3).

High–temporal and high–spatial resolution bubble kinetics in basaltic magma

Our x-ray radiography experiments provided high–temporal and high–spatial resolution bubble kinetics (growth, expansion, and coalescence) in basaltic magma during decompression, simulating magma ascent in volcanic conduits. Once nucleated, bubbles grow by the combined effect of H2O diffusion from melt to bubble and decompression-induced gas expansion (16). Incremental bubble growth rate (ΔGR = Δr/Δt, μm s−1) was calculated in the basaltic runs of the Superliq_Dec group as the incremental increase in bubble radius (Δr) over time (Δt) (table S4). We measured the growth rate of bubbles that do not coalesce. In Fig. 1A, we plotted the evolution of ΔGR, as function of pressure for each of the Superliq_Dec experiments. Each dot represents ΔGR calculated for a given bubble, and the dashed lines indicate how ΔGR evolves as function of pressure for that specific bubble. Different dashed lines indicate the evolution of ΔGR for different bubbles. We observed that, at a given decompression rate, bubble growth rate increases with decreasing pressure (megapascals) following a power-law relation∆GR=1.428·P−0.69 for a decompression rate of 0.08 MPa s−1

∆GR=2.183·P−0.685 for a decompression rate of 0.05 MPa s−1

Fig. 1. Variation of incremental bubble growth rate and bubble radius with pressure.

Plots showing (A) incremental bubble growth rate (ΔGR) and (B) bubble radius versus pressure for the Superliq_Dec experiments. In (A), the decompression rate (megapascal per second) for each run is reported in parentheses. Insets show magnification of ΔGR (A) and bubble radius (B) at low pressure (P < 3 MPa). Both (A) and (B) show the comparison between the observed bubble growth rates and radius measured from the decompression experiments with those calculated (table S5) using the numerical bubble growth model of Coumans et al. (44). Numerical results are obtained assuming different decompression rates (0.05 and 0.08 MPa s−1) and different volatile contents (1 and 2 wt % of H2O). Plotted numerical simulations have been computed using the Etna composition (table S1), a magma density of 2700 kg m−3, and assuming Nb = 1012 m−3.

As shown in Fig. 1A, it is possible to see that the ΔGR vary from ~10−1 to ~10 μm s−1 passing from ~40 to 0.1 MPa, with a major rapid increase starting from ~10 MPa. Similarly to Fig. 1A, in Fig. 1B, we plotted the evolution of bubble radius as function of pressure for each of the Superliq_Dec experiments. Also in this case, each dot represents a bubble radius calculated for a given bubble, and the dashed lines indicate how the radius of a specific bubble evolves as function of pressure. Different dashed lines indicate the evolution of radii for different bubbles. Looking at the results, the bubble radius displays a similar trend as that observed for the bubble growth rate (Fig. 1B). In particular, for bubbles grown at different times during the decompression path, we noted that both bubble radius and growth rate increase more gently for bubbles nucleated at the beginning of the decompression (P = 30 to 50 MPa; ΔGR from ~0.1 to ~0.2 μm s−1 in 10 MPa), while they increase faster with a high slope for those nucleated at the end of the decompression path (P = 0.1 MPa; ΔGR from ~2 to ~8 μm s−1 in 0.1 MPa). This can be related to the effect of both the overpressure of bubbles compared to the pressure of the melt and the oversaturation of volatiles. The effect of both bubble overpressure and volatile oversaturation on bubble growth can be noted also comparing our data with those obtained by Masotta et al. (58) for bubbles grown in a basaltic melt at constant ambient pressure. We noticed that the values they obtained (GR = ~10−1 to 10−3 μm s−1) are approximately one to three orders of magnitude lower than those calculated in our experiments once bubbles reached ambient pressure (ΔGR ~ 10 μm s−1). Masotta et al. (58) observed fast bubble growth (GR = ~10−1 to 10 μm s−1), triggered by melt degassing, shortly after nucleation (t < 20 s), followed by a nearly linear growth (GR = ~10−3 to 10−1 μm s−1) for the rest of the experiment (t > 20 s). According to Masotta et al. (58), the fast initial growth, which they did not see directly during their in situ experiments because of the opacity of the sample, is consistent with the classical formulation of bubble growth with the bubble radius proportional to the square root of time (1, 45, 47) or the logarithmic growth law (39, 41). Thus, the GR calculated in this study once bubbles reached ambient pressure (ΔGR ~ 10 μm s−1), which we measured for time intervals <20 s (table S4), are representative of the GR at the very beginning of the exsolution process as reported by Masotta et al. (58). Another interesting comparison can be done with the bubble growth rates obtained by Bai et al. (70) for a basaltic melt through 1-atm in situ degassing experiments in which bubble growth is controlled by diffusion of the volatiles from the supersaturated melt to the bubble at constant pressure. The bubble growth rates, which we extrapolated from Bubble Size Distribution data reported in Bai et al. (70), show values between 10−2 and 100 μm s−1.

We also compared our experimentally derived growth rates with those calculated using the experimentally validated numerical model of bubble growth by Coumans et al. (44). This bubble growth model is based on the mathematical formulation of Blower et al. (11), Proussevitch et al. (31), and Proussevitch and Sahagian (20), and describes bubble expansion in a viscous magma due to decompression and diffusion. This model requires information on the bubble number density (Nb), which is not accessible from our radiography experiments. We performed some numerical simulations using Coumans’ bubble growth model, implementing it with the water diffusion model for basaltic melts (equation 22) of Zhang and Ni (71) and assuming common bubble number densities for basaltic magmas (Nb = 1010 to 1013 m−3) (54). We considered two decompression rates (0.05 and 0.08 MPa s−1) and two different water content (1 and 2 wt % of H2O). The results of these numerical simulations for Nb = 1012 m−3 are plotted in Fig. 1 (blue and green, dashed and solid lines), and they show a good agreement with our observed growth rates except at low pressure (<1 MPa), where modeled growth rates exceed 10 μm s−1 (Fig. 1A, fig. S3, and table S5). However, Coumans et al. (44) reported that the numerical model overestimates bubble growth rates at high gas volume fraction (>0.4), which would explain the discrepancy with the observed values at low pressures. In addition, our growth rates at very low pressures (<1 MPa) might be underestimated because of the effect of the walls of the crucible, which exert a resistance on the melt to deform and flow as bubbles expand, resulting in a reduced expansion of bubbles.

In the Subliq_Dec runs and in the Rhyo sample, bubble growth and coalescence are hampered by the higher viscosity of the melt and not easily resolvable in the radiographic images. Consequently, we could not extrapolate any quantitative data on these processes. However, because the radiographic images show an expansion of the trapped bubbles, it was possible to obtain the incremental bubble expansion rate (ΔER = ΔA/Δt; square micrometer per second) as the incremental increase of the bubble area (ΔA) with time (Δt). For comparison, we calculated bubble ΔER also in the Superliq_Dec runs that show values approximately one to three orders of magnitude higher than those of the Subliq_Dec ones as a consequence of their different viscosity before the decompression path (Fig. 2 and table S6; see also Supplementary Text for more details).

Fig. 2. Variation of incremental expansion rate with pressure.

Plot showing incremental expansion rate (ΔER) versus pressure for (A) Subliq_Dec and (B) Superliq_Dec experiments. In (A), the liquid viscosity (log Pa·s) for each run is reported in parentheses.

In the Superliq_Dec runs, it was possible to observe in real time and identify several steps and related times leading to bubble coalescence (table S7 and Fig. 3) including the following: (1) time of contact, when two separate bubbles enter in contact; (2) time of interconnection, when two bubbles are interconnected and their films start to thin; (3) time of coalescence, when there is an open connection between the bubbles as a result of the rupture of the thinned films; and (4) time of recovery, when the coalesced bubbles recover to a spherical shape. All of these steps can be observed in detail in Fig. 3. A similar sequence for bubble coalescence has been also observed by Masotta et al. (58) during in situ high-temperature observations of bubble growth in a rhyodacitic melt and by Ohashi et al. (27) during in situ decompression experiments of viscous Newtonian analogues in a mini-desiccator box. Some textural features of bubble coalescence [i.e., bubble-melt wall thinning by bending, stretching, and dimpling (7), corresponding to steps 1 and 2 of this study] have been observed in previous ex situ experimental decompression/vesiculation studies on both basaltic (49, 50) and rhyolitic (16) compositions. As shown in this study, once two bubbles enter in contact (step 1) there is an interconnection (step 2) during which bubble walls thin until the film ruptures. During interconnection bubbles assume an “eight” shape with cusps on their walls. Because the film thickness is too thin to be resolved by radiography, we cannot see the rupture of the films in real time, but we can ascertain that it has occurred when we observe the replacement of the cusps in the “eight” shape by smooth bumps on the walls of the new coalesced bubble, and so bubbles coalescence (step 3). Once coalesced, bubbles assume an oval shape and lastly recover to a spherical one (step 4). The recovery time (i.e., the time required by coalesced bubbles to recover a spherical shape) is in the order of 1 to 3 s and results to be affected by pressure. We noted, indeed, an increase of the recovery time with decreasing pressure (Fig. 4). The recovery time (τ) for a coalesced pair of bubbles to return to a spherical shape can be compared with theoretical estimations (7, 21, 32, 72–74). These allow us to predict τ from the rheological properties of the liquid, because retraction of the common wall between two interconnected bubbles is fostered by surface tension (σ) and resisted by the effective viscosity of the liquid (μ) ([τ = (Rμ)/σ]). Theoretical results (table S7) were obtained assuming σ = 0.1 N m−1 (7, 73), μ = 148 Pa·s [table S3; calculated using the model of Giordano et al. (75)], and R as the equivalent bubble radius (i.e., the radius of an undeformed sphere of equal volume; table S7). From the comparison, we noticed that the theoretical recovery times are faster than the experimental ones with a difference up to one order of magnitude at the lowest pressures (P < 3 MPa). However, fully understanding this discrepancy requires targeted and in-depth studies that are beyond the aim of this study.

Fig. 3. Radiographic images of Superliq_Dec runs showing coalescence steps and timescales.

(1) contact between two bubbles; (2) bubble interconnection during which bubble walls thin (bubbles assume an “eight” shape with cusps on their walls); (3) bubble coalescence after the rupture of the thinned films (replacement of the cusps in the “eight” shape by smooth bumps on the walls of the new coalesced bubble); (4) bubbles assume an oval shape and lastly recover to a spherical one.

Fig. 4. Sketch reporting bubble growth and coalescence within the conduit and mechanisms of coupling (up to the last ~100 m) and decoupling between volatiles and magma.

Plots show bubble recovery time and magma time to surface versus pressure. In particular, the time required by magma to reach the surface for each pressure within the conduit using the lava fountaining simulations performed by La Spina et al. (69) for Kilauea is reported.

Coalescence and degassing were also observed in the Subliq_Dec runs (movies S1C and S3). Evidence of coalescence events can be also found in the recovered samples. We can presume that large bubbles in backscattered electron (BSE) images (fig. S4, D to G) are the product of coalescence between two or more smaller bubbles. The presence of these coalesced bubbles implies that, at the end of the decompression in the Subliq_Dec runs, the formation of single large bubbles moving toward the top of the apparatus did not occur as in the Superliq_Dec experiments. This suggests that the formation of permeable pathways in the Superliq_Dec experiments allowed gas to escape, which can be related to the presence of microlite crystals. The role of microlite crystals on coalescence and formation of permeable pathways is still a debated topic. On the one hand, the presence of microlite crystals could lead to the following: (i) an increase in melt viscosity, as their presence indicates an increased silica content of the melt; (ii) an increase in bulk viscosity, which reduces film drainage ability and consequently the rate of coalescence (13); and (iii) a physical block or impediment of bubble expansion, movement, or coalescence. On the other hand, however, some studies (76–78) observed higher permeabilities and thus a more efficient degassing in crystal-bearing than in crystal-free basaltic magmas. This would suggest that the presence of crystals, forcing the bubbles concentration in some regions of the melt rather than in others, would facilitate the contact between the bubbles and therefore the formation of permeable pathways. The presence of crystals, slowing down the recovery time because of a higher viscosity, would favor a longer opening of permeable pathways and therefore a more efficient degassing.

Influence of bubble coalescence on degassing

Volcanoes such as Etna, Stromboli, and Kilauea are generally characterized by the ascent of basaltic magma that involves different degrees of decoupled, open-system degassing, in which volatiles are able to flow faster than their originating, slowly ascending melt (69, 79–81). In the presence of more viscous melts (such as rhyolitic melts), bubbles are relatively immobile with respect to the melt (i.e., the slip velocity of bubbles is negligible compared to magma ascent velocity), and this is commonly referred to as coupled or closed-system degassing (82, 83), although in rhyolites open pathways may form, which allow degassing to take place (84, 85). The style of basaltic explosive behavior is strongly related to the ability of gas to decouple from the melt, which depends on the relative rates of ascent of melt and bubbles, establishment of percolation bubble frameworks, and the extent of bubble coalescence (80, 82, 86, 87).

While slow ascent rates, allowing gas bubbles to decouple from the melt, move, and coalesce, increase the likelihood of a Strombolian or effusive eruption (79, 88), fast ascent rates such as estimated for basaltic Plinian and sub-Plinian eruptions [with average rates ranging between 4 and 75 m s−1; (62, 89, 90)] likely restrict the time available for gas-melt decoupling during magma ascent, promoting, depending on magma viscosity, fragmentation and an explosive eruption or high-intensity fountain activity (66, 68, 69, 80, 91–94).

In this study, we measured bubble recovery time in the Superliq_Dec runs and the ascent time of Kilauea magma to investigate the role of bubble connections on degassing [i.e., both closed-system degassing (volatile exsolution and vesiculation) and open-system degassing (outgassing); (95)] of low-viscosity crystal-poor/crystal-free basaltic magmas. The quicker this process, the more likely a permeable pathway closes before connecting with the vent of the conduit, preventing gas from decoupling efficiently from the melt. For basaltic lava fountaining activity on Kilauea, La Spina et al. (69) show a maximum ascent velocity of ~60 m s−1 [average velocity ~ 15 m s−1 consistent with Ferguson et al. (96)], which means that magma would need more than 3 s to reach the surface from a depth greater than 200 m. Our obtained recovery timescales are on the order of ~3 s at low magmatic pressures, but they are even shorter at higher pressures, implying that the recovery of spherical shape at depth may occur too quickly to allow permeable pathways consisting of chains of interconnected bubbles to reach the surface and connect with the vent during fountaining. To better show this, we used results from the lava fountaining simulations performed by La Spina et al. (69) for Kilauea, and we used their model to extrapolate the time required by magma to reach the surface for each pressure. From Fig. 4 and table S8, we note that the time to reach the surface becomes greater than 3 s when the pressure is greater than ~4 MPa (~90-m depth). This implies that, for pressures higher than ~4 MPa, the formation of permeable pathways is mostly inhibited, and thus, the gas-melt system remains coupled at least up to the last 100 m of the conduit (Fig. 4). In general, our results show that lava fountaining eruptions at low-viscosity basaltic volcanoes (such as at Kilauea) are associated with rapid bubble coalescence and recovery time and a high magma ascent rate, resulting in a coupled behavior until the last ~100 m of the conduit (69, 79, 80). We speculate that some degree of decoupling may occur in the shallowest part of the conduit; otherwise, it is likely that intense fragmentation would be produced by closed-system degassing. This has been observed, for instance, in basaltic pumices from the paroxysmal events of Stromboli where bubble populations are consistent with closed-system degassing. Close to the surface, where the gas volume fraction is high, the occurrence of multiple events of coalescence at the same time may generate permeable pathways allowing some degree of open-system degassing. Such a process may contribute to near-surface open-system degassing in slow-ascending magmas, such as within a lava lake.

Another possible mechanism for gas-magma decoupling that might happen during ascent is the formation of a slug flow due to the coalescence of numerous smaller bubbles (97, 98). Our Superliq_Dec experiments show that, at low pressures, the increase in coalescence events produce ultimately a very large bubble, which indicates that the formation of a slug flow in the shallowest part of the conduit is possible. A high gas volume fraction (>0.6) (99) is required to reach the size of the conduit and develop a slug flow. For a low–volatile content magma, such as that of Kilauea, this high gas volume fraction is obtained in the shallowest part of the conduit, where most of the exsolution occurs. The formation of slug flow may occur also at depth due to accumulation of gas entrapped at some geometrical discontinuities within the plumbing system (97), but this would generate Strombolian, rather than lava fountain, activity (81, 97, 98).

Although in the Subliq_Dec runs we cannot visualize in detail the steps that lead to bubble coalescence through radiography, we can, however, obtain some insights on magma-gas coupling/decoupling. We observe that, in the presence of a more viscous magma (η ≥ 103 Pa·s), in case of a rupture of the film between bubbles, the recovery time is much longer than that obtained for the less viscous Superliq_Dec runs. This longer recovery time promotes the formation of degassing pathways that allow gas to escape. An increase in magma viscosity, indeed, can affect the eruptive behavior in several ways such as increasing the fragmentation capability due to bubble overpressure (100, 101), suppressing large bubble floatation, and increasing the capillary number and the role of shear deformation (102, 103). Our novel experiments with an in situ view show a slow bubble expansion rate and consequently a long recovery timescale at subliquidus conditions (movies S1C and S3), which supports assumptions of previous studies (72, 77, 104, 105). This result suggests that high viscosity produced by microlite crystallization restricts bubble growth and expansion and extends bubble coalescence time, promoting connected pathways between bubbles (as those visible in movie S3 from minute 15 onwards) and thus increased connectivity, which, in turn, favor outgassing (72, 77, 95, 106, 107).

Our work represents a substantial step forward in the understanding of magma and gas dynamics, even though it has limitations that result from the use of 2D radiography and inability to apply a shear stress (which are likely to affect bubble growth and coalescence). Thanks to the new IHPV apparatus presented here, we were able to capture and study the vesiculation kinetics in basaltic magmas in real time, in situ, and at pressures and temperatures compatible with those of basaltic volcanoes. Our novel x-ray transparent apparatus has proved to be an invaluable tool to capture and quantify kinetic of bubble formation (i.e., bubble growth, expansion, and coalescence) and magma dynamics (i.e., degassing and gas-magma coupling/decoupling) at syneruptive conditions. In this regard, the growth rates derived from our experiments represent a noteworthy contribution, as they confirm estimations calculated using numerical and theoretical models.

Future developments of the x-ray transparent IHPV will be dedicated to allowing fast synchrotron x-ray tomography of magmatic samples at high-pressure and high-temperature conditions, to visualize and quantify the vesiculation process (nucleation, growth, and coalescence) directly in 4D (3D space plus time) both at superliquidus and subliquidus conditions. This will allow us to further improve first the current numerical model by integrating previously unknown constraints, and then our understanding of magma behavior at pre- and syneruptive conditions and the related volcanic hazard.

MATERIALS AND METHODS

Starting material

The starting material used for our vesiculation experiments is a trachybasalt from the lower vents of the 2001 Mt. Etna eruption (62, 108–111). We used hydrous, crystal-free basaltic samples from Etna with different water contents (0.5 to 2 wt %) and a rhyolitic sample with 0.2 wt % of water. The anhydrous glassy starting material (table S1) was synthesized by melting a crushed rock sample in a Pt crucible. Melting was performed in a Nabertherm MoSi2 box furnace at 1400°C and at atmospheric pressure. The melt was left in the furnace for 4 hours to fully degas and dissolve any crystals present. The melt was then quenched in air to glass, and this procedure was repeated twice to enhance homogenization. Hydrous starting glasses with 0.5, 1, and 2 wt % H2O were obtained by melting the starting material and homogenizing it with H2O in Au80Pd20 capsules at 100 MPa and 1200°C using a Titanium Zirconium Molybdenum (TZM) cold-seal pressure vessel apparatus at the School of Earth Sciences, University of Bristol, UK. The water content of the starting materials was confirmed to be present in the glasses by Fourier Transform Infrared spectroscopy. The FTIR measurements were performed in transmission mode by a PerkinElmer Spotlight 400 spectrometer equipped with a Mercury Cadmium Telluride (MCT; or HgCdTe) array detector cooled with liquid N2 at the Department of Earth and Environmental Sciences at the University of Manchester, UK. Spectra were collected by accumulating 64 scans using a square aperture of 100 μm across with a spectrum resolution of 4 cm−1. Spectra were analyzed using Spectragryph (112). Using the density trend and Etna basalt extinction coefficient of Testemale et al. (113), the molecular H2O peak at 3550 cm−1 gives 2.05 ± 0.01 wt % for Etna 4, 1.32 ± 0.01 wt % for Etna 12, 1.48 ± 0.01 wt % for Etna 13 and 1.13 ± 0.01 wt % for Etna 15 starting materials.

In situ high-pressure, high-temperature synchrotron x-ray radiography experimental apparatus

In situ high-pressure, high-temperature experiments were performed at the x-ray tomography/radiography beamline I12-JEEP, Diamond Light Source, Harwell, UK. We used a dedicated x-ray transparent IHPV apparatus developed at Neel Institute and based on a previous one (113) combined with x-ray radiography to perform in situ vesiculation experiments under water-saturated conditions at crustal pressures. The IHPV apparatus was pressurized with He, which allowed us to precisely control the decompression rate during in situ experiments and to quantify disequilibrium in basaltic magmas after pressure perturbations. The pressurization was controlled by a pressure regulator (114). The vessel is characterized by the placement of the furnace inside the vessel (internally heated). The vessel is a thick-walled steel cylinder having both ends open. The open ends are closed by heads through which pressure, electrical lead, and thermocouple lead enter. The vessel has two sapphire windows at 180°, which allow the x-ray beam to enter the vessel, passing through the sample and reaching the camera for radiography acquisitions. Temperature was measured with a K-type thermocouple positioned close to the sample in the middle of the furnace hotspot. The K-type thermocouple measures the sample temperature with an uncertainty of ±0.5°C. The sample holder was an alumina, which is suitable for the temperature range investigated and has a low x-ray attenuation coefficient. The hydrous glass (~1.5 mm by 3 mm by 5 mm; ~22 mm3 once melted) was placed in the customized alumina crucible whose sizes are reported in fig. S1.

Experimental strategy

We combined fast x-ray synchrotron radiography with our novel IHPV apparatus to quantify bubble growth and coalescence in basaltic magmas during decompression. The experiments focused on bubble kinetics as a function of initial pressure, decompression rate, and H2O content. In all the experiments, we placed an Etna basalt with approximately 0.5 to 2 wt % of water in the sample crucible, except for the one in which we used a rhyolitic sample. We pressurized the system at first with gas (He), and then we heated up to 1180°C with a heating rate of 0.75°C s−1. At this point, we continued the experiments by keeping isothermal conditions (1180°C; Superliq_Dec experiments; table S2) or dropping the temperature to different target isothermal conditions (1050° to 1080°C; Subliq_Dec experiments; table S2) with a cooling rate of 0.75°C s−1. After that, we dropped the pressure to 0.1 MPa with a decompression rate between 0.03 and 0.08 MPa s−1 to simulate different ascent rates during basaltic eruptions, starting decompression at different initial pressures (75, 50, 30, and 20 MPa; table S2). Once reached 0.1 MPa, the temperature was dropped to the ambient one with a cooling rate of 0.75°C s−1 (fig. S2).

In situ synchrotron x-ray radiography acquisition

The x-ray radiography beamline I12-JEEP (Diamond Light Source, Harwell, UK) allowed us to perform experiments using monochromatic 53-keV x-rays, a pixel size of 6.642 μm, and a scanning time of 40 ms per frame to achieve 25 frames per second, at a sample to detector distance of 35 cm, an exposure time for a single projection of 15 ms, and an acquisition time of 988.28 s, for a total of 24,707 images. The acquisition of radiographic projections began shortly before the start of decompression and covered the entire decompression path until ambient pressure was reached.

Image processing and analysis

The radiographic images were processed and stacked using ImageJ software (115) to obtain movies (movie S1). Movies were made by importing in ImageJ the radiographic images in TIFF format as image sequences and then saved as an AVI file. The movies reported in movies S1 to S3 were edited by using DaVinci Resolve (version 18.0.4) video editing software. ImageJ was also used to measure bubble diameter and area. First of all, for all the runs, we used “set scale” to convert pixel in micrometers (1 pixel = 6.642 μm); then frames were converted in 8-bit and then adjusted by brightness/contrast to better highlight bubbles from melts. Because of the low contrast between bubbles and melt and to bubble overlays, it was not possible to use the “threshold” tool and the “tracking plugin.” For the Superliq_Dec runs bubbles diameter and area were measured by manually tracking bubbles using “oval selections” and then the “measure” tool. To better highlight edges of the bubbles the “Find edges” tool was also applied. For the Subliq_Dec runs, instead, because of the absence of spherical bubbles, bubble area was measured manually by contouring bubble edges using the “polygon selections” and then the measure tool.

Scanning electron microscope and electron microprobe analysis

BSE images were collected using a FEI Quanta 650 FEG-SEM electron microscope in the Department of Earth and Environmental Sciences, University of Manchester, UK, to analyze vesicles shapes and crystals morphologies. We used an acceleration voltage of 15 kV and a working distance of 10 mm. The starting material (glass) and the samples obtained during in situ vesiculation experiments were analyzed with a JEOL JXA-8530F field-emission electron microprobe at the Photon Science Institute, University of Manchester, UK. The operating conditions were as follows: 15-kV accelerating voltage, 10-nA beam current, and a beam diameter of 10 or 5 μm. Na and K were measured first to minimize loss by volatilization. Calibration standards were albite for Na, periclase for Mg, corundum for Al, fayalite for Fe, tephroite for Mn, apatite for P, sanidine for K, wollastonite for Ca and Si, and rutile for Ti.

Acknowledgments

We thank E. Llewellin for fruitful discussions on this work. We gratefully acknowledge support and friendship from H. Mader, who passed away in December 2022. Mader was a leader in the project which produced this work, and we will be always grateful for her warmth, creativity, and thoughtfulness.

Funding: The research leading to these results has received funding from the RCUK NERC Dis-Eqm project (NE/N018575/1) and from the UKRI FLF project 4DVOLC (MR/V023985/1). The beamtime on I12 was provided by Diamond Light Source (MG28538-1) and laboratory space by the Research Complex at Harwell. The IHPV apparatus was financially supported by French “grand emprunt” Equipex (Planex, ANR-11-EQPX-36). We thank the technical team of the Néel Institute (A. Prat, E. Lahera, W. Del Net, and S. Min) for the realization of this new autoclave and its control system. F.A. thanks the KINETICVOLC project funded by the European Union, NextGenerationEU, Ministry of University and Research, Promotion and Development Fund DM737/2021. This project has received funding from the European Union’s Horizon 2020 research and innovation program under the Marie Skłodowska-Curie grant agreement no. 101025887.

Author contributions: F.A., M.P., G.L.S., M.R.B., and J.-L.H. conceived the research project. J.-L.H. developed the original IHPV apparatus, and F.A., M.R.B., M.P., and R.A.B. worked to obtain its modified version. F.A., M.P., G.L.S., J.-L.H., M.R.B., P.D.L., S.M., and R.A. contributed to x-ray radiography experiments. F.A. collected the volcanic rocks for the starting material. R.A.B. prepared the starting material. B.B. performed radiographic image processing and image analysis. R.A.W. and B.B. performed FTIR measurements and data processing. J.F. and B.B. performed EPMA analyses and data processing. B.B. and M.P. acquired back-scattered electron images of experimental samples. B.B., M.P., F.A., and G.L.S. wrote the manuscript, with contributions from all other authors.

Competing interests: The authors declare that they have no competing interests.

Data and materials availability: All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials.

Supplementary Materials

This PDF file includes:

Supplementary Text

Figs. S1 to S7

Legends for tables S1 to S8

Legends for movies S1 to S3

References

Other Supplementary Material for this manuscript includes the following:

Tables S1 to S8

Movies S1 to S3
==== Refs
REFERENCES AND NOTES

1 R. S. J. Sparks, The dynamics of bubble formation and growth in magmas: A review and analysis. J. Volcanol. Geotherm. Res. 3 , 1–37 (1978).
2 Y. Liu, Y. Zhang, H. Behrens, Solubility of H2O in rhyolitic melts at low pressures and a new empirical model for mixed H2O–CO2 solubility in rhyolitic melts. J. Volcanol. Geotherm. Res. 143 , 219–235 (2005).
3 A. Toramaru, Measurement of bubble size distributions in vesiculated rocks with implications for quantitative estimation of eruption processes. J. Volcanol. Geotherm. Res. 43 , 71–90 (1990).
4 Y. Zhang, Z. Xu, M. Zhu, H. Wang, Silicate melt properties and volcanic eruptions. Rev. Geophys. 45 , RG4004 (2007).
5 M. T. Mangan, K. V. Cashman, S. Newman, Vesiculation of basaltic magma during eruption. Geology 21 , 157–160 (1993).
6 M. R. Burton, H. M. Mader, M. Polacci, The role of gas percolation in quiescent degassing of persistently active basaltic volcanoes. Earth Planet. Sci. Lett. 264 , 46–60 (2007).
7 J. M. Castro, A. Burgisser, C. I. Schipper, S. Mancini, Mechanisms of bubble coalescence in silicic magmas. Bull. Volcanol. 74 , 2339–2352 (2012).
8 M. Polacci, R. A. Corsaro, D. Andronico, Coupled textural and compositional characterization of basaltic scoria: Insights into the transition from Strombolian to fire fountain activity at Mount Etna Italy. Geology 34 , 201–204 (2006).
9 S. Lovejoy, H. Gaonac’h, D. Schertzer, Bubble distributions and dynamics: The expansion-coalescence equation. J. Geophys. Res. Solid Earth 109 , B11203 (2004).
10 H. R. Westrich, J. C. Eichelberger, Gas transport and bubble collapse in rhyolitic magma: An experimental approach. Bull. Volcanol. 56 , 447–458 (1994).
11 J. D. Blower, J. P. Keating, H. M. Mader, J. C. Phillips, Inferring volcanic degassing processes from vesicle size distributions. Geophys. Res. Lett. 28 , 347–350 (2001).
12 J. C. Eichelberger, C. R. Carrigan, H. R. Westrich, R. H. Price, Non-explosive silicic volcanism. Nature 323 , 598–602 (1986).
13 H. M. Gonnermann, M. Manga, The fluid mechanics inside a volcano. Annu. Rev. Fluid Mech. 39 , 321–356 (2007).
14 J. E. Gardner, Bubble coalescence in rhyolitic melts during decompression from high pressure. J. Volcanol. Geotherm. Res. 166 , 161–176 (2007).
15 T. Giachetti, H. M. Gonnermann, J. E. Gardner, A. Burgisser, S. Hajimirza, T. C. Earley, N. Truong, P. Toledo, Bubble coalescence and percolation threshold in expanding rhyolitic magma. Geochem. Geophys. Geosyst. 20 , 1054–1074 (2019).
16 C. Martel, G. Iacono-Marziano, Timescales of bubble coalescence, outgassing, and foam collapse in decompressed rhyolitic melts. Earth Planet. Sci. Lett. 412 , 173–185 (2015).
17 D. S. Martula, T. Hasegawa, D. R. Lloyd, R. T. Bonnecaze, Coalescence-induced coalescence of inviscid droplets in a viscous fluid. J. Colloid Interface Sci. 232 , 241–253 (2000).11097757
18 C. T. Nguyen, H. M. Gonnermann, Y. Chen, C. Huber, A. A. Maiorano, A. Gouldstone, J. Dufek, Film drainage and the lifetime of bubbles. Geochem. Geophys. Geosyst. 14 , 3616–3631 (2013).
19 S. Okumura, M. Nakamura, A. Tsuchiyama, Shear-induced bubble coalescence in rhyolitic melts with low vesicularity. Geophys. Res. Lett. 33 , L20316 (2006).
20 A. A. Proussevitch, D. L. Sahagian, Dynamics of coupled diffusive and decompressive bubble growth in magmatic systems. J. Geophys. Res. Solid Earth 101 , 17447–17455 (1996).
21 C. Klug, K. V. Cashman, Permeability development in vesiculating magmas: Implications for fragmentation. Bull. Volcanol. 58 , 87–100 (1996).
22 O. Navon, V. Lyakhovsky, Vesiculation processes in silicic magmas. Geol. Soc. Spec. Publ. 145 , 27–50 (1998).
23 M. Manga, H. A. Stone, Interactions between bubbles in magmas and lavas: Effects of bubble deformation. J. Volcanol. Geotherm. Res. 63 , 267–279 (1994).
24 M. Manga, H. A. Stone, Collective hydrodynamics of deformable drops and bubbles in dilute low Reynolds number suspensions. J. Fluid Mech. 300 , 231–263 (1995).
25 C. Huber, Y. Su, C. T. Nguyen, A. Parmigiani, H. M. Gonnermann, J. Dufek, A new bubble dynamics model to study bubble growth, deformation, and coalescence. J. Geophys. Res. Earth 119 , 216–239 (2014).
26 S. Mancini, L. Forestier-Coste, A. Burgisser, F. James, J. Castro, An expansion–coalescence model to track gas bubble populations in magmas. J. Volcanol. Geotherm. Res. 313 , 44–58 (2016).
27 M. Ohashi, T. Maruishi, A. Toramaru, Coalescence of growing bubbles in highly viscous liquids. Geochem. Geophys. Geosyst. 23 , e2022GC010618 (2022).
28 M. Ohashi, A. Toramaru, A. Namiki, Coalescence of two growing bubbles in a Hele–Shaw cell. Sci. Rep. 12 , 1270 (2022).35075182
29 I. Fortelný, A. Živný, J. JŮza, Coarsening of the phase structure in immiscible polymer blends. Coalescence or Ostwald ripening? J. Polym. Sci. B Polym. Phys. 37 , 181–187 (1999).
30 J. F. Larsen, M.-H. Denis, J. E. Gardner, Experimental study of bubble coalescence in rhyolitic and phonolitic melts. Geochim. Cosmochim. Acta 68 , 333–344 (2004).
31 A. A. Proussevitch, D. L. Sahagian, V. A. Kutolin, Stability of foams in silicate melts. J. Volcanol. Geotherm. Res. 59 , 161–178 (1993).
32 R. A. Herd, H. Pinkerton, Bubble coalescence in basaltic lava: Its impact on the evolution of bubble populations. J. Volcanol. Geotherm. Res. 75 , 137–157 (1997).
33 L. Karlstrom, E. M. Dunham, Excitation and resonance of acoustic-gravity waves in a column of stratified, bubbly magma. J. Fluid Mech. 797 , 431–470 (2016).
34 J. D. Paulsen, R. Carmigniani, A. Kannan, J. C. Burton, S. R. Nagel, Coalescence of bubbles and drops in an outer fluid. Nat. Commun. 5 , 3182 (2014).24458225
35 C. R. Anthony, P. M. Kamat, S. S. Thete, J. P. Munro, J. R. Lister, M. T. Harris, O. A. Basaran, Scaling laws and dynamics of bubble coalescence. Phys. Rev. Fluids 2 , 083601 (2017).
36 R. Chen, J. Zeng, H. Yu, Mechanism of damped oscillation in microbubble coalescence. Comput. Fluids 183 , 38–42 (2019).
37 R. Chen, H. W. Yu, J. Zeng, L. Zhu, General power-law temporal scaling for unequal-size microbubble coalescence. Phys. Rev. E 101 , 023106 (2020).32168553
38 J. E. Gardner, M. Hilton, M. R. Carroll, Experimental constraints on degassing of magma: Isothermal bubble growth during continuous decompression from high pressure. Earth Planet. Sci. Lett. 168 , 201–218 (1999).
39 C. Gondé, C. Martel, M. Pichavant, H. Bureau, In situ bubble vesiculation in silicic magmas. Am. Mineral. 96 , 111–124 (2011).
40 J. F. Larsen, J. E. Gardner, Experimental constraints on bubble interactions in rhyolite melts: Implications for vesicle size distributions. Earth Planet. Sci. Lett. 180 , 201–214 (2000).
41 C. Martel, H. Bureau, In situ high-pressure and high-temperature bubble growth in silicic melts. Earth Planet. Sci. Lett. 191 , 115–127 (2001).
42 O. Navon, A. Chekhmir, V. Lyakhovsky, Bubble growth in highly viscous melts: Theory, experiments, and autoexplosivity of dome lavas. Earth Planet. Sci. Lett. 160 , 763–776 (1998).
43 J. Barclay, D. S. Riley, R. S. J. Sparks, Analytical models for bubble growth during decompression of high viscosity magmas. Bull. Volcanol. 57 , 422–431 (1995).
44 J. P. Coumans, E. W. Llewellin, F. B. Wadsworth, M. C. S. Humphreys, S. A. Mathias, B. M. Yelverton, J. E. Gardner, An experimentally validated numerical model for bubble growth in magma. J. Volcanol. Geotherm. Res. 402 , 107002 (2020).
45 A. A. Proussevitch, D. L. Sahagian, Dynamics and energetics of bubble growth in magmas: Analytical formulation and numerical modeling. J. Geophys. Res. Solid Earth 103 , 18223–18251 (1998).
46 A. Toramaru, Numerical study of nucleation and growth of bubbles in viscous magmas. J. Geophys. Res. 100 , 1913–1931 (1995).
47 A. Toramaru, Vesiculation process and bubble size distributions in ascending magmas with constant velocities. J. Geophys. Res. Solid Earth 94 , 17523–17542 (1989).
48 M. Hamada, D. Laporte, N. Cluzel, K. T. Koga, T. Kawamoto, Simulating bubble number density of rhyolitic pumices from Plinian eruptions: Constraints from fast decompression experiments. Bull. Volcanol. 72 , 735–746 (2010).
49 N. Le Gall, M. Pichavant, Experimental simulation of bubble nucleation and magma ascent in basaltic systems: Implications for Stromboli volcano. Am. Mineral. 101 , 1967–1985 (2016).
50 N. Le Gall, M. Pichavant, Homogeneous bubble nucleation in H2O- and H2O-CO2-bearing basaltic melts: Results of high temperature decompression experiments. J. Volcanol. Geotherm. Res. 327 , 604–621 (2016).
51 N. G. Lensky, O. Navon, V. Lyakhovsky, Bubble growth during decompression of magma: Experimental and theoretical investigation. J. Volcanol. Geotherm. Res. 129 , 7–22 (2004).
52 M. Mangan, T. Sisson, Delayed, disequilibrium degassing in rhyolite magma: Decompression experiments and implications for explosive volcanism. Earth Planet. Sci. Lett. 183 , 441–455 (2000).
53 C. C. Mourtada-Bonnefoi, D. Laporte, Kinetics of bubble nucleation in a rhyolitic melt: An experimental study of the effect of ascent rate. Earth Planet. Sci. Lett. 218 , 521–537 (2004).
54 A. Toramaru, BND (bubble number density) decompression rate meter for explosive volcanic eruptions. J. Volcanol. Geotherm. Res. 154 , 303–316 (2006).
55 Y. Liu, Y. Zhang, Bubble growth in rhyolitic melt. Earth Planet. Sci. Lett. 181 , 251–264 (2000).
56 V. Lyakhovsky, S. Hurwitz, O. Navon, Bubble growth in rhyolitic melts: Experimental and numerical investigation. Bull. Volcanol. 58 , 19–32 (1996).
57 N. S. Bagdassarov, D. B. Dingwell, M. C. Wilding, Rhyolite magma degassing: An experimental study of melt vesiculation. Bull. Volcanol. 57 , 587–601 (1996).
58 M. Masotta, H. Ni, H. Keppler, In situ observations of bubble growth in basaltic, andesitic and rhyodacitic melts. Contrib. Mineral. Petrol. 167 , 1–14 (2014).
59 C. Gondé, D. Massare, H. Bureau, C. Martel, M. Pichavant, R. Clocchiatti, In situ study of magmatic processes: A new experimental approach. High Press. Res. 26 , 243–250 (2006).
60 D. R. Baker, F. Brun, L. Mancini, J. L. Fife, A. LaRue, C. O’Shaughnessy, R. J. Hill, M. Polacci, The importance of pore throats in controlling the permeability of magmatic foams. Bull. Volcanol. 81 , 54 (2019).
61 D. R. Baker, F. Brun, C. O’Shaughnessy, L. Mancini, J. L. Fife, M. Rivers, A four-dimensional X-ray tomographic microscopy study of bubble growth in basaltic foam. Nat. Commun. 3 , 1135 (2012).23072805
62 F. Arzilli, G. La Spina, M. R. Burton, M. Polacci, N. Le Gall, M. E. Hartley, D. Di Genova, B. Cai, N. T. Vo, E. C. Bamber, S. Nonni, R. Atwood, E. W. Llewellin, R. A. Brooker, H. M. Mader, P. D. Lee, Magma fragmentation in highly explosive basaltic eruptions induced by rapid crystallization. Nat. Geosci. 12 , 1023–1028 (2019).
63 E. C. Bamber, G. La Spina, F. Arzilli, M. de Michieli Vitturi, M. Polacci, M. E. Hartley, M. Petrelli, J. Fellowes, M. Burton, Basaltic Plinian eruptions at Las Sierras-Masaya volcano driven by cool storage of crystal-rich magmas. Commun. Earth Environ. 3 , 253 (2022).
64 E. C. Bamber, F. Arzilli, M. Polacci, M. E. Hartley, J. Fellowes, D. Di Genova, D. Chavarría, J. A. Saballos, M. R. Burton, Pre- and syn-eruptive conditions of a basaltic Plinian eruption at Masaya Volcano, Nicaragua: The Masaya Triple Layer (2.1 ka). J. Volcanol. Geotherm. Res. 392 , 106761 (2020).
65 M. Coltelli, P. Del Carlo, L. Vezzoli, Discovery of a Plinian basaltic eruption of Roman age at Etna volcano Italy. Geology 26 , 1095 (1998).
66 L. Costantini, B. F. Houghton, C. Bonadonna, Constraints on eruption dynamics of basaltic explosive activity derived from chemical and microtextural study: The example of the Fontana Lapilli Plinian eruption, Nicaragua. J. Volcanol. Geotherm. Res. 189 , 207–224 (2010).
67 B. F. Houghton, H. M. Gonnermann, Basaltic explosive volcanism: Constraints from deposits and models. Geochemistry 68 , 117–140 (2008).
68 P. Moitra, H. M. Gonnermann, B. F. Houghton, T. Giachetti, Relating vesicle shapes in pyroclasts to eruption styles. Bull. Volcanol. 75 , 691 (2013).
69 G. La Spina, F. Arzilli, E. W. Llewellin, M. R. Burton, A. B. Clarke, M. de Michieli Vitturi, M. Polacci, M. E. Hartley, D. Di Genova, H. M. Mader, Explosivity of basaltic lava fountains is controlled by magma rheology, ascent rate and outgassing. Earth Planet. Sci. Lett. 553 , 116658 (2021).
70 L. Bai, D. R. Baker, M. Rivers, Experimental study of bubble growth in Stromboli basalt melts at 1 atm. Earth Planet. Sci. Lett. 267 , 533–547 (2008).
71 Y. Zhang, H. Ni, Diffusion of H, C, and O components in silicate melts. Rev. Mineral. Geochem. 72 , 171–225 (2010).
72 A. Lindoo, An experimental study of permeability development as a function of crystal-free melt viscosity. Earth Planet. Sci. Lett. 435 , 45–54 (2016).
73 A. C. Rust, M. Manga, K. V. Cashman, Determining flow type, shear rate and shear stress in magmas from bubble shapes and orientations. J. Volcanol. Geotherm. Res. 122 , 111–132 (2003).
74 A. Toramaru, Formation of propagation pattern in two-phase flow systems with application to volcanic eruptions. Geophys. J. Int. 95 , 613–623 (1988).
75 D. Giordano, J. K. Russell, D. B. Dingwell, Viscosity of magmatic liquids: A model. Earth Planet. Sci. Lett. 271 , 123–134 (2008).
76 M. Polacci, D. R. Baker, L. Mancini, S. Favretto, R. J. Hill, Vesiculation in magmas from Stromboli and implications for normal Strombolian activity and paroxysmal explosions in basaltic systems. J. Geophys. Res. Solid Earth 114 , 10.1029/2008JB005672 (2009).
77 L. Bai, D. R. Baker, M. Polacci, R. J. Hill, In-situ degassing study on crystal-bearing Stromboli basaltic magmas: Implications for Stromboli explosions. Geophys. Res. Lett. 38 , L17309 (2011).
78 A. Theurel, M. Collombet, A. Burgisser, C. Martel, L. Arbaret, R. Champallier, Experimental evidence of primary permeability at very low gas content in crystal-rich silicic magma. Geophys. Res. Lett. 51 , e2024GL108389 (2024).
79 E. A. Parfitt, L. Wilson, Explosive volcanic eruptions-IX. The transition between Hawaiian-style lava fountaining and Strombolian explosive activity. Geophys. J. Int. 121 , 226–232 (1995).
80 J. E. Sable, B. F. Houghton, P. Del Carlo, M. Coltelli, Changing conditions of magma ascent and fragmentation during the Etna 122 BC basaltic Plinian eruption: Evidence from clast microtextures. J. Volcanol. Geotherm. Res. 158 , 333–354 (2006).
81 S. Vergniolle, C. Jaupart, Separated two-phase flow and basaltic eruptions. J. Geophys. Res. Solid Earth 91 , 12842–12860 (1986).
82 P. Papale, A. Neri, G. Macedonio, The role of magma composition and water content in explosive eruptions. J. Volcanol. Geotherm. Res. 87 , 75–93 (1998).
83 L. Wilson, R. S. J. Sparks, G. P. L. Walker, Explosive volcanic eruptions—IV. The control of magma properties and conduit geometry on eruption column behaviour. Geophys. J. Roy. Astron. Soc. 63 , 117–148 (1980).
84 J. Crozier, S. Tramontano, P. Forte, S. J. C. Oliva, H. M. Gonnermann, E. Lev, M. Manga, M. Myers, E. Rader, P. Ruprecht, H. Tuffen, R. Paisley, B. F. Houghton, T. Shea, C. I. Schipper, J. M. Castro, Outgassing through magmatic fractures enables effusive eruption of silicic magma. J. Volcanol. Geotherm. Res. 430 , 107617 (2022).
85 H. M. Gonnermann, M. Manga, Explosive volcanism may not be an inevitable consequence of magma fragmentation. Nature 426 , 432–435 (2003).14647379
86 C. Jaupart, Gas loss from magmas through conduit walls during eruption. Geol. Soc. Lond. Spec. Publ. 145 , 73–90 (1998).
87 E. A. Parfitt, A discussion of the mechanisms of explosive basaltic eruptions. J. Volcanol. Geotherm. Res. 134 , 77–107 (2004).
88 M. Polacci, D. R. Baker, L. Bai, L. Mancini, Large vesicles record pathways of degassing at basaltic volcanoes. Bull. Volcanol. 70 , 1023–1029 (2008).
89 E. C. Bamber, G. La Spina, F. Arzilli, M. Polacci, L. Mancini, M. De Michieli Vitturi, D. Andronico, R. A. Corsaro, M. R. Burton, Outgassing behaviour during highly explosive basaltic eruptions. Commun. Earth Environ. 5 , 3 (2024).
90 Y. Moussallam, E. F. Rose-Koga, K. T. Koga, E. Médard, P. Bani, J.-L. Devidal, D. Tari, Fast ascent rate during the 2017–2018 Plinian eruption of Ambae (Aoba) volcano: A petrological investigation. Contrib. Mineral. Petrol. 174 , 90 (2019).
91 K. Goepfert, J. E. Gardner, Influence of pre-eruptive storage conditions and volatile contents on explosive Plinian style eruptions of basic magma. Bull. Volcanol. 72 , 511–521 (2010).
92 P. Moitra, B. F. Houghton, Porosity-permeability relationships in crystal-rich basalts from Plinian eruptions. Bull. Volcanol. 83 , 71 (2021).
93 J. Sable, B. Houghton, C. Wilson, R. Carey, “Eruption mechanisms during the climax of the Tarawera 1886 basaltic Plinian eruption inferred from microtextural characteristics of the deposits” in Studies in Volcanology: The Legacy of George Walker, T. Thordarson, S. Self, G. Larsen, S. K. Rowland, A. Höskuldsson, Eds. (Geological Society, 2009), pp. 129–154.
94 L. A. Szramek, Mafic Plinian eruptions: Is fast ascent required? J. Geophys. Res. Solid Earth 121 , 7119–7136 (2016).
95 M. Colombier, J. Vasseur, B. F. Houghton, F. Cáceres, B. Scheu, U. Kueppers, S. Thivet, L. Gurioli, C. Montanaro, A. Soldati, A. Di Muro, D. B. Dingwell, Degassing and gas percolation in basaltic magmas. Earth Planet. Sci. Lett. 573 , 117134 (2021).
96 D. J. Ferguson, H. M. Gonnermann, P. Ruprecht, T. Plank, E. H. Hauri, B. F. Houghton, D. A. Swanson, Magma decompression rates during explosive eruptions of Kīlauea volcano, Hawaii, recorded by melt embayments. Bull. Volcanol. 78 , 71 (2016).
97 E. Del Bello, E. W. Llewellin, J. Taddeucci, P. Scarlato, S. J. Lane, An analytical model for gas overpressure in slug-driven explosions: Insights into Strombolian volcanic eruptions. J. Geophys. Res. Solid Earth 117 , B02206 (2012).
98 E. W. Llewellin, E. Del Bello, J. Taddeucci, P. Scarlato, S. J. Lane, The thickness of the falling film of liquid around a Taylor bubble. Proc. R. Soc. Math. Phys. Eng. Sci. 468 , 1041–1064 (2012).
99 D. Carbone, L. Zuccarello, A. Messina, S. Scollo, H. Rymer, Balancing bulk gas accumulation and gas output before and during lava fountaining episodes at Mt. Etna. Sci. Rep. 5 , 18049 (2016).
100 M. Cassidy, M. Manga, K. Cashman, O. Bachmann, Controls on explosive-effusive volcanic eruption styles. Nat. Commun. 9 , 2839 (2018).30026543
101 H. M. Gonnermann, Magma Fragmentation. Annu. Rev. Earth Planet. Sci. 43 , 431–458 (2015).
102 E. W. Llewellin, H. M. Mader, S. D. R. Wilson, The constitutive equation and flow dynamics of bubbly magmas. Geophys. Res. Lett. 29 , 2170 (2002).
103 E. W. Llewellin, M. Manga, Bubble suspension rheology and implications for conduit flow. J. Volcanol. Geotherm. Res. 143 , 205–217 (2005).
104 C. Oppenheimer, R. Moretti, P. R. Kyle, A. Eschenbacher, J. B. Lowenstern, R. L. Hervig, N. W. Dunbar, Mantle to surface degassing of alkalic magmas at Erebus volcano, Antarctica. Earth Planet. Sci. Lett. 306 , 261–271 (2011).
105 P. Valdivia, A. A. Marshall, B. D. Brand, M. Manga, C. Huber, Mafic explosive volcanism at Llaima Volcano: 3D x-ray microtomography reconstruction of pyroclasts to constrain shallow conduit processes. Bull. Volcanol. 84 , 2 (2022).
106 A. Burgisser, L. Chevalier, J. E. Gardner, J. M. Castro, The percolation threshold and permeability evolution of ascending magmas. Earth Planet. Sci. Lett. 470 , 37–47 (2017).
107 M. Colombier, F. B. Wadsworth, B. Scheu, J. Vasseur, K. J. Dobson, F. Cáceres, A. Allabar, F. Marone, C. M. Schlepütz, D. B. Dingwell, In situ observation of the percolation threshold in multiphase magma analogues. Bull. Volcanol. 82 , 32 (2020).32189822
108 F. Arzilli, M. Polacci, G. La Spina, N. Le Gall, E. W. Llewellin, R. A. Brooker, R. Torres-Orozco, D. Di Genova, D. A. Neave, M. E. Hartley, H. M. Mader, D. Giordano, R. Atwood, P. D. Lee, F. Heidelbach, M. R. Burton, Dendritic crystallization in hydrous basaltic magmas controls magma mobility within the Earth’s crust. Nat. Commun. 13 , 3354 (2022).35688812
109 R. A. Corsaro, L. Miraglia, M. Pompilio, Petrologic evidence of a complex plumbing system feeding the July–August 2001 eruption of Mt. Etna, Sicily, Italy. Bull. Volcanol. 69 , 401–421 (2007).
110 N. Le Gall, F. Arzilli, G. La Spina, M. Polacci, B. Cai, M. E. Hartley, N. T. Vo, R. C. Atwood, D. Di Genova, S. Nonni, E. W. Llewellin, M. R. Burton, P. D. Lee, In situ quantification of crystallisation kinetics of plagioclase and clinopyroxene in basaltic magma: Implications for lava flow. Earth Planet. Sci. Lett. 568 , 117016 (2021).
111 M. Polacci, F. Arzilli, G. La Spina, N. Le Gall, B. Cai, M. E. Hartley, D. Di Genova, N. T. Vo, S. Nonni, R. C. Atwood, E. W. Llewellin, P. D. Lee, M. R. Burton, Crystallisation in basaltic magmas revealed via in situ 4D synchrotron X-ray microtomography. Sci. Rep. 8 , 8377 (2018).29849174
112 F. Menges, Spectragryph - optical spectroscopy software, version 1.2.16 (2022); www.effemm2.de/spectragryph/.
113 D. Testemale, R. Argoud, O. Geaymond, J.-L. Hazemann, High pressure/high temperature cell for x-ray absorption and scattering techniques. Rev. Sci. Instrum. 76 , 043905 (2005).
114 R. Bruyère, A. Prat, C. Goujon, J.-L. Hazemann, A new pressure regulation device using high pressure isolation valves. J. Phys. Conf. Ser. 121 , 122003 (2008).
115 M. D. Abramoff, P. J. Magalhães, S. J. Ram, Image processing with ImageJ. Biophotonics Int. 11 , 36–42 (2004).
116 K. V. Cashman, M. T. Mangan, Physical aspects of magmatic degassing; II, Constraints on vesiculation processes from textural studies of eruptive products. Rev. Mineral. Geochem. 30 , 447–478 (1994).
117 O. Sigmarsson, D. Laporte, M. Carpentier, B. Devouard, J.-L. Devidal, J. Marti, Formation of U-depleted rhyolite from a basanite at El Hierro, Canary Islands. Contrib. Mineral. Petrol. 165 , 601–622 (2013).
118 M. S. Ghiorso, G. A. R. Gualda, An H2O–CO2 mixed fluid saturation model compatible with rhyolite-MELTS. Contrib. Mineral. Petrol. 169 , 53 (2015).
119 G. A. R. Gualda, M. S. Ghiorso, R. V. Lemons, T. L. Carley, Rhyolite-MELTS: A modified calibration of MELTS optimized for silica-rich, fluid-bearing magmatic systems. J. Petrol. 53 , 875–890 (2012).
120 A. Vona, C. Romano, D. B. Dingwell, D. Giordano, The rheology of crystal-bearing basaltic magmas from Stromboli and Etna. Geochim. Cosmochim. Acta 75 , 3214–3236 (2011).
121 A. Burgisser, J. E. Gardner, Experimental constraints on degassing and permeability in volcanic conduit flow. Bull. Volcanol. 67 , 42–56 (2004).
