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

Journal article, 2024

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.