Search

International audience; Surface deformation accompanying dike intrusions is dominated by uplift and horizontal motion directly related to the intrusions. In some cases, it includes subsidence due to associated magma reservoir deflation. When reservoir deflation is large enough, it can form, or reactivate pre-existing, caldera ring-faults. Ring-fault reactivation, however, is rarely observed during moderate-sized eruptions. On February 21st, 2015 at Ambrym volcano in Vanuatu, a basaltic dike intrusion produced more than 1 meter of co-eruptive uplift, as measured by InSAR, SAR correlation, and Multiple Aperture Interferometry (MAI). Here we show that an average of ∼40 cm of slip occurred on a normal caldera ring-fault during this moderate-sized (VEI < 3) event, which intruded a volume of ∼24 million cubic meters and erupted ∼9.3 million cubic meters of lava (DRE). Using the 3D Mixed Boundary Element Method, we explore the stress change imposed by the opening dike and the depressurizing reservoir on a passive, frictionless fault. Normal fault slip is promoted when stress is transferred from a depressurizing reservoir beneath one of Ambrym’s main craters. After estimating magma compressibility, we provide an upper-bound on the critical fraction (f = 7%) of magma extracted from the reservoir to trigger fault slip. We infer that broad basaltic calderas may form in part by hundreds of subsidence episodes no greater than a few meters, as a result of magma extraction from the reservoir during moderate- sized dike intrusions.

International audience; In December 2018, an unusually large intra- and extra-caldera eruption took place at Ambrym volcano (Vanuatu). The eruption drained the volcano’s five active lava lakes and terminated, at least momentarily, the surface activity that had been ongoing for decades to hundreds of years, sustaining the largest recorded persistent degassing on the planet. Here, we investigate the mechanisms and dynamics of this major eruption. We use major elements and volatiles in olivine and clinopyroxene-hosted melt inclusions, embayments, crystals and matrix glasses together with clinopyroxene geobarometry as well as olivine and clinopyroxene geothermometry and diffusion modelling in crystals and embayments to reconstruct the chronology and timing of the subsurface processes that accompanied the eruption. We find that the eruption began with the meeting, mingling and limited chemical mixing of mostly two magma bodies occupying similar vertical but different horizontal locations in the crust, one corresponding to the main plumbing system at Ambrym that fed the lava lakes and the other corresponding to an older, previously cutoff and more chemically evolved branch of the plumbing system. Within the primitive magma, two texturally distinct components—one microlite rich and one microlite poor—can further be identified. The 2018 eruption hence provides a detailed image of Ambrym’s complex plumbing system. Our diffusion timescales and geobarometric estimates coincide closely with geophysical observations. They point to a reconnection of the evolved magmatic branch with the main system occurring less than 10 h prior to the intra-caldera eruption and a period of 2 days for the subsequent > 30-km lateral magma transport along a deeper dike prior to submarine eruption just off the SE coast of the island with the more primitive magma reaching first followed by mingled magma containing both compositions. Magma ascent rates are estimated at 95 ± 24 m/s in the last ~ 2.5 km of ascent during the intra-caldera eruption and at 80 ± 6 m/s in the last ~ 4 km of ascent during the submarine eruption. Comparison with other lava lake draining eruptions reveals striking similarities both in terms of precursory activity, with lake level rising prior to the eruption in all cases, and in terms of plumbing system organization with the presence of peripheral magma pockets, isolated from the main magmatic system but that can be mobilized and erupted when met by dikes propagating laterally from the main system.

Emission of volcanic gas is thought to be the dominant process by which volatiles transit from the deep earth to the atmosphere. Volcanic gas emissions, remain poorly constrained, and volcanoes of Peru are entirely absent from the current global dataset. In Peru, Sabancaya and Ubinas volcanoes are by far the largest sources of volcanic gas. Here, we report the first measurements of the compositions and fluxes of volcanic gases emitted from these volcanoes. The measurements were acquired in November 2015. We determined an average SO 2 flux of 15.3 ± 2.3 kg s − 1 (1325-ton day − 1 ) at Sabancaya and of 11.4 ± 3.9 kg s − 1 (988-ton day − 1 ) at Ubinas using scanning ultraviolet spectroscopy and dual UV camera systems. In-situ Multi-GAS analyses yield molar proportions of H 2 O, CO 2 , SO 2 , H 2 S and H 2 gases of 73, 15, 10 1.15 and 0.15 mol% at Sabancaya and of 96, 2.2, 1.2 and 0.05 mol% for H 2 O, CO 2 , SO 2 and H 2 S at Ubinas. Together, these data imply cumulative fluxes for both volcanoes of 282, 30, 27, 1.2 and 0.01 kg s − 1 of H 2 O, CO 2 , SO 2 , H 2 S and H 2 respectively. Sabancaya and Ubinas volcanoes together contribute about 60% of the total CO 2 emissions from the Central Volcanic zone, and dominate by far the total revised volatile budget of the entire Central Volcanic Zone of the Andes.

Eruptive activity shapes volcanic edifices. The formation of broad caldera depressions is often associated with major collapse events, emplacing conspicuous pyroclastic deposits. However, caldera subsidence may also proceed silently by magma withdrawal at depth, more difficult to detect. Ambrym, a basaltic volcanic island, hosts a 12-km wide caldera and several intensely-degassing lava lakes confined to intra-caldera cones. Using satellite remote sensing of deformation, gas emissions and thermal anomalies, combined with seismicity and ground observations, we show that in December 2018 an intra-caldera eruption at Ambrym preceded normal faulting with >2 m of associated uplift along the eastern rift zone and 2.5 m of caldera-wide subsidence. Deformation was caused by lateral migration of >0.4 cubic kilometers of magma into the rift zone, extinguishing the lava lakes, and feeding a submarine eruption in the rift edge. Recurring rifting episodes, favored by stress induced by the D’Entrecasteaux Ridge collision against the New Hebrides arc, lead to progressive subsidence of Ambrym’s caldera and concurrent draining of the lava lakes. Although counterintuitive, convergent margin systems can induce rift zone volcanism and subsequent caldera subsidence.

Forecasting the dispersal of volcanic clouds during an eruption is of primary importance, especially for ensuring aviation safety. As volcanic emissions are characterized by rapid variations of emission rate and height, the (generally) high level of uncertainty in the emission parameters represents a critical issue that limits the robustness of volcanic cloud dispersal forecasts. An inverse modeling scheme, combining satellite observations of the volcanic cloud with a regional chemistry-transport model, allows reconstructing this source term at high temporal resolution. We demonstrate here how a progressive assimilation of freshly acquired satellite observations, via such an inverse modeling procedure, allows for delivering robust sulfur dioxide (SO2) cloud dispersal forecasts during the eruption. This approach provides a computationally cheap estimate of the expected location and mass loading of volcanic clouds, including the identification of SO2-rich parts.

Depending on the magnitude of their eruptions, volcanoes impact the atmosphere at various temporal and spatial scales. The volcanic source remains a major unknown to rigorously assess these impacts. At the scale of an eruption, the limited knowledge of source parameters, including time variations of erupted mass flux and emission profile, currently represents the greatest issue that limits the reliability of volcanic cloud forecasts. Today, a growing number of satellite and remote sensing observations of distant plumes are becoming available, bringing indirect information on these source terms. Here, we develop an inverse modelling approach combining satellite observations of the volcanic plume with an Eulerian regional chemistry-transport model (CHIMERE) to characterise the volcanic SO2 emissions during an eruptive crisis. The May 2010 eruption of Eyjafjallajökull is a perfect case study to apply this method as the volcano emitted substantial amounts of SO2 during more than a month. We take advantage of the SO2 column amounts provided by a vast set of IASI (Infrared Atmospheric Sounding Interferometer) satellite images to reconstruct retrospectively the time series of the mid-tropospheric SO2 flux emitted by the volcano with a temporal resolution of ~2 h, spanning the period from 1 to 12 May 2010. We show that no a priori knowledge on the SO2 flux is required for this reconstruction. The initialisation of chemistry-transport modelling with this reconstructed source allows for reliable simulation of the evolution of the long-lived tropospheric SO2 cloud over thousands of kilometres. Heterogeneities within the plume, which mainly result from the temporal variability of the emissions, are correctly tracked over a timescale of a week. The robustness of our approach is also demonstrated by the broad similarities between the SO2 flux history determined by this study and the ash discharge behaviour estimated by other means during the phases of high explosive activity at Eyjafjallajökull in May 2010. Finally, we show how a sequential IASI data assimilation allows for a substantial improvement in the forecasts of the location and concentration of the plume compared to an approach assuming constant flux at the source. As the SO2 flux is an important indicator of the volcanic activity, this approach is also of interest to monitor poorly instrumented volcanoes from space.

In the 15 years following the last crisis of phreatic type in 1975–77, the hydrothermal and seismic activity of La Soufriere de Guadeloupe has displayed an episodic behavior superimposed on a long-term decay, as inferred from seismicity and degassing, monitored through fumarolic emissions and Cl content in thermal springs. In particular, the evolution of the Cl content can be explained by a succession of magmatic gas injections into the hydrothermal system. This episodic supply of gas may be associated to the thermodynamic evolution of a batch of magma intruded at shallow depth around 1976. Progressive magma cooling, crystallization, gas exsolution and pressure buildup lead to cyclic failure of host rock around the intrusion, and episodic release of gas, with a decreasing frequency as cooling proceeds. This physical model allows to reproduce the pattern of degassing activity of La Soufriere, including the decreasing frequency of Cl pulses with time. Our preferred forward model of Cl pulses involves a magma intrusion volume in the range 0.01– 0.52 km 3 and a tensile strength of the surrounding rock in the range 0.1–0.6 MPa. These results reinforce the hypothesis of a shallow magma intrusion triggering the 1975–77 crisis and provide an explanation for the degassing and seismic activities observed until 1992.

International audience; The measurement of SO2 flux from volcanoes is of major importance for monitoring and hazard assessment purposes, and for evaluation of the environmental impact of volcanic emissions. We propose here a novel technique for accurate and high time resolution estimations of the gas flux. We use two wide field of view UV spectrometers capable of collecting, instantaneously, light from thin parallel cross sections of the whole gas plume, obviating the need for either traversing, scanning or imaging. It enables tracking of inho-mogeneities in the gas cloud from which accurate evaluation of the plume velocity can be made by correlation analysis. The method has been successfully applied on Mt. Erebus volcano (Antarctica). It yields estimations of the plume velocity and gas flux at unprecedented time resolution (1 Hz) and high accuracy (uncertainty of 33%). During a ∼2 h experiment on 26 December 2006, SO2 flux varied between 0.17 and 0.89 ± 0.2 kg s −1 with a vertical plume velocity varying between 1 and 2.5 ± 0.1 m s −1. These measurements provide insight into the short-term variations of the passive degassing of this volcano renowned for its active lava lake. A cyclicity in flux, ranging from about 11-24 min, is evident. We propose two physical mechanisms to explain this degassing pattern, associated to periodic supply of either gas-rich magma or gas alone into the lake. The dual-wide field of view DOAS technique promises better integration of geochemical and geophysical observations and new insights into gas and magma dynamics, as well as processes of magma storage and gas segregation at active volcanoes.

A few lava lakes, like that at Erebus volcano, Antarctica, have been continuously active for decades, reaching a steady-state. We report spectroscopic and thermal observations from Erebus that reveal remarkable, phase-locked cycles of lava lake convection and gas plume composition. We argue that the observed fluctuations in gas ratios, including the SO2/CO2 content in the plume, identify two end-member contributions to the Erebus emission: a sustained source of CO2-rich gas percolating through permeable conduit magma, and a shallower source of H2O-rich gas exsolved from magma pulses that periodically enter the lava lake. The unstable magma flow may reflect the viscosity stratification between rising and descending magma in the conduit, and the resulting oscillatory behaviour of the phonolite lake exemplifies the relative roles of closed- and open-system degassing in persistently active volcanoes.

In this paper we combine SO2/ash plume dispersion modelling, satellite and surface remote sensing observations to study the regional influence of a relatively weak volcanic eruption from Mount Etna on the optical and micro-physical properties of Mediterranean aerosols. We analyse the Mount Etna eruption episode of 25–27 October 2013. The evolution of the plume along the trajectory is investigated by means of the FLEXPART (FLEXible PARTicle dispersion model) Lagrangian dispersion model. The satellite dataset includes true colour images, retrieved values of volcanic SO2 and ash, and estimates of SO2 and ash emission rates derived from MODIS (MODerate resolution Imaging Spectroradiometer) observations, and estimates of cloud top pressure from SEVIRI (Spinning Enhanced Visible and InfraRed Imager). Surface remote sensing measurements of aerosol and SO2 made at the ENEA Station for Climate Observations (35.52° N, 12.63° E, 50 m a.s.l.) on the island of Lampedusa are used in the analysis. The combination of these different datasets suggests that SO2 and ash, despite the initial injection occurred at about 7.0 km altitude, reached altitudes around 10–12 km and influenced the aerosol size distribution at a distance more than 350 km downwind. This study indicates that even a relatively weak volcanic eruption may produce an observable effect on the aerosol properties at the regional scale. The impact of secondary sulphate particles on the aerosol size distribution at Lampedusa is discussed, and estimates of the clear sky direct aerosol radiative forcing are derived. Daily shortwave radiative forcing efficiencies are calculated with the LibRadtran model. They are estimated between −39 and −48 W m−2 AOD−1 at the top of the atmosphere, and between −66 and −49 W m−2 AOD−1, at the surface, with the variability in the estimates mainly depending on the aerosol single scattering albedo. These results suggest that sulphate particles played a large role, while the contribution by ash particles was small in the volcanic plume arriving at Lampedusa during this event.

International audience; The 1999 basaltic eruption of Shishaldin volcano (Alaska, USA) displayed both classical Strombolian activity and an explosive Subplinian plume. Strombolian activity at Shishaldin occurred in two major phases following the Subplinian activity. In this paper, we use acoustic measurements to interpret the Strombolian activity. Acoustic measurements of the two Strombolian phases show a series of explosions that are modeled by the vibration of a large overpressurised cylindrical bubble at the top of the magma column. Results show that the bubble does not burst at its maximum radius, as expected if the liquid film is stretched beyond its elasticity. But bursting occurs after one cycle of vibration, as a consequence of an instability of the air-magma interface close to the bubble minimum radius. During each Strombolian period, estimates of bubble length and overpressure are calculated. Using an alternate method based on acoustic power, we estimate gas velocity to be 30-60 m/s, in very good agreement with synthetic waveforms. Although there is some variation within these parameters, bubble length and overpressure for the first Strombolian phase are found to be c 82 F 11 m and 0.083 MPa. For the second Strombolian phase, bubble length and overpressure are estimated at 24 F 12 m and 0.15 MPa for the first 17 h after which bubble overpressure shows a constant increase, reaching a peak of 1.4 MPa, just prior to the end of the second Strombolian phase. This peak suggests that, at the time, the magma in the conduit may contain a relatively large concentration of small bubbles. Maximum total gas volume and gas fluxes at the surface are estimated to be 3.3 Â 10 7 and 2.9 Â 10 3 m 3 /s for the first phase and 1.0 Â 10 8 and 2.2 Â 10 3 m 3 /s for the second phase. This gives a mass flux of 1.2 Â 10 3 and 8.7 Â 10 2 kg/s, respectively, for the first and the second Strombolian phases.

International audience The link between seismicity and degassing is investigated during the VEI 4 eruptions of Merapi volcano (Indonesia) in October and in early November 2010. Seismicity comprised a large number and variety of earthquakes including Volcano-Tectonic events, a sustained period of Long Period Seismicity (LPS), i.e., Long-Period events (LP), Very Long Period events (VLP) and tremor. LPS seismicity is ascribed to the excitation of fluid-filled cavity resonance and inertial displacement of fluids and magma. We investigate here LPS that occurred between 17 October and 4 November 2010 to obtain insights into the volcano eruption processes which preceded the paroxysmal phase of the eruption on 4-5 November. We proceed to the moment tensor inversion of a well-recorded large VLP event during the intrusion phase on 17 October 2010, i.e., before the first explosion on 26 October. By using two simplified models (crack and pipe), we find a shallow source for this VLP event at about 1 km to the south of the summit and less than 1 km below the surface. We analyze more than 90 LP events that occurred during the multi-phase eruption (29 October-4 November). We show that most of them have a dominant frequency in the range 0.2-4 Hz. We could locate 48 of those LP events; at least 3 clusters of LP events occurred successively. We interpret these observations as generated by different fluid-filled containers in the summit area that were excited while magma rose. We also observe significant variations of the complex frequency during the course of the eruption. We discuss these changes in terms of a variable ratio of fluid to solid densities and/or by possible conduit geometry change and/or permeability of the conduit or the edifice and/or by resonance of different fluid-containers during the release of more than 0.4 Tg of SO2 and large but unknown masses of other volcanic gases. Finally, we also discuss how the major explosions of the eruption were possibly triggered by passing waves resulting from regional tectonic earthquakes on 3 and 4 November.

International audience; Merapi volcano (Indonesia) is one of the most active and hazardous volcanoes in the world. It is known for frequent small to moderate eruptions, pyroclastic flows produced by lava dome collapse, and the large population settled on and around the flanks of the volcano that is at risk. Its usual behaviour for the last decades abruptly changed in late October and early November 2010, when the volcano produced its largest and most explosive eruptions in more than a century, displacing a third of a million people, and claiming nearly 400 lives. Despite the challenges involved in forecasting this 'hundred year eruption', we show that the magnitude of precursory signals (seismicity, ground deformation, gas emissions) were proportional to the large size and intensity of the eruption. In addition and for the first time, near-real-time satellite radar imagery played an equal role with seismic, geodetic, and gas observations in monitoring eruptive activity during a major volcanic crisis. The Indonesian Center of Volcanology and Geological Hazard Mitigation (CVGHM) issued timely forecasts of the magnitude of the eruption phases, saving 10,000-20,000 lives. In addition to reporting on aspects of the crisis management, we report the first synthesis of scientific observations of the eruption. Our monitoring and petrologic data show that the 2010 eruption was fed by rapid ascent of magma from depths ranging from 5 to 30 km. Magma reached the surface with variable gas content resulting in alternating explosive and rapid effusive eruptions, and released a total of ~ 0.44 Tg of SO2. The eruptive behaviour seems also related to the seismicity along a tectonic fault more than 40 km from the volcano, highlighting both the complex stress pattern of the Merapi region of Java and the role of magmatic pressurization in activating regional faults. We suggest a dynamic triggering of the main explosions on 3 and 4 November by the passing seismic waves generated by regional earthquakes on these days.

International audience; [1] Episodic magmatic degassing has been observed at numerous volcanoes, especially those of intermediate composition. It can span timescales from years to decades. Here we propose a physical model for the degassing of a shallow magma intrusion to explain this phenomenon. The magma cools by convection, which leads to melt crystallization, volatile exsolution, and magma overpressure. When the pressure reaches a critical value, wall rocks fracture and the exsolved gas escapes. The intrusion then returns to the initial lithostatic pressure and a new cooling-crystallization-degassing cycle occurs. A series of such cycles leads to episodic degassing. The trend and timescale of the degassing process are mainly governed by magma cooling. Two degassing regimes are exhibited: an early phase with a high frequency of gas pulses and a later phase with a lower gas pulse frequency. The transition between these two regimes is caused by the viscosity increase when the magma crystallinity exceeds the crystal percolation threshold. We find that the time to this transition is dependent on magma volume, to a first approximation. Where observations are available from sustained geochemical surveillance, the model provides constraints on key aspects of the subsurface magmatic system, with estimation of the volume of an intrusion and tensile strength of the surrounding rocks. It therefore represents a relevant tool for volcanic surveillance and hazard assessment.