SO2 Flux Research Articles

Magma recharge at Manam volcano, Papua New Guinea, identified through thermal and SO2 satellite remote sensing of open-vent emissions

Manam is one of the most frequently active volcanoes in Papua New Guinea and is a top contributor to global volcanic volatile emissions due to its persistent open-vent degassing. Here, we present a multi-year time series (2018–2021) of thermal and SO2 emissions for Manam from satellite remote sensing, which we interpret in the context of open-vent feedback between magma supply, reservoir pressure, and outgassing. We classify the time series into four phases based on the varying SO2 flux and observe a transient, yet substantial, increase in time-averaged SO2 flux from background levels of ~ 0.6 to ~ 4.72 kt day−1 between March and July 2019. We also identify a transition from temporally coupled to decoupled gas and thermal emissions during this period which we explain in the context of a magma recharge event that supplied new, volatile-rich magma to the shallow plumbing system beneath Manam. We infer that the arrival of this recharge magma triggered the series of eruptions between August 2018 and March 2019. These explosive events collectively removed 0.18 km3 of degassed residual magma and signalled the onset of a renewed period of unrest that ultimately culminated in a major eruption on 28 June 2019. We quantify the magnitude of “excess” degassing at Manam after the removal of the inferred residual magma. SO2 emissions reveal that ~ 0.18 km3 of magma was supplied, but only ~ 0.08 km3 was erupted between April 2019 and December 2021. We highlight how multi-parameter remote sensing observations over months to years enable the interpretation of open-vent processes that may be missed by short-duration campaign measurements.

An investigation into Venusian atmospheric chemistry based on an open-access photochemistry-transport model at 0–112 km

Atmospheric chemistry plays a crucial role in the evolution of climate habitability on Venus. It has been widely explored by chemistry-transport models, but some characteristics are still poorly interpreted. This study is devoted to developing an open-access chemistry-transport model spanning both the middle and lower atmospheres of Venus. It provides a scheme for the structure of the chemistry, especially for the sulfur and oxygen, and investigates the influence of the cloud diffusivity and the SO2 dissolution that are adopted in the clouds. The developed model is based on the VULCAN framework and was updated with the state-of-the-art Venusian atmospheric chemistry. It includes vertical eddy diffusion retrieved recently with the Venus Express observations, and it resolves radiative transfer containing gas absorption and scattering, Mie scattering of the cloud droplets, and absorption of the unknown UV absorber. The obtained abundance profiles of SO, SO2, CO, COS, O, O2, O3, HCl, and NO are in overall agreement with the observations. The results show that the increase in cloud diffusivity has slight effects on the chemical structure. The SO2 mainly dissolves in 50–90 km and evaporates below the clouds. The rapid dissolution-release cycle is responsible for the large upward flux of SO2 at 58 km. At around 70 km, SO has a significant peak that is larger than that of previous studies by an order of magnitude, and S and SO2 also show slight increases. They are attributed to the buffering effects of liquid SO2 in the clouds. O2 is significantly eliminated by SO in this layer. We emphasize the superior regulation of the sulfur cycle on O2 at 70 km and its potential contributions to the long-standing problem of the overestimated O2 abundance.

Surface-air exchanges of H2S and SO2 in an urban wetland in eastern China

Wetlands are widely recognized as hot spots for the emission or deposition of biogenic sulfur gases, including hydrogen sulfur (H2S) and sulfur dioxide (SO2), which significantly affect air quality and climate change. With the expansion of urban wetlands, it is critical to know the roles that urban wetlands played in atmospheric H2S and SO2 budget. In this study, the surface-air exchange fluxes of H2S and SO2 were measured by the Dynamic Flux Chamber (DFC) method in a typical urban wetland in eastern China from Sep 2022 to July 2023. It was found that the urban wetland did not have the expected high H2S emission, might be caused by the relatively high pH value and low sulfate concentration in the soil. Although H2S showed emission in the daytime of spring and summer, an overall H2S flux of −0.04 kg S ha−1 yr−1 was observed throughout the year. Meanwhile, the urban wetland presented a net sink of SO2, with a deposition flux of 0.14 kg S ha−1 yr−1. The negative peaks of SO2 flux corresponded to the suddenly elevated SO2 concentration in the ambient air especially in spring and winter. Through linear fitting of SO2 flux and concentration, the concept of SO2 “compensation point” was proposed. The compensation point is the concentration level at which the observed SO2 flux equals zero. The “compensation point” changed with the season and was related to temperature and humidity. The “compensation point” in summer and autumn were larger, being 2.37 ppb and 1.40 ppb, respectively, while they were 1.07 ppb and 0.86 ppb in spring and winter respectively. Our results suggest that the urban wetland expansion may have little risk of increasing air H2S but could act as a significant sink of SO2 with high SO2 concentration in the urban region.

Experimental study and analysis of CO2 and SO2 absorption in various water-based nanofluids by response surface methodology

The absorption of acidic gases in the oil and gas industries is important due to their toxicity and corrosive effects. Recently, the application of nanofluids based on aqueous or organic solvents as absorbents has been examined by a variety of researchers. In this study, a single bubble column was exploited to study the effect of water-based nanofluids on the absorption processes of SO2 and CO2 using response surface methodology (RSM) based on Box-Behnken three-level experiment design. With this in mind, CO2 and SO2 are separately injected at the bottom of a bubble column filled with one of the nanofluids: Al2O3-water, SiO2-water, or ZnO-water for each experiment. Then, the rate of SO2 or CO2 absorption in the nanofluids has been elucidated. The effect of important parameters including the weight fraction of the nanoparticles (NPs) (0.01, 0.055, and 0.1 wt.%), gas–liquid contact time (150, 300, and 450 s), and the diameter of nozzle for gas injection (0.46, 0.57, and 0.68 mm) have been studied. Results revealed that the maximum molar flux of both gases was observed in the ZnO-water nanofluid, followed by the SiO2-water nanofluid. In addition, increasing the nanoparticle mass fraction and the bubble size causes the molar flux to rise. However, increasing the gas–liquid contact time causes the molar flux of the mentioned gases to decrease. Finally, a set of the accurate equations has been proposed to predict the molar flux of SO2 and CO2 in the various nanofluids assessed in this work.

Quantitative determination of SO2 flux from industrial chimney through machine vision with plume model verification

Accurate emission rate data are essential for assessing the environmental impact of industrial emissions and calibrating the effectiveness of desulphurization equipment. However, accurate calculation of SO2 emissions is difficult to realize due to the environmental conditions, the time-varying and spatial non-uniformity of pollution source emissions. In this paper, a Dual TV Plus machine vision algorithm is proposed, calculated and validated with the conventional Farneback optical flow algorithm in terms of industrial stack emission rates. In addition, the effect of turbulence on the inversion of SO2 emission under complex environmental conditions is evaluated by combining simulations and outdoor experiments. The results show that the Dual TV Plus algorithm performs better than the Farnebäck algorithm and the Dual TV-L1 algorithm, with 50% and 30% improvement in accuracy and 5% and 8% improvement in immunity to interference, respectively. This work is important for promoting the engineering applications of remote sensing instruments for environmental monitoring.

Automatic Estimation of Daily Volcanic Sulfur Dioxide Gas Flux From TROPOMI Satellite Observations: Application to Etna and Piton de la Fournaise

AbstractUnderstanding the dynamics of sulfur dioxide (SO2) degassing is of primary importance for tracking temporal variations in volcanic activity. Here we introduce the novel “disk method,” which aims at estimating the daily volcanic SO2 mass flux from satellite images (such as those provided by Sentinel‐5P/TROPOspheric Monitoring Instrument [TROPOMI]). The method calculates a “proto‐flux” using a regression, as a function of distance, of SO2 mass integrated in a series of nested circular domains centered on a volcano. After regression, a single multiplication by plume speed suffices to deduce the SO2 mass flux, without requiring a subsequent regression. This way, a range of plume speed and plume altitude scenarios can be easily explored. Noise level in the image is simultaneously evaluated by the regression, which allows for estimating posterior uncertainties on SO2 flux and improving the level of detection for weak sources in noisy environments. A statistical test is also introduced to automatically detect occurrences of volcanic degassing, lowering the risk of false positives. Application to multi‐year time‐series at Etna (2021) and Piton de la Fournaise (2021–2023) demonstrates (a) a reliable quantification of SO2 emissions across a broad range of degassing styles (from passive degassing to effusive or paroxysmal events), and (b) a reasonable day‐to‐day correlation between SO2 flux and seismic energy. The method is distributed as an open‐source software, and is implemented in an interactive web application within the “Volcano Space Observatory Portal,” facilitating near‐real‐time exploitation of the TROPOMI archive for both volcano monitoring and assessment of volcanic atmospheric hazards.

Contribution of emissions from the oil sands activities in Alberta, Canada to atmospheric concentration and deposition of nitrogen and sulfur species at a downwind site

Oil sands activities in the Athabasca Oil Sands Region in Alberta, Canada, are large sources of atmospheric NOx and SO2. This study investigated the impact of oil sands emissions on the atmospheric deposition of nitrogen and sulfur species at a downwind site, about 350 km from the oil sands facilities. Measurement data are from the Canadian Air and Precipitation Monitoring Network (CAPMoN) from 2015 to 2019, including ambient concentrations of HNO3, pNO3-, NO2, pNH4+, NH3, SO2, pSO42− and base cations, as well as concentrations of NO3−, SO42−, NH4+, and base cations in precipitation. Sector analysis of air mass back trajectories was conducted to distinguish measurements with different air mass origins. Median atmospheric concentrations and dry deposition fluxes of HNO3, pNO3-, NO2, pNH4+, pSO42−, and SO2 on days when the air masses came from the oil sands sector were significantly greater than those with the “Clean” sector by 34–67%, whereas the difference in NH3 concentration was not significant. Contributions of the oil sands emissions to dry deposition fluxes of these species ranged from 3.8 to 13.1%. The precipitation-weighted mean concentrations of NO3−, SO42−, and NH4+ in samples with the oil sands sector were 76 %, 65 % and 81 % greater than those with the “Clean” sector, respectively. Contributions of the oil sands emissions to wet deposition of NO3−, SO42−, and NH4+ were 12.5 ± 8.9 %, 8.7 ± 4.4 %, and 6.0 ± 3.3 %, respectively. The annual total deposition of nitrogen and sulfur were 1.9 kg-N ha−1 and 0.74 kg-S ha−1, respectively, of which 8.0 ± 3.5 % and 8.7 ± 3.6 % were from oil sands emissions. The total deposition of sulfur and nitrogen did not exceed the critical loads (CL) of acidity, but nitrogen deposition exceeded the CLs of nutrient nitrogen in the region.

Estimating emission flux of H2S from fumarolic fields using vertical sensor array system

The emission flux of volatiles from each fumarolic field in volcanic and geothermal areas can be used to evaluate the current state of magmatic activity and predict its future trends. The emission flux of SO2 has been quantified in many fumarolic fields using remote sensing techniques, such as differential optical absorption spectroscopy (DOAS). However, most of these remote sensing techniques are inapplicable to fumarolic fields emitting volatiles depleted in SO2 to which most of the geothermal fields are classified. In this study, we developed a vertical sensor array system to quantify the emission flux of H2S from each fumarolic field by integrating the cross-sectional distributions of H2S concentrations in the volcanic plume using the vertical sensor array system. In Iwo-yama of the Kirishima volcanic complex, the cross-sectional distribution of H2S concentrations was determined using the walking traverse method by moving the vertical sensor array system in the plume perpendicular to the direction of plume transport. The emission flux of SO2 (2.2 ± 0.4 ton SO2/day) was estimated from that of H2S using the walking traverse method (2.6 ± 0.5 ton H2S/day) and the molar ratio of the plume (SO2/H2S=0.45) corresponds well with that estimated optically by JMA. We concluded that the emission flux quantified using the vertical sensor array system was reliable. In the Oyunuma pond in the Kuttara volcano, the emission flux of H2S was quantified as 2.0 ton H2S/day through the fixed point method, wherein the vertical sensor array system was fixed in one point, whereas the cross sectional distribution of H2S in the plume was estimated using the natural variation in wind direction. The topography is often irregular and wind direction is variable in most fumarolic fields; thus, in general, the fixed point method should be more suitable to determine the emission flux of H2S from fumarolic fields, wherein H2S occupies a major portion of the total sulfur emission.

TROPOMI/PlumeTraj SO2 fluxes consistent with partially degassed magma supplying the 2018 Kīlauea eruption

Volcanic sulfur dioxide (SO2) emission measurements are a key element of volcano monitoring strategies and underpin our understanding of magmatic degassing impacts on air quality, aircraft, and climate. Volcanic SO2 emission monitoring is usually performed with near-field ground-based measurements, but ever-improving satellite observations increasingly offers the capacity for volcanic emission monitoring from space. Here we examine the May–August 2018 eruption at Kīlauea, Hawai'i, which produced a voluminous low-altitude gas plume. We compare SO2 emissions calculated with three approaches: ground-based networks, dissolved sulfur (S) contents and lava effusion rates, and new satellite-derived SO2 data. The high emission rates of this eruption posed a major challenge for near-field ground-based measurements, and corrections were needed to account for attenuation of the SO2 signal from the optically thick plumes. Our satellite-derived gas emissions use Sentinel-5 Precursor (S5P)/Tropospheric Monitoring Instrument (TROPOMI) observations and PlumeTraj back-trajectory analysis, deriving a total mass of 1.3 to 3.3 Mt of SO2 compared to 10.2 Mt from ground-measured fluxes (Kern et al., 2020). Our measurements deal with uncertainties such as plume aging and SO2 oxidation as well as cloud coverage which can underestimate results. The S contents also hold their own uncertainties partly due to when sampling was possible. We find agreement between satellite-derived fluxes and those that would be produced by magma which had previously lost at least 200–750 ppm (16–85%) of its initial S content during residence in the Halema'uma'u lava lake and reservoir. Our measurements allow an estimate of the efficiency of S loss from the Halema'uma'u lava lake prior to the 2018 eruption, and this implies a lava lake magma supply rate of 0.0019–0.0072 km3 per day for 2013–2018. The 2018 eruption produced a lava volume equivalent to 5 to 19 months of magma supply to Kīlauea during 3 months of eruption.

SO2 emission rates and incorporation into the air pollution dispersion forecast during the 2021 eruption of Fagradalsfjall, Iceland

During the low-effusion rate Fagradalsfjall eruption (19 March – 18 September 2021), the emission of sulfur dioxide (SO2) was frequently measured using ground-based UV spectrometers. The total SO2 emitted during the entire eruption was 970 ± 540 kt, which is only about 6% of the SO2 emitted during the similar length Holuhraun eruption (2014–2015). The eruption was divided into five phases based on visual observations, including the number of active vents and the occurrence of lava fountaining. The SO2 emission rate ranged from 44 ± 19 kg/s in Phase 2 to 85 ± 29 kg/s in Phase 5, with an average of 64 ± 34 kg/s for the entire eruption. There was notable variability in SO2 on short timescales, with measurements on 11 August 2021 ranging from 17 to 78 kg/s.SO2 flux measurements were made using scanning DOAS instruments located at different distances from and orientations relative to the eruption site augmented by traverses. Four hundred and forty-four scan and traverse measurements met quality criteria and were used, along with plume height and meteorological data, to calculate SO2 fluxes while accounting for wind-related uncertainties.A tendency for stronger SO2 flux concurrent with higher amplitude seismic tremor and the occurrence of lava fountaining was observed during Phases 4 and 5 which were characterized by intermittent crater activity including observable effusion of lava and gas release interspersed with long repose times. This tendency was used to refine the calculation of the amount of SO2 emitted during variably vigorous activity. The continuous seismic tremor time series was used to quantify how long during these eruption phases strong/weak activity was exhibited to improve the calculated SO2 flux during these Phases.The total SO2 emissions derived from field measurements align closely with results obtained by combining melt inclusion and groundmass glass analyses with lava effusion rate measurements (910 ± 230 kt SO2). Specifically, utilizing the maximum S content found in evolved melt inclusions and the least remaining S content in accompanying quenched groundmasses provides an identical result between field measurements and the petrological calculations. This suggests that the maximum SO2 release calculated from petrological estimates should be preferentially used to initialize gas dispersion models for basaltic eruptions when other measurements are lacking.During the eruption, the CALPUFF dispersion model was used to forecast ground-level exposure to SO2. The SO2 emission rates measured by DOAS were used as input for the dispersion model, with updates made when a significant change was measured. A detailed analysis of one mid-distance station over the entire eruption shows that the model performed very well at predicting the presence of volcanic SO2 when it was measured. However, it frequently predicted the presence of SO2 that was not measured and the concentrations forecasted had no correlation with the concentrations measured. Various approaches to improve the model forecast were tested, including updating plume height and SO2 flux source terms based on measurements. These approaches did not unambiguously improve the model performance but suggest that improvements might be achieved in more-polluted conditions.

Tracing Magma Migration at Mt. Etna Volcano during 2006–2020, Coupling Remote Sensing of Crater Gas Emissions and Ground Measurement of Soil Gases

The geochemical monitoring of volcanic activity today relies largely on remote sensing, but the combination of this approach together with soil gas monitoring, using the appropriate parameters, is still not widely used. The main purpose of this study was to correlate data from crater gas emissions with flank emissions of soil gases at Mt. Etna volcano from June 2006 to December 2020. Crater SO2 fluxes were measured from fixed stations around the volcano using the DOAS technique and applying a modeled clear-sky spectrum. The SO2/HCl ratio in the crater plume was measured with the OP-FTIR technique from a transportable instrument, using the sun as an IR source. Soil CO2 efflux coupled with the 220Rn/222Rn activity ratio in soil gases (named SGDI) were measured at a fixed monitoring site on the east flank of Etna. All signals acquired were subject both to spectral analysis and to filtering of the periodic signals discovered. All filtered signals revealed changes that were nicely correlated both with other geophysical signals and with volcanic eruptions during the study period. Time lags between parameters were explained in terms of different modes of magma migration and storage inside the volcano before eruptions. A comprehensive dynamic degassing model is presented that allows for a better understanding of magma dynamics in an open-conduit volcano.

Exceptionally low mercury concentrations and fluxes from the 2021 and 2022 eruptions of Fagradalsfjall volcano, Iceland

Mercury (Hg) is naturally released by volcanoes and geothermal systems, but the global flux from these natural sources is highly uncertain due to a lack of direct measurements and uncertainties with upscaling Hg/SO2 mass ratios to estimate Hg fluxes. The 2021 and 2022 eruptions of Fagradalsfjall volcano, southwest Iceland, provided an opportunity to measure Hg concentrations and fluxes from a hotspot/rift system using modern analytical techniques. We measured gaseous Hg and SO2 concentrations in the volcanic plume by near-source drone-based sampling and simultaneous downwind ground-based sampling. Mean Hg/SO2 was an order of magnitude higher at the downwind locations relative to near-source data. This was attributed to the elevated local background Hg at ground level (4.0 ng m−3) likely due to emissions from outgassing lava fields. The background-corrected plume Hg/SO2 mass ratio (5.6 × 10−8) therefore appeared conservative from the near-source to several hundred meters distant, which has important implications for the upscaling of volcanic Hg fluxes based on SO2 measurements. Using this ratio and the total SO2 flux from both eruptions, we estimate the total mass of gaseous Hg released from the 2021 and 2022 Fagradalsfjall eruptions was 46 ± 33 kg, equivalent to a flux of 0.23 ± 0.17 kg d−1. This is the lowest Hg flux estimate in the literature for active open-conduit volcanoes, which range from 0.6 to 12 kg d−1 for other hotspot/rift volcanoes, and 0.5–110 kg d−1 for arc volcanoes. Our results suggest that Icelandic volcanic systems are fed from an especially Hg-poor mantle. Furthermore, we demonstrate that the aerial near-source plume Hg measurement is feasible with a drone-based active sampling configuration that captures all gaseous and particulate Hg species, and recommend this as the preferred method for quantifying volcanic Hg emissions going forward.

Excess degassing drives long-term volcanic unrest at Nevado del Ruiz

This study combines volcanic gas compositions, SO2 flux and satellite thermal data collected at Nevado del Ruiz between 2018 and 2021. We find the Nevado del Ruiz plume to have exhibited relatively steady, high CO2 compositions (avg. CO2/ST ratios of 5.4 ± 1.9) throughout. Our degassing models support that the CO2/ST ratio variability derives from volatile exsolution from andesitic magma stored in the 1–4 km depth range. Separate ascent of CO2-rich gas bubbles through shallow (< 1 km depth), viscous, conduit resident magma causes the observed excess degassing. We infer that degassing of ~ 974 mm3 of shallow (1–4 km) stored magma has sourced the elevated SO2 degassing recorded during 2018–2021 (average flux ~ 1548 t/d). Of this, only < 1 mm3 of magma have been erupted through dome extrusion, highlighting a large imbalance between erupted and degassed magma. Escalating deep CO2 gas flushing, combined with the disruption of passive degassing, through sudden accumulation and pressurization of bubbles due to lithostatic pressure, may accelerate volcanic unrest and eventually lead to a major eruption.

A satellite chronology of plumes from the April 2021 eruption of La Soufrière, St Vincent

Abstract. Satellite instruments play a valuable role in detecting, monitoring and characterising emissions of ash and gas into the atmosphere during volcanic eruptions. This study uses two satellite instruments, the Infrared Atmospheric Sounding Interferometer (IASI) and the Advanced Baseline Imager (ABI), to examine the plumes of ash and sulfur dioxide (SO2) from the April 2021 eruption of La Soufrière, St Vincent. The frequent ABI data have been used to construct a 14 d chronology of a series of explosive events at La Soufrière, which is then complemented by measurements of SO2 from IASI, which is able to track the plume as it is transported around the globe. A minimum of 35 eruptive events were identified using true, false and brightness temperature difference maps produced with the ABI data. The high temporal resolution images were used to identify the approximate start and end times, as well as the duration and characteristics of each event. From this analysis, four distinct phases within the 14 d eruption have been defined, each consisting of multiple explosive events with similar characteristics: (1) an initial explosive event, (2) a sustained event lasting over 9 h, (3) a pulsatory phase with 25 explosive events in a 65.3 h period and (4) a waning sequence of explosive events. It is likely that the multiple explosive events during the April 2021 eruption contributed to the highly complex plume structure that can be seen in the IASI measurements of the SO2 column amounts and heights. The bulk of the SO2 from the first three phases of the eruption was transported eastwards, which based on the wind direction at the volcano implies that the SO2 was largely in the upper troposphere. Some of the SO2 was carried to the south and west of the volcano, suggesting a smaller emission of the gas into the stratosphere, there being a shift in wind direction around the height of the tropopause. The retrieved SO2 heights show that the plume had multiple layers but was largely concentrated between 13 and 19 km, with the majority of the SO2 being located in the upper troposphere and around the height of the tropopause, with some emission into the stratosphere. An average e-folding time of 6.07±4.74 d was computed based on the IASI SO2 results: similar to other tropical eruptions of this magnitude and height. The SO2 was trackable for several weeks after the eruption and is shown to have circulated the globe, with parts of it reaching as far as 45∘ S and 45∘ N. Using the IASI SO2 measurements, a time series of the total SO2 mass loading was produced, with this peaking on 13 April (descending orbits) at 0.31±0.09 Tg. Converting these mass values to a temporally varying SO2 flux demonstrated that the greatest emission occurred on 10 April with that measurement incorporating SO2 from the second phase of the eruption (sustained emission) and the beginning of the pulsatory phase. The SO2 flux is then shown to fall during the later stages of the eruption: suggesting a reduction in eruptive energy, something also reflected in ash height estimates obtained with the ABI instrument. A total SO2 emission of 0.63±0.5 Tg of SO2 has been derived, although due to limitations associated with the retrieval, particularly in the first few days after the eruption began, this, the retrieved column amounts and the total SO2 mass on each day should be considered minimum estimates. There are a number of similarities between the 1979 and 2021 eruptions at La Soufrière, with both eruptions consisting of a series of explosive events with varied heights and including some emission into the stratosphere. These similarities highlight the importance of in-depth investigations into eruptions and the valuable contribution of satellite data for this purpose; as these studies aid in learning about a volcano's behaviour, which may allow for better preparation for future eruptive activity.

Exceptional eruptive CO2 emissions from intra-plate alkaline magmatism in the Canary volcanic archipelago

Alkaline mafic magmas forming intra-plate oceanic islands are believed to be strongly enriched in CO2 due to low-degree partial melting of enriched mantle sources. However, until now, such CO2 enhancement has not been verified by measuring CO2 degassing during a subaerial eruption. Here, we provide evidence of highly CO2-rich gas emissions during the 86-day 2021 Tajogaite eruption of Cumbre Vieja volcano on La Palma Island, in the Canary archipelago. Our results reveal sustained high plume CO2/SO2 ratios, which, when combined with SO2 fluxes, melt inclusion volatile contents and magma production rates at explosive and effusive vents, imply a magmatic CO2 content of 4.5 ± 1.5 wt%. The amount of CO2 released during the 2021 eruptive activity was 28 ± 14 Mt CO2. Extrapolating to the volume of alkaline mafic magmas forming La Palma alone (estimated as 4000 km3 erupted over 11 Ma), we infer a maximum CO2 emission into the ocean and atmosphere of 1016 moles of CO2, equivalent to 20% of the eruptive CO2 emissions from a large igneous province eruption, suggesting that the formation of the Canary volcanic archipelago produced a CO2 emission of similar magnitude as a large igneous province.

Sulfur origin and flux variations in fumarolic fluids of Vulcano Island, Italy

A sharp increase in volatiles, especially SO2 fluxes from the solfataric plume and diffuse CO2 from the soils of the La Fossa crater area, started in June 2021, and subsequently from the Levante Bay area, suggests renewed unrest at Vulcano Island, Italy. This event has encouraged monitoring activities and stimulated new research activities aimed at understanding the recent evolution of the volcanic system. In this study, the chemical and isotopic composition of fumaroles, thermal waters, and soil gases from the main degassing areas of Vulcano Island with a special focus on sulfur isotopes, are used to investigate the fluid transfer mechanism inside the volcano. Sulfur is one of the most abundant volatile elements present in magmas and volcanic fluids from the La Fossa crater, where it mostly occurs as SO2 and H2S at variable relative concentrations depending on oxygen fugacity and temperature. The isotope composition and the chemical ratio of sulfur species depict a complex hydrothermal-magmatic system. In addition, we utilize the installed SO2 monitoring network that measures the total outgassing of SO2 with the UV-scanning DOAS technique. The SO2 fluxes from the La Fossa crater fumaroles, coupled with the SO2/CO2 and SO2/H2O ratios, were measured to evaluate the total mass of fluids emitted by the shallow plumbing system and its relationship with the status of volcanic activity. Combining the whole chemical composition of fumaroles analyzed with a discrete, direct sampling of high-temperature fumaroles located on the crater summit, the output of discharged water vapor has been estimated (5,768 t·d−1). On the basis of the water output, we estimated the total thermal energy dissipated by the crater during the last enhanced degassing activity (167 MW). This strong and sharp increase in energy observed during the current crisis confirms the long-growing trend in terms of mass and energy recorded in recent decades, which has brought the surface system of Vulcano Island to a critical level that has never been recorded since the last eruptive event of 1888–91.

Environmental and Volcanic Implications of Volatile Output in the Atmosphere of Vulcano Island Detected Using SO2 Plume (2021–23)

The volatiles released by the volcanic structures of the world contribute to natural environmental pollution both during the passive and active degassing stages. The Island of Vulcano is characterized by solfataric degassing mainly localized in the summit part (Fossa crater) and in the peripheral part in the Levante Bay. The normal solfataric degassing (high-temperature fumarolic area of the summit and boiling fluids emitted in the Levante Bay area), established after the last explosive eruption of 1888–90, is periodically interrupted by geochemical crises characterized by anomalous degassing that are attributable to increased volcanic inputs, which determine a sharp increase in the degassing rate. In this work, we have used the data acquired from the INGV (Istituto Nazionale di Geofisica e Vulcanologia) geochemical monitoring networks to identify, evaluate, and monitor the geochemical variations of the extensive parameters, such as the SO2 flux from the volcanic plume (solfataric cloud) and the CO2 flux from the soil in the summit area outside the fumaroles areas. The increase in the flux of volatiles started in June–July 2021 and reached its maximum in November of the same year. In particular, the mean monthly flux of SO2 plume of 22 tons day−1 (t d−1) and of CO2 from the soil of 1570 grams per square meter per day (g m2 d−1) increased during this event up to 89 t d−1 and 11,596 g m2 d−1, respectively, in November 2021. The average annual baseline value of SO2 output was estimated at 7700 t d−1 during normal solfataric activity. Instead, this outgassing increased to 18,000 and 24,000 t d−1 in 2021 and 2022, respectively, indicating that the system is still in an anomalous phase of outgassing and shows no signs of returning to the pre-crisis baseline values. In fact, in the first quarter of 2023, the SO2 output shows average values comparable to those emitted in 2022. Finally, the dispersion maps of SO2 on the island of Vulcano have been produced and have indicated that the areas close to the fumarolic source are characterized by concentrations of SO2 in the atmosphere higher than those permitted by European legislation (40 μg m−3 for 24 h of exposition) on human health.

A SO2 flux study of the Etna volcano 2020–2021 paroxysmal sequences

The persistent open-vent degassing of Mt. Etna is often punctuated by months-long paroxysmal sequences characterized by episodes of violent Strombolian to lava fountaining activity. Understanding these gas-fueled transitions from quiescence to eruption requires routine measurement of gas fluxes. Here, we report SO2 flux measurements, obtained from a permanent UV camera system, collected over a two-year-long period spanning two paroxysmal sequences of Etna’s New South East Crater (NSEC) in December 2020/April 2021 and May/October 2021. In both cases, SO2 flux increased from ≤ 3250 Mg/day during “ordinary” activity to ≥ 4200 Mg/day. We interpret these distinct SO2 degassing regimes in light of seismic and thermal observations and drawing on numerical simulations of sulfur degassing constrained by parental melt sulfur contents in Etna’s hawaiites. We find that initiation of a paroxysmal sequence results from an approximate doubling of the time-averaged rate of magma supply (and degassing) above the sulfur exsolution level (∼150 MPa pressure), to &amp;gt;4 m3/s. This corroborates recent models that argue for the triggering of paroxysmal sequences by escalating supply of volatile-rich magma to a reservoir ∼3–4 km below the summit region. The non-stationary nature of magma flow and volcanic degassing we identify highlights the need for sustained surveillance to characterize long-term atmospheric budgets of volcanic volatiles.

Temporal Variability in Gas Emissions at Bagana Volcano Revealed by Aerial, Ground, and Satellite Observations

AbstractBagana is a remote, highly active volcano, located on Bougainville Island in southeastern Papua New Guinea. The volcano has exhibited sustained and prodigious sulfur dioxide gas emissions in recent decades, accompanied by frequent episodes of lava extrusion. The remote location of Bagana and its persistent activity have made it a valuable case study for satellite observations of active volcanism. This remoteness has also left many features of Bagana relatively unexplored. Here, we present the first measurements of volcanic gas composition, achieved by unoccupied aerial system (UAS) flights through the volcano's summit plume, and a payload comprising a miniaturized MultiGAS. We combine our measurements of the molar CO2/SO2 ratio in the plume with coincident remote sensing measurements (ground‐ and satellite‐based) of SO2 emission rate to compute the first estimate of CO2 flux at Bagana. We report low SO2 and CO2 fluxes at Bagana from our fieldwork in September 2019, ∼320 ± 76 td−1 and ∼320 ± 84 td−1, respectively, which we attribute to the volcano's low level of activity at the time of our visit. We use satellite observations to demonstrate that Bagana's activity and emissions behavior are highly variable and advance the argument that such variability is likely an inherent feature of many volcanoes worldwide and yet is inadequately captured by our extant volcanic gas inventories, which are often biased to sporadic measurements. We argue that there is great value in the use of UAS combined with MultiGAS‐type instruments for remote monitoring of gas emissions from other inaccessible volcanoes.

Temporal evolution of the magmatic system at Nevado del Ruiz Volcano (Colombia) inferred from long-period seismic events in 2003–2020

Long-period (LP) seismic events have been broadly assumed to be generated by oscillations of crack fluid-filled resonators. We investigated the waveforms and spectral characteristics of LP seismic events recorded by five stations at Nevado del Ruiz Volcano, Colombia. The observed event types of LP events exhibited decaying harmonic oscillations with either monochromatic or polychromatic spectral features (NLP events) or events characterized by non-decaying harmonic oscillations (BLP events). The NLP events occurred over four periods: September 2003 (period I), January 2004–October 2004 (period II), October 2010–March 2012 (period III), and October 2019–April 2020 (period IV). Their source locations were determined to be within the first kilometer of the conduit placed below the Arenas crater. Their locations involved a distinct resonator for each of the four periods but a single resonator for all the events during a single period. Assuming a crack geometry at their sources and based on their locations, we considered a dusty gas (mixture of ash and water vapor) filling a crack resonator. We estimated the crack model parameters to investigate the temporal processes generating variations in frequency and quality factor Q. Then, the estimated parameters indicate that NLP events were triggered by magma fragmentation. We found that events from period III preceded phreatomagmatic eruptions, increased ash and gas emissions and SO2 flux, and increased overall seismic activity. Three years after this period III, a dome was emplaced at the base of the Arenas crater. Besides, events from period IV also preceded by ~3 years a significant increase in the overall seismic activity, which led to a change in alert activity at the volcano. We argue that BLP events were driven by resonances provoked by the unstable degassing through the uppermost magma plug's open cracks (tuffisite veins).In contrast, NLP events occurred by overpressure released through fractures along the conduit's walls (shear fractures). The shear fractures at the conduit wall were in a region of high strain rate where a decrease in confining pressure triggered fragmentation, resulting in dusty gas and, subsequently, a resonance of the crack. NLP and BLP events were in a narrow relationship with the ascent of vesiculated magma to the surface. Our results show that the LP events help diagnose the state of magmatic fluids and understand the eruption processes of Nevado del Ruiz Volcano. LP events are not always followed by explosive eruptions, but they do indicate significant overpressure developing in the shallow conduit.