SO2 Clouds Research Articles

Improved estimation of volcanic SO<sub>2</sub> injections from satellite retrievals and Lagrangian transport simulations: the 2019 Raikoke eruption

Abstract. Monitoring and modeling of volcanic plumes are important for understanding the impact of volcanic activity on climate and for practical concerns, such as aviation safety or public health. Here, we apply the Lagrangian transport model Massive-Parallel Trajectory Calculations (MPTRAC) to estimate the SO2 injections into the upper troposphere and lower stratosphere by the eruption of the Raikoke volcano (48.29∘ N, 153.25∘ E) in June 2019 and its subsequent long-range transport and dispersion. First, we used SO2 retrievals from the AIRS (Atmospheric Infrared Sounder) and TROPOMI (TROPOspheric Monitoring Instrument) satellite instruments together with a backward trajectory approach to estimate the altitude-resolved SO2 injection time series. Second, we applied a scaling factor to the initial estimate of the SO2 mass and added an exponential decay to simulate the time evolution of the total SO2 mass. By comparing the estimated SO2 mass and the mass from TROPOMI retrievals, we show that the volcano injected 2.1 ± 0.2 Tg SO2, and the e-folding lifetime of the SO2 was about 13 to 17 d. The reconstructed SO2 injection time series are consistent between using the AIRS nighttime and the TROPOMI daytime products. Further, we compared forward transport simulations that were initialized by AIRS and TROPOMI SO2 products with a constant SO2 injection rate. The results show that the modeled SO2 change, driven by chemical reactions, captures the SO2 mass variations from TROPOMI retrievals. In addition, the forward simulations reproduce the SO2 distributions in the first ∼10 d after the eruption. However, diffusion in the forward simulations is too strong to capture the internal structure of the SO2 clouds, which is further quantified in the simulation of the compact SO2 cloud from late July to early August. Our study demonstrates the potential of using combined nadir satellite retrievals and Lagrangian transport simulations to further improve SO2 time- and height-resolved injection estimates of volcanic eruptions.

Volcanic SO<sub>2</sub> layer height by TROPOMI/S5P: evaluation against IASI/MetOp and CALIOP/CALIPSO observations

Abstract. Volcanic eruptions eject large amounts of ash and trace gases such as sulfur dioxide (SO2) into the atmosphere. A significant difficulty in mitigating the impact of volcanic SO2 clouds on air traffic safety is that these gas emissions can be rapidly transported over long distances. The use of space-borne instruments enables the global monitoring of volcanic SO2 emissions in an economical and risk-free manner. Within the European Space Agency (ESA) Sentinel-5p+ Innovation project, the S5P SO2 layer height (S5P+I: SO2LH) activities led to the improvements of the retrieval algorithm and generation of the corresponding near real-time S5P SO2 LH products. These are currently operationally provided, in near real-time, by the German Aerospace Center (DLR) within the framework of the Innovative Products for Analyses of Atmospheric Composition (INPULS) project. The main aim of this paper is to present its extensive verification, accomplished within the S5P+I: SO2LH project, over major recent volcanic eruptions, against collocated space-borne measurements from the IASI/Metop and CALIOP/CALIPSO instruments as well as assess its impact on the forecasts provided by the Copernicus Atmospheric Monitoring Service (CAMS). The mean difference between S5P and IASI observations for the Raikoke 2019, the Nishinoshima 2020 and the La Soufrière-St Vincent 2021 eruptive periods is ∼ 0.5 ± 3 km, while for the Taal 2020 eruption, a larger difference was found, between 3 ± 3 km and 4 ± 3 km. The comparison of the daily mean SO2 LH further demonstrates the capabilities of this near real-time product, with slopes between 0.8 and 1 and correlation coefficients ranging between 0.6 and 0.8. Comparisons between the S5P SO2 LH and the CALIOP/CALIPSO ash plumes revealed an expected bias at −2.5 ± 2 km, considering that the injected SO2 and ash plume locations do not always coincide over an eruption. Furthermore, the CAMS assimilation of the S5P SO2 LH product led to much improved model output against the non-assimilated IASI LH, with a mean difference of 1.5 ± 2 km, compared to the original CAMS analysis, and improved the geographical spread of the Raikoke volcanic plume following the eruptive days.

Ensemble-Based Forecast of Volcanic Clouds Using FALL3D-8.1

Operational forecasting of volcanic ash and SO2 clouds is challenging due to the large uncertainties that typically exist on the eruption source term and the mass removal mechanisms occurring downwind. Current operational forecast systems build on single-run deterministic scenarios that do not account for model input uncertainties and their propagation in time during transport. An ensemble-based forecast strategy has been implemented in the FALL3D-8.1 atmospheric dispersal model to configure, execute, and post-process an arbitrary number of ensemble members in a parallel workflow. In addition to intra-member model domain decomposition, a set of inter-member communicators defines a higher level of code parallelism to enable future incorporation of model data assimilation cycles. Two types of standard products are automatically generated by the ensemble post-process task. On one hand, deterministic forecast products result from some combination of the ensemble members (e.g., ensemble mean, ensemble median, etc.) with an associated quantification of forecast uncertainty given by the ensemble spread. On the other hand, probabilistic products can also be built based on the percentage of members that verify a certain threshold condition. The novel aspect of FALL3D-8.1 is the automatisation of the ensemble-based workflow, including an eventual model validation. To this purpose, novel categorical forecast diagnostic metrics, originally defined in deterministic forecast contexts, are generalised here to probabilistic forecasts in order to have a unique set of skill scores valid to both deterministic and probabilistic forecast contexts. Ensemble-based deterministic and probabilistic approaches are compared using different types of observation datasets (satellite cloud detection and retrieval and deposit thickness observations) for the July 2018 Ambae eruption in the Vanuatu archipelago and the April 2015 Calbuco eruption in Chile. Both ensemble-based approaches outperform single-run simulations in all categorical metrics but no clear conclusion can be extracted on which is the best option between these two.

Tracking aerosols and SO<sub>2</sub> clouds from the Raikoke eruption: 3D view from satellite observations

Abstract. The 21 June 2019 eruption of the Raikoke volcano (Kuril Islands, Russia; 48∘ N, 153∘ E) produced significant amounts of volcanic aerosols (sulfate and ash) and sulfur dioxide (SO2) gas that penetrated into the lower stratosphere. The dispersed SO2 and sulfate aerosols in the stratosphere were still detectable by multiple satellite sensors for many months after the eruption. For this study of SO2 and aerosol clouds we use data obtained from two of the Ozone Mapping and Profiler Suite sensors on the Suomi National Polar-orbiting Partnership satellite: total column SO2 from the Nadir Mapper and aerosol extinction profiles from the Limb Profiler as well as other satellite data sets. We evaluated the limb viewing geometry effect (the “arch effect”) in the retrieval of the LP standard aerosol extinction product at 674 nm. It was shown that the amount of SO2 decreases with a characteristic period of 8–18 d and the peak of stratospheric aerosol optical depth recorded at a wavelength of 674 nm lags the initial peak of SO2 mass by 1.5 months. Using satellite observations and a trajectory model, we examined the dynamics of an unusual atmospheric feature that was observed, a stratospheric coherent circular cloud of SO2 and aerosol from 18 July to 22 September 2019.

Day–Night Monitoring of Volcanic SO2 and Ash Clouds for Aviation Avoidance at Northern Polar Latitudes

We describe NASA’s Applied Sciences Disasters Program, which is a collaborative project between the Direct Readout Laboratory (DRL), ozone processing team, Jet Propulsion Laboratory, Geographic Information Network of Alaska (GINA), and Finnish Meteorological Institute (FMI), to expedite the processing and delivery of direct readout (DR) volcanic ash and sulfur dioxide (SO2) satellite data. We developed low-latency quantitative retrievals of SO2 column density from the solar backscattered ultraviolet (UV) measurements using the Ozone Mapping and Profiler Suite (OMPS) spectrometers as well as the thermal infrared (TIR) SO2 and ash indices using Visible Infrared Imaging Radiometer Suite (VIIRS) instruments, all flying aboard US polar-orbiting meteorological satellites. The VIIRS TIR indices were developed to address the critical need for nighttime coverage over northern polar regions. Our UV and TIR SO2 and ash software packages were designed for the DRL’s International Planetary Observation Processing Package (IPOPP); IPOPP runs operationally at GINA and FMI stations in Fairbanks, Alaska, and Sodankylä, Finland. The data are produced within 30 min of satellite overpasses and are distributed to the Alaska Volcano Observatory and Anchorage Volcanic Ash Advisory Center. FMI receives DR data from GINA and posts composite Arctic maps for ozone, volcanic SO2, and UV aerosol index (UVAI, proxy for ash or smoke) on its public website and provides DR data to EUMETCast users. The IPOPP-based software packages are available through DRL to a broad DR user community worldwide.

The 2019 Raikoke volcanic eruption – Part 1: Dispersion model simulations and satellite retrievals of volcanic sulfur dioxide

Abstract. Volcanic eruptions can cause significant disruption to society, and numerical models are crucial for forecasting the dispersion of erupted material. Here we assess the skill and limitations of the Met Office's Numerical Atmospheric-dispersion Modelling Environment (NAME) in simulating the dispersion of the sulfur dioxide (SO2) cloud from the 21–22 June 2019 eruption of the Raikoke volcano (48.3∘ N, 153.2∘ E). The eruption emitted around 1.5±0.2 Tg of SO2, which represents the largest volcanic emission of SO2 into the stratosphere since the 2011 Nabro eruption. We simulate the temporal evolution of the volcanic SO2 cloud across the Northern Hemisphere (NH) and compare our model simulations to high-resolution SO2 measurements from the TROPOspheric Monitoring Instrument (TROPOMI) and the Infrared Atmospheric Sounding Interferometer (IASI) satellite SO2 products. We show that NAME accurately simulates the observed location and horizontal extent of the SO2 cloud during the first 2–3 weeks after the eruption but is unable, in its standard configuration, to capture the extent and precise location of the highest magnitude vertical column density (VCD) regions within the observed volcanic cloud. Using the structure–amplitude–location (SAL) score and the fractional skill score (FSS) as metrics for model skill, NAME shows skill in simulating the horizontal extent of the cloud for 12–17 d after the eruption where VCDs of SO2 (in Dobson units, DU) are above 1 DU. For SO2 VCDs above 20 DU, which are predominantly observed as small-scale features within the SO2 cloud, the model shows skill on the order of 2–4 d only. The lower skill for these high-SO2-VCD regions is partly explained by the model-simulated SO2 cloud in NAME being too diffuse compared to TROPOMI retrievals. Reducing the standard horizontal diffusion parameters used in NAME by a factor of 4 results in a slightly increased model skill during the first 5 d of the simulation, but on longer timescales the simulated SO2 cloud remains too diffuse when compared to TROPOMI measurements. The skill of NAME to simulate high SO2 VCDs and the temporal evolution of the NH-mean SO2 mass burden is dominated by the fraction of SO2 mass emitted into the lower stratosphere, which is uncertain for the 2019 Raikoke eruption. When emitting 0.9–1.1 Tg of SO2 into the lower stratosphere (11–18 km) and 0.4–0.7 Tg into the upper troposphere (8–11 km), the NAME simulations show a similar peak in SO2 mass burden to that derived from TROPOMI (1.4–1.6 Tg of SO2) with an average SO2 e-folding time of 14–15 d in the NH. Our work illustrates how the synergy between high-resolution satellite retrievals and dispersion models can identify potential limitations of dispersion models like NAME, which will ultimately help to improve dispersion modelling efforts of volcanic SO2 clouds.

A multi-sensor satellite-based archive of the largest SO<sub>2</sub> volcanic eruptions since 2006

Abstract. We present a multi-sensor archive collecting spatial and temporal information about volcanic SO2 clouds generated by the 11 largest eruptions of this century. The detection and monitoring of volcanic clouds are an important topic for aviation management, climate issues and weather forecasts. Several studies focusing on single eruptive events exist, but no archive available at the moment combines quantitative data from as many instruments. We archived and collocated the SO2 vertical column density estimations from three different satellite instruments (AIRS, IASI and GOME-2), atmospheric parameters as vertical profiles from the Global Navigation Satellite Systems (GNSS) Radio Occultations (RO), and the cloud-top height and aerosol type from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP). Additionally, we provide information about the cloud-top height from three different algorithms and the atmospheric anomaly due to the presence of the cloud. The dataset is gathering 206 d of SO2 data, collocated with 44 180 backscatter profiles and 64 764 radio occultation profiles. The modular structure of the archive allows an easy collocation of the datasets according to the users' needs, and the cross-comparison of the datasets shows different consistency of the parameters estimated with different sensors and algorithms, according to the sensitivity and resolution of the instruments. The data described here are published with a DOI at https://doi.org/10.5880/fidgeo.2020.016 (Tournigand et al., 2020a).

Extension of the WRF-Chem volcanic emission preprocessor to integrate complex source terms and evaluation for different emission scenarios of the Grimsvötn 2011 eruption

Abstract. Volcanic eruptions may generate volcanic ash and sulfur dioxide (SO2) plumes with strong temporal and vertical variations. When simulating these changing volcanic plumes and the afar dispersion of emissions, it is important to provide the best available information on the temporal and vertical emission distribution during the eruption. The volcanic emission preprocessor of the chemical transport model WRF-Chem has been extended to allow the integration of detailed temporally and vertically resolved input data from volcanic eruptions. The new emission preprocessor is tested and evaluated for the eruption of the Grimsvötn volcano in Iceland 2011. The initial ash plumes of the Grimsvötn eruption differed significantly from the SO2 plumes, posing challenges to simulate plume dynamics within existing modelling environments: observations of the Grimsvötn plumes revealed strong vertical wind shear that led to different transport directions of the respective ash and SO2 clouds. Three source terms, each of them based on different assumptions and observational data, are applied in the model simulations. The emission scenarios range from (i) a simple approach, which assumes constant emission fluxes and a predefined vertical emission profile, to (ii) a more complex approach, which integrates temporarily varying observed plume-top heights and estimated emissions based on them, to (iii) the most complex method that calculates temporal and vertical variability of the emission fluxes based on satellite observations and inversion techniques. Comparisons between model results and independent observations from satellites, lidar, and surface air quality measurements reveal the best performance of the most complex source term.

Fast ascent rate during the 2017–2018 Plinian eruption of Ambae (Aoba) volcano: a petrological investigation

In September 2017, after more than a hundred years of quiescence, Ambae (Aoba), Vanuatu’s largest volcano, entered a new phase of eruptive activity, triggering the evacuation of the island’s 11,000 inhabitants resulting in the largest volcanic disaster in the country’s history. Three subsequent eruptive phases in November 2017, March 2018, and July 2018 expelled some of the largest tropospheric and stratospheric SO2 clouds observed in the last decade. Here, we investigate the mechanisms and dynamics of this eruption. We use major elements, trace elements, and volatiles in olivine and clinopyroxene hosted melt inclusions, embayments, crystals, and matrix glasses together with clinopyroxene geobarometry and olivine, plagioclase, and clinopyroxene geothermometry to reconstruct the physical and chemical evolution of the magma, as it ascends to the surface. Volatile elements in melt inclusions and geobarometry data suggest that the magma originated from depth of ~ 14 km before residing at shallow (~ 0.5 to 3 km) levels. Magma ascent to the surface was likely facilitated by shallow phreatic eruptions that opened a pathway for magma to ascend. Succeeding eruptive phases are characterised by increasingly primitive compositions with evidence of small amounts of mixing having taken place. Mg–Fe exchange diffusion modelling yields olivine residence times in the magma chamber ranging from a few days to a year prior to eruption. Diffusion modelling of volatiles along embayments (melt channels) from the first two phases of activity and microlite number density suggests rapid magma ascent in the range of 15–270 km/h, 4–75 m/s (decompression rates of 0.1 to ~ 2 MPa/s) corresponding to a short travel time between the top of the shallow reservoir and the surface of less than 2 min.

First Observations of Volcanic Eruption Clouds From the L1 Earth‐Sun Lagrange Point by DSCOVR/EPIC

AbstractVolcanic sulfur dioxide (SO2) emissions have been measured by ultraviolet sensors on polar‐orbiting satellites for several decades but with limited temporal resolution. This precludes studies of key processes believed to occur in young (~1–3 hr old) volcanic clouds. In 2015, the launch of the Earth Polychromatic Imaging Camera (EPIC) aboard the Deep Space Climate Observatory (DSCOVR) provided an opportunity for novel observations of volcanic eruption clouds from the first Earth‐Sun Lagrange point (L1). The L1 vantage point provides continuous observations of the sunlit Earth, offering up to eight or nine observations of volcanic SO2 clouds in the DSCOVR/EPIC field of view at ~1‐hr intervals. Here we demonstrate DSCOVR/EPIC's sensitivity to volcanic SO2 using several volcanic eruptions from the tropics to midlatitudes. The hourly cadence of DSCOVR/EPIC observations permits more timely measurements of volcanic SO2 emissions, improved trajectory modeling, and novel analyses of the temporal evolution of volcanic clouds.

Intercomparison of SO2 camera systems for imaging volcanic gas plumes

SO2 camera systems are increasingly being used to image volcanic gas plumes. The ability to derive SO2 emission rates directly from the acquired imagery at high time resolution allows volcanic process studies that incorporate other high time-resolution datasets. Though the general principles behind the SO2 camera have remained the same for a number of years, recent advances in CCD technology and an improved understanding of the physics behind the measurements have driven a continuous evolution of the camera systems. Here we present an intercomparison of seven different SO2 cameras. In the first part of the experiment, the various technical designs are compared and the advantages and drawbacks of individual design options are considered. Though the ideal design was found to be dependent on the specific application, a number of general recommendations are made. Next, a time series of images recorded by all instruments at Stromboli Volcano (Italy) is compared. All instruments were easily able to capture SO2 clouds emitted from the summit vents. Quantitative comparison of the SO2 load in an individual cloud yielded an intra-instrument precision of about 12%. From the imagery, emission rates were then derived according to each group's standard retrieval process. A daily average SO2 emission rate of 61±10t/d was calculated. Due to differences in spatial integration methods and plume velocity determination, the time-dependent progression of SO2 emissions varied significantly among the individual systems. However, integration over distinct degassing events yielded comparable SO2 masses. Based on the intercomparison data, we find an approximate 1-sigma precision of 20% for the emission rates derived from the various SO2 cameras. Though it may still be improved in the future, this is currently within the typical accuracy of the measurement and is considered sufficient for most applications.

Circumpolar transport of a volcanic cloud from Hekla (Iceland)

Hekla volcano (Iceland) erupted on 17 August 1980 and emplaced a sulfur dioxide (SO2) cloud into the north polar stratosphere at a maximum altitude of ∼15 km. The SO2 is tracked using satellite data from the ultraviolet (UV) Nimbus‐7 Total Ozone Mapping Spectrometer (N7/TOMS) and the infrared (IR) High‐resolution Infrared Radiation Sounder (HIRS/2) on the NOAA TIROS Operational Vertical Sounder (TOVS) platform. The eruption emitted ∼0.5–0.7 Tg of SO2, which later split into three distinct clouds, one of which circled the North Pole at the perimeter of an atypically persistent Arctic cyclone for six days, impacting airspace on three continents. Separate clouds drifted across eastern Russia, Alaska, and Canada. Maximum SO2 columns derived from TOMS and HIRS/2 accurately define the volcanic cloud's path and fit trajectories produced by the NOAA Hybrid Single‐Particle Lagrangian Integrated Trajectory (HYSPLIT) model, providing confidence in the model. When combined with the SO2 measurements, trajectory altitudes derived from HYSPLIT provide robust estimates of the altitudes of the SO2 clouds (8–15 km), which would be elusive using either the satellite data or the trajectory model in isolation. Near‐coincident, spectrally discrete UV and IR retrievals are compared in the volcanic cloud and indicate good agreement between TOMS and HIRS/2 SO2 columns for pixels with similar viewing geometry. Hekla eruptions, which follow a pattern of early explosive venting of volcanic gases with significant stratospheric injection, could play a role in promoting Arctic ozone loss, depending on the phase of the North Atlantic Oscillation during the eruption.

Impact of atmospheric water vapor on the thermal infrared remote sensing of volcanic sulfur dioxide emissions: A case study from the Pu'u ‘O’ vent of Kilauea Volcano, Hawaii

The December 18, 1999, launch of NASA's Terra satellite put two multispectral thermal infrared imaging instruments into Earth orbit. Experiments with airborne instruments have demonstrated that the data from such instruments can be used to detect volcanic SO2 plumes and clouds. However, one of the greatest challenges that will confront efforts to monitor volcanic SO2 emissions from space is the need to characterize the local atmosphere. In this paper we evaluate the sensitivity of the SO2 retrieval procedure to our knowledge of the local atmospheric conditions. We compare SO2 retrievals obtained with distant (radiosonde) and local (Fourier transform infrared (FTIR) soundings) atmospheric measurements and find that the relative difference is typically ±25%. For ground temperature retrievals the relative difference is ±1.5%. These results lead us to conclude that while local measurements of atmospheric conditions are preferable, useful retrievals can be obtained using atmospheric measurements from distant sites. In addition, we find very good agreement between SO2 and ground temperature retrievals obtained from thermal infrared imagery and FTIR soundings.

Detection of volcanic ash clouds from Nimbus 7/total ozone mapping spectrometer

Measured radiances from the Version 7 reprocessing of the Nimbus 7/total ozone mapping spectrometer (TOMS) 340‐ and 380‐nm channels are used to detect absorbing particulates injected into the atmosphere after the El Chichon eruption on April 4, 1982. It is shown that while the single‐channel reflectivity determined from the 380‐nm channel is able to detect clouds and haze composed of nonabsorbing aerosols, the spectral contrast between the 340‐ and 380‐nm channels is sensitive to absorbing particulates such as volcanic ash, desert dust, or smoke from biomass burning. In this paper the spectral contrast between these two channels is used to detect the volcanic ash injection into the atmosphere and to track its evolution for several days. The movement of the ash clouds is shown to be consistent with the motions expected from the National Centers for Environmental Prediction (NCEP)‐derived balanced wind fields in the troposphere and lower stratosphere. The movement of the volcanic SO2 cloud detected from TOMS data was also in agreement with the NCEP wind at higher altitudes of up to 100–10 mbar. The vertical wind shear in the neighborhood of the eruption site resulted in a clear separation of the ash and SO2 clouds. The location and movement of the ash cloud are consistent with information obtained by the advanced very high resolution radiometer (AVHRR) instrument on board the NOAA 7 satellite and to ground reports of ash fall.

Oxidation of volcanic SO2: A sink for stratospheric OH and H2O

The oxidation of volcanic SO2 to sulphate in the stratosphere is investigated for an anomalously large volcanic sulphur loading (∼200 Mt) similar to that generated by the Tambora eruption in 1815. Model calculations suggest that stratospheric OH levels can be severely reduced by the gas‐phase oxidation of SO2. One implication of this is that dense volcanic SO2 clouds may last substantially longer than previously thought. SO2 oxidation is also found to lead to significant stratospheric dehydration because approximatively three molecules of H2O are consumed for each molecule of SO2 converted to sulphate. Therefore it is possible to view major volcanic eruptions as a potential sink for stratospheric H2O.

Transport of Cerro Hudson SO2 clouds

The Cerro Hudson volcano in southern Chile (45.92°S, 73.0°W) emitted large ash and sulfur dioxide clouds on August 12–15, following several days of minor activity [Global Volcanism Network Bulletin, 1991]. The SO2 clouds were observed using (preliminary) near real‐time data from the Total Ozone Mapping Spectrometer (TOMS) as they encircled the south polar region. The injection of SO2 into the stratosphere has essentially created a gigantic chemical tracer that could provide new insights into the wind patterns and seasonal circulation around the Antarctic region.around the Antarctic region. The TOMS instrument, on board the National Aeronautic and Space Administration's Nimbus 7 satellite, measures the ratio of backscattered Earth radiance to incoming solar irradiance in the ultraviolet spectrum. Although originally designed to measure ozone, it was later discovered that the TOMS instrument could also detect and quantify SO2 [Krueger, 1985]. After this discovery, measurements from TOMS were examined for SO2 emissions for all recorded volcanic eruptions since Nimbus‐7 was launched in October 1978, and current data are analyzed as new eruptions occur. The satellite is in a polar, Sun‐synchronous orbit so that it crosses the equator at local noon and observes the whole sunlit Earth in approximately 14 orbits each day. Total column amounts of SO2 are determined that represent the amount of gas affecting the reflection of ultraviolet light through a column of the atmosphere from the satellite to the reflecting surface, Earth, given in terms of milli atmospheres centimeter (1000 milli atm cm = a gas layer 1‐cm thick at STP). The mass of SO2 is calculated by integrating over the cloud area to obtain a volume, then converting to tons.