Volcanic air pollution and human health: recent advances and future directions

Volcanic air pollution from both explosive and effusive activity can affect large populations as far as thousands of kilometers away from the source, for days to decades or even centuries. Here, we summarize key advances and prospects in the assessment of health hazards, effects, risk, and management. Recent advances include standardized ash assessment methods to characterize the multiple physicochemical characteristics that might influence toxicity; the rise of community-based air quality monitoring networks using low-cost gas and particulate sensors; the development of forecasting methods for ground-level concentrations and associated public advisories; the development of risk and impact assessment methods to explore health consequences of future eruptions; and the development of evidence-based, locally specific measures for health protection. However, it remains problematic that the health effects of many major and sometimes long-duration eruptions near large populations have gone completely unmonitored. Similarly, effects of prolonged degassing on exposed populations have received very little attention relative to explosive eruptions. Furthermore, very few studies have longitudinally followed populations chronically exposed to volcanic emissions; thus, knowledge gaps remain about whether chronic exposures can trigger development of potentially fatal diseases. Instigating such studies will be facilitated by continued co-development of standardized protocols, supporting local study teams and procuring equipment, funding, and ethical permissions. Relationship building between visiting researchers and host country academic, observatory, and agency partners is vital and can, in turn, support the effective communication of health impacts of volcanic air pollution to populations, health practitioners, and emergency managers.

Methodology to obtain highly resolved SO₂ vertical profiles for representation of volcanic emissions in climate models

In this study we describe a methodology to create high-vertical-resolution SO2 profiles from volcanic emissions. We demonstrate the method's performance for the volcanic clouds following the eruption of Sarychev in June 2009. The resulting profiles are based on a combination of satellite SO2 and aerosol retrievals together with trajectory modelling. We use satellite-based measurements, namely lidar backscattering profiles from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) satellite instrument, to create vertical profiles for SO2 swaths from the Atmospheric Infrared Sounder (AIRS) aboard the Aqua satellite. Vertical profiles are created by transporting the air containing volcanic aerosol seen in CALIOP observations using the FLEXible PARTicle dispersion model (FLEXPART) while preserving the high vertical resolution using the potential temperatures from the MERRA-2 (Modern-Era Retrospective analysis for Research and Application) meteorological data for the original CALIOP swaths. For the Sarychev eruption, air tracers from 75 CALIOP swaths within 9 d after the eruption are transported forwards and backwards and then combined at a point in time when AIRS swaths cover the complete volcanic SO2 cloud. Our method creates vertical distributions for column density observations of SO2 for individual AIRS swaths, using height information from multiple CALIOP swaths. The resulting dataset gives insight into the height distribution in the different sub-clouds of SO2 within the stratosphere. We have compiled a gridded high-vertical-resolution SO2 inventory that can be used in Earth system models, with a vertical resolution of 1 K in potential temperature, 61 ± 56 m, or 1.8 ± 2.9 mbar.

Volcanic origin of the mercury anomalies at the Cretaceous-Paleogene transition of Bidart, France

The timing and mechanisms of the climatic and environmental perturbations induced by the emplacement of the Deccan Traps large igneous province (India) and their contribution to the Cretaceous-Paleogene (K-Pg) mass extinction are still debated. In many marine sediment archives, mercury (Hg) enrichments straddling the K-Pg boundary have been interpreted as the signature of Deccan Traps volcanism, but Hg may also have been derived from the Chicxulub (Mexico) impact. We investigated the Hg isotope composition, as well as the behavior of iridium (Ir) and other trace elements, in K-Pg sediments from the Bidart section in southwest France. Above the K-Pg boundary, Ir content gradually decreases to background values in the Danian carbonates, which is interpreted to indicate the erosion and redistribution of Ir-rich fallouts. No significant enrichment in Ir and W, or Zn and Cu, is observed just below the K-Pg boundary, excluding the hypothesis of downward remobilization of Hg from the boundary clay layer. Positive Δ199Hg and slightly negative values in the upper Maastrichtian and lower part of the early Danian are consistent with the signature of sediments supplied by atmospheric Hg2+ deposition and volcanic emissions. Up section, large shifts to strongly negative mass-dependent fractionation values (δ202Hg) result from the remobilization of Hg formerly sourced by the impactor or by a mixture of different sources including biomass burning, volcanic eruption, and asteroid impact, requiring further investigation. Our results provide additional support for the interpretation that the largest eruptions of the Deccan Traps began just before, and encompassed, the K-Pg boundary and therefore may have contributed to the K-Pg mass extinction.

Validity of Geochemical Signatures of Abiotic Hydrocarbon Gases on Earth

Abiotic synthesis was hypothesized for the occurrences of hydrocarbon gases with atypical molecular and isotopic compositions. This work provides biotic interpretations on these atypical compositions: 1) microbial CH4 oxidation and CO2-sourced methanogenesis may enrich 13C in hydrothermal CH4; 2) microbial hydrocarbon generation using serpentinization-derived H2 may deplete deuterium in hydrocarbons; 3) three processes may cause isotopic reversal with the carbon number in biotic hydrocarbons: i. the decrease of kinetic isotope effect (KIE) with the increase of carbon number during alkane biodegradation, ii. inverse KIE during the thermal decomposition of higher alkanes, and iii. isotopic fractionation during gas diffusion in rock samples; 4) random scission of long alkyl chains may form the exponential distribution of alkane abundance with respect to carbon number (“Schulz-Flory distribution”); 5) isotopic compositions are often not equilibrated; even if they are, the equilibrium temperatures are not necessarily the same as the temperature of hydrocarbon generation. Case studies demonstrate that previously proposed abiotic hydrocarbon gases in continental serpentinite-hosted seepages, continental or oceanic hydrothermal vents, volcanic emissions, gas fields in volcanic reservoirs, and fluid inclusions in alkaline or granitic rocks were formed with various biotic processes. The occurrence of abiotic hydrocarbon gas with decisive evidence is limited to fluid inclusions in some mantle minerals.Supplementary material: Appendices 1 and 2 and tables S1 and S2 are available at https://doi.org/10.6084/m9.figshare.c.5660278.

The CAMS volcanic forecasting system utilizing near-real time data assimilation of S5P/TROPOMI SO₂ retrievals

The Copernicus Atmosphere Monitoring Service (CAMS), operated by the European Centre for Medium-Range Weather Forecasts on behalf of the European Commission, provides daily analyses and 5-day forecasts of atmospheric composition, including forecasts of volcanic sulphur dioxide (SO2) in near-real time. CAMS currently assimilates total column SO2 retrievals from the GOME-2 instruments on MetOp-B and -C and the TROPOMI instrument on Sentinel-5P which give information about the location and strength of volcanic plumes. However, the operational TROPOMI and GOME-2 retrievals do not provide any information about the height of the volcanic plumes and therefore some prior assumptions need to be made in the CAMS data assimilation system about where to place the resulting SO2 increments in the vertical. In the current operational CAMS configuration, the SO2 increments are placed in the mid-troposphere, around 550 hPa or 5 km. While this gives good results for the majority of volcanic emissions, it will clearly be wrong for eruptions that inject SO2 at very different altitudes, in particular exceptional events where part of the SO2 plume reaches the stratosphere. A new algorithm, developed by DLR for GOME-2 and TROPOMI and optimized in the frame of the ESA-funded Sentinel-5P Innovation–SO2 Layer Height Project, the Full-Physics Inverse Learning Machine (FP_ILM) algorithm, retrieves SO2 layer height from TROPOMI in NRT in addition to the SO2 column. CAMS is testing the assimilation of these data, making use of the NRT layer height information to place the SO2 increments at a retrieved altitude. Assimilation tests with the TROPOMI SO2 layer height data for the Raikoke eruption in June 2019 show that the resulting CAMS SO2 plume heights agree better with IASI plume height retrievals than operational CAMS runs without the TROPOMI SO2 layer height information and that making use of the additional layer height information leads to improved SO2 forecasts than when using the operational CAMS configuration. By assimilating the SO2 layer height data the CAMS system can predict the overall location of the Raikoke SO2 plume up to 5 days in advance for about 20 days after the initial eruption.

Efficient release of bromine by super-eruptions

Bromine is a key halogen element in the quantification of volcanic volatiles, but analytical difficulties in measuring its very low abundances have prevented progress in understanding its behavior and its role in volcanic emissions. We present a new data set of bromine, chlorine, and fluorine concentrations in melt inclusions and matrix glasses for two rhyolitic super-eruptions from the Toledo and Valles calderas, New Mexico, USA. We show that before eruption, Br and Cl were efficiently partitioned from the gas-saturated magma into a separate fluid phase, and we calculate the mass of halogens in the fluid phase. We further demonstrate that syn-eruptive magma degassing was negligible during the super-eruptions, so that the main source of halogen emissions must have been the fluid phase. If the fluid phase were erupted, the large mass of Br and Cl could have severely impacted the atmospheric chemistry upon eruption.

Modeling study of the impact of SO₂ volcanic passive emissions on the tropospheric sulfur budget

Well constrained volcanic emissions inventories in chemistry transport models are necessary to study the impacts induced by these sources on the tropospheric sulfur composition and on sulfur species concentrations and depositions at the surface. In this paper, the changes induced by the update of the volcanic sulfur emissions inventory are studied using the global chemistry transport model MOCAGE (MOdèle de Chimie Atmosphérique à Grande Échelle). Unlike the previous inventory (Andres and Kasgnoc, 1998), the updated one (Carn et al., 2016, 2017) uses more accurate information and includes contributions from both passive degassing and eruptive emissions. Eruptions are provided as daily total amounts of sulfur dioxide (SO2) emitted by volcanoes in the Carn et al. (2016, 2017) inventories, and degassing emissions are provided as annual averages with the related mean annual uncertainties of those emissions by volcano. Information on plume altitudes is also available and has been used in the model. We chose to analyze the year 2013, for which only a negligible amount of eruptive volcanic SO2 emissions is reported, allowing us to focus the study on the impact of passive degassing emissions on the tropospheric sulfur budget. An evaluation against the Ozone Monitoring Instrument (OMI) SO2 total column and MODIS (Moderate-Resolution Imaging Spectroradiometer) aerosol optical depth (AOD) observations shows the improvements of the model results with the updated inventory. Because the global volcanic SO2 flux changes from 13 Tg yr−1 in Andres and Kasgnoc (1998) to 23.6 Tg yr−1 in Carn et al. (2016, 2017), significant differences appear in the global sulfur budget, mainly in the free troposphere and in the tropics. Even though volcanic SO2 emissions represent 15 % of the total annual sulfur emissions, the volcanic contribution to the tropospheric sulfate aerosol burden is 25 %, which is due to the higher altitude of emissions from volcanoes. Moreover, a sensitivity study on passive degassing emissions, using the annual uncertainties of emissions per volcano, also confirmed the nonlinear link between tropospheric sulfur species content with respect to volcanic SO2 emissions. This study highlights the need for accurate estimates of volcanic sources in chemistry transport models in order to properly simulate tropospheric sulfur species.

Abundance and Structure of the Zooplankton Community During a Post-eruptive Process: The Case of the Submarine Volcano Tagoro (El Hierro; Canary Islands), 2013-2018

The mesozooplankton community was analyzed over a 6-year period (2013-2018) during the post-eruptive stage of the submarine volcano Tagoro, located south of the island of El Hierro (Canary Archipelago, Spain). Nine cruises from March 2013 to March 2018 were carried out in two different seasons, spring (March-April) and autumn (October). A high-resolution study was carried out across the main cones of Tagoro volcano, as well as a large number of reference stations surrounding El Hierro (unaffected by the volcano). The zooplankton community at the reference stations showed a high similarity with more than 85% of the variation in abundance and composition attributable to seasonal differences. Moreover, our data showed an increase in zooplankton abundance in waters affected by the volcano with a higher presence of non-calanoid copepods and a decline in the diversity of the copepod community, indicating that volcanic inputs have a significant effect on these organisms. Fourteen different zooplankton groups were found but copepods were dominant (79%) with 59 genera and 170 species identified. Despite the high species number, less than 30 presented a larger abundance than 1%. Oncaea and Clausocalanus were the most abundant genera followed by Oithona and Paracalanus (60%). Nine species dominated (>2%): O. media, O. plumifera, and O. setigera among the non-calanoids and M. clausi, P. nanus, P. parvus, C. furcatus, C. arcuicornis, and N. minor among the calanoids. After the initial low abundance of the copepods as a consequence of the eruption, an increase was observed in the last years of the study, where besides the small Paracalanus and Clausocalanus, the Cyclopoids seem to have a good adaptive strategy to the new water conditions. The increase in zooplankton abundance and the decline in the copepod diversity in the area affected by the volcano indicate that important changes in the composition of the zooplankton community have occurred. The effect of the volcanic emissions on the different copepods was more evident in spring when the water was cooler and the mixing layer was deeper. Further and longer research is recommended to monitor the zooplankton community in the natural laboratory of the Tagoro submarine volcano.