# Bubble connectivity in experimentally-sheared crystal-bearing silicic melts

Camille Daffos, Caroline Martel, Laurent Arbaret, Rémi Champallier

► To cite this version:

Camille Daffos, Caroline Martel, Laurent Arbaret, Rémi Champallier. Bubble connectivity in experimentally-sheared crystal-bearing silicic melts. Comptes Rendus. Géoscience, 2024, 356 (S1), pp.1-18. 10.5802/crgeos.214. insu-04267660

HAL Id: insu-04267660

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

Submitted on 2 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

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

Distributed under a Creative Commons Attribution 4.0 International License
# Comptes Rendus
## Géoscience
### Sciences de la Planète

Camille Daffos, Caroline Martel, Laurent Arbaret and Rémi Champallier

**Bubble connectivity in experimentally-sheared crystal-bearing silicic melts**

Published online: 12 June 2023

https://doi.org/10.5802/crgeos.214

**Part of Special Issue:** Magma degassing and its impact on the Earth's atmosphere: from magma oceans to lava lakes

**Guest editors:** Manuel Moreira (Institut des Sciences de la Terre d'Orléans Université d'Orléans-CNRS-BRGM 1a rue de la Férollerie 45071 Orléans France), Bruno Scaillet (Institut des Sciences de la Terre d'Orléans Université d'Orléans-CNRS-BRGM 1a rue de la Férollerie 45071 Orléans France) and Clive Oppenheimer (Department of Geography, University of Cambridge, Downing Place, Cambridge CB2 3EN, UK)

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

This article is licensed under the
CREATIVE COMMONS ATTRIBUTION 4.0 INTERNATIONAL LICENSE.
http://creativecommons.org/licenses/by/4.0/

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

Les Comptes Rendus. Géoscience — Sciences de la Planète sont membres du
Centre Mersenne pour l'édition scientifique ouverte
www.centre-mersenne.org
e-ISSN : 1778-7025
et al., 2021, Popa et al., 2021] and emplacement conditions [e.g. flow and conduit properties, rates of decompression and shear; Cassidy et al., 2018].

An ascending volatile-saturated magma exsolves gases that are accommodated through bubble nucleation or growth of pre-existing bubbles in case of an initial volatile phase in excess. In highly viscous magmas (viscosity >10⁶ Pa·s), bubbles are usually trapped in the melt and cannot easily migrate independently of the liquid to outgas via separated flow. If physical and time conditions permit, bubbles connect by coalescing and can develop gas channels that ultimately make the magma permeable to gases, thus allowing outgassing through fractured dome rocks or conduit walls. Rates of gas escape in silica-rich melts are difficult to assess because they depend on intrinsic parameters, such as the bubble characteristics [e.g., vesicularity, Klug and Cashman, 1996; pore aperture size, Wright et al., 2009, Degruyter et al., 2010; bubble number density; Bain et al., 2019; bubble network, Burgisser et al., 2017], and external parameters related to magma emplacement kinetics. Among these external parameters, decompression rate has been demonstrated to affect bubble connectivity and permeability. The studies dealing with decompression-induced bubbles developing isotropically highlighted permeability increase at vesicularities mostly >50 vol% [Burgisser and Gardner, 2005, Takeuchi et al., 2005, Okumura et al., 2012, Martel and Iacono-Marziano, 2015]. Thereafter, shear has been demonstrated to drastically affect bubble connectivity and permeability, by promoting bubble coalescence [Okumura et al., 2008] and outgassing [Okumura et al., 2009, Shields et al., 2016], but also by alternately favouring and preventing connectivity where deformation resulted in magma compaction [Gonnermann et al., 2017]. If magma viscosity and/or shear rate are high enough, permeability development may proceed through magma fracturing [Stasiuk et al., 1996, Tuffen et al., 2003, Castro et al., 2012, Gaunt et al., 2014, Kushnir et al., 2017].

Many studies have been dedicated to crystal-free or -poor melts, but magmas can be highly crystallized, which affects the development of bubble connectivity. Estimates of the percolation threshold commonly range from 30 vol% [Saar and Manga, 1999, Blower, 2001] to 80 vol% vesicularity [Westrich and Eichelberger, 1994, Takeuchi et al., 2008, Wright et al., 2009], with the highest values characterizing

rapidly-decompressed crystal-free magmas. Okumura et al. [2012] experimentally decompressed rhyolitic melts containing 30 and 50 vol% corundum crystals, and concluded that gas permeability remained low until vesicularity reached >68 vol%. Lindoo et al. [2017] and deGraffenried et al. [2019] decompressed mafic and rhyolitic magmas, respectively, and reported a decrease of the percolation threshold from >60 to <50 vol% in the presence of more than about 20 vol% crystals. Parmigiani et al. [2017] numerically modelled the outgassing efficiency of a magmatic volatile phase in crystal-rich (mush-like) magmas, highlighting gas permeability from very low porosities (~10 vol%) via a mechanism of channelling by which crystals build sustainable channels for gas percolation. Such a permeability development and outgassing at low vesicularity and high crystal content could explain the low porosity values inferred to occur in the volcanic conduit prior to Vulcanian eruptions [Collombet et al., 2021].

In crystal-bearing silicic magmas, the bubble network mostly rearranges under decompression in the central part of the volcanic conduit and simple-shear deformation at conduit edges or in lava domes where the rearrangement of the stressed crystal framework is facilitated by strain localization. To understand the processes governing the development of degassing pathways in crystal-rich silicic magmas, few studies have been dedicated to deform three-phase (melt, bubbles, crystals) magmas in simple shear using high-temperature and high-pressure deformation rigs of Paterson type. Laumonier et al. [2011] deformed a haplotonalitic melt containing 50 vol% plagioclase crystals and 11 vol% porosity, at temperature (T) of 600 °C, confining pressure (P) of 200 MPa, bulk finite shear strains (γ) of 1.3 and shear rates (γr) from 3 × 10⁻⁵ to 1 × 10⁻³ s⁻¹. They highlighted gas accumulation in local microstructures caused by shear-induced crystal fabric (local γ up to ~9). Shields et al. [2014] sheared haplogranitic melts spanning crystal contents from 0 to 42 vol% and CO₂-bubble contents from 12 to 36 vol%, at T < 600 °C, P of 150–200 MPa, γ up to 10, and γr from 1 × 10⁻⁴ to 5 × 10⁻⁴ s⁻¹. They did not observe bubble coalescence and outgassing took place via sample-wide fracturing. These results of sample permeability reached by fracturing agreed with those obtained by Kushnir et al. [2017] using crystal-free haplogranitic melts with ~15 vol% argon bubbles deformed at magmatic T
of 880 °C, *P* of 60 MPa, γ up to 7, and γ<sub>r</sub> of 1 × 10<sup>-4</sup> to 8 × 10<sup>-4</sup> s<sup>-1</sup>. Pistone et al. [2012] sheared haplogranitic melts with 24 to 65 vol% quartz crystals and 9–12 vol% CO<sub>2</sub> bubbles, at *T* of 450–750 °C, *P* of 200 MPa, γ up to 8, and γ<sub>r</sub> from 5 × 10<sup>-6</sup> to 4 × 10<sup>-3</sup> s<sup>-1</sup>. In the crystal-poor (24–44 vol%) samples, the authors reported bubble stretching with increasing strain, facilitating coalescence (for γ > 2 and for increasing γ<sub>r</sub>) and gas channel formation, eventually leading to outgassing along the walls of the sample container. In the crystal-rich (55–65 vol%) samples, however, bubble stretching was localized in shear bands where the crystal framework (including crystal breakage) prevented gas loss.

These experimental studies provided information on the deformation of three-phase magmas at *P* of ~200 MPa, *T* mostly close to the glass transition, for bubble contents mostly <15 vol%, and for γ up to ~10. We aim at extending these studies to the investigation of the effect of a crystal network on bubble connectivity and gas percolation in bubbly crystal-bearing magmas emplacing at shallow depth, such as in upper volcanic conduits or lava domes, from which explosive eruptions may be triggered. We performed simple shear experiments in a Paterson rig, using three-phase magmas consisting of 0 to 50 bulk vol% plagioclase crystals and 15–30 bulk vol% H<sub>2</sub>O bubbles (20–45 vol% recalculated on the melt phase) in a haplogranitic melt, under magmatic *T* of 650–750 °C, confining *P* of 50 MPa, γ up to ~2, and γ<sub>r</sub> of 2 × 10<sup>-4</sup> s<sup>-1</sup>. We characterized the sample microstructures and bubble connectivity as a function of crystal content and shear strain, and have discussed implications for the natural magmas and eruptive dynamics.

**2. Methods**

**2.1. *Experimental methods***

The experiments consisted of three main phases: (i) synthesis of crystal-bearing hydrated glasses, (ii) decompression-induced H₂O-bubble formation, and (iii) deformation in simple shear of the three-phase magmas.

2.1.1. *Synthesis of crystal-bearing hydrated glasses*

The method of preparation is detailed in Supple-
mentary Section 1.1.1 and a summary is given below.

The anhydrous starting glass is a haplogranite (HPG8; composition in wt%: 78.6 SiO₂, 12.5 Al₂O₃, 4.6 Na₂O, 4.2 K₂O). The HPG8 glass was chosen because (i) it is rheologically well-characterized [Hess and Dingwell, 1996], (ii) its eutectic composition allows experiments at relatively low *T* without crystallization, and (iii) it simulates late-stage crystallization of evolved liquids such as rhyolitic melts. The HPG8 hydrations were performed in internally-heated pressure vessels (at the Institut des Sciences de la Terre d'Orléans, ISTO) in order to obtain glasses with 5.0 and 9.6 wt% H₂O, to which different amounts of plagioclase crystals (size of 50–90 µm) were manually added. Plagioclase crystals were chosen because they represent the main crystalline phase of silicic magmas, as phenocrysts or microcrysts. The crystal contents recalculated on the melt are 0, 21, 50, and 70 vol%, leading to bulk crystal volume (Φ<sub>cBulk_3D</sub>, recalculated on a bubble-bearing basis) of about 0, 10, 30, and 50 vol%.

2.1.2. *Synthesis of a bubble-bearing magmas*

The powder mixtures of plagioclase crystals and hydrated HPG8 glass were decompressed from 300 to 50 MPa at 850 °C in the Paterson gas-medium apparatus [Paterson and Olgaard, 2000; Australian Scientific Instruments Pty Ltd, at ISTO; Figure 1a], in order to trigger a bubble nucleation event of supercritical H₂O fluid (hereafter referred as to gas bubbles). The strategy was to obtain bulk amounts of bubbles around 20–40 vol%, ideally not connected to each other, in order to check whether further deformation will promote gas connectivity. The experimental details are given in Supplementary Section 1.1.2.

2.1.3. *Torsion experiments*

Keeping the experimental setup as it was after
decompression, deformation in right-lateral simple
shear was carried out at *P* of 50 MPa, *T* of 650 or
750 °C, γ<sub>r</sub> of 2 × 10<sup>-4</sup> s<sup>-1</sup>, and γ from 0 to 2. Such γ are
rather small compared to those observed at natural
conduit margins or dome bases, but the aim was to
determine the gas percolation thresholds occurring
during the initiation of deformation. The experimen-
tal details and calculations confirming limited relax-
ation of the bubbles during quenching are given in
Supplementary Section 1.1.3.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Experimental device and strategy. (a) Scheme of the Paterson gas-medium press, (b) Photography of the sample assembly, and (c) Pressure-temperature-time path showing the isothermal step (850 °C) of decompression-induced (300 to 50 MPa) bubble nucleation, followed by an isobaric (50 MPa) and isothermal (750 or 650 °C) step of sample torsion (γ < 2) at a rate of 2 × 10⁻⁴ s⁻¹.</figcaption></figure>

## 2.2. Analytical methods

After experiment, the sample was cut following a section parallel to the shear plane and exposing the longitudinal tangential surfaces (maximum shear strain) for microstructural analyses using a scanning electron microscope (SEM; Merlin Compact

Zeiss at ISTO). The whole section was imaged at high resolution using the SEM, and the images were segmented using the SPO2003 image-processing software [Launeau and Robin, 1996, Launeau and Cruden, 1998], in order to determine the two-dimensional (2D) contents, sizes, and ori-
entations of the bubbles and the crystals, as detailed in Supplementary Section 1.2.1. A core sample bored perpendicularly to the shear plane was used for three-dimensional (3D) analyses using X-ray computed tomography (XCT; Phoenix NanoTOM at ISTO). The 3D bulk porosity ($\Phi_{\text{bulk 3D}}$) and the amount of bubbles connected to the sample borders ($\Phi_{\text{connect 3D}}$) were determined by segmenting the SCT images using commercial software, such as VGStudio Max and Blob3D, as detailed in Supplementary Section 1.2.2. The plagioclase crystals were hardly discernible due to their low contrast of density with the glass and were not subjected to 3D analyses.

The glass H₂O contents of the bubble-bearing samples were checked following the “by-difference” method using an electron microprobe (EMP; Cameca SX Five at ISTO). As some of the residual glasses were unexpectedly rich in H₂O, they were double-checked by micro-Raman spectroscopy (ISTO), as detailed in Supplementary Section 1.2.3.

## 2.3. Calculations

The calculations mainly concerned the gas fraction and the bulk viscosity of the samples. The expected gas fraction was determined using the equation of Jaupart and Allègre [1991], which calculates the gas volume ($\alpha_{\text{melt}}$) resulting from a closed-system degassing (in that gas bubbles remain in contact with the melt), as described in Supplementary Section 1.3.1.

The samples were not viscous enough at 750 °C to allow viscosity measurements in the Paterson press, so we estimated the bulk viscosity ($\eta_{\text{bulk}}$) from Hess and Dingwell [1996] for the hydrated melt and using the equation of Truby et al. [2015] for the three-phase suspensions, as detailed in Supplementary Section 1.3.2.

# 3. Results

Four series with different bulk crystal contents ($\Phi_{\text{bulk 3D}}$) of 0, 10, 30, 50 vol%, hereafter referred as to serie0, serie10, serie30, and serie50, respectively, were decompressed in the Paterson press for bubble formation and subsequent deformation in simple shear ($P = 50$ MPa, $T = 650–750$ °C, $\gamma_r = 2 \times 10^{-4}$ s⁻¹) to finite deformations $\gamma = 0$ (no deformation), $\gamma = 0.6$, and $\gamma = 1.3–1.8$ (no run for serie10). The experimental conditions of the deformed

bubble- and crystal-bearing magmas are given in Table 1. The analytical results are reported in Table 2 and detailed below.

## 3.1. Crystal-free series (serie0)

The SEM pictures of the samples from serie0 show relatively homogeneous spatial distributions of the bubbles (Figure 2a–c), in agreement with the absence of strain localization expected from the Newtonian behaviour of pure melts. Melt porosity, $\Phi_{\text{melt 2D}}$, is 31±3, to 24±2 and 30±4 area%, at $\gamma = 0, 0.6$, and 1.3, respectively (Figure 3a), in relatively good agreement with the porosities measured from the SCT images, $\Phi_{\text{bulk 3D}}$, which are 30±3, 21±3, and 15±3 vol%, respectively. The content of bubbles connected to the sample outside measured from the SCT images, $\Phi_{\text{connect 3D}}$, represents about half of the $\Phi_{\text{bulk 3D}}$ (i.e., 49±3, 52±3, and 50±3 vol%, at $\gamma = 0, 0.6$, and 1.3, respectively). The 2D bubble number densities measured from the SEM images, BND_2D, are of $10^{5.6–5.7} \pm 0.1$ m⁻² (Figure 4a) and those determined using the SCT, BND_3D, are $10^{10.5–10.6}$ m⁻³. In the undeformed sample, 95% of the bubbles have equivalent diameters of 20–40 µm, representing about half of the total vesicularity. The other half volume is made of bubbles with up to 300 µm diameter (Supplementary Figure 3a). Bubble eccentricity, $R_b$, increases from 1.3±0.1 to 1.9±0.2 and 2.1±0.3 with increasing $\gamma$, with a slight preferential orientation of the bubbles, $\theta_b$ of ~22±6° (Figure 4b). Glass H₂O contents measured following the EMP by-difference method are 3.6±0.7 wt% at $\gamma = 0$ and 4.1±0.7 wt% at $\gamma = 1.3$, giving bulk viscosities of $10^{7.3}$ and $10^{7.0}$ Pa·s, respectively.

## 3.2. Series with 10% crystals (serie10)

Serie10 consists of only a $\gamma = 0$ (undeformed) sample and a $\gamma = 0.6$ sample, since higher shear strains inevitably conducted to capsule tearing. The deformed sample shows incipient shear zones marked by discrete local reorganisation of the bubbles and crystals (Figure 2e). $\Phi_{\text{melt 2D}}$ slightly decreases from 31±3 area% at $\gamma = 0$ to 27±3 area% at $\gamma = 0.6$ (Figure 3a). Recalculating bulk porosities (i.e. including crystals) give $\Phi_{\text{bulk 2D}}$ of 29±3 area% at $\gamma = 0$ and 24±3 area% at $\gamma = 0.6$, in good agreement with the SCT $\Phi_{\text{bulk 3D}}$ of 28±3 vol% at $\gamma = 0$ and 18±3 vol% at $\gamma =$
<table><caption>**Table 1.** Experimental conditions</caption><thead><tr><th rowspan="3">Run #</th><th rowspan="3">Series #</th><th rowspan="3">Φ<sub>c_melt_3D</sub><sup>a</sup><br>(vol%)</th><th rowspan="3">H<sub>2</sub>O<sub>i</sub><sup>b</sup><br>(wt%)</th><th colspan="5">Pre-deformation bubble forming</th><th rowspan="3">Shear at P-T<br>(γ<sub>r</sub> = 2 × 10<sup>-4</sup> s<sup>-1</sup>)</th></tr><tr><th colspan="2">Decompression (from<br>300 MPa and 850 °C)<sup>c</sup></th><th colspan="3">Cooling at P</th></tr><tr><th>P<br>(MPa)</th><th>τ<sub>DP</sub><sup>d</sup><br>(min)</th><th>Dwell1<br>(min)</th><th>T<br>(°C)</th><th>Dwell2<br>(min)</th><th>τ<sub>def</sub><sup>d</sup><br>(min)</th><th>γ</th></tr></thead><tbody><tr><td>Pp17</td><td>serie0</td><td>0</td><td>5.0</td><td>50</td><td>19</td><td>45</td><td>650</td><td>15</td><td>-</td><td>0</td></tr><tr><td>Pp16</td><td>serie0</td><td>0</td><td>5.0</td><td>50</td><td>18</td><td>45</td><td>650</td><td>15</td><td>60</td><td>0.6</td></tr><tr><td>Pp15</td><td>serie0</td><td>0</td><td>5.0</td><td>50</td><td>18</td><td>45</td><td>650</td><td>15</td><td>120</td><td>1.3</td></tr><tr><td>Pp5</td><td>serie10</td><td>21</td><td>9.6</td><td>50</td><td>16</td><td>45</td><td>750</td><td>240</td><td>-</td><td>0</td></tr><tr><td>Pp7</td><td>serie10</td><td>21</td><td>9.6</td><td>50</td><td>18</td><td>45</td><td>750</td><td>15</td><td>60</td><td>0.6</td></tr><tr><td>Pp6</td><td>serie30</td><td>50</td><td>9.6</td><td>50</td><td>9</td><td>45</td><td>750</td><td>15</td><td>-</td><td>0</td></tr><tr><td>Pp8</td><td>serie30</td><td>50</td><td>9.6</td><td>50</td><td>19</td><td>45</td><td>750</td><td>15</td><td>60</td><td>0.6</td></tr><tr><td>Pp9</td><td>serie30</td><td>50</td><td>9.6</td><td>50</td><td>16</td><td>45</td><td>750</td><td>15</td><td>164</td><td>1.8</td></tr><tr><td>Pp12</td><td>serie50</td><td>70</td><td>9.6</td><td>50</td><td>18</td><td>45</td><td>750</td><td>15</td><td>-</td><td>0</td></tr><tr><td>Pp11</td><td>serie50</td><td>70</td><td>9.6</td><td>50</td><td>20</td><td>45</td><td>750</td><td>15</td><td>62</td><td>0.6</td></tr><tr><td>Pp10</td><td>serie50</td><td>70</td><td>9.6</td><td>50</td><td>17</td><td>45</td><td>750</td><td>15</td><td>145</td><td>1.5</td></tr></tbody></table>

<sup>a</sup>Φ<sub>c_melt_3D</sub> is the content of plagioclase crystals manually added to the hydrated haplogranitic (HPG8) glass powder.

<sup>b</sup>H<sub>2</sub>O<sub>i</sub> is the H<sub>2</sub>O content of the hydrated glasses that were used for the Paterson experiments, calculated after Newman and Lowenstern [2002] at 950 °C and 170 MPa (for the crystal-free series) or 400 MPa (for the crystal-bearing series).

<sup>c</sup>Decompression follows a 15-min dwell at 300 MPa and 850 °C.

<sup>d</sup>τ<sub>DP</sub><sup>c</sup> and τ<sub>def</sub><sup>c</sup> are the duration of the decompressions and deformations, respectively.

0.6. BND_2D is 10<sup>5.6±0.2</sup> m<sup>-2</sup> at γ = 0 and 10<sup>5.2±0.2</sup> m<sup>-2</sup> at γ = 0.6 (Figure 4a), and BND_3D is 10<sup>9.2-9.3</sup> m<sup>-3</sup>. In the undeformed sample, about 80% of the bubbles have diameters of 20–40 µm, but their volumes only represent less than 10% of total vesicularity. More than one third of the vesicularity is occupied by bubbles with equivalent diameters of 220–240 µm (Supplementary Figure 3b). The bubbles have R<sub>b</sub> of 1.0 ± 0.1 at γ = 0 and 1.3 ± 0.1 at γ = 0.6, showing no preferential orientation (θ<sub>b</sub> of 20 ± 43°) in the unsheared sample (Figure 4b). The glass H<sub>2</sub>O contents are of 6.1 ± 0.7 wt% at γ = 0 and 5.3 ± 0.7 wt% at γ = 0.6. The crystal contents, Φ<sub>c bulk 2D</sub>, are of 9 ± 1 area% at γ = 0 and 12 ± 2 area% at γ = 0.6. The crystal number densities, CND, are 10<sup>4.8±0.2</sup> m<sup>-2</sup> at γ = 0 and 10<sup>5.1±0.3</sup> m<sup>-2</sup> at γ = 0.6 (Figure 4a). The crystals have eccentricities, R<sub>c</sub>, that are not varying from 1.1 ± 0.1 with γ, and their main orientations, θ<sub>c</sub>, is 15 ± 33° at γ = 0 and 18 ± 3° at γ = 0.6 (Figure 4b).

### 3.3. Series with 30% crystals (serie30)

Serie30 shows shear zones of reorganised bubbles and crystals, evidencing an intensification of strain localisation from γ = 0.6 to 1.8 (Figure 2g,h). Φ<sub>b melt 2D</sub> are 20 ± 3, 26 ± 3, 26 ± 12 area% at γ = 0, 0.6, and 1.8, respectively (Figure 3a). Recalculating bulk porosities gives Φ<sub>b bulk 2D</sub> of 15 ± 3, 20 ± 3, 23 ± 12 area% at γ = 0, 0.6, and 1.8, respectively. The XCT Φ<sub>b bulk 3D</sub> gives 22 ± 3, 25 ± 3 and 15 ± 3 vol% at γ = 0, 0.6, and 1.8, respectively. BND_2D is 10<sup>5.6-5.8±0.1</sup> m<sup>-2</sup> (Figure 4a) and BND_3D is 10<sup>9.3</sup> m<sup>-3</sup> at γ = 0 and 0.6, and 10<sup>9.9</sup> m<sup>-3</sup> at γ = 1.8. The undeformed sample shows more than 90% of bubbles with 20-µm diameters, the volume of which represents one third of total vesicularity. The rest of the vesicularity is mostly occupied by bubbles with equivalent diameters from 40 to 140 µm, with some larger bubbles up to 240 µm diameter (Supplementary Figure 3f). The bubbles have R<sub>b</sub> of 1.0 ± 0.1 at any γ, showing
<figure><img src="image_5.png" /><figcaption>**Figure 2.** SEM pictures of the samples cut following the sketch shown in Supplementary Section 1.2.1. Increasing shear strain from $\gamma = 0$ to 2 (from left to right) and crystal content from 0 to 50 vol% (from top to down) promotes the development of non-spherical bubbles and heterogeneities in bubble spatial distribution. All main pictures have a yellow scale bar of 4 mm and subpanels are enlargements (with different coloured scale bars). The red arrows in (i) and (k) point to crystal fragmentations.</figcaption></figure>

no preferential orientations ($\theta_b$ of 33–38°; Figure 4b). Glass H₂O content was only analysed at $\gamma = 1.8$, giving 3.0 ± 0.7 wt%. $\Phi_{c\_bulk\_2D}$ is 28 ± 1 area% at $\gamma = 0$, 29 ± 4 area% at $\gamma = 0.6$, and 33 ± 4 area% at $\gamma = 1.8$. The CND is $10^{5.5-5.6\pm0.2}$ m⁻² at any $\gamma$ (Figure 4a). $R_c$ are not varying from 1.1 ± 0.1 with $\gamma$. $\theta_c$ are 44 ± 31°, 58 ± 38°, and 26 ± 13° in the samples at $\gamma = 0, 0.6$, and 1.8, respectively (Figure 4b).

### 3.4. Series with 50% crystals (serie50)

Serie50 shows strong features of strain localization in the deformed samples, with alternation of strongly bubble-depleted zones and channelized

bubble zones (Figure 2j,k). $\Phi_{b\_melt\_2D}$ increases from 21 ± 2, 23 ± 10, 44 ± 17 area% with $\gamma$ increasing from 0, to 0.6 and 1.5 (Figure 3a). Bulk $\Phi_{b\_bulk\_2D}$ gives 14 ± 2, 16 ± 10, and 23 ± 17 area% at $\gamma = 0, 0.6$, and 1.5, respectively, in good agreement with the ACT $\Phi_{b\_bulk\_3D}$ of 19 ± 3, 16 ± 3 and 25 ± 3 vol%. $BND\_2D$ is $10^{6.1\pm0.1}$ m⁻² at $\gamma = 0$ and 0.6, and $10^{5.9\pm0.3}$ m⁻² at $\gamma = 1.5$ (Figure 4a). $BND\_3D$ is $10^{10.0-10.2}$ m⁻³ at $\gamma = 0$ and 0.6, and $10^{9.5}$ m⁻³ at $\gamma = 1.5$. In the undeformed sample, 95% of the bubbles have diameters of 20 µm, representing about 25% of the total vesicularity. Another third of the volume is occupied by bubbles with equivalent diameters of 240–260 µm (Supplementary Figure 3i). The bubbles have $R_b$ of
<figure><img src="image_6.png" /></figure>

**Figure 3.** Bubble formation. (a) Measured porosity versus porosity calculated for a H₂O degassing from the initial content (5.0 wt% for the crystal-free samples and 9.6 wt% for the others) to the final content measured in the residual glasses, after the equation of Jaupart and Allègre [1991]; the two arrows on the X-axis give the gas content calculated for an equilibrium H₂O exsolution (i.e., final content is H₂O solubility at 50 MPa); the vertical error bars are not statistical uncertainties but span the whole range of the measured values; (b) Difference between the H₂O contents measured in the residual glasses (circles and stars using electron microprobe and Raman spectrometry, respectively; technical details in Supplementary Section 1) and H₂O solubility at final pressure (dashed box), showing incomplete H₂O exsolution in some samples.

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

**Figure 4.** Bubble and crystal characteristics. (a) Number density and (b) Preferential orientations (0° is parallel to shear direction), with the dashed red line giving the theoretical angles between particle major axis and shear direction calculated for solid particles with eccentricity of 1.1, using the formulation of Fernandez et al. [1983]; the X-axis values have been slightly shifted for clarity; the vertical bars are not statistical uncertainties but span the whole range of the measured values.

1.0 ± 0.1 at any $\gamma$, not showing significant preferential orientations ($\theta_b$ from -13 to +33°, with standard deviations >35°; Figure 4b). The glass H₂O contents are 4.9 ± 0.7 wt% at $\gamma$ = 0 and 2.7 ± 0.7 wt% at $\gamma$ = 1.5. $\Phi_{c\_bulk\_2D}$ is of 39–48 area% whatever $\gamma$. The CND increases from $10^{5.5} \pm 0.3$ m⁻² at $\gamma$ = 0 to $10^{5.8} \pm 0.1$ m⁻² at $\gamma$ = 0.6 and $10^{6.1} \pm 0.2$ m⁻² at $\gamma$ = 1.5 (Figure 4a). $R_c$ are not varying from 1.1 ± 0.1 with $\gamma$. $\theta_c$ are 40 ± 27°, 27 ± 8°, and 28 ± 2° in the samples at $\gamma$ = 0, 0.6, and 1.5, respectively (Figure 4b).
<table><caption>Table 2. Analytical results</caption><thead><tr><th rowspan="2">Run #</th><th rowspan="2">Series #</th><th rowspan="2">γ</th><th colspan="2">Glass</th><th colspan="2">Crystals</th><th colspan="2">Bubbles</th><th colspan="2">Viscosity<sup>c</sup></th></tr><tr><th>Φ<sub>g,bulk_2D</sub> (calc) (vol% ± 3)</th><th>H<sub>2</sub>O<sub>calc</sub> (wt%)</th><th>Φ<sub>c,bulk_2D</sub> (area%) (wt%)</th><th>Log CND<sub>2D</sub> (m<sup>-2</sup>)</th><th>R<sub>c</sub><sup>b</sup></th><th>θ<sub>b</sub><sup>b</sup> (°)</th><th>Φ<sub>b,bulk_3D</sub> (vol% ± 3)</th><th>Φ<sub>b,connect_3D</sub> (%) ± 3</th><th>Φ<sub>b,bulk_2D</sub> (area%)</th><th>Φ<sub>b,metl_2D</sub> (calc) (area%)</th><th>Log BND<sub>3D</sub> (m<sup>-2</sup>)</th><th>R<sub>b</sub><sup>b</sup></th><th>θ<sub>b</sub><sup>b</sup> (°)</th><th>α<sub>metl</sub><sup>d</sup> (vol%) with H<sub>2</sub>O<sub>meas</sub> (wt%)</th><th>Log η<sub>metl</sub> (Pa·s)</th><th>Log η<sub>bulk</sub> (Pa·s)</th></tr></thead><tbody><tr><td>Pp17</td><td>ser10</td><td>0</td><td>69</td><td>3.2</td><td>3.6</td><td>0</td><td>n.a.</td><td>n.a.</td><td>n.a.</td><td>30(28)</td><td>49</td><td>31 ± 3</td><td>31 ± 3</td><td>5.6 ± 0.1</td><td>10.6</td><td>1.3 ± 0.1</td><td>-10 ± 21</td><td>22</td><td>7.1</td><td>7.3</td></tr><tr><td>Pp16</td><td>ser10</td><td>0.6</td><td>76</td><td>3.2</td><td>n.a.</td><td>0</td><td>n.a.</td><td>n.a.</td><td>n.a.</td><td>21(20)</td><td>52</td><td>24 ± 2</td><td>24 ± 2</td><td>5.7 ± 0.1</td><td>10.6</td><td>1.9 ± 0.2</td><td>+22 ± 6</td><td>27</td><td>n.a.</td><td>n.a.</td></tr><tr><td>Pp15</td><td>ser10</td><td>1.3</td><td>79</td><td>3.2</td><td>4.1</td><td>0</td><td>n.a.</td><td>n.a.</td><td>n.a.</td><td>15(21)</td><td>50</td><td>30 ± 4</td><td>30 ± 4</td><td>5.7 ± 0.1</td><td>10.5</td><td>2.0 ± 0.3</td><td>+23 ± 6</td><td>27</td><td>15</td><td>6.8</td><td>7.0</td></tr><tr><td>Pp5</td><td>ser10</td><td>0</td><td>63</td><td>2.8</td><td>6.1(5.1)</td><td>9 ± 1</td><td>4.8 ± 0.2</td><td>1.1 ± 0.1</td><td>+15 ± 33</td><td>28(29)</td><td>59</td><td>29 ± 3</td><td>31 ± 3</td><td>5.6 ± 0.2</td><td>9.3</td><td>1.0 ± 0.1</td><td>+20 ± 43</td><td>60</td><td>41</td><td>4.9</td><td>5.2</td></tr><tr><td>Pp7</td><td>ser10</td><td>0.6</td><td>65</td><td>2.8</td><td>5.3(5.9)</td><td>12 ± 2</td><td>5.1 ± 0.3</td><td>1.1 ± 0.1</td><td>+18 ± 3</td><td>18(23)</td><td>64</td><td>24 ± 3</td><td>27 ± 3</td><td>5.2 ± 0.2</td><td>9.2</td><td>1.3 ± 0.1</td><td>+35 ± 4</td><td>60</td><td>46</td><td>5.2</td><td>5.4</td></tr><tr><td>Pp6</td><td>ser10</td><td>0</td><td>63</td><td>2.8</td><td>n.a.(4.3)</td><td>28 ± 1</td><td>5.5 ± 0.2</td><td>1.1 ± 0.1</td><td>+44 ± 31</td><td>22(22)</td><td>37</td><td>15 ± 3</td><td>20 ± 3</td><td>5.8 ± 0.1</td><td>9.3</td><td>1.0 ± 0.1</td><td>+33 ± 17</td><td>60</td><td>n.a.</td><td>n.a.</td><td>n.a.</td></tr><tr><td>Pp8</td><td>ser10</td><td>0.6</td><td>57</td><td>2.8</td><td>n.a.</td><td>29 ± 4</td><td>5.5 ± 0.1</td><td>1.1 ± 0.1</td><td>+58 ± 38</td><td>25(27)</td><td>83</td><td>20 ± 3</td><td>26 ± 3</td><td>5.6 ± 0.1</td><td>9.3</td><td>1.0 ± 0.1</td><td>+38 ± 22</td><td>60</td><td>n.a.</td><td>n.a.</td><td>n.a.</td></tr><tr><td>Pp9</td><td>ser10</td><td>1.8</td><td>62</td><td>2.8</td><td>3.0</td><td>33 ± 4</td><td>5.6 ± 0.2</td><td>1.1 ± 0.1</td><td>+26 ± 13</td><td>15(10)</td><td>80</td><td>23 ± 12</td><td>26 ± 12</td><td>5.7 ± 0.1</td><td>9.9</td><td>1.0 ± 0.1</td><td>+38 ± 22</td><td>60</td><td>57</td><td>6.2</td><td>6.9</td></tr><tr><td>Pp12</td><td>ser10</td><td>0</td><td>53</td><td>2.8</td><td>4.9</td><td>48 ± 3</td><td>5.5 ± 0.3</td><td>1.1 ± 0.1</td><td>+40 ± 27</td><td>19(19)</td><td>87</td><td>14 ± 2</td><td>21 ± 2</td><td>6.1 ± 0.1</td><td>10.2</td><td>1.0 ± 0.1</td><td>-13 ± 42</td><td>60</td><td>49</td><td>5.3</td><td>6.5</td></tr><tr><td>Pp11</td><td>ser10</td><td>0.6</td><td>54</td><td>2.8</td><td>n.a.</td><td>43 ± 4</td><td>5.8 ± 0.1</td><td>1.1 ± 0.1</td><td>+27 ± 8</td><td>16(15)</td><td>76</td><td>16 ± 10</td><td>23 ± 10</td><td>6.1 ± 0.1</td><td>10.0</td><td>1.0 ± 0.1</td><td>+33 ± 36</td><td>60</td><td>n.a.</td><td>n.a.</td><td>n.a.</td></tr><tr><td>Pp10</td><td>ser10</td><td>1.5</td><td>49</td><td>2.8</td><td>2.7</td><td>39 ± 7</td><td>6.1 ± 0.2</td><td>1.1 ± 0.1</td><td>+28 ± 2</td><td>25(27)</td><td>99</td><td>39 ± 17</td><td>44 ± 17</td><td>5.9 ± 0.3</td><td>9.5</td><td>1.0 ± 0.1</td><td>-1 ± 36</td><td>60</td><td>58</td><td>6.4</td><td>7.4</td></tr></tbody></table>

n.a. indicates not analyzed.

ªH<sub>2</sub>O<sub>meas</sub> is the H<sub>2</sub>O solubility at 650 or 750 °C and 50 MPa calculated after Newman and Lowenstein [2002]; H<sub>2</sub>O<sub>meas</sub> is the glass H<sub>2</sub>O content measured by EMP (analytical uncertainty of ±0.7 wt%) and by Raman spectroscopy in bracket (average of 2 to 3 analyses with statistical error from ±0.2 to ±0.7 wt%).

ᵇθ<sub>b</sub> and θ<sub>b</sub> are the orientations to the horizontal of the bubbles and crystals, respectively; R gives the bubble mean eccentricity as the ratio of the long axes to the short axes.

ᶜΦ<sub>b,bulk_3D</sub> gives the bulk porosity measured on cubic samples of 1400-μm side volume and at smaller scale (600-μm side) in brackets; Φ<sub>b,connect_3D</sub> gives the percentage of bubbles connected to the sample outside, as defined in Supplementary Section I.

ᵈα<sub>metl</sub> is calculated after Jaupart and Allegre [1991] for a degassing from the H<sub>2</sub>O content of the pre-hydrated glass (H<sub>2</sub>O, in Table 1) to either H<sub>2</sub>O<sub>calc</sub> or H<sub>2</sub>O<sub>meas</sub>.

ᵉη<sub>metl</sub> is the viscosity of the hydrated melt calculated at 650 or 750 °C after Hess and Dingwell [1996] using H<sub>2</sub>O<sub>meas</sub>; η<sub>bulk</sub> is the viscosity of the three-phase (crystal, bubble, melt) suspension calculated after Truby et al. [2015] using η<sub>metl</sub>. Φ<sub>c,bulk_2D</sub>, and Φ<sub>b,metl_2D</sub> (see text and Supplementary Section 1 for details).
# 4. Discussion

## 4.1. Bubble nucleation and growth

Exsolving H₂O from 9.6 to 2.8 wt% [solubility at 50 MPa and 750 °C; Newman and Lowenstern, 2002] corresponds to a calculated vesicularity of 58 vol% [Jaupart and Allègre, 1991; Supplementary Section 1 (Equation (1))]. Yet the average porosities measured in the decompressed samples are mostly below 35 vol% (Figure 3a). We propose two explanations to account for this deviation; incomplete H₂O exsolution and partial loss of bubbles.

### 4.1.1. Incomplete H₂O exsolution

Checking the H₂O contents of the residual glasses in the undeformed samples using both electron microprobe and Raman spectroscopy, revealed contents in the crystal-bearing samples higher by 1.5 to 3.5 wt% than expected for an equilibrium H₂O exsolution (Figure 3b). This suggests incomplete H₂O exsolution during the ~1.5 h ($\tau_{DP} + \text{Dwell2} + \text{Dwell3}$ in Table 1) allocated to decompression-induced bubble nucleation. Yet, all series show equilibrium H₂O exsolution (~3.0 wt% H₂O) at $\gamma$ of 1.5–1.8, that is, for another couple of hours spent at 50 MPa, thus confirming that H₂O exsolution under the chosen experimental conditions requires more than 1 h to reach equilibrium. This strongly disagrees with the conclusions of Martel and Schmidt [2003] suggesting H₂O exsolution timescales of less than 1 min in microlite-bearing rhyolitic melts at 860 °C, but agrees with complete bubble expansion requiring several hours. Therefore, we infer that in our experiments, the timescale of bubble expansion controlled H₂O exsolution from the melt. In contrast, the undeformed crystal-free samples degassed to equilibrium, which we rather attribute to their lower initial H₂O content (5.0 instead of 9.6 wt%) rather than the absence of crystals, since serie10 which contains 10 vol% crystals, shows the highest H₂O exsolution disequilibrium. Further investigation would be required to validate the positive correlation between the initial H₂O content of the melt and the timescale of bubble expansion, which in turn controls the H₂O content of the residual melt.

### 4.1.2. Partial loss of bubbles

Knowing that H₂O exsolution was not completed in the H₂O-rich samples, the gas fractions were

calculated using the final H₂O contents measured in the residual glasses instead of using H₂O solubility at 50 MPa. The results show that the measured porosities are lower by 10 to 30 vol% than expected for an equilibrium degassing (Figure 3a), suggesting that these samples lost a part of their bubble cargo. One explanation could be that the melt partially outgassed by bubble connection at some point during bubble growth, and partially collapsed via compaction. Yet the released gas pockets cannot escape the sealed capsule and were considered in the measurement of total vesicularity. Therefore, we infer that these gas pockets were heterogeneously distributed in the capsule corners and were not systematically visible on our SEM pictures (for instance, a large gas pocket is visible below the capsule roof in Figure 2c), so that they may miss from some porosity budget. This hypothesis is supported by the lack of clear relationships between the bubble deficiency and parameters such as H₂O content, crystal content, or finite deformation.

### 4.1.3. Pre-deformation porosity

The undeformed samples show melt vesicularities of 19–30 vol% (Figure 3a) and bubble number densities not varying by more than 0.5 log unit (Figure 4a), which suggests a comparable and reproducible process of bubble formation. Within this 0.5 log unit, however, BND increases with crystal content, as experimentally observed by Cáceres et al. [2021]. The bubble number densities of $10^{9.3-10.6} m^{-3}$ measured from the XCT images are lower by 1 log unit compared to those predicted in crystal-free magmas by the formulation of Toramaru [2006] for decompression rates of 0.1–1.0 MPa/s (i.e., BNDs of $10^{10.5-11.5} m^{-3}$), highlighting a possible role of the crystals in controlling nucleation sites for bubbles, as already pointed out by Cáceres et al. [2021]. In the undeformed samples from all series, more than half (and up to 95% in serie0 and serie50) of the bubbles have diameters of about 20 µm. These small bubbles occupy about one third of total porosity in all series but serie0. In serie0, the volume of the small bubble population nearly disappeared at the expense of larger (220–240 µm diameter) bubble volume (Supplementary Figure 3a). In serie50, the population of 20-µm diameter bubbles coexists with some 240–260 µm diameter bubbles (Figure 2i) representing one third of total porosity (Supplementary Figure 3i).
The homogeneous spatial distribution of the small 20–40 µm bubbles (Figure 2a) confirms a process of homogeneous nucleation [Mourtada-Bonnefoi and Laporte, 2002], in agreement with the observations in natural and experimental products that plagioclase is not a wetting phase for nucleating H₂O bubbles [Navon et al., 1998, Cáceres et al., 2021]. The population of 200–300 µm bubbles may have three possible origins: coalescence of the 20–40 µm bubbles, a separate nucleation event, or pre-decompression bubbles. A separate nucleation event or bubble coalescence would lead to more or less similar bubble size distributions in the different series, which is not the case (Supplementary Figure 3a,d,f,i). Therefore, we argue in favour of bubbles that formed before decompression, as often observed in experiments starting with powdered material that is able to trap H₂O and air in the intergranular voids before melting [i.e., the “hydration bubbles” in Gardner et al., 1999]. The number and size of these bubbles would then be affected by the crystal content in the glass powder. The undeformed samples of serie30 and serie50 show that bubbles started channelling along the crystals (Figure 2f,i). According to the model of Burgisser et al. [2017], crystal-free melts containing 20–30 vol% bubbles, number density of 10⁹⁻¹¹ m⁻³, and size of 200–300 µm, would be impermeable, even with 50–70% connected porosity. In contrast, gas permeability measurements in natural crystal-bearing samples with connected porosities from 50 to 80% give permeabilities from ~10⁻¹⁴ to ~10² m² [Bain et al., 2019]. Therefore, in absence of permeability measurements in our samples, it is not clear whether the gas percolation thresholds were met before deformation.

### 4.1.4. Bulk viscosity

The bulk viscosities calculated using the formulation of Truby et al. [2015] range from 10⁵.²⁻⁵.⁴ Pa·s for the series containing 10 vol% crystals at 750 °C, 10⁶.⁵⁻⁷.⁴ Pa·s for the series containing 30 and 50 vol% crystals at 750 °C, to 10⁷.¹⁻⁷.³ Pa·s for the crystal-free series at 650 °C (Table 2). The regime of low capillary number (Ca < 1) calculated for our sample deformation conditions (see Supplementary Section 1.3.2 for details) results in considering the 20–30 bulk vol% bubbles as rigid particles that increase the hydrated melt viscosity by about 0.2 log unit (e.g. serie0 in Table 2). Adding crystals further increases bulk viscosity, from ~0.3 log unit for 10 vol% added crystals to

~0.7 log unit for 30 vol% crystals, and 1.0–1.2 log unit for 40–50 vol% crystals, which agrees with most of the common models of particle suspension rheology [see Mader et al., 2013].

## 4.2. Effect of the crystal network on porosity development under small shear ($\gamma < 2$)

### 4.2.1. Crystalline fabric

The crystals were nearly isotropic in shape ($R_c = 1.1 \pm 0.1$; Table 2), suggesting that the cutting of the labradorite block using the Selfrag apparatus did not permit to reproduce the elongated plagioclase crystals commonly observed in the natural magmas. This small anisotropy likely limited the development of strong crystalline fabrics in the sheared samples. Nevertheless, the crystals show main orientations that agree with those calculated as a function of $\gamma$ by Fernandez et al. [1983] for solid particles of $R_c = 1.1$ (Figure 4b). Yet, the high standard deviations associated with these angles, especially at $\gamma = 0$ and 0.6, suggest that the crystalline fabric is very weak, with the exception of serie50 at $\gamma = 0.6$ and 1.5, for which the crystalline fabrics are well marked (standard deviations <10°). The serie30 and serie50 samples show crystal breakage features, evidenced by the appearance of a population of small crystals (Figure 2i,k) that are also visible in the 2D crystal size distributions (Supplementary Figure 4). Crystals breakage was also observed in crystal-rich samples deformed under low strains by Forien et al. [2011], likely resulting from intense stress localization at some grain contact.

### 4.2.2. Bubble spatial distribution

The melt porosity of 20–30 area% is maintained in whatever crystal content and $\gamma$, with the exception of one sample (Figure 5a), which agrees with the results of Pistone et al. [2012] obtained on their samples sheared at T high enough to prevent bubble loss through shear fracturing. Yet in details, increasing both bulk crystal content $\ge$30 vol% and finite strain, increases heterogeneities in bubble spatial distribution, as reflected by porosity variations up to $\pm$20% between bubble-poor and bubble-rich zones, with the variation being positively correlated to $\gamma$ (Figure 5a). For these samples, the gas accumulated in local microstructures caused by the shear-induced crystal fabric (Figure 2h,k), in agreement with the observations of Laumonier et al. [2011]. In samples with
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Effects of the crystals on (a) melt porosity and (b) connected porosity; the solid line indicates the linear fit of the data but one point ($\gamma = 0$) from serie30 with an anomalously low connected porosity; the dashed lines contour the evolution of the connected porosity as a function of bulk crystallinity (0–50 vol%) for bulk vesicularities of 20–30 vol% and shear strain of $0 < \gamma < 2$.</figcaption></figure>

bulk crystal contents <30 vol%, however, our results suggest that the spatial distribution of the bubbles remain homogeneous (variations below ±5 vol%; Figure 5a).

### 4.2.3. Connectivity development

At $\gamma = 0$, bubble connectivity increases with crystal content (Figure 6), reaching ~87% connected porosity in the sample with bulk crystal contents of 50 vol% (Figure 5b; one exception is the undeformed sample of serie30 with only 37% connected porosity, for which we have no explanation). These results agree with those of Cáceres et al. [2021] who observed that phenocryst-bearing magma noticeably showed a higher degree of bubble coalescence, but disagree

with those of Okumura et al. [2012] who reported no clear effect of crystallinity (30–50 vol% recalculated with respect to the melt volume) on the degree of bubble coalescence and connectivity until melt vesicularities reached ~68 vol%. This disagreement may be explained by differences in bubble texture. Indeed, their crystal-rich samples with vesicularity ~20–30 vol% (i.e. those decompressed to 20 MPa; number densities of ~$10^{12}$ m⁻³) do not show bubble channelling features, such as those observed in the undeformed samples of our serie50 (Figure 2). This could also result from different bulk viscosities, different bubble to crystal size ratios, or unimodal versus polymodal bubble size distributions (i.e. the ~200–300 µm bubble population in our samples may have played a key role in early connectivity development).

In our samples, increasing bubble connectivity with increasing crystal content before shearing ($\gamma = 0$) highlights the role of the crystal network on the decompression-induced formation of bubbles. Firstly, the crystals possibly influence bubble nucleation, by increasing bubble number density by up to 0.5 log unit with increasing crystal content (see log BND at $\gamma = 0$ in Figure 4a). We speculate that physical processes at the origin of more bubble nuclei in crystal-richer magmas could involve modified gas-melt interfacial tensions or excess free energy. After nucleation, bubble growth must have been partly dictated by the crystal network, forcing bubble coalescence and channelling, eventually approaching full bubble connection in serie50 (i.e. >80% connectivity; Figures 5b and 6i). These observations confirm the results of the numerical modelling of Parmigiani et al. [2017] predicting maximum outgassing efficiency via 40–50% connected porosity in crystal-rich (40–70 vol%) magmas.

Increasing $\gamma$ from 0 to 2 does not significantly influence the amount of connected porosity (Figure 5b), thus highlighting the large predominance of crystal content over shear strain up to $\gamma = 2$ in favouring bubble connections. Fitting all the data but one (the low-connected sample at $\gamma = 0$ from serie30) shows a positive linear correlation between connected porosity and bulk crystal content ($\gamma = 0.84x + 52$; Figure 5b), thus valid for silicic magmas with 20–30 vol% bulk vesicularity and sheared at low strain ($\gamma < 2$) with shear rates of $10^{-4}$ s⁻¹ under a low confining pressure of 50 MPa. The rather poor
<figure><img src="image_9.png" /><figcaption>**Figure 6.**onde (3 mm diameters) for increasing shear strain from
the core samples (3 mm diameters) for increasing shear strain from
$\gamma = 0$ to 2 (from left to right) and crystal content from 0 to 50 vol% (from top to down). The left blocks were
processed using VGStudioMax, with the melt and crystals in yellow, the isolated bubbles in red and the
bubbles connected to the sample outside in grey. $\Phi_{b\_bulk\_3D}$ and $\Phi_{b\_connect\_3D}$ give the 3D percentage of
bulk and connected bubbles, as defined in Supplementary Section 1. The upper right insets show details
of the bubbles processed using Blob3D, with another colour code that show one colour per isolated
bubble or identified network of connected bubbles. The black line in (e) highlights an alignment of
bubbles connected to the sample outside, which denotes a post-experimental fracture.</figcaption></figure>

coefficient of determination, R² = 0.84, mostly reflects a large variability of the connected porosity at high crystallinity (±15% for crystal content ~45 vol%) with respect to crystal-free samples (±2%; Figure 5b).

### 4.3. Implications for natural volcanic systems

The experimental samples are silica-rich melts containing 20–30 vol% bubbles and 0 to 50 bulk vol%

plagioclase crystals that cover a range from natural moderately-porous obsidians to crystal-rich (mush-like) rhyolitic magmas. The 50–90 µm crystal size typically scales microlites to microphenocrysts [e.g. Hammer et al., 1999, Martel and Poussineau, 2007]. The low confining pressure of 50 MPa and the low shear strains ($\gamma < 2$) applied to the samples simulate a shear initiation at shallow depth (about 2 km deep, considering rock density of 2400 kg/m³), such
<figure><img src="image_10.png" /></figure>

**Figure 7.** Schematic representation of bubble connectivity in magmas with (a) 0–10 vol% crystals or (b) 30–50 vol% crystals, deformed under shear at magmatic temperatures (650–750 °C), moderate shear-strain rates (2 × 10⁻⁴ s⁻¹) and low pressures (50 MPa, that is at level of upper conduit to lava dome interior). The gas bubbles in crystal-poor magmas with 20–30 vol% porosity hardly connect under strain below γ = 2, so that outgassing requires fracture opening under high strains or strain rates, likely only reachable at conduit margins [Kushnir et al., 2017]. In contrary, a dense crystalline fabric favours bubble connectivity under strains as low as γ = 2 and bulk porosity of 20–30 vol%, so that outgassing may be an early process upon magma emplacement. Crystal fragmentation commonly occurs in crystal-rich magmas and heightens with increasing crystal content and strain. If temperature decreases and/or strain increases as expected at conduit margins, fractures may open and assist outgassing [Pistone et al., 2012, Shields et al., 2014].

as in a volcanic upper conduit. Overall, these conditions may be representative of Volcanian events or dome-related eruptions, during which a moderately-porous magma degasses and crystallizes microcrysts at shallow depths, eventually emplacing effusively (as domes or plugs) or erupting explosively (as surges), depending on outgassing efficiency. In particular, vesicularities of 20–30 vol% are representative of those measured in lava domes and blast-generated pyroclastic density currents, such as the ones generated in 1350 AD (P1) and 1902 at Montagne Pelée, Martinique [Martel et al., 2000] and in 1980 at Mount St Helens, USA [Hoblitt and Harmon, 1993].

The implications of our results to volcanic systems are schematized in Figure 7 and detailed

below. Our results suggest that H₂O bubbles nucleating in undeformed aphyric or crystal-poor (0–10 vol%) magmas by decompression at rates of 0.1–1.0 MPa/s (i.e. bubble number densities around 10⁹⁻¹¹ m⁻³) can connect to 50% (Figure 6a,d), but the magma could still be impermeable mainly due to the moderate bulk porosity of 20–30 vol% [Burgisser et al., 2017]. Such moderately-porous crystal-poor melts may reach permeability via a process of fracturing at high strains (γ > 2) along the conduit margins, as experimentally demonstrated by Kushnir et al. [2017]. Adding crystals drastically increases bubble connectivity (>70%; Figure 5b), by favouring channelling even at low shear strains (Figure 2h,k). Channelling creates gas permeability of
the order of $10^{-10}$ m$^2$, which is two to five orders of magnitude higher than the other permeability processes at the equivalent bubble content [Collombet et al., 2021]. Therefore, magmas with matrix crystallinity >30 vol% and vesicularity of 20–30 vol%, as commonly observed in lava domes [Boudon et al., 2015], may show bubble connectivity >70% (Figure 5b) and gas permeabilities >$10^{-12}$ m$^2$ [Bain et al., 2019]. Our data suggest that bubble connectivity in low-porous (20–30 vol%) silica-rich magmas slightly sheared ($\gamma < 2$; strain rate of $10^{-4}$ s$^{-1}$) at low confining pressure (50 MPa; about 2 km deep) increases linearly with crystal content (following $y = 0.84x + 52$, with y being the connected porosity in vol% and x being the bulk crystal content in area%; Figure 5b). In natural samples, connectivity would be enhanced by unequal bubble sizes that favour bubble coalescence [Lin and Lin, 2009], microwRKTEs made of crystals with high aspect ratios [Lindoo et al., 2017], and polydisperse crystal size distributions involving phenocrysts and microlites [Cáceres et al., 2021]. Therefore, in comparison to crystal-poor magmas, crystal-rich melts are able to outgas much more easily (at lower $\gamma$) for comparable temperature, pressure, and strain rate, confirming the effusive mode of crystal-rich lava-dome emplacement and the unnecessary to invoke dominant outgassing mechanisms other than channelling [Collombet et al., 2021]. From a crystallinity point of view and regardless of melt H$_2$O content and magma storage pressure, the gas-permeable nature of the 30–50 vol% crystal-bearing rhyolitic melts investigated in this study agrees with the eruptions of low magnitude reported by [Takeuchi et al., 2021] and the crystallinity-controlled effusive eruptions reported by Popa et al. [2021]. Yet, the present experiments highlight the major control of an exsolved volatile phase (20–30 vol%) and low shear ($\gamma < 2$) on the permeability development, with consequences on eruptive style that may deviate from those expected from the sole consideration of crystallinity and melt H$_2$O content. Since low-sheared crystal-rich magmas at upper conduit level is permeable to gases, this poses the question of the mechanisms behind dome-related explosive eruptions, such as Vulcanian or Pelean explosions. External factors can be invoked, such as depressurization due to unloading by gravitational collapse of a part of the dome [Voight et al., 1981]. Among internal factors of self-explosivity, one

could invoke processes capable of building up overpressure in gas bubbles from an initially permeable magma, such as a second boiling event in response to late-stage extensive microlite crystallization [Sparks, 1997] and gas volume reduction by silica-phase deposition in bubbles [Boudon et al., 2015]. Another hypothesis for self-explosivity is permeability barriers and pressure accumulation via vertical rheological and lithological gradients in the magma column [e.g., Voight et al., 1999].

## 5. Conclusions

By deforming silicic melts with 20–30 vol% H$_2$O bubbles and 0 to 50 vol% crystals at magmatic temperatures (650–750 °C) and low confining pressures (50 MPa), we highlighted that:

* Torsion experiments using three-phase (bubbles, crystals, and melt) suspensions are valuable analogues to shear initiation in moderately-porous aphyric to crystal-bearing silicic magmas at shallow depth, such as in lava domes or upper conduits.

* Bubble connectivity linearly increases with crystal content.

* Increasing strain from $\gamma = 0$ to 2 in crystal-bearing magmas drastically increases the heterogeneity in bubble spatial distribution, by concentrating pores in shear zones.

* Crystal-rich magmas develop gas permeability via a process of channelling in shear zones, even at moderate porosities (20–30 vol%) and low strains ($\gamma < 2$).

* Natural crystal-rich magmas stalling in upper conduits or lava domes are thus likely permeable, so that violent dome explosions into devastating surges require additional mechanisms of gas pressurisation in shallow crystal-rich magmas.

## Conflicts of interest

Authors have no conflict of interest to declare.

## Acknowledgements

The authors greatly thank I. di Carlo, S. Erdmann, and P. Benoist for assistance with the SEM and EMP;
M. Hatton for the analyses using the elemental analyser; A. Slodczyk and F. Faranda for the analyses using Raman spectroscopy; P. Penhoud for the technical assistance with the X-ray microtomograph; and M. Beaulieu for the crystal block grinding using the Selfrag device at BRGM. This work is part of CD PhD thesis that was funded by the French MNRT grant. The study costs were provided by the CNRS-INSU TelluS program, the Agence Nationale de la Recherche (ANR-19-CE31-0007; LA), the EQUIPEX PLANEX project (ANR-11-EQPX-0036; B. Scaillet), and the LABEX VOLTAIRE project (ANR-10-LABX-100-01; B. Scaillet). The authors greatly thank B. J. Andrews and two anonymous reviewers for their helpful comments on a previous version of the manuscript, as well as R. Cioni and an anonymous reviewer for their comments on the present manuscript. We would also like to thank B. Scaillet for the editorial handling.

## Supplementary data

Supporting information for this article is available on the journal's website under https://doi.org/10.5802/crgeos.214 or from the author.

## References

* Bain, A. A., Lamur, A., Kendrick, J. E., Lavallée, Y., Calder, E. S., Cortés, J. A., Butler, I. B., and Cortés, G. P. (2019). Constraints on the porosity, permeability and porous micro-structure of highly-crystalline andesitic magma during plug formation. *J. Volcanol. Geotherm. Res.*, 379, 72–89.

* Blower, J. (2001). Factors controlling permeability-porosity relationships in magma. *Bull. Volcanol.*, 63(7), 497–504.

* Boudon, G., Balcone-Boissard, H., Villemant, B., and Morgan, D. J. (2015). What factors control superficial lava dome explosivity? *Sci. Rep.*, 5, article no. 14551.

* Burgisser, A., Chevalier, L., Gardner, J. E., and Castro, J. M. (2017). The percolation threshold and permeability evolution of ascending magmas. *Earth Planet. Sci. Lett.*, 470, 37–47.

* Burgisser, A. and Gardner, J. E. (2005). Experimental constraints on degassing and permeability in volcanic conduit flow. *Bull. Volcanol.*, 67, 42–56.

* Cáceres, F., Scheu, B., Colombier, M., Hess, K.-U., Feisel, Y., Ruthensteiner, B., and Dingwell, D. B. (2021). The roles of microlites and phenocrysts during degassing of silicic magma. *Earth Planet. Sci. Lett.*, 577, article no. 117264.

* Cassidy, M., Manga, M., Cashman, K., and Bachmann, O. (2018). Controls on explosive-effusive volcanic eruption styles. *Nat. Commun.*, 9, article no. 2839.

* Castro, J. M., Cordonnier, B., Tuffen, H., Tobin, M. J., Puskar, L., Martin, M. C., and Bechtel, H. A. (2012). The role of melt-fracture degassing in defusing explosive rhyolite eruptions at volcano Chaiten. *Earth Planet. Sci. Lett.*, 333, 63–69.

* Collombet, M., Burgisser, A., Colombier, M., and Gaunt, E. (2021). Evidence for deep gas loss in open volcanic systems. *Bull. Volcanol.*, 83(7), 1–16.

* deGraffenried, R. L., Larsen, J. F., Graham, N. A., and Cashman, K. V. (2019). The influence of phenocrysts on degassing in crystal-bearing magmas with rhyolitic groundmass melts. *Geophys. Res. Lett.*, 46, 5127–5136.

* Degruyter, W., Bachmann, O., and Burgisser, A. (2010). Controls on magma permeability in the volcanic conduit during the climactic phase of the Kos Plateau Tuff eruption (Aegean Arc). *Bull. Volcanol.*, 72, 63–74.

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

* Fernandez, A. N., Febesse, J. L., and Mezure, J. F. (1983). Theoretical and experimental study of fabrics developed by different shaped markers in two-dimensional simple shear. *Bull. Soc. Géol. France*, 7-XXV(3), 319–326.

* Forien, M., Arbaret, L., Burgisser, A., and Champalier, R. (2011). Experimental constraints on shear-induced crystal breakage in magmas. *J. Geophys. Res.*, 116, article no. B08217.

* Gardner, J. E., Hilton, M., and Carroll, M. R. (1999). Experimental constraints on degassing of magma: isothermal bubble growth during continuous decompression from high pressure. *Earth Planet. Sci. Lett.*, 168, 201–218.

* Gaunt, H. E., Sammonds, P. R., Meredith, P. G., Smith, R., and Pallister, J. S. (2014). Pathways for degassing during the lava dome eruption of Mount St. Helens 2004–2008. *Geology*, 42(11), 947–950.

* Gonnermann, H. M., Giachetti, T., Fliedner, C.,
Nguyen, C. T., Houghton, B. F., Crozier, J. A., and Carey, R. J. (2017). Permeability during magma expansion and compaction. *J. Geophys. Res., Solid Earth*, 122, 9825–9848.

Hammer, J. E., Cashman, K. V., Hoblitt, R. P., and Newman, S. (1999). Degassing and microlite crystallization during pre-climactic events of the 1991 eruption of Mt Pinatubo, Philippines. *Bull. Volcanol.*, 60, 355–380.

Hess, K.-U. and Dingwell, D. B. (1996). Viscosities of hydrous leucogranitic melts: non-Arrhenian model. *Am. Mineral.*, 81, 1297–1300.

Hoblitt, R. P. and Harmon, R. (1993). Bimodal density distribution of cryptodome dacite from the 1980 eruption of Mount St Helens, Washington. *Bull. Volcanol.*, 55, 421–438.

Jaupart, C. and Allègre, C. J. (1991). Gas content, eruption rate and instabilities of eruption regime in silicic volcanoes. *Earth Planet. Sci. Lett.*, 102, 413–429.

Klug, C. and Cashman, K. V. (1996). Permeability development in vesiculating magmas: implications for fragmentation. *Bull. Volcanol.*, 58, 87–100.

Kushnir, A. R. L., Martel, C., Champallier, R., and Arbaret, L. (2017). In situ confirmation of permeability development in shearing bubble-bearing melts and implications for volcanic outgassing. *Earth Planet. Sci. Lett.*, 458, 315–326.

Laumonier, M., Arbaret, L., Burgisser, A., and Champallier, R. (2011). Porosity redistribution enhanced by strain localization in crystal-rich magmas. *Geology*, 39(8), 715–718.

Launeau, P. and Cruden, A. R. (1998). Magmatic fabric acquisition mechanisms in a syenite: results of a combined anisotropy of magnetic susceptibility and image analysis study. *J. Geophys. Res.*, 103, 5067–5089.

Launeau, P. and Robin, P.-Y. (1996). Fabric analysis using the intercept method. *Tectonophysics*, 267, 91–119.

Lin, T.-J. and Lin, G.-M. (2009). Mechanisms of inline coalescence of two-unequal bubbles in a non-Newtonian fluid. *Chem. Eng. J.*, 155, 750–756.

Lindoo, A., Larsen, J. F., Cashman, K. V., and Oppenheimer, J. (2017). Crystal controls on permeability development and degassing in basaltic andesite magma. *Geology*, 45, 831–834.

Mader, H. M., Llewellin, E. W., and Mueller, S. P. (2013). The rheology of two-phase magmas: a review and analysis. *J. Volcanol. Geotherm. Res.*, 257, 135–158.

Martel, C., Bourdier, J.-L., Pichavant, M., and Traineau, H. (2000). Textures, water content and degassing of silicic andesites from recent plinian and dome-forming eruptions at Mount Pelée volcano (Martinique, Lesser Antilles arc). *J. Volcanol. Geotherm. Res.*, 96, 191–206.

Martel, C. and Iacono-Marziano, G. (2015). Timescales of bubble coalescence, outgassing, and foam collapse in decompressed rhyolitic melts. *Earth Planet. Sci. Lett.*, 412, 173–185.

Martel, C. and Poussineau, S. (2007). Diversity of eruptive style inferred from the microlites of Mt. Pelée andesite (Martinique, Lesser Antilles). *J. Volcanol. Geotherm. Res.*, 166, 233–254.

Martel, C. and Schmidt, B. C. (2003). Decompression experiments as an insight into magma ascent rates. *Contrib. Mineral. Petrol.*, 144, 397–415.

Mourtada-Bonnefoi, C. and Laporte, D. (2002). Homogeneous bubble nucleation in rhyolitic magmas: an experimental study of the effect of H₂O and CO₂. *J. Geophys. Res.*, 107(B4), article no. 2066.

Navon, O., Chekhmir, A., and Lyakhovsky, V. (1998). Bubble growth in highly viscous melts: theory, experiments, and autoexplosivity of dome lavas. *Earth Planet. Sci. Lett.*, 160, 763–776.

Newman, S. and Lowenstern, J. B. (2002). VOLATILE-CALC: a silicate melt-H₂O-CO₂ solution model written in Visual Basic for Excel. *Comput. Geosci.*, 28(5), 597–604.

Okumura, S., Nakamura, M., Nakano, T., Uesugi, K., and Tsuchiyama, A. (2012). Experimental constraints on permeable gas transport in crystalline silicic magmas. *Contrib. Mineral. Petrol.*, 164(3), 493–504.

Okumura, S., Nakamura, M., Takeuchi, S., Tsuchiyama, A., Nakano, T., and Uesugi, K. (2009). Magma deformation may induce non-explosive volcanism via degassing through bubble networks. *Earth Planet. Sci. Lett.*, 281(3–4), 267–274.

Okumura, S., Nakamura, M., Tsuchiyama, A., Nakano, T., and Uesugi, K. (2008). Evolution of bubble microstructure in sheared rhyolite: formation of a channel-like bubble network. *J. Geophys. Res., Solid Earth*, 113, article no. B07208.

Parmigiani, A., Degruyter, W., Leclaire, S., Huber, C., and Bachmann, O. (2017). The mechanics of shal-
low magma reservoir outgassing. *Geochem. Geo-phys. Geosyst.*, 18, 2887–2905.

Paterson, M. S. and Olgaard, D. L. (2000). Rock de-formation tests to large shear strains in torsion. *J. Struct. Geol.*, 22, 1341–1358.

Pistone, M., Caricchi, L., Ulmer, P., Burlini, L., Ardia, P., Reusser, E., Marone, F., and Arbaret, L. (2012). Deformation experiments of bubble-and crystal-bearing magmas: rheological and microstructural analysis. *J. Geophys. Res., Solid Earth*, 117, article no. B05208.

Popa, R. G., Bachmann, O., and Huber, C. (2021). Explosive or effusive style of volcanic eruption determined by magma storage conditions. *Nat. Geosci.*, 14, 781–786.

Saar, M. O. and Manga, M. (1999). Permeability-porosity relationship in vesicular basalts. *Geophys. Res. Lett.*, 26(1), 111–114.

Shields, J. K., Mader, H. M., Caricchi, L., Tuffen, H., Mueller, S., Pistone, M., and Baumgartner, L. (2016). Unravelling textural heterogeneity in ob-sidian: shear-induced outgassing in the Rocche Rosse flow. *J. Volcanol. Geotherm. Res.*, 310, 137–158.

Shields, J. K., Mader, H. M., Pistone, M., Caricchi, L., Floess, D., and Putlitz, B. (2014). Strain-induced outgassing of three-phase magmas during simple shear. *J. Geophys. Res., Solid Earth*, 119(9), 6936–6957.

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

Stasiuk, M. V., Barclay, J., Carroll, M. R., Jaupart, C., Ratte, J. C., Sparks, R. S. J., and Tait, S. R. (1996). Degassing during magma ascent in the Mule Creek vent (USA). *Bull. Volcanol.*, 58(2–3), 117–130.

Takeuchi, S., Nakashima, S., and Tomiya, A. (2008). Permeability measurements of natural and experimental volcanic materials with a simple permeameter: toward an understanding of magmatic degassing processes. *J. Volcanol. Geotherm. Res.*, 177(2), 329–339.

Takeuchi, S., Nakashima, S., Tomiya, A., and Shinohara, H. (2005). Experimental constraints on the low gas permeability of vesicular magma dur-

ing decompression. *Geophys. Res. Lett.*, 32, article no. L10312.

Takeuchi, S., Toshida, K., Miura, D., Ito, H., and Uesawa, S. (2021). Relationships between magmatic properties and eruption magnitude of explosive eruptions at Japanese arc volcanoes during the last one hundred thousand years. *J. Volcanol. Geotherm. Res.*, 419, article no. 107345.

Toramaru, A. (2006). BND (bubble number density) decompression rate meter for explosive volcanic eruptions. *J. Volcanol. Geotherm. Res.*, 154, 303–316.

Truby, J. M., Mueller, S. P. E., Llewellin, W., and Mader, H. M. (2015). The rheology of three-phase suspensions at low bubble capillary number. *Proc. R. Soc. A*, 471, article no. 20140557.

Tuffen, H., Dingwell, D. B., and Pinkerton, H. (2003). Repeated fracture and healing of silicic magma generate flow banding and earthquakes? *Geology*, 31(12), 1089–1092.

Voight, B., Glicken, H., Janda, R. J., and Douglass, P. M. (1981). Catastrophic rockslide avalanche of May 18. *U.S. Geol. Surv. Prof. Paper*, 1250, 347–377.

Voight, B., Sparks, R. S. J., Miller, A. D., Stewart, R. C., Hoblitt, R. P., Clarke, A., Ewart, J., Aspinall, W. P., Baptie, B., Calder, E. S., Cole, P., Druitt, T. H., Hartford, C., Herd, R. A., Jackson, P., Lejeune, A. M., Lockhart, A. B., Loughlin, S. C., Luckett, R., Lynch, L., Norton, G. E., Robertson, R., Watson, I. M., Watts, R., and Young, S. R. (1999). Magma flow instability and cyclic activity at Soufriere Hills volcano, Montserrat, British West Indies. *Science*, 283(5405), 1138–1142.

Westrich, H. R. and Eichelberger, J. C. (1994). Gas transport and bubble collapse in rhyolitic magma: an experimental approach. *Bull. Volcanol.*, 56(6–7), 447–458.

Woods, A. W. and Koyaguchi, T. (1994). Transitions between explosive and effusive eruptions of silicic magmas. *Nature*, 370(6491), 641–644.

Wright, H. M. N., Cashman, K. V., Gottesfeld, E. H., and Roberts, J. J. (2009). Pore structure of volcanic clasts: Measurements of permeability and electrical conductivity. *Earth Planet. Sci. Lett.*, 280(1–4), 93–104.