# Modest volcanic SO2 emissions from the Indonesian archipelago

Philipson Bani, Clive Oppenheimer, Vitchko Tsanev, Bruno Scaillet, Sofyan Primulyana, Ugan Boyson Saing, Hilma Alfianti, Mita Marlia

► To cite this version:

Philipson Bani, Clive Oppenheimer, Vitchko Tsanev, Bruno Scaillet, Sofyan Primulyana, et al., Modest volcanic SO2 emissions from the Indonesian archipelago. Nature Communications, 2022, 13 (1), pp.3366. 10.1038/s41467-022-31043-7. insu-03696475

HAL Id: insu-03696475

https://insu.hal.science/insu-03696475

Submitted on 16 Jun 2022

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

ARTICLE

https://doi.org/10.1038/s41467-022-31043-7

OPEN

<figure><img src="image_3.png" /></figure>

Check for updates

# Modest volcanic SO₂ emissions from the Indonesian archipelago

Philipson Bani<fnref n="1" />,²✉, Clive Oppenheimer<fnref n="1" />,³, Vitchko Tsanev³, Bruno Scaillet<fnref n="1" />,⁴, Sofyan Primulyana⁵, Ugan Boyson Saing⁵, Hilma Alfianti⁵ & Mita Marlia⁵

Indonesia hosts the largest number of active volcanoes, several of which are renowned for climate-changing historical eruptions. This pedigree might suggest a substantial fraction of global volcanic sulfur emissions from Indonesia and are intrinsically driven by sulfur-rich magmas. However, a paucity of observations has hampered evaluation of these points—many volcanoes have hitherto not been subject to emissions measurements. Here we report new gas measurements from Indonesian volcanoes. The combined SO₂ output amounts to 1.15 ± 0.48 Tg/yr. We estimate an additional time-averaged SO₂ yield of 0.12-0.54 Tg/yr for explosive eruptions, indicating a total SO₂ inventory of 1.27-1.69 Tg/yr for Indonesian. This is comparatively modest—individual volcanoes such as Etna have sustained higher fluxes. To understand this paradox, we compare the geodynamic, petrologic, magma dynamical and shallow magmatic-hydrothermal processes that influence the sulfur transfer to the atmosphere. Results reinforce the idea that sulfur-rich eruptions reflect long-term accumulation of volatiles in the reservoirs.

<footnote n="1">¹Laboratoire Magmas et Volcans, Université Blaise Pascal-CNRS-IRD, OPGC, 63170 Aubière, France. ²Centre IRD de la Nouvelle-Calédonie, 101, Promenade Roger Laroque, BP A5, 98848 Nouméa Cedex, Nouvelle-Calédonie, France. ³Department of Geography, University of Cambridge, Downing Place, Cambridge CB2 3EN, UK. ⁴Institut des sciences de la Terre d'Orléans, Université d'Orléans-CNRS-BRGM, 1a rue de la Férollerie, 45071 Orléans, France. ⁵Center for Volcanology and Geological Hazard Mitigation, Jl. Diponegoro No. 57, Bandung 40122, Indonesia. ✉email: philipson.bani@ird.fr</footnote>
While not the most abundant species in volcanic gases, sulfur dioxide ($SO_2$) is the easiest to measure remotely with the aim of deriving a flux. This owes principally to its absorption of ultraviolet (UV) light, enabling daytime spectroscopic measurements from the ground, air and space¹. Measurements of $SO_2$ flux are a cornerstone of volcano monitoring and contribute to the understanding volcanic degassing. They permit the calculation of fluxes of other volcanic gas species ($X$) from measurements of their ratios to sulfur dioxide ($X/SO_2$) and underpin global inventories of volcanic gas emissions to the atmosphere. In this respect, $SO_2$ is particularly important given its roles in atmospheric chemistry and radiation².

Advances in satellite remote sensing of $SO_2$ in both ultraviolet and infrared wavebands are adding to our knowledge³,⁴ and the proliferation of compact UV spectrometers and cameras⁵⁻⁷, is enabling measurements at less accessible volcanoes. However, the compilation of global inventories still faces numerous challenges, including temporal and spatial data gaps, measurement uncertainties, the presence of multiple sulfur species (including S, $H_2S$ and $H_2SO_4$) in volcanic emissions, and the challenges of processing large volumes of data.

One notable lacuna in $SO_2$ inventories is the Indonesian archipelago. According to Siebert et al.⁸, there are 78 historically active volcanoes in Indonesia, i.e., those with at least one historically-recorded eruption. But such a definition finds its limit in Indonesia where the documentary record is incomplete and traditional knowledge lost or not fully integrated into scientific records. A more complete inventory of the Indonesian volcanoes can be found in the “Badan Geologi” database⁹ which lists 126 active volcanoes (including six submarine edifices). They are subdivided into 77 type-A volcanoes, which have experienced at least one increase in magmatic and/or phreatic activity since 1600; 29 type-B volcanoes with solfataric and/or fumarolic manifestations but no eruption since 1600; and 20 type-C, which are solfataras and/or fumarole fields lacking a defined volcanic edifice (Fig. 1 and Table 2). These volcanoes fall within four distinct arcs: Sunda, Banda, Sangihe, and Halmahara (Fig. 1).

In the first global compilations of the volcanic $SO_2$ budget, the Indonesian contribution was unspecified¹⁰. Subsequently, over four decades, new observations have furnished estimates of annual $SO_2$ inventories for Indonesia (Table 1). These have varied considerably, beginning with the work of Le Guern¹¹, who estimated 0.07 Tg/yr (representing just 0.15% of the global volcanic $SO_2$ budget) compared with a more recent estimate of Carn et al.³ of 2.2 Tg $SO_2$/yr (representing 9.5% of the global volcanic budget). Note that these and intervening works have also reported disparate figures for the global total (Table 1), which is unsurprising given the different datasets (and their timespans) and methods employed. Despite developments in $SO_2$ sensing, hitherto only a fraction of the more than 100 Indonesian volcanoes classified as active has been subject to $SO_2$ flux measurement campaigns.

Here we present a new inventory of volcanic $SO_2$ emissions for Indonesia based on portable ground-based remote sensing instruments and systematic program of fieldwork observations. We focus our efforts on the subaerial type-A volcanoes with passive degassing, which we consider based on field observation to be the main volcanic degassing sources in Indonesia (Table 2 and Fig. 1). We use the term ‘passive’ to refer to the style of gas emission so as to distinguish it from larger, sporadic explosive emissions, though the term can encompass a wide range of sources from magmatic to fumarolic. We evaluate the factors influencing the variations identified between volcanoes and between sub-regions of Indonesia and consider the total $SO_2$ emission rate for the archipelago in the global context.

## Results

**$SO_2$ emission budget.** Of the 73 aerial type-A volcanoes across Indonesia, we conducted measurements at 47 (Fig. 1 and Table 3), including 12 that exhibit negligible $SO_2$ release. Of the remaining 26 volcanoes that were not visited, 20 are either inactive or exhibit negligible $SO_2$ emission, according to local observatory reports and available data and images (https://vsi.esdm.go.id/). There were six volcanoes known for persistent degassing that we did not visit: Banda Api, Serua, Batu Tara, Sangeang Api, Rinjani, and Arjuno Welirang (Fig. 2). However, satellite observations provide some constraints for Batu Tara, Rinjani and Sangeang Api³.

The Sunda arc makes the largest contribution, with a collective daily output of 1313 ± 539 Mg (Fig. 2). Sinabung, Kawah Ijen, Slamet, Anak Krakatau and Bromo volcanoes are the strongest $SO_2$ sources of the arc at 275 ± 25 Mg/d, 238 ± 194 Mg/d, 206 ± 66 Mg/d, 190 ± 77 Mg/d, and 166 ± 2 Mg/d, respectively. Moderate to small $SO_2$ emission rates have been reported for Rinjani (74 ± 65 Mg/d), Sangeang Api (71 ± 75 Mg/d)³, Semeru (48 ± 22 Mg/d)¹², Merapi (20 ± 7 Mg/d), Kerinci (9.8 ± 4.3 Mg/d) and Kaba (9.0 ± 3.1 Mg/d) (Table 3). These 11 volcanoes are the main degassing sources of the Sunda arc. Measurements have also quantified minor to negligible $SO_2$ emission rates (0.2–2.6 Mg/d, Table 3) for five other volcanoes, namely Marapi, Tangkuban Parahau, Papandayan, Talang and Guntur. In total there are 16 volcanic $SO_2$ degassing sources across the Sunda arc. Two of them, Dempo and Kelut, host crater lakes that trap condensable gases, limiting their atmospheric contribution. The remaining 19 type A volcanoes out of a total 37 volcanoes across the Sunda are either quiescent (non-emitters) or characterized by low temperature solfataras and/or fumaroles (Table 3), except Ajurno Welirang, which sustains a persistent but minor degassing.

The Banda arc has the lowest passive $SO_2$ degassing budget with a total daily output of 330 ± 175 Mg. Ili Lewotolo, Sirung, and Lewotobi Perempuan are the main sources but with moderate emission rates corresponding to 75 ± 40 Mg/d, 48 ± 22 Mg/d, 15 ± 10 Mg/d. Batu Tara and Rokatenda with 102 ± 51 and 60 Mg/d, respectively³, are among the main $SO_2$ degassing sources of the Banda arc. Note that the figure for Batu Tara was obtained by subtracting our measurement for Ili Lewotolo from the reported combined flux for both volcanoes³. Lesser $SO_2$ emission rates are found for Iya (8 ± 6 Mg/d), Wurlali, (8 ± 6 Mg/d) and Ebulobo (6 ± 3 Mg/d). We obtain negligible fluxes for Egon (3 ± 2 Mg/d), Kelimutu (2.0 ± 0.7 Mg/d), Lewotobi Lakilaki (2.0 ± 0.7 Mg/d), and Ili Werung (1.0 ± 0.8 Mg/d) (Table 3). Two other volcanoes of the Banda arc, including Serua, and Banda Api, were not visited and therefore their $SO_2$ emissions remain unknown. However, based on information from the local observatories, degassing strength of Banda Api is comparable to that of Wurlali, and exceeds that of Serua. The lack of measurements is thus unlikely to bias significantly our arc-scale flux estimate. The $SO_2$ degassing associated with low temperature soltaras and fumaroles from the remaining eight type A volcanoes of this arc is negligible (Table 3).

The Sangihe arc hosts three volcanoes with relatively strong $SO_2$ emission rates, including Soputan (376 ± 100 Mg/d), Lokon (117 ± 10 Mg/d) and Karangetang (120 ± 55 Mg/d). Awu the northernmost volcano of the arc emits 13 ± 5 Mg $SO_2$/d. The $SO_2$ contribution from the five remaining volcanoes is negligible (Table 4). Hence, with a total $SO_2$ degassing budget of 626 ± 170 Mg per day, the Sangihe arc constitutes a notable arc-scale volcanic degassing source to the atmosphere. Note, however, that Soputan clearly stands out as the strongest source, representing 60% of the total for this arc.

Finally, the $SO_2$ emission rate from the Halmahera arc amounts to 897 ± 437 Mg/d with more than 90% of this flux
<figure><img src="image_4.png" /><figcaption>**Fig. 1** Indonesian active volcanoes. The distribution of the 126 active volcanoes across the archipelago of Indonesia, including 120 aerial and six known submarine edifices (not shown on the map). 77 are classified as Type-A (red triangles), 29 as type-B (yellow squares) and 20 as type-C (green circles) The volcanoes visited in this work are highlighted in red-bold-italic.</figcaption></figure>

accounted for by Dukono which emits 819 ± 394 Mg/d. The other volcanoes are low to moderate sources, with 59 ± 32 Mg/d from Ibu, 16 ± 10 Mg/d from Gamalama and 3.4 ± 1.0 Mg/d from Gamkonora. On Kie Besi, only a small fumarole is present and its SO₂ emission is considered negligible.

Based on the above results, and as summarized in Table 3, the total daily SO₂ emission passively released into the atmosphere from the entire Indonesian archipelago is 3200 ± 1300 Mg/d, equivalent to 1.15 ± 0.48 Tg SO₂ yr⁻¹ (Table 3). We emphasise that this figure is representative of the periods of observations and must be viewed cautiously but we believe it gives a useful guide to the scale of emissions at the scale of the entire archipelago.

**Principal point sources.** Our ranking of Indonesian volcanic sources of passive degassing is shown in Table 3. Dukono (Halmahera arc) is the strongest, representing more than a quarter (26%) of the total. Soputan (Sangihe arc), Sinabung and Kawah Ijen (Sunda arc) are also notable, representing, 12%, 9% and 8% of the total, respectively. These four volcanoes alone constitute around half of the total inventory. They are sustained by different magma compositions, i.e., basaltic andesite to andesite on Sinabung¹³, basaltic to dacite on Kawah Ijen, basaltic on Soputan¹⁴, and andesite to trachyandesite on Dukono¹⁵. Six other volcanoes exhibit moderate rates of SO₂ emission, including Slamet (206 Mg/d), Anak Krakatau (190 Mg/d), Bromo (166 Mg/d), Karangetang (120 Mg/d), Lokon (117 Mg/d) and Batu Tara (102 Mg/d). Thus, ten volcanoes contribute 82% of the total passive volcanic SO₂ emission budget of Indonesia. There are five other volcanoes with modest SO₂ fluxes,

between 50 and 100 Mg/d that together represent 11% of the budget, seven with SO₂ emission rates between 10 and 50 Mg/d, representing 5% of the budget, and finally 14 volcanoes whose SO₂ degassing is below 10 Mg/d.

**Arc scale variations.** Our new SO₂ inventory reveals substantial variations in SO₂ output between the arcs of the Indonesian archipelago. The 3000-km-long Sunda arc, with 37 type A volcanoes, is the largest SO₂ source at 0.48 ± 0.20 Tg/yr representing 41% of the Indonesian total. The 2000-km-long Banda arc, in contrast, contributes just 0.12 ± 0.06 Tg/yr representing only 10% of the total, despite hosting 24 Type A volcanoes. The 600-km-long Sangihe and 500-km-long Halmahera arcs are stronger sources (0.23 ± 0.06 Tg and 0.33 ± 0.16 Tg SO₂ yr⁻¹, respectively) despite their shorter extents and comparatively few volcanoes (eleven and five, respectively). It is thus plausible that the geodynamic contexts play a key role in the SO₂ emission budget of Indonesia; in terms of SO₂ emission rate per km of arc, the Halmahera arc is the strongest source with an output of 655 Mg SO₂ per km yr⁻¹ followed by the Sangihe, Sunda and Banda arcs with 380, 160 and 60 Mg SO₂ per km yr⁻¹, respectively.

**Passive and explosive degassing.** Our total Indonesian SO₂ inventory of 1.15 ± 0.48 Tg/yr based on ground-based and airborne surveys of persistent volcanic degassing across the archipelago is half the estimated emission from the 20 Indonesian volcanoes reported in ref. ³, derived from satellite remote sensing measurements. But this latter approach focused on a different
<table><caption>Table 1 Reported global volcanic SO₂ inventories and contribution from Indonesian volcanoes.</caption><thead><tr><td>Authors</td><td>Total volcanic SO₂ (Tg/yr)</td><td>Method(s)</td><td>Contribution from Indonesian volcanoes</td></tr></thead><tbody><tr><td>Le Guern, 1982<sup>11</sup></td><td>50.0</td><td>Using Correlation Spectrometer (COSPEC) data from ref. 50, and extrapolating to a larger number of active volcanoes in different geodynamic provinces.</td><td>0.073 Tg yr<sup>-1</sup> from Indonesia (Merapi volcano was considered as the main degassing source).</td></tr><tr><td>Spiro et al., 1992<sup>51</sup></td><td>19.2</td><td>Based on plume size, following ref. 52 and referring to volcanism in 1964-1972 and 1980: 28% of the SO₂ emission budget from passive degassing, assuming 61% from eruption to troposphere and 11% to the stratosphere.</td><td>0.41 Tg yr<sup>-1</sup> attributed to Indonesian volcanoes.</td></tr><tr><td>Andres and Kasgnoc, 1998<sup>53</sup></td><td>13.4</td><td>From a compilation of S fluxes in 214 published references, personal comm., and conference presentations. Electronic mail messages were sent to the VOLCANO list for data discussion with volcanologists and atmospheric scientists. Two categories were distinguished: continuously (49) and sporadically erupting (25) volcanoes.</td><td>Indonesian volcanoes contributed to only 0.10 Tg yr<sup>-1</sup> (four volcanoes were considered, including Merapi, Tangkuban parahu, Bromo and Slamet).</td></tr><tr><td>Halmer et al., 2002<sup>54</sup></td><td>15.0-21.0</td><td>Considering the SO₂ emissions of 50 volcanoes recorded by the Total Ozone Mapping Spectrometer (TOMS) and COSPEC, then extrapolated to 310 unmeasured volcanoes based on the VEI-SO₂ relationship, magma composition, tectonic setting and the state of activity.</td><td>2.1-2.6 Tg yr<sup>-1</sup> attributed to Indonesian subduction zone.</td></tr><tr><td>Shinohara, 2013<sup>55</sup></td><td>19.8</td><td>Based on a literature review: 76 persistently degassing volcanoes release an estimated 18.5 Tg/yr of SO₂ and the time-averaged annual SO₂ flux from explosive eruption (1.3 Tg/yr) is obtained based on VEI-SO₂ emission correspondence.</td><td>0.1 Tg from Indonesia, four volcanoes were considered: Merapi, Tangkubanparahu, Slamet and Bromo.</td></tr><tr><td>Carn et al., 2017<sup>3</sup></td><td>23.0</td><td>Based on OMI data spanning 2005-2015 and focused on passive degassing from 91 volcanoes worldwide.</td><td>2.2 Tg from Indonesia. 20 volcanoes were considered: Dukono, Bromo-Semeru, Lewotolo-Batu Tara, Ijen-Raung, Sirung, Sinabung, Karangetang, Krakatau, Kerinci, Slamet, Lokon, Ebulobo, Rinjani, Sangeang api, Paluweh Marapi and Merapi (from the highest to the lowest SO₂ emission).</td></tr></tbody></table>

study period (2005–2015) and integrated the SO₂ contribution from explosive events. Between 2010 and 2020, there were 110 eruptive episodes across the Indonesian archipelago reported in the Global Volcanism Program (https://volcano.si.edu/) and Bandan Geologi (https://vsi.esdm.go.id/). Most of these were minor to moderate in scale (VEI < 3) and their SO₂ emissions were, mostly, not captured by satellite sensors. The year 2014 was the most active year with 14 eruptions while in 2020 there were only six eruptive events reported. The mean value is ten eruptions per year (Table 4 and Fig. 3). The Sunda arc dominates this list with 60 eruptions at 16 different volcanoes. Sinabung and Anak Krakatau were the most active with nine and eight events, respectively. Kerinci and Marapi were also notably active with, respectively, six and seven eruptions, whilst Merapi, Semeru, and Sangeang Api experienced four eruptions each. The Banda arc produced 13 eruptions at seven different volcanoes over the last decade. Batu Tara was the most active with four eruptions. For the Sangihe arc, 19 eruptions were reported at three different volcanoes, including Karangetang, the most active with ten eruptions, Soputan with six eruptions and Lokon with three eruptions. Eighteen eruptive events were recorded for the Halmahera arc. Dukono was the most active volcano with 11 eruptions, followed by Gamalama with six eruptions.

To estimate the SO₂ contribution from these explosive events, we first used a formulation¹⁶ relating volcanic explosivity index (VEI) and SO₂ yield:

$$ \log10(\text{SO}_2, \text{Tg}) = 0.71\text{VEI} - 3.15 \quad (1) $$

We took VEI values reported in the Global Volcanism Program (https://volcano.si.edu/). This indicates a total eruptive SO₂

output over the 2010–2020 period of 5.99 ± 0.31 Tg with annual totals varying between 0.25 Tg and 1.03 Tg (Table 4, Fig. 3). The Sunda arc released 3.7 ± 0.3 Tg, representing 62% of the total, with the main contributions from Sinabung, Merapi, Kelut, Anak Krakatau, and Agung. For the Banda arc we estimate 0.40 ± 0.05 Tg SO₂ (0.04 Tg/yr) accounting for 7% of the budget, mostly contributed by Rokatenda, Lewotolo and Batu Tara. For the Sangihe arc, Soputan, Karangetang and Lokon volcanoes were the only sources with an estimated combined yield of 0.75 ± 0.05 Tg (0.07 Tg/yr) of SO₂ or 12% of the total. Lastly, the Halmahera arc released 1.12 ± 0.01 Tg (0.10 Tg/yr) of SO₂ through eruptions, representing 19% of the total. Dukono was the main contributor, accounting for 83% of the arc's output (Table 4).

We also analyzed available satellite data (http://SO2.gsfc.nasa.gov/) for the SO₂ mass over Indonesia between 2010 and 2020. Out of the 110 eruptions, 71 were captured by the satellite (64%) and the corresponding SO₂ mass amounts to a total of 1.31 ± 0.18 Tg, with a mean annual value of 0.12 ± 0.04 Tg. The Sunda and Halmahera arcs are the main contributors representing, respectively, 81% (1.07 ± 0.18 Tg) and 11% (0.14 ± 0.01 Tg) of the total SO₂. The main contributors from the Sunda arc are Kelut (44%), owing to its 2014 eruption¹⁷, Merapi (14%), primarily related to its 2010 event¹⁸, and Sinabung (12%), which has experienced episodic dome growth and collapse since 2010 (ref. ¹⁹). For the Halmahera arc, 99% of the arc contribution is from Dukono, reflecting its continuous eruptive activity. For the Sangihe arc, the total SO₂ yield amounts to 0.03 ± 0.005 Tg, and Soputan, with its recurrent eruptive activity, is the main source, representing 67% of the arc contribution. Finally, the Banda arc contribution corresponds to 0.068 ± 0.019 Tg and is dominated by the 2020
<table><caption>Table 2 The distribution of the Indonesian active volcanoes per type (A, B or C), region and arcs.</caption><thead><tr><th></th><th>Type-A</th><th>Type-B</th><th>Type-C</th><th></th><th>Type-A</th><th>Type-B</th><th>Type-C</th></tr></thead><tbody><tr><td rowspan="13">Sunda Arc</td><td><strong>Sumatra</strong></td><td>Bur Ni Geureudong</td><td>Pulau Weh</td><td rowspan="13">Banda Arc</td><td><strong>Flores-Lembata-Pantar</strong></td><td>Ili Muda</td><td>Waisano</td></tr><tr><td>Seulawah Agam</td><td>Sibayak</td><td>Gayolesten</td><td>Anak Ranaka</td><td>Labalekan</td><td>Poco Leak</td></tr><tr><td>Peut Sague</td><td>Pusuk Bukit</td><td>Halatobi Tarutung</td><td>Inelika</td><td>Yersey (sub. volc.)</td><td>Ndetu Napi</td></tr><tr><td>Bur Ni Telong</td><td>Bual Buali</td><td>Halubelu</td><td>Inerie</td><td></td><td>Riang Kotang</td></tr><tr><td>Sinabung</td><td>Talak Mau</td><td>Marga Bayur</td><td>Ebulobo</td><td></td><td></td></tr><tr><td>Sorik Marapi</td><td>Kunyit</td><td>Permatang Bata</td><td>Iya</td><td></td><td></td></tr><tr><td>Tadikat</td><td>Bemerang Beriti</td><td></td><td>Kelimutu (Paluweh)</td><td></td><td></td></tr><tr><td>Marapi</td><td>Bukit Daun</td><td></td><td>Rokatenda</td><td></td><td></td></tr><tr><td>Talang</td><td>Lumut Balai</td><td></td><td>Batu Tara</td><td></td><td></td></tr><tr><td>Kirinci</td><td>Sikicau Belirang</td><td></td><td>Egon</td><td></td><td></td></tr><tr><td>Sumbing</td><td>Rajabasa</td><td></td><td>Lewotobi Laki-Laki</td><td></td><td></td></tr><tr><td>Kaba</td><td></td><td></td><td>Lewotobi Perempuan</td><td></td><td></td></tr><tr><td>Dempo</td><td></td><td></td><td>Lereboleng</td><td></td><td></td></tr><tr><td>Anak Krakatau</td><td></td><td></td><td>Ili Boleng</td><td></td><td></td></tr><tr><td><strong>Java</strong></td><td></td><td></td><td>Ili Werung</td><td></td><td></td></tr><tr><td>Salak</td><td>Karang</td><td>Kiaraberes Gagak</td><td>Ili Lewotolo</td><td></td><td></td></tr><tr><td>Gede</td><td>Pulosar</td><td>Perbakti</td><td>Sirung</td><td></td><td></td></tr><tr><td>Tangkuban Parahu</td><td>Patuha</td><td>Kawah Manuk</td><td>Hobal (sub. volc.)</td><td></td><td></td></tr><tr><td>Papandayan</td><td>Wayang Windu</td><td>Kawah Kamojang</td><td><strong>South Maluku</strong></td><td></td><td></td></tr><tr><td>Guntur</td><td>Talaga Bodas</td><td>Kawah Karaha</td><td>Wetar</td><td>Manuk</td><td></td></tr><tr><td>Galunggung</td><td>Ungaran</td><td></td><td>Nieuwerkerk (sub. volc.)</td><td>Emp. China (sub. volc.)</td><td></td></tr><tr><td></td><td></td><td></td><td>Wurlali</td><td></td><td></td></tr><tr><td>Ceremai</td><td>Merbabu</td><td></td><td>Teon</td><td></td><td></td></tr><tr><td>Slamet</td><td>Lawu</td><td></td><td>Nila</td><td></td><td></td></tr><tr><td>Dieng</td><td>Wilis</td><td></td><td>Serua</td><td></td><td></td></tr><tr><td>Sundoro</td><td>Lyang Argapura</td><td></td><td>Banda Api</td><td></td><td></td></tr><tr><td>Sumbing</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Merapi</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Kelut</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ajurno Welirang</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Semeru</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Bromo</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Lamongan</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Raung</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Kawah Ijen</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td><strong>Bali-Lombok-Sumbawa</strong></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Batur</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Agung</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Rinjani</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Tambora</td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Sangeang Api</td><td></td><td></td><td></td><td></td><td></td></tr></tbody></table>

eruption of Ili Lewotolo, representing 92% of the arc contribution. These figures obtained from satellite data are lower than those calculated from the VEIs (Fig. 3 and Table 4), except for Kelut and Ili Lewotolo, where the estimates converge. Despite the discrepancy, particularly in the lower SO₂ estimates, both approaches highlight the Sunda and Halmahera arcs as the main SO₂ contributors to explosive emissions.

Combining the figures we derive for passive and explosive SO₂ degassing yields a total source between 1.27 and 1.69 Tg SO₂ /yr for the Indonesia archipelago. The lower and the higher range integrate respectively the mean annual figures from satellite data and from the VEIs (Table 4). About 10% to 30% of total SO₂ release was sustained by larger, sporadic explosive emissions. This Indonesia total budget corresponds to 3–7% of the global volcanic SO₂ emission budget based on estimates in ref. 3 and 20 (23–33 Tg/yr) and is comparable to the total emissions from Japan²¹

(Fig. 4), although much less when the degassing budget is normalized by arc length.

## Discussion

This work constitutes the first near-comprehensive SO₂ emission survey across the Indonesian archipelago. We estimate a passive degassing flux of 1.15 Tg SO₂/yr for Indonesia in the period of 2010–2019. This represents the cumulative emission from twenty volcanoes, the four strongest sources being Dukono, Soputan, Sinabung and Kawah Ijen, which together represent 54% of the total emission budget, while Slamet, Anak Krakatau, Bromo, Karangetang, and Lokon account for another 25% of the total. More modest sources include Batu Tara, Ili Lewotolo, Rinjani, Sangeang Api, Rokatenda, and Ibu, representing 14% of the total. Seven minor sources, Semeru, Slamet, Merapi, Gamalama,
<table><caption>Table 3 Combined SO₂ flux results for the Indonesian volcanoes.</caption><thead><tr><th>Volcano Name</th><th>long/lat</th><th>Rank</th><th>Typical degassing status</th><th>Average plume height (m)</th><th>Mean SO₂ flux (Mg/d)</th><th>error (Mg/d)</th><th>Method/ Technique</th><th>Source/ (measurement date)</th></tr></thead><tbody><tr><td>Sunda Arc<br>Sinabung</td><td>98.392E/<br>3.170 N</td><td>3</td><td>emission from<br>lava dome</td><td>2600-2800</td><td>275</td><td>25</td><td>DOAS<br>scanner<br>(NOVAC)</td><td>Primulyana et al. 2018<br>(2010-2016)</td></tr><tr><td>Marapi</td><td>100.474E/<br>0.380 S</td><td>29</td><td>Minor degassing<br>from the main<br>central crater</td><td>2800-2900</td><td>2.6</td><td>1.5</td><td>DOAS<br>Walking<br>traverses</td><td>This work (25/10/2014)</td></tr><tr><td>Talang</td><td>100.681E/<br>0.979 S</td><td>35</td><td>Minor degassing<br>from fumaroles<br>formed along a NE-<br>SW fracture<br>transecting<br>the summit</td><td>2400-2500</td><td>0.3</td><td>0.04</td><td>DOAS<br>Walking<br>traverses</td><td>This work (27/06/<br>2012)</td></tr><tr><td>Kerinci</td><td>101.264E/<br>1.697 S</td><td>22</td><td>Emission from a lava<br>dome in the crater<br>and frequent<br>eruptions</td><td>3200-4300</td><td>9.8</td><td>4.3</td><td>DOAS<br>Walking<br>traverses</td><td>This work (03/05/<br>2012)</td></tr><tr><td>Kaba</td><td>102.615E/<br>3.522 S</td><td>23</td><td>Degassing of the<br>extended fumarole<br>zone in the crater</td><td>1700-1800</td><td>9.0</td><td>3.1</td><td>DOAS<br>Walking<br>traverses</td><td>This work (02/06/<br>2015)</td></tr><tr><td>Anak Krakatau</td><td>105.423E/<br>6.102 S</td><td>6</td><td>Outgassing through<br>the summit crater<br>filled by a lava flow</td><td>400-500</td><td>190</td><td>77</td><td>DOAS<br>Airborne<br>traverses</td><td>This work (02/04/<br>2013)</td></tr><tr><td>Tangkuban Parahu</td><td>107.600E/<br>6.770 S</td><td>32</td><td>Minor gas release<br>from the main crater</td><td>1850-1900</td><td>1.8</td><td>0.4</td><td>DOAS<br>Walking<br>traverses</td><td>This work (06/09/<br>2012)</td></tr><tr><td>Papandayan</td><td>107.730E/<br>7.320 S</td><td>33</td><td>Hydrothermal<br>dominated gas<br>through three<br>primary</td><td>2250-2300</td><td>1.4</td><td>0.8</td><td>DOAS<br>Walking<br>traverses</td><td>Bani et al. 2013 (18/06/<br>2011)</td></tr><tr><td>Guntur</td><td>107.841E/<br>7.143 S</td><td>36</td><td>Gas emissions from<br>solfatara</td><td>2100-2150</td><td>0.2</td><td>0.1</td><td>DOAS<br>Walking<br>traverses</td><td>This work (21/10/2012)</td></tr><tr><td>Slamet</td><td>109.208E/<br>7.242 S</td><td>5</td><td>Persistent degassing<br>from the summit<br>crater.</td><td>3300-3400</td><td>206</td><td>66</td><td>OMI</td><td>Carn et al. 2017<br>(2005-2015)</td></tr><tr><td>Merapi</td><td>110.446E/<br>7.540 S</td><td>18</td><td>Degassing from<br>lava dome</td><td>2800-2900</td><td>20</td><td>7</td><td>DOAS<br>scanner</td><td>This work (May-<br>Jun.-Jul. 2016)</td></tr><tr><td>Arijuno Welirang</td><td>112.575E/<br>7.733 S</td><td></td><td>Degassing associated<br>with vigorous<br>fumaroles—not<br>measured</td><td></td><td>?</td><td></td><td></td><td></td></tr><tr><td>Semeru</td><td>112.922E/<br>8.108 S</td><td>16</td><td>Continuous eruptive<br>activity with<br>intermittent<br>strong events</td><td>3500-3700</td><td>48</td><td>22</td><td>UV-Camera</td><td>Smekens et al. 2015<br>(16-22/05/2013: 31/<br>05-03/06/2013)</td></tr><tr><td>Bromo</td><td>112.950E/<br>7.942 S</td><td>7</td><td>Degassing via<br>main crater</td><td>2600-2700</td><td>166</td><td>2</td><td>UV-Camera</td><td>Airuppa et al. 2015<br>(20-21/09/2014)</td></tr><tr><td>Kawah Ijen</td><td>114.242E/<br>8.058 S</td><td>4</td><td>Degassing via a<br>solfatara provoked by<br>mining activity</td><td>2300-2600</td><td>238</td><td>194</td><td>UV-Camera</td><td>This work (12/05/2015)</td></tr></tbody></table>
<table><caption>Table 3 (continued)</caption><thead><tr><th>Volcano Name</th><th>long/lat</th><th>Rank</th><th>Typical degassing status</th><th>Average plume height (m)</th><th>Mean SO<sub>2</sub> flux (Mg/d)</th><th>error (Mg/d)</th><th>Method/ Technique</th><th>Source/ (measurement date)</th></tr></thead><tbody><tr><td>Rinjani</td><td>116.470E/ 8.420 S</td><td>12</td><td>Degassing from intracaldera cone</td><td>2300-2600</td><td>74</td><td>65</td><td>OMI</td><td>Carn et al. 2017 (2005-2015)</td></tr><tr><td>Sangeang Api</td><td>119.070E/ 8.200 S</td><td>13</td><td>Intermittent steam releases</td><td>1700-1800</td><td>71</td><td>75</td><td>OMI</td><td>Carn et al. 2017 (2005-2015)</td></tr><tr><td><strong>Sanghe Arc</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Soputan</td><td>124.737E/ 1.112 N</td><td>2</td><td>Sustained degassing at the summit and frequent eruptions</td><td>1800-2000</td><td>376</td><td>100</td><td>UV-camera</td><td>This work (21/07/2014)</td></tr><tr><td>Lokon</td><td>124.792E/ 1.358 N</td><td>9</td><td>Sustained degassing and frequent eruptions</td><td>1200-1500</td><td>117</td><td>10</td><td>UV-camera</td><td>This work (19/07/2014)</td></tr><tr><td>Karangetang</td><td>125.407E/ 2.781 S</td><td>8</td><td>Degassing associated with lava dome</td><td>1700-1900</td><td>120</td><td>55</td><td>DOAS scanner</td><td>This work (24/07/ 2015)</td></tr><tr><td>Awu</td><td>125.447E/ 3.689 N</td><td>21</td><td>Degassing via crater wall and from small intracrater lava dome</td><td>900-1100</td><td>13</td><td>5</td><td>DOAS scanner</td><td>This work (27/07/ 2015)</td></tr><tr><td><strong>Banda Arc</strong></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Ebulobo</td><td>121.191E/ 8.817 S</td><td>26</td><td>Degassing via summit fumarolic activity</td><td>2100-2300</td><td>6</td><td>3</td><td>UV-camera</td><td>This work (02/10/ 2014)</td></tr><tr><td>Iya</td><td>121.641E/ 8.891 S</td><td>24</td><td>Degassing via fumarole activity</td><td>600-650</td><td>8</td><td>6</td><td>DOAS Walking traverses</td><td>This work (01/10/2014)</td></tr><tr><td>Kelimutu</td><td>121.820E/ 8.770 S</td><td>30</td><td>Negligible SO<sub>2</sub> emission through crater lake</td><td>1550-1600</td><td>2.0</td><td>0.7</td><td>DOAS walking traverse</td><td>This work (30/09/ 2014)</td></tr><tr><td>Rokatenda (Paluweh)</td><td>121.708E/ 8.320 S</td><td>14</td><td>Minor gas release via lava dome</td><td>700-1000</td><td>60</td><td>32</td><td>OMI</td><td>Carn et al. 2017 (2005-2015)</td></tr><tr><td>Egon</td><td>122.455E/ 8.676 S</td><td>28</td><td>Gas release via fumarolic activity at the summit</td><td>1600-1700</td><td>3</td><td>2</td><td>UV-camera</td><td>This work (13/05/2013)</td></tr><tr><td>Lewotobi Lakilaki</td><td>122.768E/ 8.537 S</td><td>31</td><td>Minor gas release from summit fumaroles</td><td>1350-1400</td><td>2.0</td><td>0.7</td><td>UV-camera</td><td>This work (12/05/2013)</td></tr><tr><td>Lewotobi Perempuan</td><td>122.781E/ 8.551 S</td><td>20</td><td>Degassing from a small intracrater dome</td><td>1500-1600</td><td>15</td><td>10</td><td>UV-camera</td><td>This work (12/05/2013)</td></tr><tr><td>Ili Lewotolo</td><td>123.508E/ 8.274 S</td><td>11</td><td>Gas release from the crater and fumaroles</td><td>1400-1500</td><td>75</td><td>40</td><td>UV-camera</td><td>This work (07/05/ 2013)</td></tr><tr><td>Ili Werung</td><td>123.573E/ 8.532 S</td><td>34</td><td>Minor gas releases from fumaroles on the flank and crater wall</td><td>550-580</td><td>1</td><td>0.8</td><td>UV-camera</td><td>This work (10/05/2013)</td></tr><tr><td>Batu Tara</td><td>123.585E/ 7.791 S</td><td>10</td><td>Degassing from the main crater Eruptive activity: 2012-2015</td><td>750-800</td><td>102</td><td>51</td><td>OMI</td><td>Carn et al. 2017 (2005-2015)</td></tr><tr><td>Sirung</td><td>124.130E/ 8.508 S</td><td>17</td><td>Degassing via secondary craters.</td><td>600-700</td><td>48</td><td>22</td><td>DOAS scanner</td><td>This work (13/08/2015)</td></tr></tbody></table>
<table><caption>Table 3 (continued)</caption><thead><tr><th>Volcano Name</th><th>long/lat</th><th>Rank</th><th>Typical degassing status</th><th>Average plume height (m)</th><th>Mean SO<sub>2</sub> flux (Mg/d)</th><th>error (Mg/d)</th><th>Method/ Technique</th><th>Source/ (measurement date)</th></tr></thead><tbody><tr><td>Wurlali</td><td>128.678E/ 7.125 S</td><td>25</td><td>Main crater hosts a crater lake<br>Gas release from numerous solfataro zones</td><td>700-800</td><td>8</td><td>6</td><td>UV-Camera</td><td>This work (21/10/2019)</td></tr><tr><td>Serua</td><td>130.017E/ 6.312 S</td><td></td><td>Minor degassing from the summit</td><td></td><td>?</td><td></td><td></td><td></td></tr><tr><td>Banda Api</td><td>129.881E/ 4.523 S</td><td></td><td>Minor degassing from the summit</td><td></td><td>?</td><td></td><td></td><td></td></tr><tr><td><strong>Halmahera Arc</strong></td><td></td><td></td><td></td><td>Total SO<sub>2</sub> for Banda arc</td><td></td><td></td><td></td><td></td></tr><tr><td>Gamalama</td><td>127.330E/ 0.800 N</td><td>19</td><td>Degassing via a large fracture at the summit</td><td>1650-1700</td><td>16</td><td>10</td><td>DOAS scanner</td><td>This work (27/07/ 2014)</td></tr><tr><td>Gamkonora</td><td>127.530E/ 1.380 N</td><td>27</td><td>Persistent minor degassing</td><td>1350-1400</td><td>3.4</td><td>1.0</td><td>DOAS scanner</td><td>This work (24/08/ 2018)</td></tr><tr><td>Ibu</td><td>127.630E/ 1.488 N</td><td>15</td><td>Degassing associated with dome growth and explosions</td><td>1200-1300</td><td>59</td><td>32</td><td>DOAS scanner</td><td>This work (25/09/ 2018)</td></tr><tr><td>Dukono</td><td>127.894E/ 1.693 N</td><td>1</td><td>Continuous eruptive activity, with variable intensity</td><td>1300-1400</td><td>819</td><td>394</td><td>DOAS scanner</td><td>This work (12/07/2015)</td></tr><tr><td></td><td></td><td></td><td></td><td>Total SO<sub>2</sub> for Halmahera arc</td><td></td><td></td><td></td><td></td></tr><tr><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr></tbody><tfoot><tr><td colspan="9">TOTAL FOR INDONESIA ARCHIPELAGO = 3166 ± 1321 Mg/d (1.15 ± 0.48 Tg/yr)</td><td></td></tr><tr><td colspan="9">Question mark (?) indicates persistently degassing volcanoes not yet measured.</td><td></td></tr></tfoot></table>
<figure><img src="image_5.png" /><figcaption>**Fig. 2** The main volcanic degassing points of Indonesia. The SO₂ emission rates across the four volcanic arcs of Indonesia highlight the Sunda arc as the largest SO₂ contributor and Dukono is the strongest individual source. The question marks (?) denote the unmeasured sources and the error bars correspond to standard deviation.</figcaption></figure>

<table><caption>**Table 4** Number of eruptive events per year and the corresponding SO₂ release per arc for 2010-2020 period for both from satellite and VEI results.</caption><thead><tr><th></th><th>2010</th><th>2011</th><th>2012</th><th>2013</th><th>2014</th><th>2015</th><th>2016</th><th>2017</th><th>2018</th><th>2019</th><th>2020</th><th>Total events/arc</th><th>Mean number of events/yr</th></tr></thead><tbody><tr><td>Number of eruptive events</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>Sunda</td><td>5</td><td>4</td><td>2</td><td>6</td><td>9</td><td>6</td><td>3</td><td>8</td><td>8</td><td>6</td><td>3</td><td>60</td><td>5</td></tr><tr><td>Banda</td><td>0</td><td>0</td><td>4</td><td>3</td><td>2</td><td>2</td><td>0</td><td>1</td><td>0</td><td>0</td><td>1</td><td>13</td><td>1</td></tr><tr><td>Sangihe</td><td>1</td><td>3</td><td>3</td><td>1</td><td>1</td><td>3</td><td>2</td><td>1</td><td>2</td><td>1</td><td>1</td><td>19</td><td>2</td></tr><tr><td>Halmahera</td><td>1</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>1</td><td>2</td><td>1</td><td>1</td><td>18</td><td>2</td></tr><tr><td>Total events/year</td><td>7</td><td>9</td><td>11</td><td>12</td><td>14</td><td>13</td><td>7</td><td>11</td><td>12</td><td>8</td><td>6</td><td>110</td><td>10</td></tr><tr><td>SO₂ emission per arc (Tg)</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td>Total SO₂ (Tg)</td><td>Mean Tg/yr</td></tr><tr><td>Sunda</td><td>0.172</td><td>0.050</td><td>0.009</td><td>0.004</td><td>0.574</td><td>0.006</td><td>0.047</td><td>0.000</td><td>0.021</td><td>0.028</td><td>0.002</td><td>1.069 ± 0.184</td><td>0.097 ± 0.011</td></tr><tr><td>Banda</td><td>0.641</td><td>0.044</td><td>0.022</td><td>0.082</td><td>0.824</td><td>0.235</td><td>0.056</td><td>0.151</td><td>0.821</td><td>0.722</td><td>0.118</td><td>3.715 ± 0.337</td><td>0.337 ± 0.021</td></tr><tr><td>Sangihe</td><td>0.000</td><td>0.000</td><td>0.001</td><td>0.005</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.000</td><td>0.062</td><td>0.068 ± 0.019</td><td>0.006 ± 0.003</td></tr><tr><td>Halmahera</td><td>0.000</td><td>0.010</td><td>0.003</td><td>0.000</td><td>0.000</td><td>0.001</td><td>0.000</td><td>0.000</td><td>0.014</td><td>0.005</td><td>0.000</td><td>0.397 ± 0.05</td><td>0.036 ± 0.008</td></tr><tr><td>Total Tg/yr</td><td>0.095</td><td>0.133</td><td>0.133</td><td>0.018</td><td>0.004</td><td>0.118</td><td>0.099</td><td>0.004</td><td>0.114</td><td>0.017</td><td>0.019</td><td>0.755 ± 0.055</td><td>0.069 ± 0.018</td></tr><tr><td></td><td>0.095</td><td>0.114</td><td>0.063</td><td>0.047</td><td>0.099</td><td>0.080</td><td>0.203</td><td>0.104</td><td>0.055</td><td>0.248</td><td>0.407</td><td>0.143 ± 0.011</td><td>0.013 ± 0.003</td></tr><tr><td></td><td>0.233</td><td>0.121</td><td>0.075</td><td>0.057</td><td>0.673</td><td>0.088</td><td>0.250</td><td>0.104</td><td>0.090</td><td>0.281</td><td>0.473</td><td>1.121 ± 0.008</td><td>0.102 ± 0.033</td></tr><tr><td></td><td>0.832</td><td>0.291</td><td>0.390</td><td>0.302</td><td>0.979</td><td>0.489</td><td>0.254</td><td>0.253</td><td>1.034</td><td>0.835</td><td>0.327</td><td>5.988 ± 0.310</td><td>0.119 ± 0.045</td></tr></tbody></table>

Underlined values are those obtained from VEIs.

wotobi Perempuan and Awu account for the remainder. We have also estimated explosive SO₂ emissions from Indonesia for the period 2010–2020 based on a simple scaling from reported VEI values and satellite records. The mean annual explosive-SO₂ obtained range between 0.12 and 0.54 Tg/yr, 63–81% of which is associated with the Sunda arc (0.10–0.34 Tg/yr), 5–7% (0.04–0.07 Tg/yr) the Banda arc, 3–13% (0.03–0.07 Tg/yr) the Sangihe arc and 11–18% (0.10–0.14 Tg/yr) the Halmahera arc. Combining the calculated passive and explosive SO₂ data suggests a total volcanic

SO₂ yield for the Indonesian archipelago of 1.27–1.69 Tg/yr. We consider this a representative figure, acknowledging that it is based on very limited temporal sampling of the volcanoes in question.

Our SO₂ inventory indicates a surprisingly modest SO₂ emission budget for Indonesian volcanoes, considering the 6000 km extent of the archipelago, four distinct volcanic arcs, 126 active volcanoes, and on the order of ten larger eruptions per year. For comparison, Ambrym (Vanuatu) and Kilauea (Hawaii) volcanoes
<figure><img src="image_6.png" /><figcaption>**Fig. 3** The explosive SO₂ released per volcano over the period of 2010-2020. The names of the volcanoes that erupted over the decade are grouped by arc. The SO₂ mass per volcano obtained from satellite data are displayed on the left column whilst the right column shows the SO₂ amount obtained from the VEIs. The 0 Tg correspond to undetected eruptive emission by satellite sensors. The color code differentiates the years of observation and the height corresponds to the amount of SO₂ released per year. The number of eruptions per volcano is provided above each SO₂ mass value on the left column. Note that Dukono exhibits a continuous eruptive manifestation but only the largest event with ash fall on the nearby cities are considered.</figcaption></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 4** The new SO₂ flux results compared to other estimates. A Estimates of the global volcanic SO₂ inventory that include contributions from Indonesian volcanoes, highlighted by gray square with the corresponding values. B The Indonesian SO₂ emission budget compared with other arcs (data from ref. ²⁰). The annual SO₂ emission per km of each arc²¹ are shown for comparison. C The SO₂ emission budgets from the four Indonesian arcs. D Strength of passive SO₂ emissions by altitude (in 500-m bins) from the observations reported in Table 3.</figcaption></figure>

alone have passively released more SO₂ into the atmosphere: 2.7 Tg/yr and 1.8 Tg/yr, respectively³. Several individual eruptions of the last 15 years also released comparable or higher SO₂ amounts compared with the annual Indonesian output, including Kasa-tochi (2.7 Tg) in 2008, Sarychev Peak (1.2 Tg) in 2009, Eyjafjal-lajökull (1.2 Tg) in 2010, and Nabro (4.5 Tg) in 2011¹⁶. This modest SO₂ emission budget also contrasts with the picture of renowned climate-changing Indonesian eruptions, including Agung 1963 (ref. ²²), Tambora 1815 (ref. ²³), Krakatau 1883 (ref. ²⁴), and Samalas 1257 (ref. ²⁵). More recently, the Galung-gung eruption of 1982–1983 (Java) yielded 2.5 Tg of SO₂ (ref. ²⁶).

While the overall SO₂ budget is unremarkable, the bulk of the emissions are into the free troposphere (Fig. 4D) likely to extend timescales of atmospheric processing and deposition of sulfur²⁷. Previous studies have highlighted the contribution of sulfur deposition from volcanic plume to sulfur emissions from peat fires²⁸.

Factors controlling sulfur output are numerous and include deep source characteristics and chemical processes occurring during magma storage and transfer through the crust. Hydro-thermal scavenging and scrubbing of sulfur from magmatic-hydrothermal discharges is often invoked as a process for sulfur
<figure><img src="image_8.png" /><figcaption>Fig. 5 Sulfur content in fluids and melt inclusions versus SO₂ fluxes.</figcaption></figure>

**A** Variation of the sulfur content of volcanic fluids ($XS_{\text{Total}} = XSO_2 + XH_2S$) with the measured SO₂ flux along the Indonesian arc. References for fluid compositions are from Allard et al. 1981⁵⁶ (Krakatau); Poorter et al. 1989⁵⁷ (Lowotolo); Giggenbach et al. 2001⁵⁸ (Merapi, Tangkuban Parahu, Papandayan); Clor et al. 2005⁵⁹ (Soputan); Aiuppa et al. 2015⁶⁰ (Bromo); Gunawan et al. 2017⁶¹ (Kawah Ijen); Bani et al. 2017⁶² (Sirung); Bani et al. 2017¹⁵ (Dukono); Saing et al. 2020⁶³ (Gamkonora); Kunrat et al. 2020⁶⁴ (Gamalama); Bani et al. 2020⁶⁵ (Awu). **B** The record of sulfur in melt inclusions along the Indonesian arc in respect to the main degassing sources. **C** Relationships between the maximum sulfur content analysed in melt inclusions (MI) and the measured SO₂ flux. Note, we did not find melt inclusion values for data points with zero SO₂ flux and conversely. References for MI are: Vidal et al., 2016²⁵ (Rinjani); Mandeville et al., 1996⁶⁶ (Krakatau); Bani et al., 2017¹⁵ (Dukono); Preece et al., 2014⁶⁷ (Merapi); de Hoog et al. 2001⁶⁸ (Guntur, Ili Boleng); Vigouroux et al. 2012⁶⁹ (Kawah Ijen, Galunggung, Tambora); Self and King, 1996²² (Agung); Self et al. 2004⁷⁰ (Tambora), Kunrat, 2017⁴¹ (Soputan).

depletion in volcanic fluids, significantly altering the magmatic signature²⁹. Acidic crater lakes, which are numerous in Indonesia, are perhaps the most obvious manifestation of such processes³⁰. Substantial sulfur deposits are known to be sequestered by volcanic lake systems and, conceivably, variations in climatic conditions, notably rainfall, across the archipelago could play a role in volcanic emissions to the atmosphere.

The subaerial sulfur output will depend on initial gas composition, the flow path, gas-wall rock heat transfer, and the effective water to rock ratio, all parameters that are difficult to constrain and which vary greatly between volcanoes. However, were hydrothermal scavenging and scrubbing leading mechanisms, one would expect flux strength to correlate with the concentration of sulfur in fluids. Instead, for the few volcanic centres for which, in addition to SO₂ flux, there are constraints on H₂O, CO₂, H₂S + SO₂ species in fumaroles, we find no correlation between the mass fraction of sulfur in the fluid, $XS_{\text{tot}} (=XSO_2 + XH_2S)$ and SO₂ flux (Fig. 5A), which presumably is a function of both the initial volatile content of the magma and degassing conditions. Dukono's gas, in particular, does not have higher sulfur than other Indonesian fumaroles (the case of Krakatau needs to be confirmed by more measurements). This indicates that SO₂ flux need not reflect particular enrichment/depletion in sulfur of the emitted gas, and implicates the role of degassing vigour during the time interval considered (which can scale with conduit radius and presence of an open vent to the atmosphere). From this perspective, it is worth emphasising that to establish robust links between volcanic degassing and processes at depth, it requires comprehensive measurement of gas composition.

While hydrothermal sequestration of sulfur is likely to play a significant role in modulating subaerial emissions, we consider also whether arc scale differences in SO₂ emissions across the Indonesian archipelago might reflect geodynamic or source controls, as proposed for CO₂ in arc magmas worldwide³¹. The amount of SO₂ released per km of arc per year reveals the Halmahera arc as the strongest SO₂ source, followed by the Sangihe arc. The magmatic sources of these two arcs are sustained by the double subduction of the Molucca Sea plate that deepens to the west beneath the Sangihe arc, and to the east under the Halmahera arc¹⁵. Geochemistry of lavas sampled along these arcs indicates enriched magma sources in fluid-mobile elements and notable sediment contributions³²,³³, which may play a role in subaerial sulfur budgets. The steepening of the subducted slab, the downward force from the Philippine Sea plate, and the westward motion of the continental fragments along the Sorong fault could have promoted fluid fluxes into the mantle wedge along the Halmahera arc¹⁵. To a first approximation, this peculiar geodynamic context may explain elevated SO₂ fluxes at both Halmahera and Sangihe arcs simply because of enhanced magmatic activity.

In contrast, the Banda arc stretches 2000 km but exhibits a remarkably low SO₂ emission, the weakest in our inventory. The arc is also characterized by anomalously low ³He/⁴He ratios³⁴ reflecting the arc collision with the Australian continental block and subduction of continental material that ultimately supplies less sulfur to the mantle wedge, compared to subduction of oceanic plate³⁵.

The Sunda arc is the largest SO₂ source, representing 43–48% of the total, however, its annual SO₂ emission per km of arc is modest compared with the Halmahera and Sangihe arcs, and with other arcs worldwide²⁰. The magma source beneath the Sunda arc is sustained by subduction of the Indo-Australian plate. However, while deep sea drilling has revealed a 1400 m sediment column in front of Sumatra, 300 m in front of Java, and 500 m in front of Sumbawa³⁶, less than 15% of these sediment columns is subducted³⁷. The mass transfer along the Sunda arc is dominated by an active frontal accretionary prism that strongly limits sediment subduction. Each year only 2.6 × 10⁷ m³ of sediment is subducted beneath the Sunda arc, compared with the 1.8 × 10⁸ m³/yr available, given the average subduction speed of 6.7 cm/yr. Furthermore, the sediment input from the Sunda arc is mostly trapped in the forearc basin and does not reach the trench. Given that subducted sediment can strongly contribute to the volcanic sulfur budget, it is possible that this active accretional
prism plays a key role in modulating the SO₂ emission budget of the Sunda arc by limiting the mass transfer of sediment-derived sulfur into the mantle wedge. However, as shown in the next section, such a variability in sediment contribution is not evident in variable sulfur abundance in magmas sustaining Indonesian arc volcanism.

Arc volcanoes are typically supplied by reservoirs in the shallow crust, which are in turn fed by basaltic melts rising from the mantle wedge and carrying an imprint of slab volatiles. Because they are hot, these mafic magmas have a higher sulfur carrying capacity than cooler silicic magmas, hence any volcano erupting mafic magmas should generally be associated with stronger sulfur emissions (though of course arc-scale variations in mafic melt sulfur content are possible). Detailed petrological studies of Indonesian volcanoes remain scarce, and only a few have had their sulfur content characterized via analysis of melt inclusions (MI). Figure 5B draws on these studies and shows the highest sulfur contents measured in MI along the Indonesian arcs. Drawing rough relationships from these data, basaltic MI have 2000–3500 ppm S, andesitic MI 800–500 ppm S, and rhyodacitic 200–300 ppm S, with no obvious geographical trends along arc being apparent, assuming that primary sulfur contents in the melt inclusions are little affected by post entrapment processes.

The most evident feature is that Dukono, the strongest SO₂ source we identify, has basaltic melt inclusions with rather low S content (1000 ppm) relative to other, currently, weaker SO₂ sources, such as Rinjani or Kawah Ijen, which have MI with sulfur in excess of 2000 ppm. Similarly, Soputan emits half the SO₂ of Dukono, yet has basaltic MI with much higher sulfur content (3500 ppm). This is illustrated in Fig. 5C, which shows that a broad positive correlation between SO₂ flux and the maximum sulfur content, Smax, in MI. Only Dukono departs significantly from the trend indicated by other centres. This again argues for strong decoupling between the fertility of the immediate source of magma degassing (the crustal reservoir) and its ultimate surface manifestation. Such a behavior may simply reflect conduit dynamics, such as convection, which is strongly dependent on conduit radius³⁸⁻⁴⁰, and which can sustain strong SO₂ degassing of an otherwise comparatively sulfur-poor reservoir. Alternatively, a more sulfur-rich magma may give rise to a low sulfur output simply because of a low magma influx that cannot sustain conduit convection. A critical parameter is the conduit radius, R, since magma (and gas) flux scales with R⁴ (ref. ⁴⁰) such that small variations in conduit size can result in large fluctuation in sulfur flux.

The composition and temperature of the magma supplied to the conduit will also be important, owing to their influence on viscosity and rheology. Hotter and fluid material will promote not only higher rates of magma overturn in the conduit but also more efficient degassing of slowly diffusing species, such as sulfur, in silicate melts. From this perspective, systems lying above the dashed line toward Soputan, whose high SO₂ flux may be considered as directly related to its basaltic magma source¹⁴ and high sulfur content in MI⁴¹ (Fig. 5C), may reflect more pronounced scavenging by the aquifer/hydrothermal system overlying magma reservoirs, limiting sulfur emissions to the atmosphere. These local controls will be superimposed on any deeper source signatures and may even obliterate them, as exemplified here by Dukono.

The SO₂ emission budget of the Indonesian archipelago thus reflects a complex interplay between deep (geodynamic) factors that control primary magma compositions and their availability along the arc and superficial processes such as hydrothermal scavenging and conduit dynamics. Reservoirs regularly supplied by fresh hot magma may promote sustained and vigorous degassing via conduit convection, leading to strong SO₂ outputs

even when sulfur complements in the melt are comparatively poor. In other words, the vigour of convection may largely compensate for, or even offset, any deep source deficiency in sulfur. The relatively low SO₂ output for the Indonesian archipelago documented here may appear in stark contrast with the record of explosive eruptions at several Indonesian volcanoes and their recognized global climate impacts (e.g. Tambora, Krakatau, Rinjani). These events, however, essentially reflect the long term accumulation of magmas and volatiles in closed crustal reservoirs, which cool and fractionate with little volatile loss, a process that differs from the persistent or passive degassing operating at open-conduit systems such as those we document here. An in depth knowledge of the petrology of volcanic products, and a robust characterization of emanating fluids, are both required if sound connections between the plumbing system and degassing are to be established at any active volcano.

## Methods

**Passive ultraviolet spectrometers.** We used two techniques, USB-controlled ultraviolet spectrometers and Differential Optical Absorption Spectroscopy (DOAS)⁴² and ultraviolet cameras (UV-cam)⁶. The passive ultraviolet spectrometers were either carried beneath the plume on a moving platform⁴³ or located in a fixed position and attached to scanning optics⁷. Both approaches yield the SO₂ profile across the plume. The spectrometers used were the Ocean Optics USB2000 (280–400 nm, 0.5 nm FWHM resolution), USB4000 (290–440 nm, and 0.3 nm FWHM), and USB2000 + (290–440 nm and 0.5 nm). For traverse measurements, the spectrometer was connected via an optic fibre bundle to a vertically pointed telescope of 8 mrad FOV (Field Of View). The location of each recorded spectrum was obtained using a continuously recording GPS unit. The DOAS traverse setup requires no additional power supply since the spectrometer is powered by the laptop. We operated the equipment onboard a light aircraft, from a 4WD vehicle, and on foot.

For the scanning observations, we used a rotating window that accepts light from selected directions across the plume. The light that transits through the window is redirected to an embedded telescope by a 45° optical prism, then transmitted to the spectrometer via optical fibre. The rotating window was attached to a stepper motor controlled by the laptop via a microcontroller. The system was designed to perform a 180° scanning angle with a minimum step angle of 1.8°. The scanner required an external 12 V power supply. The scanning setup could be readily operated by one person. Spectra were acquired using Jscript executed by DOASIS software⁴⁴. The script used in both traverse and scanning allowed optimization of the signal-to-noise ratio by automatically adjusting exposure time and numbers of co-added spectra⁴⁵. This was particularly useful for scanning, given the change of light intensity with scan angle. Both traverse and stationary recording were carried out at distance varying between few tens of meters from the craters to around 5 km downwind, depending on the access difficulties, the plume size and the volcanic activity.

SO₂ column amounts (ppm m) were retrieved using standard DOAS calibration and analysis procedures outlined in ref. ⁴³. Reference spectra included in the nonlinear fit were obtained by convolving high-resolution SO₂ and O₃ cross-sections with the instrument line shape. A Fraunhofer reference spectrum and ring spectrum, calculated in DOASIS, were also included in the fit. The optimum fitting windows were selected where they provided a near-random fit residual with minimum deviation. The total SO₂ column amount across the plume was then multiplied by the estimated plume speed to obtain the SO₂ flux. The plume velocities were measured mainly using videography and handheld anemometers, except in the case of airborne measurements where the plume speed was obtained by flying along and against the plume axis.

**Ultraviolet cameras.** The imaging setup consisted of two Apogee Alta U260 UV cameras. Each was coupled to a Pentax B2528-UV lens, with a focal length of 25 mm allowing a full angle FOV of around 24°. Immediately in front of each lens, a 10 nm (FWHM) bandpass filter was placed, one filter was centered at 310 nm (Asahi Spectra XBPA310) where SO₂ absorbs and the other at 330 nm (Asahi Spectra XBPA330) outside the SO₂ absorption region. Image acquisition and processing were achieved using Vulcamera⁴⁶. For each pixel the optical depth (OD) was obtained according to the following equation:

$$OD = -\ln\left[\frac{(PA - DA)}{(CA - DA}\right] + \frac{PB - DB}{(CB - DB)}\right]$$ where A and B represent the camera with the 310 nm and 330 nm filters respectively, and P, C, and D represent plume, dark and clear images. To correlate the OD values with the SO₂ slant column densities (SCDs), four calibration cells with known amounts of SO₂ (94, 189, 475, 982 ppm.m) were used. Calibration images were acquired at the beginning of measurements and repeated with long series of measurements. The UV-cam was generally positioned with a view perpendicular to the plume transport direction, at distance between <1 km and 6 km depending on plume size. Plume
speeds were derived during data processing by following plume structures between two fixed lines perpendicular to plume transport direction⁴⁶.

**Uncertainties.** Uncertainties in DOAS and UV-cam SO₂ flux measurements are discussed in many past works, including the following ref. ⁴⁷. The dominant error in the retrieved SO₂ column amount is induced by the variability of light intensity and the distance between the plume and instruments. With increasing distance, light that has not traveled through the plume can contribute significantly to the signal. This leads to light dilution of the plume signal that can easily cause more than 50% underestimate in SO₂ emission rate⁴⁸. To reduce this effect in the UV-cam measurement, we deployed the system during clear sky conditions, at distances <6 km, and performed calibration every hour during long series of measurements. UV-cam measurements were performed mainly in the late morning before the clouds started to form, generally at 9–11 am.

For the DOAS measurement, we compensate for light intensity changes using an artificial constant dark, calculated from each recorded spectrum, in the 'UV blind' region (below 290 nm). Such corrections account for dark spectrum, offset and stray light. We estimate that the error in the column amount contributes ~0.01 to the squared variation coefficient of the total flux. We also assumed that the plume and transport direction is homogeneous and in a straight line since it is difficult to rigorously assess in this work. We, therefore, performed flux calculations for direction φ, φ±3, and φ±6. The mean contribution to the square variation of the total flux is in the order of a thousandth. These errors are negligible in comparison to uncertainties in the plume speed that resulted from the complexity of wind field around volcanoes and frequent variations in both time and space. The plume transport speed relative error is conservatively assumed to be about 30–35%, which is towards the higher end of the range of past estimates⁴⁹. These errors are applied to each traverse and profile then the mean value is calculated for each series of measurement with the corresponding standard deviation. The global estimates for the arc is the sum of the mean values.

## Data availability

The data that support the plots within this article are available from the corresponding author upon reasonable request.

Received: 9 August 2021; Accepted: 30 May 2022;
Published online: 11 June 2022

## References

1. Oppenheimer, C. 3.04—Volcanic Degassing. In *Treatise on Geochemistry* (eds. Holland, H. D. & Turekian, K. K.) 123–166 (Pergamon, 2003). https://doi.org/10.1016/B0-08-043751-6/03020-6.

2. Grainger, R. G. & Highwood, E. J. Changes in stratospheric composition, chemistry, radiation and climate caused by volcanic eruptions. *Geol. Soc. Lond. Spec. Publ.* **213**, 329–347 (2003).

3. Carn, S. A., Fioletov, V. E., McLinden, C. A., Li, C. & Krotkov, N. A. A decade of global volcanic SO 2 emissions measured from space. *Sci. Rep.* **7**, 44095 (2017).

4. Clarisse, L. et al. Retrieval of sulphur dioxide from the infrared atmospheric sounding interferometer (IASI). *Atmos. Meas. Tech.* **5**, 581–594 (2012).

5. McGonigle, A. J. S. Volcano remote sensing with ground-based spectroscopy. *Philos. Trans. R. Soc. Math. Phys. Eng. Sci.* **363**, 2915–2929 (2005).

6. Mori, T. & Burton, M. The SO2 camera: a simple, fast and cheap method for ground-based imaging of SO2 in volcanic plumes. *Geophys. Res. Lett.* **33** (2006).

7. Galle, B. et al. Network for Observation of Volcanic and Atmospheric Change (NOVAC)—A global network for volcanic gas monitoring: Network layout and instrument description. *J. Geophys. Res. Atmospheres* **115** (2010).

8. Siebert, L., Cottrell, E., Venzke, E. & Andrews, B. Chapter 12—Earth's Volcanoes and Their Eruptions: An Overview. In *The Encyclopedia of Volcanoes* (Second Edition) (ed. Sigurdsson, H.) 239–255 (Academic Press, 2015). https://doi.org/10.1016/B978-0-12-385938-9.00012-2.

9. Geologi, I. B. *Data dasar gunung api Indonesia*. (Kementerian Energi dan Sumber Daya Mineral, Badan Geologi, 2011).

10. Bani, P., Surono, Hendrasto, M., Gunawan, H. & Primulyana, S. Sulfur dioxide emissions from Papandayan and Bromo, two Indonesian volcanoes. *Nat. Hazards Earth Syst. Sci.* **13**, 2399–2407 (2013).

11. Le Guern, F. Les débits de CO2 et de SO2 volcaniques dans l'atmosphère. *Bull. Volcanol.* **45**, 197–202 (1982).

12. Smekens, J.-F., Clarke, A. B., Burton, M. R., Harijoko, A. & Wibowo, H. E. SO2 emissions at Semeru volcano, Indonesia: Characterization and quantification of persistent and periodic explosive activity. *J. Volcanol. Geotherm. Res.* **300**, 121–128 (2015).

13. Nakada, S. et al. Growth process of the lava dome/flow complex at Sinabung Volcano during 2013–2016. *J. Volcanol. Geotherm. Res.* **382**, 120–136 (2019).

14. Kushendratno et al. Recent explosive eruptions and volcano hazards at Soputan volcano—a basalt stratovolcano in north Sulawesi, Indonesia. *Bull. Volcanol.* **74**, 1581–1609 (2012).

15. Bani, P. et al. Dukono, the predominant source of volcanic degassing in Indonesia, sustained by a depleted Indian-MORB. *Bull. Volcanol.* **80**, 5 (2017).

16. Carn, S. A., Clarisse, L. & Prata, A. J. Multi-decadal satellite measurements of global volcanic degassing. *J. Volcanol. Geotherm. Res.* **311**, 99–134 (2016).

17. Goode, L. R., Handley, H. K., Cronin, S. J. & Abdurrachman, M. Insights into eruption dynamics from the 2014 pyroclastic deposits of Kelut volcano, Java, Indonesia, and implications for future hazards. *J. Volcanol. Geotherm. Res.* **382**, 6–23 (2019).

18. Surono et al. The 2010 explosive eruption of Java's Merapi volcano—A '100-year' event. *J. Volcanol. Geotherm. Res.* **241–242**, 121–135 (2012).

19. Gunawan, H. et al. Overview of the eruptions of Sinabung Volcano, 2010 and 2013–present and details of the 2013 phreatomagmatic phase. *J. Volcanol. Geotherm. Res.* **382**, 103–119 (2019).

20. Fischer, T. P. et al. The emissions of CO 2 and other volatiles from the world's subaerial volcanoes. *Sci. Rep.* **9**, 18716 (2019).

21. Hilton, D. R., Fischer, T. P. & Marty, B. Noble gases and volatile recycling at subduction zones. *Rev. Mineral. Geochem.* **47**, 319–370 (2002).

22. Self, S. & King, A. J. Petrology and sulfur and chlorine emissions of the 1963 eruption of Gunung Agung, Bali, Indonesia. *Bull. Volcanol.* **58**, 263–285 (1996).

23. Oppenheimer, C. Climatic, environmental and human consequences of the largest known historic eruption: Tambora volcano (Indonesia) 1815. *Prog. Phys. Geogr. Earth Environ.* **27**, 230–259 (2003).

24. Self, S. & Rampino, M. R. The 1883 eruption of Krakatau. *Nature* **294**, 699–704 (1981).

25. Vidal, C. M. et al. The 1257 Samalas eruption (Lombok, Indonesia): the single greatest stratospheric gas release of the Common Era. *Sci. Rep.* **6**, 34868 (2016).

26. Bluth, G. J. S. et al. Evaluation of sulfur dioxide emissions from explosive volcanism: the 1982–1983 eruptions of Galunggung, Java, Indonesia. *J. Volcanol. Geotherm. Res.* **63**, 243–256 (1994).

27. Pfeffer, M. A., Langmann, B. & Graf, H.-F. Atmospheric transport and deposition of Indonesian volcanic emissions. *Atmos. Chem. Phys.* **6**, 2525–2537 (2006).

28. Langmann, B. & Graf, H. F. Indonesian smoke aerosols from peat fires and the contribution from volcanic sulfur emissions. *Geophys. Res. Lett.* **30** (2003).

29. Bani, P. et al. Remarkable geochemical changes and degassing at Vouïi crater lake, Ambae volcano, Vanuatu. *J. Volcanol. Geotherm. Res.* **188**, 347–357 (2009).

30. Delmelle, P. & Bernard, A. The remarkable chemistry of sulfur in hyper-acid crater lakes: A Scientific Tribute to Bokuichiro Takanو and Minoru Kusakabe. In *Volcanic Lakes* (eds. Rouwet, D., Christenson, B., Tassi, F. & Vandemeulebrouck, J.) 239–259 (Springer, 2015). https://doi.org/10.1007/978-3-642-36833-2_10.

31. Aiuppa, A., Fischer, T. P., Plank, T., Robidoux, P. & Di Napoli, R. Along-arc, inter-arc and arc-to-arc variations in volcanic gas CO2/ST ratios reveal dual source of carbon in arc volcanism. *Earth-Sci. Rev.* **168**, 24–47 (2017).

32. Hanyu, T. et al. Across- and along-arc geochemical variations of lava chemistry in the Sangihe arc: Various fluid and melt slab fluxes in response to slab temperature. *Geochem. Geophys. Geosystems* **13** (2012).

33. Bani, P. et al. Heterogeneity of volatile sources along the Halmahera arc, Indonesia. *J. Volcanol. Geotherm. Res.* **107342** https://doi.org/10.1016/j.jvolgeores.2021.107342 (2021).

34. Poorter, R. P. E., Varekamp, J. C., Poreda, R. J., Van Bergen, M. J. & Kreulen, R. Chemical and isotopic compositions of volcanic gases from the east Sunda and Banda arcs, Indonesia. *Geochem. Cosmochim. Acta* **55**, 3795–3807 (1991).

35. Li, J.-L. et al. Uncovering and quantifying the subduction zone sulfur cycle from the slab perspective. *Nat. Commun.* **11**, 514 (2020).

36. Plank, T. & Langmuir, C. H. The chemical composition of subducting sediment and its consequences for the crust and mantle. *Chem. Geol.* **145**, 325–394 (1998).

37. Kopp, H. et al. The Java margin revisited: Evidence for subduction erosion off Java. *Earth Planet. Sci. Lett.* **242**, 130–142 (2006).

38. Kazahaya, K., Shinohara, H. & Saito, G. Excessive degassing of Izu-Oshima volcano: magma convection in a conduit. *Bull. Volcanol.* **56**, 207–216 (1994).

39. Shinohara, H. & Tanaka, H. K. M. Conduit magma convection of a rhyolitic magma: Constraints from cosmic-ray muon radiography of Iwodake, Sumatra-Iwojima volcano, Japan. *Earth Planet. Sci. Lett.* **349–350**, 87–97 (2012).

40. Stevenson, D. S. & Blake, S. Modelling the dynamics and thermodynamics of volcanic degassing. *Bull. Volcanol.* **60**, 307–317 (1998).

41. Kunrat, S. L. Soputan Volcano, Indonesia: Petrological Systematics of Volatiles and Magmas and their Bearing on Explosive Eruptions of a Basalt Volcano. (Portland State University, 2017).
42. Platt, U. & Stutz, J. Differential Absorption Spectroscopy. In *Differential Optical Absorption Spectroscopy: Principles and Applications* (eds. Platt, U. & Stutz, J.) 135–174 (Springer, 2008). https://doi.org/10.1007/978-3-540-75776-4_6.

43. Bani, P. et al. First measurement of the volcanic gas output from Anak Krakatau, Indonesia. *J. Volcanol. Geotherm. Res.* **302**, 237–241 (2015).

44. Kraus, S. DOASIS: a framework design for DOAS. In (2006).

45. Tsanev, V. A collection of JScripts for retrieval of gas column amounts using DOAS methodology (2008).

46. Tamburello, G., Kantzas, E. P., McGonigle, A. J. S. & Aiuppa, A. Vulcamera: a program for measuring volcanic SO2 using UV cameras. *Ann. Geophys.* **54** (2011).

47. Kern, C. et al. Radiative transfer corrections for accurate spectroscopic measurements of volcanic gas emissions. *Bull. Volcanol.* **72**, 233–247 (2010).

48. Varnam, M. et al. Quantifying light dilution in ultraviolet spectroscopic measurements of volcanic SO2 using dual-band modeling. *Front. Earth Sci.* **8** (2020).

49. Stoiber, R. E., Malinconico, L. L. & Williams, S. N. Use of the correlation spectrometer at volcanoes. In *Forecasting Volcanic Events* 425–444 (Elsevier, 1983).

50. Stoiber, R. E. & Jepsen, A. Sulfur dioxide contributions to the atmosphere by volcanoes. *Science* **182**, 577–578 (1973).

51. Spiro, P. A., Jacob, D. J. & Logan, J. A. Global inventory of sulfur emissions with 1°×1° resolution. *J. Geophys. Res. Atmospheres* **97**, 6023–6036 (1992).

52. Stoiber, R. E., Williams, S. N. & Huebert, B. Annual contribution of sulfur dioxide to the atmosphere by volcanoes. *J. Volcanol. Geotherm. Res.* **33**, 1–8 (1987).

53. Andres, R. J. & Kasgnoc, A. D. A time-averaged inventory of subaerial volcanic sulfur emissions. *J. Geophys. Res. Atmospheres* **103**, 25251–25261 (1998).

54. Halmer, M. M., Schmincke, H.-U. & Graf, H.-F. The annual volcanic gas input into the atmosphere, in particular into the stratosphere: a global data set for the past 100 years. *J. Volcanol. Geotherm. Res.* **115**, 511–528 (2002).

55. Shinohara, H. Volatile flux from subduction zone volcanoes: Insights from a detailed evaluation of the fluxes from volcanoes in Japan. *J. Volcanol. Geotherm. Res.* **268**, 46–63 (2013).

56. Allard, P., Jehanno, C. & Sabroux, J.-C. Composition chimique et isotopique des produits gazeux et solides de l'activité éruptive du Krakatau (Indonésie) pendant la période 1978–1980. *Comptes Rendus Académie Sci. — Ser. IIA—* *Earth Planet. Sci.* **634**, 1095–1098 (1981).

57. Poorter, R. P. E. et al. Geochemistry of hot springs and fumarolic gases from the Banda Arc. *Neth. J. Sea Res.* **24**, 323–331 (1989).

58. Giggenbach, W. F. et al. Evaluation of results from the fourth and fifth IAVCEI field workshops on volcanic gases, Volcano island, Italy and Java, Indonesia. *J. Volcanol. Geotherm. Res.* **108**, 157–172 (2001).

59. Clor, L. E., Fischer, T. P., Hilton, D. R., Sharp, Z. D. & Hartono, U. Volatile and N isotope chemistry of the Molucca Sea collision zone: Tracing source components along the Sangihe Arc, Indonesia. *Geochem. Geophys. Geosystems* **6** (2005).

60. Aiuppa, A. et al. First determination of magma-derived gas emissions from Bromo volcano, eastern Java (Indonesia). *J. Volcanol. Geotherm. Res.* **304**, 206–213 (2015).

61. Gunawan, H. et al. New insights into Kawah Ijen's volcanic system from the wet volcano workshop experiment. *Geol. Soc. Lond. Spec. Publ.* **437**, 35–56 (2017).

62. Bani, P. et al. First study of the heat and gas budget for Sirung volcano, Indonesia. *Bull. Volcanol.* **79** (2017).

63. Saing, U. B. et al. First characterization of Gamkonora gas emission, North Maluku, East Indonesia. *Bull. Volcanol.* **82**, 37 (2020).

64. Kunrat, S. et al. First gas and thermal measurements at the frequently erupting Gamalama volcano (Indonesia) reveal a hydrothermally dominated magmatic system. *J. Volcanol. Geotherm. Res.* **407**, 107096 (2020).

65. Bani, P. et al. Elevated CO2 Emissions during Magmatic-Hydrothermal Degassing at Awu Volcano, Sangihe Arc, Indonesia. *Geosciences* **10**, 470 (2020).

66. Mandeville, C. W., Carey, S. & Sigurdsson, H. Magma mixing, fractional crystallization and volatile degassing during the 1883 eruption of Krakatau volcano, Indonesia. *J. Volcanol. Geotherm. Res.* **74**, 243–274 (1996).

67. Preece, K. et al. Pre- and syn-eruptive degassing and crystallisation processes of the 2010 and 2006 eruptions of Merapi volcano, Indonesia. *Contrib. Mineral. Petrol.* **168**, 1061 (2014).

68. de Hoog, J. C. M., Taylor, B. E. & van Bergen, M. J. Sulfur isotope systematics of basaltic lavas from Indonesia: implications for the sulfur cycle in subduction zones. *Earth Planet. Sci. Lett.* **189**, 237–252 (2001).

69. Vigouroux, N. et al. The sources of volatile and fluid-mobile elements in the Sunda arc: A melt inclusion study from Kawah Ijen and Tambora volcanoes, Indonesia. *Geochem. Geophys. Geosystems* **13** (2012).

70. Self, S., Gertisser, R., Thordarson, T., Rampino, M. R. & Wolff, J. A. Magma volume, volatile emissions, and stratospheric aerosols from the 1815 eruption of Tambora. *Geophys. Res. Lett.* **31** (2004).

## Acknowledgements

This work was carried out under the collaboration between IRD (Institut de Recherche pour le Développement) and CVGHM (Center for Volcanology and Geological Hazard Mitigation). Field expeditions were supported by JEAI-COMMISSION and ANR-DOMERAPI projects. B.S. activity was supported by LabEx VOLTAIRE (LABX-100-01). This is Laboratory of Excellence Clervolc contribution n° 542.

## Author contributions

The study was designed by P.B. and C.O. Data processing protocols and instrument calibrations were supervised by V.T. Field measurements were carried out by P.B., S.P., U.B.S., H.A. and M.M. Data analyses were performed by P.B. The paper was primarily written by P.B., C.O. and B.S. with input from all authors.

## Competing interests

The authors have no competing interests.

## Additional information

**Supplementary information** The online version contains supplementary material available at https://doi.org/10.1038/s41467-022-31043-7.

**Correspondence** and requests for materials should be addressed to Philipson Bani.

**Peer review information** *Nature Communications* thanks Simon Carn, Emma Liu and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

**Reprints and permission information** is available at http://www.nature.com/reprints

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_9.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022