Muon imaging has recently emerged as a powerful method to complement standard geophysical tools in the understanding of the Earth's subsurface. Muon measurements yield a “radiography” of the average ...density along the muon path, allowing to image large volumes of a geological body from a single observation point. Here we jointly invert muon data from three simultaneous telescope acquisitions together with gravity data to estimate the three‐dimensional density structure of the La Soufrière de Guadeloupe lava dome. Our unique data set allows us to achieve an unprecedented spatial resolution with this novel technique. The retrieved density model reveals an extensive, low‐density anomaly where the most active part of the volcanic hydrothermal system is located, supporting previous studies that indicate this region as the most likely to be involved in a partial edifice collapse.
Key Points
We simultaneously image the La Soufrière de Guadeloupe lava dome with three muon telescopes aiming at the volcano from different locations
We jointly invert the multitelescope muon data with gravity data to retrieve the three‐dimensional density distribution in the dome
We find a low‐density, mechanically weak region where the hydrothermal system is most active, suggestive of an edifice collapse hazard
Catastrophic collapses of the flanks of stratovolcanoes constitute a major hazard threatening numerous lives in many countries. Although many such collapses occurred following the ascent of magma to ...the surface, many are not associated with magmatic reawakening but are triggered by a combination of forcing agents such as pore-fluid pressurization and/or mechanical weakening of the volcanic edifice often located above a low-strength detachment plane. The volume of altered rock available for collapse, the dynamics of the hydrothermal fluid reservoir and the geometry of incipient collapse failure planes are key parameters for edifice stability analysis and modelling that remain essentially hidden to current volcano monitoring techniques. Here we derive a high-resolution, three-dimensional electrical conductivity model of the La Soufrière de Guadeloupe volcano from extensive electrical tomography data. We identify several highly conductive regions in the lava dome that are associated to fluid saturated host-rock and preferential flow of highly acid hot fluids within the dome. We interpret this model together with the existing wealth of geological and geochemical data on the volcano to demonstrate the influence of the hydrothermal system dynamics on the hazards associated to collapse-prone altered volcanic edifices.
SUMMARY
Density muon radiography is a new method to determine the average density of geological bodies by measuring the attenuation produced by rocks on the flux of cosmic muons. We present such ...density radiographies obtained for the Soufrière of Guadeloupe lava dome, both in the north–south and east–west planes. These radiographies reveal the highly heterogeneous density structure of the volcano, with low‐density regions corresponding to recognized hydrothermally altered areas. The main structures observed in the density radiographies correlate with anomalies in electrical resistivity cross‐sections and a density model obtained from gravity data.
Imaging geological structures through cosmic muon radiography is a newly developed technique which shows a great potential in volcanology. Here we demonstrate that muon radiography permits to detect ...and characterize mass movements in shallow hydrothermal systems of low-energy active volcanoes like the La Soufrière lava dome. We present an experiment conducted on this volcano during the Summer 2014 and bring evidence that very important density changes occurred in three domains of the lava dome. Depending on their position and on the medium porosity the volumes of these domains vary from 1 × 10(6) m(3) to 7 × 10(6) m(3). However, the total mass budget remains approximately constant : two domains show a mass loss (Δm∈ -0.8;-0.4 × 10(9) kg) and the third one a mass gain (Δm∈ 1.5; 2.5 × 10(9) kg). We attribute the negative mass changes to the formation of steam in shallow hydrothermal reservoir previously partly filled with liquid water. This coincides with the emergence of new fumaroles on top of the volcano. The positive mass change is synchronized with the negative mass changes indicating that liquid water probably flowed from the two reservoirs invaded by steam toward the third reservoir.
Polar ice core records attest to a colossal volcanic eruption that took place ca. A.D. 1257 or 1258, most probably in the tropics. Estimates based on sulfate deposition in these records suggest that ...it yielded the largest volcanic sulfur release to the stratosphere of the past 7,000 y. Tree rings, medieval chronicles, and computational models corroborate the expected worldwide atmospheric and climatic effects of this eruption. However, until now there has been no convincing candidate for the mid-13th century “mystery eruption.” Drawing upon compelling evidence from stratigraphic and geomorphic data, physical volcanology, radiocarbon dating, tephra geochemistry, and chronicles, we argue the source of this long-sought eruption is the Samalas volcano, adjacent to Mount Rinjani on Lombok Island, Indonesia. At least 40 km ³ (dense-rock equivalent) of tephra were deposited and the eruption column reached an altitude of up to 43 km. Three principal pumice fallout deposits mantle the region and thick pyroclastic flow deposits are found at the coast, 25 km from source. With an estimated magnitude of 7, this event ranks among the largest Holocene explosive eruptions. Radiocarbon dates on charcoal are consistent with a mid-13th century eruption. In addition, glass geochemistry of the associated pumice deposits matches that of shards found in both Arctic and Antarctic ice cores, providing compelling evidence to link the prominent A.D. 1258/1259 ice core sulfate spike to Samalas. We further constrain the timing of the mystery eruption based on tephra dispersal and historical records, suggesting it occurred between May and October A.D. 1257.
We have used a three-dimensional, non-equilibrium multiphase flow numerical model to simulate subplinian eruption scenarios at La Soufrière de Guadeloupe (Lesser Antilles, France). Initial and ...boundary conditions for computer simulations were set on the basis of independent estimates of eruption source parameters (i.e. mass eruption rate, volatile content, temperature, grain size distribution) from a field reconstruction of the 1530 CE subplinian eruption. This event is here taken as a reference scenario for hazard assessment at La Soufrière de Guadeloupe. A parametric study on eruption source parameters allowed us to quantify their influence on the simulated dynamics and, in particular, the increase of the percentage of column collapse and pyroclastic density current (PDC) intensity, at constant mass eruption rate, with variable vent diameter. Numerical results enabled us to quantify the effects of the proximal morphology on distributing the collapsing mass around the volcano and into deep and long valleys and to estimate the areas invaded by PDCs, their associated temperature and dynamic pressure. Significant impact (temperature > 300 °C and dynamic pressure > 1 kPa) in the inhabited region around the volcano is expected for fully collapsing conditions and mass eruption rates > 2 × 10
7
kg/s. We thus combine this spatial distribution of temperature and dynamic pressure with an objective consideration of model-related uncertainty to produce preliminary PDC hazard maps for the reference scenario. In such a representation, we identify three areas of varying degree of susceptibility to invasion by PDCs—very likely to be invaded (and highly impacted), susceptible to invasion (and moderately impacted), and unlikely to be invaded (or marginally impacted). The study also raises some key questions about the use of deterministic scenario simulations for hazard assessment, where probability distributions and uncertainties are difficult to estimate. Use of high-performance computing techniques will in part allow us to overcome such difficulties, but the problem remains open in a scientific context where validation of numerical models is still, necessarily, an incomplete and ongoing process. Nevertheless, our findings provide an important contribution to the quantitative assessment of volcanic hazard and risk at La Soufrière de Guadeloupe particularly in the context of the current unrest of the volcano and the need to prepare for a possible future reawakening of the volcano that could culminate in a magmatic explosive eruption.
In May 2018, a seismically quiet region of the Indian Ocean awoke. More than 130 magnitude 4+ earthquakes were recorded in the first month, including a MW 5.9 event on May 15th, 2018. This seismic ...activity was later identified as being related to an exceptional underwater volcanic eruption offshore Mayotte island, which had emitted more than 6.5 km3 of lava by the time of writing. To better constrain the geodynamic processes responsible for the seismic and volcanic activity, a new network of ocean-bottom seismometers and land stations has been deployed around the seismically active region since February 2019. We present here an improved 1D velocity model for the active area and relocations of manually-picked earthquakes using this new model. The best-constrained events image detailed structures within two clusters of seismic activity east of Mayotte. The westernmost, proximal cluster, close to Mayotte's Petite-Terre island, has a “donut” shape horizontally and an “hourglass” shape in depth. The events distribution suggests the presence of a magma reservoir at around 27 km depth, with earthquakes focused along its sides, and a collapsing system underneath, related to the drainage of another, deeper magma storage zone. The distal cluster, focused 30–50 km offshore of Petite-Terre island, highlights the propagation of a dike between 45 and 25 km depth, aligned towards the new volcanic activity on the seafloor. We interpret this cluster as the fluid pathway towards the new volcano and nearby active seafloor lava fields. The improved velocity model also permits more robust daily monitoring of the seismicity using land stations, allowing local authorities to better assess seismic and volcanic hazards and to communicate them to the island's population.
•We develop a new local 1D velocity model for Mayotte.•The locations of the detected earthquakes are significantly improved.•The earthquakes' distribution images the active structures of the magmatic system.•Daily monitoring is now more robust, helping hazard assessments and communication.
La Soufrière de Guadeloupe volcano is characterized by seismo-volcanic activity dominantly linked to an active hydrothermal system. This microseismicity is shallow, mainly triggered in swarms and ...characterized by repeating earthquakes. Four recurrent families of Volcano-Tectonic (VT) repeaters have been identified and the main repeater accounts for 80% of detections. Stacking repeating seismic waveforms produces MASTER events with a high Signal-to-Noise Ratio (SNR) and thanks to the deployment of a temporary seismic nodal array, their absolute locations can be better constrained. Because we observe systematic positive residuals of P-wave arrival times at almost all stations, we investigate potential bias in the velocity model in the shallow part of the dome. By increasing the P-wave velocity and decreasing the VP/VS ratio from 1.8 to 1.69, we reduce these residuals and improve the local velocity model. We use this new velocity model to locate each VT event relatively to is own MASTER hypocenter. This procedure offers a new image of the hydrothermal seismic activity that we locate under the acid lake of the summit Tarissan crater, along a sub-vertical conduit. We also define a linear relationship between the logarithm of the peak amplitude of seismic events and their duration to obtain a pseudo local magnitude. The approach enables to accurately and automatically estimate the magnitude during the event detection. We show that after the occurrence in April 2018 of the largest earthquake (Mlv 4.1) since the last phreatic eruption in 1976–1977, the swarm frequency and amplitude increased significantly, and a new family of VT repeater emerges, accounting for up to 14% of detections. We show a quiescence in the main repeaters activity in the month following the April 2018 earthquake, suggesting a significant stress release within the volcanic system, which is likely related to dynamic stress changes rather than static stress variations caused by the earthquake. Moreover, we detect the emergence of a new family of repeaters 3 months after the event, located ∼100 m above the main repeater family, that could be related to a fluid pressure increase induced by the regional hydrological cycle, superimposed to the internal forcing driven by the hydrothermal system. Finally, using a statistical approach, we highlight seasonal periodicities in the number of events and the released seismic moment at La Soufrière de Guadeloupe, with a dominant peak of seismic activity in October–November and a second, less significant, peak of in April.
•Increased location accuracy with the implementation of a relative location procedure•Microseismicity extends towards the surface following the April 2018 earthquake.•April 2018 earthquake resulted in increased shallow seismicity and energy release.•Evidence of seasonality in the shallow seismo-volcanic activity
Signals of volcanic unrest have been used successfully to provide insights into the timing, magnitude, intensity, and style of future eruptions. However, in order to provide context for the ...subsequent activity, analysis of past eruptions is required. This provides useful information in order to understand processes of magma genesis, storage, evolution and ascent which lead to the onset of eruptions. Here, we examine basaltic-andesitic to andesitic deposits from La Soufrière de Guadeloupe Holocene eruptions, covering a range of explosive eruption styles, ages and magnitudes. Our work is timely given unrest at this system has increased over the last 25 years, with a potential eruption capable of directly impacting up to 80,000 people in Southern Basse-Terre and potentially thousands more indirectly on a regional scale. We report on the geochemistry of pre-eruptive magmas using detailed analyses of glass (melt inclusions and groundmass glass) from four Holocene explosive eruptions: 1657 Cal. CE (Vulcanian, VEI 2), 1530 Cal. CE (sub-Plinian, VEI 3), 1010 Cal. CE (Plinian, VEI 4), and 5680 Cal. BCE (Plinian, VEI 4). Major element concentrations vs SiO2 in whole rock (WR), groundmass glass (GM) and melt inclusions (MI) show a strong linear trend. MIs reveal a relatively homogenous melt composition from the first to the most recent eruptions, ranging from 63.6–78.7 wt% SiO2. Volatiles, including H2O (2.3–4.4 wt%), CO2 (35–866 ppm) and sulphur (30–202 ppm), are also consistent across the various eruptions. The major element and volatile compositional homogeneity across the eruptions indicates that composition and volatiles do not have a direct control on eruption explosivity at this system. Instead, we find differences in ascent rate, groundmass glass viscosity and microlite volume percentage, indicating that explosive eruptive style at La Soufrière is controlled by a combination of ascent rate and top-down controls affecting rock strength, stress distribution and the development of fluid overpressure. Rapid ascent in the absence of top-down controls (processes with a cause external to the magma but affecting the plumbing system) will result in explosive eruptions driven from the bottom-up (internal to magma dynamic response with varying pressure and temperature, e.g., 1010 Cal. CE in the case of very rapid ascent or 1657 Cal. CE in the case of rapid ascent). However, we also highlight the importance of top-down controls, such as conduit sealing which can promote the onset of explosive eruptions, even in the case of slow magma ascent (e.g., 5680 Cal. BCE). External effects (including ingress of water and rapid edifice unloading) can also favour explosive eruptions with flank collapses involved in some scenarios (e.g., 1530 Cal. CE). The multiple controls on explosive eruption style make this system more hazardous and complex to model and monitor. In order to improve early-warning system efficiency, forecast models, eruption scenario crisis response and long-term risk reduction planning, we stress that internal processes such as fracture and host-rock sealing (fluid pore pressure) as well as external processes such as water moving into the system and the mechanical stability of the edifice should be monitored and modelled closely.
•La Soufrière de Guadeloupe explosive eruption styles vary from Vulcanian to Plinian (VEI 2 to VEI 4).•The mush storage zone and the magma composition have remained stable across the eruptions studied•Differences in magma ascent rate result in explosive eruption style variations.•Top-down process such as conduit sealing can also lead to explosive eruptions.•Multiple controls on explosive eruption style make this system hazardous.
Volcanoes are unstable structures that deform laterally and frequently experience mass wasting events. Hydrothermal alteration is often invoked as a mechanism that contributes significantly to ...volcano instability. We present a study that combines laboratory experiments, geophysical data, and large‐scale numerical modeling to better understand the influence of alteration on volcano stability, using La Soufrière de Guadeloupe (Eastern Caribbean) as a case study. Laboratory experiments on variably altered (advanced argillic alteration) blocks show that uniaxial compressive strength, Young's modulus, and cohesion decrease as a function of increasing alteration, but that the internal friction angle does not change systematically. Simplified volcano cross sections were prepared (a homogenous volcano, a volcano containing the alteration zone identified by a recent electrical survey, and a volcano with an artificially enlarged area of alteration) and mechanical properties were assigned to zones corresponding to unaltered and altered rock. Numerical modeling performed on these cross sections, using a hydro‐thermo‐mechanical modeling code, show (a) the importance of using upscaled values in large‐scale models and (b) that alteration significantly increases volcano deformation and collapse volume. Finally, we combined published muon tomography data with our laboratory data to create a 3D strength map, exposing a low‐strength zone beneath the southern flank of the volcano coincident with the hydrothermal system. We conclude that hydrothermal alteration decreases volcano stability and thus expedites volcano spreading and increases the likelihood of mass wasting events and associated volcanic hazards. Hydrothermal alteration, and its evolution, should therefore be monitored at active volcanoes worldwide.
Plain Language Summary
The rocks forming a volcanic edifice can be altered by circulating hydrothermal fluids. This alteration can influence the physical and mechanical properties of these rocks, which could jeopardize volcano stability. The stability of a volcanic edifice is an important consideration in volcanic hazards and risk assessments due to the potentially dire consequences of partial volcanic flank collapse. Using a combination of experimental data, geophysical data, and modeling, and La Soufrière de Guadeloupe (Eastern Caribbean, France) as a case study, we find that hydrothermal alteration decreases volcano stability and thus promotes volcano instability and associated volcanic hazards. As a result, we conclude that hydrothermal alteration, and its evolution, should be monitored at active volcanoes worldwide.
Key Points
Laboratory experiments show that hydrothermal alteration reduces the strength of volcanic rock from La Soufrière
Numerical modeling shows that hydrothermal alteration significantly increases volcano deformation and collapse volume
We provide a 3D strength map of La Soufrière that exposes a low‐strength zone coincident with the hydrothermal system