Syn-Eruptive Conditions of the AD 1530 Sub-Plinian Eruption of La Soufrière of Guadeloupe (Lesser Antilles)

A prerequisite in refining volcanic hazard at explosive volcanoes is a better quantification of the timescales of the syn-eruptive processes, such as magma degassing and crystallization prior to eruption. To this aim, new data on the matrices (microlites, residual glasses, and bubbles) of pumice, scoria, and dense clasts erupted during the AD 1530 andesitic eruption of La Soufrière of Guadeloupe are combined with published data from phase-equilibrium and kinetic experiments, in order to estimate pressures of microlite crystallization and magma ascent rates. From the timescale data, we infer that the AD 1530 eruption started with phreatomagmatic explosions tapping magmas that ascended during about 1 month (decompression rate of ∼50 Pa/s) from the coldest parts of the reservoir (∼825°C and a 74 wt% SiO2 melt). These magmas continuously crystallized microlites (∼25 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, quartz, and apatite), but did not outgas likely due to sealed conduit margins. The flank collapse (preexisting “cold” edifice) that followed the phreatomagmatic phase triggered a sub-Plinian eruption that progressively tapped the hotter main reservoir (∼875°C and 71 wt% SiO2 interstitial melt), emitting banded and homogeneous pumice. The banded pumice did not significantly outgas and mostly lack decompression-induced microlites, suggesting short ascent durations of the order of 0.5–1 day (decompression rates of 1,400–4,000 Pa/s). The following Strombolian phase emitted dark scoria that did not significantly outgas and only crystallized rare microlites, suggesting ascent duration of the order of 2 days (decompression rates of ∼550 Pa/s). The terminal lava dome growth involved fully outgassed magmas ascended during more than 1 month, giving time for microlite crystallization (∼40 vol% plagioclase, orthopyroxene, clinopyroxene, magnetite, and cristobalite). The detection of any shallow new magmatic intrusion is crucial, since it can trigger a sequence of conduit processes leading to an eruption marked by a succession of different and disastrous eruptive styles, following a scenario similar to the AD 1530 eruption. Overall, we provide a combined approach of petrological, geochemical, and experimental data that may be used to infer ascent conditions and rates at other volcanic systems.


INTRODUCTION
Subduction-arc volcanoes are highly hazardous because they commonly alternate effusive with explosive eruptions, and magmatic with non-magmatic (phreatic) activity. In the Lesser Antilles arc, the historical volcanic crises had major societal and economic impacts: the sadly notorious 1902 dome-forming eruption of Montagne Pelée (Martinique) killed 30,000 persons and devastated the town of Saint Pierre (Lacroix, 1904), the 1976 explosive phreatic event of La Soufrière (Guadeloupe) led to a controversial evacuation of 73,000 inhabitants from the surrounding areas (Hincks et al., 2014), and the 1995 on-going eruption of Soufrière Hills (Montserrat) has been responsible for the relocation of nine-tenths of the island inhabitants and brought a severe socio-economic crisis (Hicks and Few, 2015). These crisis situations are partly ascribable to our failure to forecast the date and the type of the eruption, mainly coming from the difficulty to unequivocally identify precursory signals of impending eruption in slow-moving viscous magmas. Recent studies on felsic magmas have proposed timescales between the last magmatic recharge/ mixing in the reservoir and the eruption (e.g., Metcalfe et al., 2021;Ostorero et al., 2021), which could be a step forward in hazard assessment, should this last deep magmatic injection be detected. However, from the study of four magmatic eruptions of La Soufrière (but not the AD 1530 eruption), including Plinian, Vulcanian, and Strombolian eruptive styles, Metcalfe et al. (2021) deduced that the recharge-to-eruption timescales were not related to the eruption style. This highlights the need to decouple the timescales of the magmatic processes in the reservoir from the timescales of the ascent processes in the volcanic conduit, when aiming at investigating the relationships between the syn-eruptive processes and the eruption style.
The processes controlling eruptive dynamics mostly take place in the volcanic conduit during magma ascent. It has been demonstrated that the eruption dynamic directly depends on the syn-eruptive degassing process: overpressurized gas bubbles trapped in the magma favor explosive eruptions, whereas gas permeability development and outgassing (through permeable or fractured conduits) allow effusive eruptions (Eichelberger et al., 1986;Jaupart and Allègre, 1991;Gonnermann and Manga, 2013). The key parameter controlling the degassing kinetics is magma ascent rate, with slow ascent rates giving time for outgassing whereas rapid ascent rates prevent gas escape (e.g., Woods and Koyaguchi, 1994). As a result, andesitic to rhyolitic magmas ascending at rates of the order of m/s may lead to highly explosive eruptions whereas those ascending at rates of the order of mm to cm/s may lead to lava dome growths (Rutherford and Gardner, 2000).
The decompression of a volatile-saturated magma forces the volatiles to exsolve (following the gas solubility laws), forming gas bubbles. Because the liquidus T of a hydrous magma generally increases with decreasing melt volatile content, ascending magmas are expected to crystallize. The strong disequilibrium that comes along with rapid changes in volatile content generates nucleation-driven crystallization, i.e., lot of small crystals designated as microlites (crystals <50-100 µm in length, as opposed to the larger crystals that grew at near-equilibrium in the reservoir and referred as phenocrysts). Based on experimental studies, microlite textures (number density, size, shape) in silicic melts have been used to back-track crystallization kinetics in natural samples (e.g., Hammer and Rutherford, 2002;Martel and Schmidt, 2003) and decompression rate-meters have been proposed from number densities of bubbles and microlites (Toramaru, 2006;Toramaru et al., 2008). This highlights that combining experimental data (from equilibrium and kinetic experiments) and natural data on sample matrices, provides valuable information on the degassing and crystallization kinetics that decipher eruptive dynamics, as demonstrated for instance for Montagne Pelée, Martinique, recent eruptions (Martel, 2012), Mount Saint Helens, United States, 1980-1986eruption (Geschwind and Rutherford, 1995Riker et al., 2015), Mount Unzen, Japan, 1991-1995 eruption (Cichy et al., 2011), and Rabaul, Papua New Guinea, recent eruptions (Bernard and Bouvet de Maisonneuve, 2020). Therefore, our ability to forecast future eruptions and assess hazards also relies on our ability to measure in-situ magma ascent rates (e.g., via geophysical methods) and to relate them to eruptive scenarios (established via volcanological, petrological, and numerical models).
La Soufrière of Guadeloupe is a well-monitored active volcano of the Lesser Antilles subduction arc ( Figure 1A). Since its construction 0.25 Ma ago, La Soufrière has shown a large diversity of eruptive activity, alternating magmatic to non-magmatic events , which generates a multitude of eruptive scenarios making hazard assessment highly complex . The present lava dome ( Figure 1B) was formed during the magmatic eruption of AD 1530. Since this eruption, a small Vulcanian eruption (Legendre, 2012) and six phreatic eruptions have occurred. Since the last phreatic eruption in 1976-1977, La Soufrière volcano has become intensively monitored with a network of pluridisciplinary methods that are implemented by the Guadeloupe Volcanological and Seismological Observatory (OVSG-IPGP). Since 1992, the OVSG has recorded a systematic and progressive increase in shallow depth low-energy seismicity and slow temperature rise of the thermal springs close to the lava dome (OVSG-IPGP website). This unrest situation has boosted studies aiming at understanding and forecasting future eruptions at La Soufrière. For instance, recent geophysical studies of the internal structure of the lava dome have warned about its instability due to the volume of hydrothermally-altered rocks and the likelihood of a future major flank collapse (Rosas-Carbajal et al., 2016). Because the last 400 years of La Soufrière activity have not involved juvenile magma, little attention has been given to magmatic eruptions. Yet, during the 1976-77 crisis, geophysical and geochemical evidence supported the emplacement of a smallvolume intrusion of a viscous magma body that stopped within a few kilometers below the surface (Feuillard et al., 1983;Villemant et al., 2005Villemant et al., , 2014Boichu et al., 2011). This underscores the importance of integrating magmatic eruptions in future eruptive scenarios of La Soufrière. The AD 1530 eruption is currently viewed as a model for a future magmatic reactivation of the volcano . Because it involved a partial flank collapse, a sub-Plinian phase (characterized by a Volcano Explosivity Index VEI of 2-3 and an unsteady sustained convective column of 9-12 km high) with pumice-and-scoria fallout and flows, lava fountaining, and a terminal lava dome growth , this eruption is one of the most complex and dreaded in terms of hazard assessment Spence et al., 2008). In the light of the present volcanic unrest, for which a magmatic origin cannot be clearly excluded, there is an urgent need of knowledge of the eruptive mechanisms that may lead to an explosive magmatic event at La Soufrière.
To this aim, the present study constrains the magma ascent conditions and timescales that lead to the AD 1530 eruption of La Soufrière, by characterizing the textures and compositions of the matrix of the erupted products, and relate them to published experimental and numerical data, in order to assess the pressures (depths) of magma degassing, crystallization, and fragmentation, and the ascent durations of the magma.

Chronology of the Recent Activity
The composite "Grande Découverte-La Soufrière" volcanic complex, hereafter referred as "La Soufrière" was built during three main phases: Grande Découverte (0.25 Ma-42 ka), Carmichaël (42-11.5 ka), and La Soufrière phase that started 11.5 ka ago (Boudon et al., 1988(Boudon et al., , 1989. During the ongoing La Soufrière phase, the volcanic activity was roughly centered on the area where the current Soufrière lava dome is located. A series of ten flank collapses of low volume occurred during this phase, generating debris avalanches that flowed on the southwestern flank of the volcano, where the towns of Basse Terre and Saint Claude are now settled. Most of these debris avalanche reached the sea, 10 km away from the summit Boudon et al., 2007;Legendre, 2012). Several dome-forming and small sub-Plinian eruptions occurred during this period (Legendre, 2012), including the last magmatic AD 1530 eruption at the origin of the Soufrière lava dome (see description below). Some deposits dated to AD 1657 have been interpreted as a small explosive Vulcanian eruption (Legendre, 2012), which would be the last magmatic event of La Soufrière. The historical period is characterized by six phreatic eruptions (AD 1690(AD , 1797(AD -98, 1812(AD , 1836(AD -37, 1956(AD , and 1976 taking place from fractures and vents on La Soufrière lava dome (Barabé and Jolivet, 1958;Feuillard et al., 1983;Komorowski et al., 2005). The last phreatic eruption in 1976-1977 generated an important seismic crisis, 26 phreatic FIGURE 1 | (A) Location map of the Lesser Antilles, Guadeloupe archipelago, and La Soufrière volcanic complex, after Boudon et al. (2008), (B) Photo of La Soufrière present dome and (C) Typical stratigraphy section of the AD 1530 eruption, combining the sequence deposits at the sites of the Carbet Springs and Galion River , showing the origin of the studied samples (from Boudon et al., 2008;Pichavant et al., 2018).

Eruption Dynamics
The AD 1530 eruption represents the last major magmatic event of La Soufrière and likely the only magmatic eruption during the last 1,400 years. It was first described by Vincent et al. (1979) and Semet et al. (1981), and revised by Boudon et al. (2008). Following a small phreatic phase, the magma reached the surface and generated a phreatomagmatic explosion followed by a flank collapse of the summit part of the volcano, the debris avalanche of which flowed on the southwestern flank of the volcano and reached the sea. Subsequently to the flank collapse, the volcanic conduit was decompressed down to the magmatic chamber, leading to an explosive sub-Plinian phase emitting pumice fallout and column-collapse pyroclastic density currents channeled into valleys. The column reached a height between 9 and 12 km, producing about 3.9 × 106 m 3 DRE of juvenile products . The emitted products were first grey-white followed by white and blackand-white banded pumices. This event was followed by nearvent pyroclastic scoria fountains and ended with the growth of the present-day lava dome .

Pre-Eruptive Conditions
The evolution of the magma reservoir prior to eruption was investigated by Pichavant et al. (2018), through a petrological and experimental study aiming at constraining both, the storage conditions of the resident magma and the characteristics of the intruding magma. These authors concluded that the main andesitic magma contained ∼46-60 vol% phenocrysts in equilibrium with a rhyolitic melt bearing 5.0 to 5.5 wt% H 2 O, and was stored at a temperature (T) of 875 ± 25°C and pressure (P) of 175 ± 20 MPa (about 7 km deep, considering a rock density of 2.4 g/cm 3 ). These pre-eruptive H 2 O contents agree with the highest H 2 O contents of ∼5.5 wt% measured in glass inclusions trapped in plagioclase and orthopyroxene phenocrysts. There is evidence from melt inclusions and residual glasses that T was locally as cold as 825 ± 25°C. The residual glasses are indicators of mixing between the resident and intruding melts, showing strong compositional heterogeneities from 64 to >76 wt% SiO 2 .

Deep Mafic Intrusions
Based on the wide compositional range of whole rocks and textural plus mineralogical evidence, Semet et al. (1981) first proposed that La Soufrière andesitic reservoir was remobilized by the intrusion of a deep mafic magma, with an about 40-60% proportion of arrival on resident magma. Pichavant et al. (2018) experimentally constrained the end-member magma composition to ∼50 wt% SiO 2 for the basaltic intrusion and 62 wt% SiO 2 for the resident andesite. Although long timescales of reservoir replenishment of about 35 ka have been proposed using U-Th-Ra disequilibria methods (Touboul et al., 2007), Pichavant et al. (2018) highlighted much shorter timescales of the order of tens of hours, based on diffusion modelling in orthopyroxenes (following the method given in Solaro et al., 2020).

Parental Melts to Ascent Crystallization
In order to reconstruct the ascent-crystallization paths of the magmas in the AD 1530 sequence, one needs to infer the chemical composition of the melt from which the microlites crystallized. Ideally, this parental melt is the one in chemical equilibrium with the phenocryst rims in the storage region prior to magma ascent, and may be preserved as glass inclusions in these crystals. Yet, eruptions are often triggered by an intruding more mafic magma, whose melt mixes with the resident one, so that the melt from which microlites crystallize upon ascent may be more mafic than expected from the initial resident melt. Using the AD 1530 fused whole rocks as starting material (i.e., chemically including any late injection of mafic melt), the composition of the pre-eruptive interstitial melt was determined to be 71 ± 1 wt% SiO 2 from the phase-equilibrium experiments (Pichavant et al., 2018). Yet, the glass inclusions trapped in phenocrysts of plagioclase and orthopyroxene suggest melt SiO 2 contents of 74 ± 1 wt% (recalculated on anhydrous basis; Pichavant et al., 2018). Therefore, melts parental to microlites could span a SiO 2 range from ∼71 to ∼74 wt%, after and before mafic melt injection, respectively. As for the P-T conditions of ascent crystallization, it is commonly accepted that decompressioninduced microlites crystallize nearly isothermally, that is, for the AD 1530 eruption, at pre-eruptive T of ∼875 or 825°C and P lower than the pre-eruptive P of 175 MPa.

Natural Samples
The studied samples come from the lava dome ( Figure 1B) and the deposit sequences of the Carbet Springs O1215 site and Galion River site . Both deposit sequences have been combined to provide a typical sequence for the AD 1530 deposits ( Figure 1C). The samples encompass: three white pumice clasts, including SG7B from the phreatomagmatic unit and SG5A and SG6A from the base of the sub-Plinian unit; two banded pumice samples from the upper part of the sub-Plinian unit: SG11A, in which SG11A_white and SG11A_dark refer to the white and dark layer, respectively, and O1215Eb that represents the white layer of a banded pumice; one dark scoria (SG3) from the Strombolian phase, and a dense block from the AD 1530 lava dome (1227A). The upper sub-Plinian O1215Eb pumice and the 1227A dome clast were sampled by Boudon et al. (2008). The Strombolian SG3 dark scoria, the lower sub-Plinian SG5A and SG6A pumice samples, the phreatomagmatic SG7B pumice, and the upper sub-Plinian SG11A pumice were sampled by Pichavant et al. (2018). These samples are basaltic andesites to andesites, with bulk-rock 56.1 to 61.7 wt% SiO 2 , 2.7 to 4.2 wt% MgO, and 3.4 to 4.3 wt% alkalis Pichavant et al., 2018).

Sample Textures
The petrological studies were conducted on thin sections made from the samples, using a scanning electron microscopy (SEM; Merlin Compact ZEISS at ISTO). They served first to identify Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 glass areas and microlite phases favorable to chemical analysis and, second, to select representative zones for the textural analysis of the matrices.

Chemical Compositions
Crystal and glass chemistries were obtained using either an electron microprobe (EMP; Cameca SX Five at ISTO) or by energy dispersive spectroscopy (EDS; Bruker Quantax XFlash6 30 mm 2 ) on the SEM. The glass H 2 O contents were measured using either the EMP following a modified "by-difference" method or a Secondary Ion Mass Spectrometer (SIMS; IMS 1280 HR2 spectrometer at the Centre de Recherches Pétrographiques et Géochimiques, CRPG, Nancy). The analytical conditions for these techniques are detailed in Supplementary Material. Importantly for the present study is that the CO 2 contents are below 50 ppm in the residual glasses and mostly below 250 ppm in the glass inclusions. Therefore, we disregarded CO 2 as a component of the fluid phase and we considered a pure H 2 O fluid composition.

Silica-Polymorph Identification
Silica polymorphs (quartz and cristobalite) were identified by Raman spectroscopy (Renishaw InVia spectrometer at a wavelength of 514 nm and a maximum power of 100 mW; BRGM-ISTO), following the method described in Martel et al. (2021).

Textural Analysis
Selected SEM images were processed using the SPO software (Launeau and Robin, 1996), following the method described in Martel and Poussineau (2007). The selected images were only matrix (glass, microlites, and bubbles), i.e., avoiding phenocrysts. Bubble area percentage (Φ b ), plagioclase microlite area percentage (Φ Plag ), and area number density of plagioclase microlites (N A ) were determined as follows: Φ b (bubble area/matrix area) x 100; note that Φ b was measured on SEM images nesting scales from about 200 × 300 μm to 1,000 × 2000 µm for each sample, which may be slightly underestimated compared to bulk vesicularities performed at macroscopic scale, due to presence of possible fractures or larger bubbles. Φ Plag (plagioclase microlite area/bubble-free matrix area) x 100. N A number of plagioclase microlites/bubble-free matrix area.

Quench Pressure
The H 2 O content of residual glasses can be used to estimate the P at which the magma was quenched, using H 2 O solubility models and neglecting post-quenching alteration processes, such as re-or dehydration of the glasses. The quench P were calculated using the rhyolite solubility model of Newman and Lowenstern (2002) at T of 875 and 825°C (experimentally-determined main pre-eruptive T and coldest pre-eruptive T of the AD 1530 magma, respectively; Pichavant et al., 2018).
Quench P may also be back-calculated from the gas fraction equation of Jaupart and Allègre (1991), which calculates the gas volume fraction (α) resulting from isothermal degassing of a H 2 O-saturated melt from initial to final (quench) P, assuming a closed-system degassing (in that gas bubbles remain in contact with the melt). Back-calculating the quench P implies that the gas fraction becomes the vesicularity measured in the natural samples (Φ b /100) and that the H 2 O content of the residual melt follows a solubility law for silicic melts. The details of the calculation are given in Supplementary Material.

Crystallization Pressure
Because plagioclase composition is highly sensitive to melt H 2 O content, natural microlites can give information about the H 2 O content of the melt from which they crystallize and, by inference, on their crystallization P assuming melt H 2 O saturation and unaltered melt composition due to other processes. This approach requires experimental calibrations of plagioclase composition for melt and P-T conditions appropriate to the case study. It is also possible to use experiments to calibrate the plagioclase fraction and the residual melt SiO 2 content as a function of the P-T conditions. The selection of the experiments and the methodology to calibrate microlite crystallization in the case of the AD 1530 eruption are detailed in Supplementary Material. The results show that the crystallization P can be determined using two sets of equations: one for melts with 71 wt% SiO 2 and another one for melts with 74 wt% SiO 2.

RESULTS
Selected SEM images of the samples are shown in Figure 2 and textural and compositional analyses are summarized in Table 1.

H 2 O Content of the Residual Glasses
The H 2 O contents of the residual glasses determined by EMP and SIMS are reported in Table 1 and displayed in Figure 6. The EMP and SIMS values are in good agreement within respective analytical errors, i.e., not differing by more than 0.5 wt%. One

Quench Pressure
The quench P calculated using the H 2 O solubility model of Newman and Lowenstern (2002) are reported in Figure 7; note that calculating P at 825 and 875°C does not differ by more than ∼5 (Figure 7). Yet for the 1227A dome sample, the degassing obviously did not occur in closed-system, so that the quench P of 130 MPa can be discarded.

Crystallization Pressure
The calculated equilibrium P (using Eqs 2-6 in Supplementary Material) are shown in Figure 7 for p ≤ 200 MPa. p > 200 MPa, i.e., higher than the pre-eruptive P, were calculated for all samples, and mostly concern An 65-80 microlites that likely grew from melt compositions other than rhyolitic. The phreatomagmatic SG7B pumice shows crystallization P calculated from microlite composition from 5 to 115 MPa (71 wt% SiO 2 melt) and from 15 to 160 MPa (74 wt% SiO 2 melt). The crystallization P calculated from microlite crystallinity (i.e., plagioclase content) is 25 ± 5 MPa. If the residual glass SiO 2 equations are used, P are <10 and <65 MPa, respectively with the 71 and 74 wt% SiO 2 melts. In the lower sub-Plinian SG6 pumice, the P calculated from plagioclase composition range from ∼100 to >175 MPa for a 71 wt% SiO 2 melt and is > 135 MPa for a 74 wt% SiO 2 melt. From microlite crystallinity, P is 130 ± 25 MPa. From glass SiO 2 , P is ∼55-80 MPa and ∼110-130 MPa when calculated from a 71 and 74 wt% SiO 2 melt, respectively. In the lower sub-Plinian SG5A pumice, most of the P calculated from microlite composition are >125 MPa, but rare P of ∼75 MPa are also obtained. P calculated from crystallinity are 50 ± 10 MPa. P deduced from the residual glass SiO 2 are ∼100 and ∼160 MPa, when calculated from a 71 and 74 wt% SiO 2 melt, respectively. In the upper sub-Plinian O1215Eb, SG11A_white, and SG11A_dark pumice samples,  Table 1; the curves represent equilibrium crystallization starting from melts with 71 and 74 wt% SiO 2 (plain and dashed lines, respectively), as calculated from combining Eq. 3 with Eq. 4 and Eq. 6 in Supplementary Material. For comparison, the blue triangles represent the November 5, 2010 sub-Plinian event of Merapi, Indonesia (from Preece et al., 2014), and the red diamonds show the P1, P2, P3, and P4 Plinian eruptions of Montagne Pelée, Martinique (Martel and Poussineau, 2007).
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 8 plagioclase microlite compositions are too calcic and no p < 175 MPa are calculated. In the O1215Eb pumice, the P is 140 ± 5 MPa from microlite crystallinity. From residual glass SiO 2 , P are of 70-120 MPa and 125-180 MPa, respectively with the 71 and 74 wt% SiO 2 starting melts. In the Strombolian SG3 dark scoria, most of the P calculated from microlite composition are >125 MPa, but rare P of ∼75 MPa are also obtained. P calculated from crystallinity are 60 ± 22 MPa. P deduced from the residual glass SiO 2 are 100-135 MPa and 150-200 MPa, for a 71 and 74 wt% SiO 2 melt, respectively. In the 1227A dome sample, P calculated from microlite compositions range from 25 to >200 MPa, is of 20 ± 5 MPa from crystallinity, and of ∼50 MPa and <5 MPa, when calculated from residual glass SiO 2 starting from a 71 and 74 wt% SiO 2 melt, respectively.

Refining Starting Melt Composition
To better fix the crystallization conditions of the microlites, it is useful to look for further constraints on the melt composition (i.e., either 71 or 74 wt% SiO 2 ). The phreatomagmatic SG7B matrix contains quartz, which requires SiO 2 -saturated melts, i.e., SiO 2 contents > ∼79 wt% at ∼50 MPa (Martel et al., 2021). From the H 2 O-SiO 2 contents of the residual glasses, quartz crystallization in SG7B would not be possible at 875°C starting from a 71 wt% SiO 2 melt, whereas it may be possible at 825°C if the starting melt contains 74 wt% SiO 2 ( Figure 6). In contrast, the cristobalite-bearing 1227A dome sample is consistent with an origin by equilibrium crystallization from a 71 wt% SiO 2 starting melt at 875°C (Figure 6). In the upper sub-Plinian SG11A banded pumice, the residual glass compositions differ between the white and dark bands ( Figure 5). Yet, crystallization from a melt with 74 wt% SiO 2 in SG11A_white and 71 wt% SiO 2 in SG11A_dark yields overlapping P of residual glass formation in the 130-175 MPa range (Figure 7), validating crystallization from a 74 wt% SiO 2 melt in the white bands whereas from a 71 wt% SiO 2 melt in the dark bands. The lower sub-Plinian SG6A pumice and the upper sub-Plinian O1215Eb pumice show residual glasses that are slightly more silicic than expected from equilibrium crystallization from a 71 wt% SiO 2 melt ( Figure 3C). In addition, both samples contain apatite microlites, which were also found in Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 10 the phreatomagmatic SG7B pumice for which the initial melt has been determined to contain 74 wt% SiO 2 . Therefore, the microlites in the SG6A and O1215Eb pumice samples likely crystallized from a 74 wt% SiO 2 melt. For the lower sub-Plinian SG5A pumice and the Strombolian SG3 dark scoria, the choice of the starting melt composition is less clear. Yet, the SG3 sample does not contain apatite, is dark in colour and has a bulk-rock composition (∼58 wt% SiO 2 and ∼3.8 wt% MgO) close to that of the lava dome, comes from a Strombolian event typically involving melts more mafic than for sub-Plinian events, and chronologically just preceded the terminal dome growth, so that it is reasonable to assume a pre-eruptive melt composition similar to that of the dome, i.e., with 71 wt% SiO 2 . As for the SG5A pumice, both SiO 2 contents, 71 and 74 wt%, will be considered.

Degassing Regime
Comparing the vesicularities measured in the samples to the calculated gas fraction for a closed-system degassing (from Eq. 1 in Supplementary Material) provides information on the magma degassing regime. Figure 8 shows that the phreatomagmatic SG7B pumice, the lower sub-Plinian SG6A pumice, the upper sub-Plinian O1215Eb pumice, and the Strombolian SG3 dark scoria have measured vesicularities close to the theoretical ones, suggesting that they probably degassed in closed system. The large variability in the calculated vesicularities for the phreatomagmatic SG7B pumice result from the high variability in the measured H 2 O contents of the residual glasses (Table 1). In particular, locally high H 2 O contents (∼3.0 wt%) lead to theoretical vesicularities (∼20 vol%) much lower than actually measured (∼60 vol%), which possibly result from rehydration of the residual glasses in the external H 2 O-rich context that characterizes phreatomagmatism. Contrary to these samples, the lower sub-Plinian SG5A homogeneous pumice and the upper sub-Plinian SG11A banded pumice seem to have partially lost bubbles, which could result from local permeability development and partial outgassing in an opensystem degassing regime. For the SG5A pumice, the H 2 O contents of the residual glasses are lower than for the other pumice (0.7 vs >1.1 wt%; Table 1), whereas the measured FIGURE 6 | H 2 O versus SiO 2 content of the matrix glasses. The H 2 O contents are determined by the by-difference method using the electron microprobe, except where "(sims)" is mentioned which refers to H 2 O determination by ion probe (as given in Table 1). Pressures on the right Y-axis are H 2 O solubility pressures calculated after Newman and Lowernstern (2002). The purple open circles are the glass inclusion data given in Pichavant et al. (2018). The crossed circles are experiments from Pichavant et al. (2018) and the colored box gives the experimentally-determined pre-eruptive conditions. The curves show the trends for equilibrium crystallization starting from a 71 wt% SiO 2 melt (plain curve; after Eq. 4 in Supplementary Material) and from a 74 wt% SiO 2 melt (dashed curve; after Eq. 6 in Supplementary Material). The dashed lines show the silica polymorph boundaries in rhyolitic melts as defined in Martel et al. (2021).
The 1227A lava dome clast shows a much lower vesicularity (<10 vol%) than calculated (>90 vol%) from its low residual glass H 2 O content (0.4 wt%), highlighting a strong outgassing, as expected from the slow ascent rates commonly prevailing during lava dome emplacements. Martel and Iacono-Marziano (2015) demonstrated that the timescales of degassing, i.e. of bubble nucleation, growth, coalescence, and outgassing, primarily depend on the viscosity of the vesiculated magma under very low shear. The formulation of the bubble-melt viscosity used by these authors was set for crystal-free melts, so that the degassing timescales have to be taken with caution when dealing with crystal-bearing magmas. The bubble-melt viscosity was calculated from the steady simpleshear flow equation of Llewellin et al. (2002), calculating the viscosity of the hydrated residual glass at 875 and 825°C after Hess and Dingwell (1996) and using a capillary number (Ca) from 1.0 to 2.5 (characterizing little deformed bubbles). The calculations give bulk viscosities of 10 5.1-5.3 Pa s at 875°C and 10 5.6-5.9 Pa s at 825°C, and do not significantly differ for Ca varying from 1.0 to 2.5. Since all pumice samples show textural evidence of bubble coalescence, the ascent durations must have been longer than 1 h (Martel and Iacono-Marziano, 2015; their Figure 6). The magmas of the phreatomagmatic SG7B pumice, the lower sub-Plinian SG6A pumice, the upper sub-Plinian O1215Eb pumice, and the Strombolian SG3 dark scoria, did not outgas significantly, so that maximum ascent durations could be 10-100 h (0.5-4 days). The partially-outgassed lower sub-Plinian SG5A pumice and upper sub-Plinian SG11A banded pumice have calculated bulk viscosity from 10 5.2 Pa s for T of 875°C to 10 5.9 Pa s at 825°C, which suggest ascent durations >100 h. The 1227A lava dome sample that nearly fully outgassed has a calculated bulk viscosity of 10 5.2 and 10 5.8 Pa s at 875 and 825°C, respectively, suggesting emplacement duration >1,000 h (>40 days).

Ascent Crystallization
Ascent crystallization of microlites can be characterized by two main criteria. Firstly, microlite compositions should be compatible with crystallization from the pre-eruptive melt composition, i.e., with compositions more evolved than coexisting phenocryst rims (assuming an ideal case of phenocrysts no more growing upon ascent). Secondly, the proportion of microlites grown during ascent cannot be higher than expected for decompression rates permitting equilibrium crystallization; higher microlite contents would reflect inherited crystals. Based on these criteria, only two samples (the phreatomagmatic SG7B pumice and the 1227A dome clast) show major ascent-crystallization of microlites, three of them (the sub-Plinian O1215Eb, SG6A, and SG11A pumice samples) Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 appear devoid of decompression-induced crystallization, and two others (the sub-Plinian SG5A pumice and the Strombolian SG3 scoria) show scarce ascent-crystallization, as detailed below. The phreatomagmatic SG7B pumice fulfills the criteria for ascent-crystallization: the microlite compositions are more evolved than those crystallizing under pre-eruptive conditions ( Figure 4) and the microlite content agrees with equilibrium crystallization at 825°C from a 74 wt% SiO 2 starting melt ( Figure 3C). The continuous range of plagioclase microlite compositions suggests steady-state crystallization from reservoir to P of ∼15 MPa (Figure 7); note that ascentcrystallization also took place as overgrowths on phenocrysts, as evidenced by some ∼ An 28-40 plagioclase and En 53-54 orthopyroxene phenocryst rims ( Figures 4A-C). The good agreement between crystallization P deduced from compositions and contents of plagioclase microlites and those calculated from the residual glass compositions (Figure 7) suggests a near-equilibrium crystallization process. Moreover, these crystallization P agree with the degassing P calculated from vesicularity and H 2 O contents of the residual glasses (Figure 7), so that similar conditions are established for degassing and crystallization in the conduit.
For the 1227A dome sample, the P calculated from the compositions of the plagioclase microlites suggest that microlites continuously crystallized during ascent (Figure 7). Orthopyroxene and clinopyroxene microlite compositions ( Figures 4C,D) confirm continuous crystallization during ascent. Degassing P retrieved from residual glass H 2 O contents are <30 MPa, reaching p < 5 MPa using residual glass compositions, highlighting a crystallization process close to equilibrium, although likely partly hampered by the slow crystallization kinetics prevailing at low P. Yet, lava domes commonly cool upon emplacement, thus driving plagioclase to more albitic compositions than expected from isothermal decompression, so that the crystallization P calculated for this sample have to be taken as minimum values. The Mt 55-58 titanomagnetite microlites ( Figure 4B) show about 10 mol% Mt less than both the microlites in the other samples and the dome phenocrysts (Supplementary Material). This confirms late-emplacement cooling of the AD 1530 lava dome, rather than evidencing more reduced conditions that appear unlikely since the dome cooled in air.
The lower sub-Plinian SG6A pumice and the upper sub-Plinian O1215Eb, SG11A_white, and SG11A_dark pumice, all show microlite compositions similar or less evolved than the preeruptive ones (Figure 4), with calculated p > 130 MPa (Figure 7). Therefore, most of microlites must have crystallized at reservoir depth prior to eruption, without significant decompressioninduced crystallization. The strong contrast between these crystallization P and the degassing P (≤50 MPa) deduced from vesicularity and residual glass H 2 O contents (Figure 7) shows that crystallization and degassing were largely decoupled. Since degassing is a process much faster than crystallization (e.g., Martel and Schmidt, 2003), strong kinetic crystallization delays are implied. In both bands of the SG11A pumice, the orthopyroxene microlites are more magnesian than the phenocryst rims ( Figure 4C), likely reflecting crystallization from a more mafic melt. The SG11A_dark sample additionally contains Wo 44-47 En 41-46 clinopyroxene microlites, which cover the compositions of the phenocrysts, from common Wo 42 (∼1 wt % Al 2 O 3 ) augites to Wo 46 (∼6 wt% Al 2 O 3 ) that can be considered as xenocrysts inherited from more mafic compositions (Pichavant et al., 2018). Also, SG11A_dark has a strongly higher microlite content than expected from equilibrium crystallization ( Figure 3C). Overall, these observations argue in favor of mostly inherited microlites that crystallized from a more mafic melt intruding the reservoir, as described by Martel et al. (2006). Considering that SG6, O1215Eb, and SG11A_white pumice did not crystallize microlites during ascent, i.e., all microlites are inherited, suggests about 5-10 vol% inherited plagioclase microlites (Φ plag values in Table 1) for the magmas coming from the coldest part of the reservoir.
The lower sub-Plinian SG5A pumice and the Strombolian SG3 dark scoria mostly crystallized plagioclase microlites at P of ∼125-200 MPa, i.e., similar to the pre-eruptive P, and some rare microlites at P of 50-75 MPa (Figure 7). The P deduced from the residual melt compositions are similar to the preeruptive P when calculated starting with a 74 wt% SiO 2 melt (SG5A) and ∼100-130 MPa when calculated with a 71 wt% SiO 2 melt (SG5A and SG3). Therefore, the amount of plagioclase microlites grown at 50-75 MPa was likely too low to drive the  Table 1) versus the vesicularity calculated after Eq. 1 in Supplementary Material (Jaupart and Allègre, 1991) for a closedsystem degassing for which the gas bubbles stay in the magma. The 1:1 line (enlarged to 10% uncertainty) delimits a closed-system degassing, for which measured and calculated vesicularities compare, to an open-system degassing regime, for which bubbles are lost (outgassing) so that the measured vesicularities are lower than calculated for a closed-system degassing. The error bars give the statistical errors on the vesicularity measurements (X-axis) and the variability of the calculations resulting from the range of residual glass H 2 O contents (Y-axis), as given in Table 1.
Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 residual glass compositions to those expected at 50-70 MPa. This may be illustrated by the results of the phase-equilibrium experiments, showing that increasing the plagioclase content by 5 vol% from 75 to 50 MPa (Supplementary Material) does not increase the glass SiO 2 contents, which is capped at 79-80 wt % SiO 2 (Supplementary Material). The sub-Plinian SG5A sample contains orthopyroxene and clinopyroxene microlites that are compositionally less evolved than expected from decompression-induced crystallization at p < 175 MPa, either from 71 or 74 wt% SiO 2 melts (Figure 4). The Strombolian SG3 dark scoria contains orthopyroxene, clinopyroxene, and titanomagnetite microlites, whose compositions are on both sides, either mafic or felsic, of the pre-eruptive compositions ( Figure 4). Additionally, the microlite contents are higher than expected from equilibrium crystallization starting from a 71 wt% SiO 2 melt ( Figure 3C), suggesting a significant amount of inherited microlites. Therefore, most of the microlites in these two samples are inherited from mixing processes at reservoir level, with only a very few fraction crystallizing during magma ascent. Neglecting ascent-crystallization in the cases of the Strombolian SG3 scoria, the upper sub-Plinian SG11A_dark pumice, and the lower sub-Plinian SG5A pumice, may also suggest about 20-27 vol% inherited microlites from the Φ plag values in Table 1.

Magma Decompression Rate
Bubble number density has been demonstrated to correlate with decompression rate in rhyolitic melts (Toramaru, 2006). Yet, the AD 1530 magmas obviously underwent too much bubble coalescence (e.g., Figure 2) to make it possible to restore the original bubble number densities from the pumice samples. Alternatively, the decompression rate can be estimated from the microlite number density, in the case of a single event of crystal nucleation driven by melt H 2 O exsolution (Toramaru et al., 2008). Among the samples that nucleated plagioclase microlites during ascent, the presence of inherited microlites (predominant in SG5A), around which a part of the decompression-induced crystallization occurred (zoned microlites), means that ascentcrystallization partly occurred as growth around preexisting crystals instead of nucleation of new crystals. This disagrees with the requirements of the decompression-rate meter of Toramaru et al. (2008), implying that the calculated data must be treated with caution. The calculations are reported in Figure 9, and converted below into ascent durations assuming a magma travel from the ∼7 km deep reservoir to subsurface. The phreatomagmatic SG7B pumice could have ascended at rates of ∼400 Pa/s (ascent velocities of ∼0.002 m/s; duration of ∼30 days), the lower sub-Plinian SG5A pumice at rates from 40 to 200 Pa/s (velocities between ∼0.003 and 0.015 m/s and duration from 30 to 2 days, calculated for 74 and 71 wt% SiO 2 melts, respectively), and the Strombolian SG3 dark scoria at rates of ∼300 Pa/s (velocity of ∼0.015 m/s; duration of ∼2 days). The decompression rate of 900 Pa/s (velocity of ∼0.05 m/s; duration of 1.5 days) calculated for the 1227A dome sample is unrealistically fast, likely biased by continuous microlite nucleation upon ascent and by nonisothermal crystallization (late microlites formed upon cooling), so that it will not be further considered.
In comparison with the timescales of degassing deduced from the viscosity of a vesiculated magma (from Martel and Iacono-Marziano, 2015), the timescales of magma ascent calculated from the microlite content (from Toramaru et al., 2008) are about the same order of magnitude for the Strombolian SG3 dark scoria and the sub-Plinian SG5A pumice when calculated using a 71 wt% SiO 2 melt (∼2 days versus maximum 0.5-4 days) and much longer for the phreatomagmatic SG7B pumice and the sub-Plinian SG5A pumice when calculated using a 74 wt% SiO 2 melt (∼30 days versus maximum 0.5-4 days). Concerning the sub-Plinian SG5A pumice, this could argue for a starting interstitial melt closer to 71 than 74 wt% SiO 2 , i.e., tapping the main reservoir. As for the phreatomagmatic SG7B pumice, we conclude that the microlite number density could actually reflect a magma ascent duration of ∼30 days, whereas the shorter duration of 0.5-4 days attributed to degassing could reflect outgassing difficulties, as explained below.

The AD 1530 Eruptive Sequence
The AD 1530 eruptive sequence is revisited below using temporal constraints from microlite crystallization and degassing, as illustrated in Figure 10.

The Basal Phreatomagmatic Explosions
The sequence started with phreatic and phreatomagmatic explosions likely triggered by the shallow emplacement of a magma intrusion that overheated the surrounding hydrothermal system. This early magma was probably the FIGURE 9 | Decompression rates calculated from the plagioclase number density using the formulation of Toramaru et al. (2008;their Eq. 10). Plain and open symbols are calculation from a 71 and 74 wt% SiO 2 melt, respectively. Volume number density are calculated as N A 2/3 following Armienti (2008), initial H 2 O content is 5.4 wt%, and decompression rate is converted into ascent velocity using a 1.7 g/cm 3 bulk density for the vesiculated magma. Note that the decompression rate calculated for the dome magma is unrealistically high, due to multiple events of microlite nucleation (see text).
"SG7B magma", coming from the coldest parts of the reservoir (∼825°C) and pushed upwards by a more mafic magma at depth, as evidenced by the Ca-rich cores of the plagioclase microlites (∼An 75 ; Figure 4). During ascent, the magma continuously crystallized microlites of plagioclase, pyroxene, clinopyroxene, apatite, and quartz microlites under near-equilibrium conditions ( Figure 3C). Near-equilibrium crystallization of the plagioclase microlites is confirmed by their low number density and rather large (20-100 µm) euhedral shapes, which reflect a growthdominated crystallization under small undercoolings of < ∼40°C (Hammer and Rutherford, 2002;Mollard et al., 2012). The microlite content suggests magma ascent at low velocity of ∼0.002 m/s (Figure 9; from Toramaru et al., 2008), i.e., in about 1 month from reservoir to quench P of ∼50 MPa (decompression rate of ∼50 Pa/s). Interestingly, the magma did not outgas (Figure 8), which suggest a low permeability of the conduit walls or the presence of a magma plug, including the possibility of a small dome or cryptodome, that prevented bubble loss.

The Sub-Plinian Phase
The flank collapse that followed the phreatomagmatic explosions likely depressurized the underlying magma and triggered the sub-Plinian eruption. The white homogeneous pumice (SG6A, SG5A) and the white layers of the banded pumice (O1215Eb, SG11A_white) share close bulk-rock compositions (59-62 wt% SiO 2 ; Figure 5) with the phreatomagmatic SG7B pumice, so that they likely tapped the coldest parts of the reservoir. The magma mixing features in the banded pumice (e.g., upper sub-Plinian SG11A_dark) suggest that the main reservoir (at 875°C) was progressively tapped. Most of these sub-Plinian pumice (i.e., SG5A, SG6A, O1215Eb) did not crystallize upon ascent. This means that the decompression rates of ∼50 Pa/s deduced from the microlite number density in the SG5A sample using the FIGURE 10 | Magma ascent conditions of the AD 1530 eruption of la Soufrière of Guadeloupe. The eruption was likely triggered by the injection of a basaltic magma (green arrow) into the residing andesitic magma (red and blue shaded parts of the reservoir). The first magma that emplaced in AD 1530 likely tapped the coldest parts of the main reservoir (in blue; 825°C, interstitial melt with 74 wt% SiO 2 and 5.5 wt% H 2 O) and ascended during about 1 month up to ∼2 km deep. During ascent, the magma had time to crystallize microlites but did not outgas, likely due to sealed conduit walls. Heating and pressurizing the surrounding hydrothermal system triggered a phreatomagmatic explosion, that destabilized the summit of the volcano, generating a flank collapse of the lave dome. This triggered a sub-Plinian event that progressively tapped the main reservoir (in red; 875°C, interstitial melt with 71 wt% SiO 2 and 5.5 wt% H 2 O), emitting homogeneous and banded pumice. The magma ascent durations are estimated to 0.5-4 days. The following short Strombolian event emitted dark scoria that ascended in 1-2 days before fragmentation at ∼1 km deep, and ended by the growth of the present lava dome whose magma ascended in more than 1 month.
formulation of Toramaru et al. (2008) (Figure 9) may not be valid, because these are microlites inherited from reservoir processes. Decompression durations preventing crystallization are of the order of <1 day, according to decompression experiments performed on compositionally close melts (75-76 wt% SiO 2 ) and T (850-860°C; Martel and Schmidt, 2003;Martel, 2012), but should be > ∼10 h for the bubbles to reach complete expansion (Martel and Schmidt, 2003), as required for a closed-system equilibrium degassing (Figure 8).
Such an estimated duration of ∼0.5-1 day converts into a decompression rate of 1,750-3,500 Pa/s (considering a decompression of 150 MPa), and agrees with the maximum ascent duration of 0.5-4 days deduced from the degassing timescale (Martel and Iacono-Marziano, 2015). The "SG5A and SG11A magmas" started to outgas (Figure 8), suggesting local permeability of the conduit walls, and slower decompression rates of ∼270-1700 Pa/s (1-4 days, from Martel and Iacono-Marziano 2015).

The Strombolian Phase
The eruption went on by extruding a basaltic andesitic magma (56-59 wt% SiO 2 ) as lava fountaining, from which the SG3 dark scoria was collected. At pre-eruptive conditions of 875°C and 175 MPa, these magmas have melt compositions (in equilibrium with the phenocrysts) of ∼71 wt% SiO 2 , maybe somewhat less.

The Terminal Dome Growth
The eruption ended with the growth of the present lava dome. The "1227A" dome sample represents the most mafic magma (∼56 wt% SiO 2 and 4.2 wt% MgO) of the eruptive sequence, with pre-eruptive residual melt compositions of ≤71 wt% SiO 2 at 875°C and 175 MPa. The ascent duration of the "1227A" magma is not well constrained: it was long enough to allow steady-state microlite crystallization ( Figure 7) and full outgassing (Figure 8), the latter argument suggesting ascent duration of more than ∼40 days (from Martel and Iacono-Marziano, 2015; decompression rate <50 Pa/s).

Comparison With Sub-Plinian Eruptions From Other Volcanic Systems
Explosive eruptions, including Plinian, sub-Plinian, and Vulcanian events, and involving dacitic to rhyolitic magmas, commonly show syneruptive decompression rates from 100 to ∼3,000 Pa/s (ascent rates of 0.01-0.1 m/s), based on petrological, experimental, and monitoring data (Cassidy et al., 2018). These rates of decompression agree with the decompression rates of 270-3,500 Pa/s that we propose for the AD 1530 sub-Plinian event of La Soufrière. Nevertheless, refinements can be proposed by restricting the comparison to sub-Plinian events for which magma ascent rates are known, such as the 1996 explosion of Soufrière Hills, Montserrat, and the 2010 eruption of Merapi, Indonesia. The September 17th, 1996, short-lived sub-Plinian event (VEI 3) of Soufrière Hills volcano (Robertson et al., 1998) produced homogeneous pumices. From the absence of reaction rims around the amphibole phenocrysts in the pumice samples, Devine et al. (1998) estimated magma ascent durations >0.01 m/s, i.e., < 4 days. These ascent rates are comparable to the ascent durations <4 days proposed for the AD 1530 sub-Plinian pumice samples of La Soufrière.
The 2010 eruption of Merapi volcano started on October 26th by a phreatomagmatic explosion, followed by 7 days of dome growth and collapse, a sub-Plinian event (VEI 4) on November 5th, and terminated by a new dome growth (Surono et al., 2012;Komorowski et al., 2013). This eruptive sequence compares to the AD 1530 sequence of La Soufrière. The 2010 Merapi sub-Plinian event produced white pumice and grey scoria. The matrices of the white pumice have some similarities with those of the AD 1530 pumice of Soufrière: they show limited outgassing and tabular feldspar microlites with number densities of 10 9.8 m −2 and contents of 16-32 vol% (Preece et al., 2016) (Figure 3). Although the 2010 Merapi magma is more mafic (bulk-rock of ∼55-56 wt% SiO 2 and residual glasses of ∼68 wt% SiO 2 for the white pumice; Preece et al., 2014) and was stored at higher temperature (925-950°C; Erdmann et al., 2016) than for the AD 1530 La Soufrière magma, such textural similarities indicate comparable ascent conditions. The Merapi sub-Plinian event itself lasted 2 h and the preceding dome emplaced during 7 days, thus bracketing the ascent timescales of Merapi magma between 2 h for the pumice to 7 days for the denser clasts. These rates are in the same order of magnitude than those proposed for the AD 1530 sub-Plinian event (0.5-4 days; Figure 10).

Volcanic Hazard in Guadeloupe
Volcanic eruptions often start with a phreatic event that uncorks the conduit, and whether a magmatic eruption follows, it is still not predictable. The trigger of a magmatic eruption, and the possible transitions between effusive and explosive events, relies on interrelated and competing magmatic properties and processes, such as magma viscosity, degassing, conduit geometry and permeability, which ultimately control magma ascent rate (Cassidy et al., 2018). However, as only slight variations in the involved parameters may force a change in the eruptive style, we need to understand the factors controlling the conduit processes.
The study of the phreatomagmatic SG7B pumice from the AD 1530 deposit suggests that the early magma emplaced in a month or so at depth of ∼2 km, likely with little surface outgassing since bubbles did not escape the melt at depth. Such an intrusion could have been similar to that inferred during the 1976-77 phreatic eruption (Feuillard et al., 1983;Villemant et al., 2014). At the present day, a magma intrusion at depth could be detected via geophysical (seismic and ground deformation) methods, depending on the magma volume involved. Such geophysical signals of magma intrusion combined with an absence of surface degassing signals may Frontiers in Earth Science | www.frontiersin.org November 2021 | Volume 9 | Article 686342 be interpreted as gas pressurization at depth and potential precursors for an impending explosion. In AD 1530, the phreatomagmatic explosion likely destabilized the summit area of the volcano, generating a flank collapse that depressurized the underlying magma, leading to a sub-Plinian eruption. Such a scenario is nowadays all the more probable considering that the highly fractured and hydrothermalized summit lava dome is unstable (Rosas-Carbajal et al., 2016). In AD 1530, the collapse of the eruptive column during the sub-Plinian phase generated pyroclastic density currents channeled through two or three valleys, but the growth of the lava dome was not accompanied by destruction into pyroclastic density currents, as frequently occurred during the 1995-2010 Soufrière Hills eruption (Sparks and Young, 2002). Yet, superficial laterally directed explosions can also occur, particularly at the beginning of lava dome growth (e.g., historical eruptions at Montagne Pelée; Boudon and Balcone-Boissard, 2021). Such lava dome-related eruptions, which commonly depend on magma ascent rate (Martel, 2012) and lava dome size (Taisne and Jaupart, 2008), have to be considered as likely hazards in future eruptive scenarios of La Soufrière. Therefore, geophysical and geochemical signals of shallow intrusion of new magma are of prime importance at La Soufrière volcano, because the produced heat could enhance the alteration of the lava-dome rocks eventually generating cascade processes including lava-dome and/or flank collapse and magmatic eruptions. Such a dramatic eruptive scenario threatens up to 70,000 people in southern Basse-Terre island of the region.

CONCLUSION
The textures and compositions of the matrices from samples emitted from various eruptive styles during the AD 1530 eruption of La Soufrière of Guadeloupe were compared to experimental samples from phase equilibria and decompression experiments from the literature, in order to provide timescales of syn-eruptive crystallization and degassing. We highlighted that: • The early magma ascended from the reservoir to about 2 km deep in about 1 month. This intrusion heated the surrounding hydrothermal system, which likely lead to a flank collapse a following phreatomagmatic eruption. • Dome removal decompressed the underlying magma generating magma ascent from the reservoir within 0.5-4 days, leading to a sub-Plinian event emitting homogeneous and banded pumice. • The sub-Plinian event ended by a Strombolian phase, during which magma ascended in 1-2 days, followed by a terminal dome growth involving magma ascent during more than 1 month.
• Understanding the conduit processes that lead to the AD 1530 eruption is crucial for the hazard evaluation of a future volcanic crisis at La Soufriere, especially in a context of a present altered and instable lava dome that could collapse and trigger a major explosive event.

DATA AVAILABILITY STATEMENT
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

AUTHOR CONTRIBUTIONS
CM, MP, HB-B, and GB contributed to the conception of the work. CM and MP acquired the data. CM organized the database and wrote the first draft of the manuscript. All authors contributed to manuscript revision, read, and approved the submitted version.