Abstract
The underwater Hunga Tonga-Hunga Ha-apai volcano erupted in the early hours of 15th January 2022, and injected volcanic gases and aerosols to over 50 km altitude. Here we synthesise ...satellite, ground-based, in situ and radiosonde observations of the eruption to investigate the strength of the stratospheric aerosol and water vapour perturbations in the initial weeks after the eruption and we quantify the net radiative impact across the two species using offline radiative transfer modelling. We find that the Hunga Tonga-Hunga Ha-apai eruption produced the largest global perturbation of stratospheric aerosols since the Pinatubo eruption in 1991 and the largest perturbation of stratospheric water vapour observed in the satellite era. Immediately after the eruption, water vapour radiative cooling dominated the local stratospheric heating/cooling rates, while at the top-of-the-atmosphere and surface, volcanic aerosol cooling dominated the radiative forcing. However, after two weeks, due to dispersion/dilution, water vapour heating started to dominate the top-of-the-atmosphere radiative forcing, leading to a net warming of the climate system.
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.
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 preexisting, caldera ring‐faults. Ring‐fault reactivation, however, is rarely observed during moderate‐sized eruptions. On February 21, 2015 at Ambrym volcano in Vanuatu, a basaltic dike intrusion produced more than 1 m of coeruptive uplift, as measured by InSAR, synthetic aperture radar correlation, and Multiple Aperture Interferometry. 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 × 106 m3 and erupted ∼9.3 × 106 m3 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.
Plain Language Summary
Many volcanoes feature large depressions, called calderas. Calderas form when enough magma leaves a deep reservoir, and the solid rock lid above this reservoir can no longer support its own weight. Caldera faults, or cracks surrounding the reservoir which extend from the reservoir to Earth's surface, form as the lid collapses. Ambrym volcano (Vanuatu) has a 12‐km wide caldera, and researchers propose it formed during an explosive eruption 2,000 years ago. However, in 2018, Ambrym's caldera sunk along caldera faults during a nonexplosive eruption. This observation questions whether an explosive eruption was necessary to form Ambrym's caldera in the first place. Furthermore, in February 2015, an eruption 10 times smaller than in 2018 also caused the ground to sink along caldera faults. Utilizing ground motion data obtained from satellite radar systems to model magma reservoir outflow and fault displacement, we conclude that, in 2015, the ground sank along caldera faults. This sinking is explained by the removal of as little as <7% of the stored magma from the reservoir. We therefore propose that Ambrym's wide caldera may have formed as a result of many frequently occurring medium‐sized eruptions. This challenges the thought that wide calderas mainly form as a result of large eruptions.
Key Points
Ground displacement at Ambrym in February 2015 was caused by a dike intrusion, deflating reservoir, and normal slip on a caldera ring‐fault
Extracting at most 7% of the magma from Ambrym's reservoir suffices to reactivate the caldera ring‐faults
Normal slip along Ambrym's ring‐fault can occur during moderate‐sized eruptions, resulting in subsidence and further caldera development
•Persistent volcanic degassing explains steady subsidence at Ambrym from 2015–2017.•Subsidence ends in 2017, as magma influx balances depressurization from degassing.•We estimate a magma influx of ...∼0.16 km3 over 14 months before the 2018 eruption.
Despite being a widespread and common process, the impact of passive volcanic degassing on the pressurization state of a magma reservoir is not well understood. If mass loss due to gas emissions results in reservoir depressurization and surface subsidence, the pressure difference between a shallow reservoir and deep magma source may result in magma recharge and eventually trigger an eruption. It is therefore important to determine how a simplified reservoir-conduit system responds to such degassing processes. Here we use an extreme example of persistent volcanic degassing—Ambrym—as a case study to relate sulphur dioxide mass flux with reservoir depressurization and edifice-scale subsidence, both measured from satellite-based remote sensing observations. A geodetic inversion of surface displacements measured with Interferometric Synthetic Aperture Radar modeled using the Boundary Element Method provides bounds on the reservoir pressure change during an episode of subsidence at Ambrym from 2015 to 2017. These results are input into a lumped parameter theoretical model developed by Girona et al. (2014), and the free parameters (e.g., reservoir size and conduit radius) are systematically explored. We find that the 2015–2017 subsidence episode is consistent with pressure decreasing at a rate of −5.2 to −2.0 MPa year−1 in a reservoir at ∼2 km b.s.l., as a result of passive degassing. The subsidence episode is observed to end abruptly in October 2017, and no significant deformation is detected in the 14 months leading up to a rift zone intrusion and submarine eruption in December 2018, despite substantial degassing. We explain this lack of pre-eruptive deformation by an influx of ∼0.16 km3 of magma into a shallow (<2 km b.s.l.) reservoir that counterbalances the depressurization caused by degassing. This recharge volume is comparable with the volume of magma subsequently extracted from Ambrym's reservoir in December 2018. We conclude that at some open-vent passively degassing volcanoes, deflation caused by degassing may reduce or even cancel any inflation signal caused by magma influx. Nonetheless, detection of pre-eruptive recharge can be achieved by monitoring changes in the long-term deformation rate.
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 ≈82±11 m and 0.083 MPa. For the second Strombolian phase, bubble length and overpressure are estimated at 24±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.
Volcanic sulfate aerosols play a key role in air quality and climate. However, the rate of oxidation of sulfur dioxide (SO2) precursor gas to sulfate aerosols (SO42-) in volcanic plumes is poorly ...known, especially in the troposphere. Here we determine the chemical speciation as well as the intensity and temporal persistence of the impact on air quality of sulfate aerosols from the 2014–2015 Holuhraun flood lava eruption of Icelandic volcano Bárðarbunga. To do so, we jointly analyse a set of SO2 observations from satellite (OMPS and IASI) and ground-level measurements from air quality monitoring stations together with high temporal resolution mass spectrometry measurements of an Aerosol Chemical Speciation Monitor (ACSM) performed far from the volcanic source. We explore month/year long ACSM data in France from stations in contrasting environments, close and far from industrial sulfur-rich activities. We demonstrate that volcanic sulfate aerosols exhibit a distinct chemical signature in urban/rural conditions, with NO3:SO4 mass concentration ratios lower than for non-volcanic background aerosols. These results are supported by thermodynamic simulations of aerosol composition, using the ISORROPIA II model, which show that ammonium sulfate aerosols are preferentially formed at a high concentration of sulfate, leading to a decrease in the production of particulate ammonium nitrate. Such a chemical signature is however more difficult to identify at heavily polluted industrial sites due to a high level of background noise in sulfur. Nevertheless, aged volcanic sulfates can be distinguished from freshly emitted industrial sulfates according to their contrasting degree of anion neutralization. Combining AERONET (AErosol RObotic NETwork) sunphotometric data with ACSM observations, we also show a long persistence over weeks of pollution in volcanic sulfate aerosols, while SO2 pollution disappears in a few days at most. Finally, gathering 6-month long datasets from 27 sulfur monitoring stations of the EMEP (European Monitoring and Evaluation Programme) network allows us to demonstrate a much broader large-scale European pollution, in both SO2 and SO4, associated with the Holuhraun eruption, from Scandinavia to France. While widespread SO2 anomalies, with ground-level mass concentrations far exceeding background values, almost entirely result from the volcanic source, the origin of sulfate aerosols is more complex. Using a multi-site concentration-weighted trajectory analysis, emissions from the Holuhraun eruption are shown to be one of the main sources of SO4 at all EMEP sites across Europe and can be distinguished from anthropogenic emissions from eastern Europe but also from Great Britain. A wide variability in SO2:SO4 mass concentration ratios, ranging from 0.8 to 8.0, is shown at several stations geographically dispersed at thousands of kilometres from the eruption site. Despite this apparent spatial complexity, we demonstrate that these mass oxidation ratios can be explained by a simple linear dependency on the age of the plume, with a SO2-to-SO4 oxidation rate of 0.23 h−1. Most current studies generally focus on SO2, an unambiguous and more readily measured marker of the volcanic plume. However, the long persistence of the chemical fingerprint of volcanic sulfate aerosols at continental scale, as shown for the Holuhraun eruption here, casts light on the impact of tropospheric eruptions and passive degassing activities on air quality, health, atmospheric chemistry and climate.
In this paper we combine SO2 and ash plume dispersion modelling with 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 FLEXible PARTicle Lagrangian dispersion (FLEXPART) model. The satellite data set includes true colour images, retrieved values of volcanic SO2 and ash, 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 data sets suggests that SO2 and ash, despite the initial injection at about 7.0 km altitude, reached altitudes around 10–12 km and influenced the column average aerosol particle size distribution at a distance of 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 sulfate 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, i.e. radiative forcing per unit AOD (aerosol optical depth), 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 sulfate particles played a large role in the transported plume composition and radiative forcing, while the contribution by ash particles was small in the volcanic plume arriving at Lampedusa during this event.
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.
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.