Supplemental Material: Eruption dynamics leading to a volcanic thunderstorm—The January 2020 eruption of Taal volcano, Philippines

Expanded methodology and calculations, plume modeling, and details of lightning data and photographs.<br>

Supplemental Material: Eruption dynamics leading to a volcanic thunderstorm—The January 2020 eruption of Taal volcano, Philippines

Expanded methodology and calculations, plume modeling, and details of lightning data and photographs.<br>

Tectonic controls on geomorphology and spatial distribution of monogenetic volcanoes in the Central Southern Volcanic Zone of the Andes (Argentina)

Monogenetic volcanoes are among the most common volcanic landforms on Earth. The morphology and distribution of small volcanoes can provide important information about eruption dynamics and tectonics. The Southern Volcanic Zone of the Andes (CSVZ) comprises one of the most active magmatic regions on Earth. Characterized by the presence of polygenetic volcanoes and calderas in a complex tectonic setting, this region also hosts hundreds of small, back-arc monogenetic volcanoes. In this contribution, we apply a Geographic Information System (GIS) that combines imagery data and digital elevation models to establish the first comprehensive dataset of monogenetic volcanoes in the CSVZ (38° to 40° S), exploring their eruption dynamics and relationship to tectonic and structural processes. Combining spatial analysis and geomorphological observations, we identify the presence of 356 monogenetic volcanoes distributed into nine clusters, now grouped in the Zapala Volcanic Field (ZVF). The ZVF is marked by the predominance of cinder cones (80%) followed by phreatomagmatic volcanoes (20%), suggesting some influence of external water in the eruption dynamics. Generally, monogenetic vents present a clear association with local and regional lineaments, suggesting a strong structural control on the occurrence of the monogenetic deposits. The higher vent densities are observed in the southern Loncopué Though, an important extensional feature related to tearing of the subducted Nazca plate underneath the South American Plate. Morphometric parameters of cinder cones indicate variable stress orientations in the CSVZ that possibly result from the oblique tectonics in the region. From north to south, the maximum principal stress rotates from NE-SW to E-W and becomes progressively less constrained as it distances from the current magmatic arc. Based on the relative ages, we map the evolution of monogenetic volcanism through time. Our results suggest a waning in the monogenetic activity in ZVF over time. When compared to monogenetic fields in the Central Andes, the ZVF is marked by higher vent densities and number of phreatomagmatic landforms, with the absence of lava domes. This ultimately reflects the contrasting crustal structure and climate conditions of these two regions.

External surface water influence on explosive eruption dynamics, with implications for stratospheric sulfur delivery and volcano-climate feedback

Explosive volcanic eruptions can inject sulfur dioxide (SO2) into the stratosphere to form aerosol particles that modify Earth’s radiation balance and drive surface cooling. Eruptions involving interactions with shallow layers (< 500 m) of surface water and ice modify the eruption dynamics that govern the delivery of SO2 to the stratosphere. External surface water potentially controls the evolution of explosive eruptions in two ways that are poorly understood: (1) by modulating the hydrostatic pressure within the conduit and at the vent, and (2) through the ingestion and mixing of external water, which governs fine ash production as well as eruption column buoyancy flux. To make progress, we couple one-dimensional models of magma flow in the conduit and atmospheric column rise through a novel ”magma-water interaction” model that simulates the occurrence, extent and consequences of water entrainment depending on the depth of a surface water layer. We explore the effects of hydrostatic pressure on magma ascent in the conduit and gas decompression at the vent, and the conditions for which water entrainment drives fine ash production by quench fragmentation, eruption column collapse, or outright failure of the jet to breach the water surface. We show that the efficiency of water entrainment into the jet is the predominant control on jet behavior. For an increase in water depth of 50 to 100 m, the critical magma mass eruption rate required for eruption columns to reach the tropopause increases by an order of magnitude. Finally, we estimate that enhanced emission of fine ash leads to up to a 2-fold increase in the mass flux of particles < 125 microns to spreading umbrella clouds, together with up to a 10-fold increase in water mass flux, conditions that can enhance the removal of SO2 via chemical scavenging and ash sedimentation. Overall, compared to purely magmatic eruptions, we suggest that hydrovolcanic eruptions will be characterized by a reduced delivery of SO2 to the stratosphere. Our results thus suggest the possibility of an unrecognized volcano-climate feedback mechanism arising from modification of volcanic climate forcing by direct interaction of erupting magma with varying distributions of water and ice at the Earth’s surface.

Rapid eruptive transitions from low to high intensity explosions and effusive activity: insights from textural analysis of a small-volume trachytic eruption, Ascension Island, South Atlantic

Proximal deposits of small-volume trachytic eruptions are an under-studied record of eruption dynamics despite being common across a range of settings. The 59 ± 4 ka Echo Canyon deposits, Ascension Island, resulted from a small-volume explosive-effusive trachytic eruption. Variations in juvenile clast texture reveal changes in ascent dynamics and transitions in eruption style. Five dominant textural types are identified within the pumice lapilli population. Early Strombolian-Vulcanian eruption phases are typified by macro- and micro-vesicular equant clast types. Sheared clasts are most abundant at the eruption peak, transitioning to dense clasts in later phases due to shear-induced coalescence, outgassing and vesicle collapse. Melt densification and outgassing via tuffisite veins increased plume density, contributing to partial column collapse and the explosive-effusive transition. Bulk vesicularity distributions indicate a shift in dominant fragmentation mechanism during the eruption, from early-stage bubble interference and rupture to late-stage transient fragmentation, with a transient peak of Plinian activity. Dome and lava groundmass crystallinities of up to 70% indicate near-complete degassing during effusive phases, followed by shallow over pressurisation and a final less explosive phase. We provide textural evidence for high-intensity explosive phases and rapid transitions in eruptive style during small-volume trachytic eruptions and consider the impact of trachytic melt compositions on underlying dynamics of these short-lived, explosive events. This analysis demonstrates the value of detailed stratigraphy in understanding critical changes in eruption dynamics and the timescales over which they may occur which is of particular value in anticipating future eruptions of this type.

Online treatment of eruption dynamics improves the volcanic ash and SO₂ dispersion forecast: case of the Raikoke 2019 eruption

In June 2019, the Raikoke volcano, Kuril Islands, emitted 0.4–1.8 × 109 kg of very fine ash and 1–2 × 109 kg of SO2 up to 14 km into the atmosphere. The eruption was characterized by several phases or puffs of different duration and eruption heights. Resolving such complex eruption dynamics is required for precise volcanic plume dispersion forecasts. To address this issue, we coupled the atmospheric model system ICON-ART (ICOsahedral Nonhydrostatic – Aerosols and Reactive Trace gases) with the 1-D plume model FPlume to calculate the eruption source parameters (ESPs) online. The main inputs are the plume heights for the different eruption phases that are geometrically derived from satellite data. An empirical relationship is used to derive the amount of very fine ash (particles < 32 µm), which is relevant for long range transport in the atmosphere. On the first day after the onset of the eruption, the modeled ash loading agrees very well with the ash loading estimated from AHI (Advanced Himawari Imager) observations due to the resolution of the eruption phases and the online treatment of the ESPs. In later hours, aerosol dynamical processes (nucleation, condensation, coagulation) explain the loss of ash in the atmosphere in agreement with the observations. However, a direct comparison is partly hampered by water and ice clouds overlapping the ash cloud in the observations. We compared 6-hourly means of model and AHI data with respect to the structure, amplitude, and location (SAL-method) to further validate the simulated dispersion of SO2 and ash. In the beginning, the structure and amplitude values differed largely because the dense ash cloud leads to an underestimation of the SO2 amount in the satellite data. On the second and third day, the SAL values are close to zero for all parameters indicating a very good agreement of model and observations. Furthermore, we found a separation of the ash and SO2 plume after one day due to particle sedimentation, chemistry, and aerosol-radiation interaction. The results confirm that coupling the atmospheric model system and plume model enables detailed treatment of the plume dynamics (phases and ESPs) and leads to significant improvement of the ash and SO2 dispersion forecast. This approach can benefit the operational forecast of ash and SO2 especially in case of complex and non-continuous volcanic eruptions like the Raikoke 2019.

Using apatite records of volatile budget and magma ascent rates to investigate eruption dynamics

&lt;p&gt;One of the major challenges faced by volcanologists to investigate controls on eruption dynamics is to quantify both pre-eruptive volatile budgets and timescales of magma ascent. Indeed, petrological investigations of the two parameters usually rely on different methods/analytical techniques that are not always applicable/accessible. Recent studies have shown that the abundance and zoning pattern of F, Cl, and OH in apatite can be used to determine both pre-eruptive volatile budget and magma degassing rates that can, under some conditions, be related to magma ascent rates ([1],[2]).&lt;/p&gt;&lt;p&gt;Here we apply the two methods to apatite in the Rabaul 2006 eruption deposits (Papua-New-Guinea). This was a VEI-4 eruption and occurred in three main phases: (1) a sub-plinian onset followed 12h after its start by (2) a mixed strombolian-effusive phase, which subsequently evolved into (3) discrete vulcanian explosions. We sampled deposits of the three phases: (1) pumices, (2) fragments of lava flow, and (3) fragments of cow-pad bombs.&lt;/p&gt;&lt;p&gt;We calculated pre-eruptive water contents using apatite included in clinopyroxene as they keep a better record of reservoir conditions from the time of entrapment. We found that the magma that fed the sub-plinian phase contained the highest water content of about 2 wt.%, while magmas that fed the lava flow and the vulcanian phase were drier, with 0.2 to 0.5 wt.% less H&lt;sub&gt;2&lt;/sub&gt;O. X-ray maps acquired with an EPMA show that only apatite crystals in the groundmass of the vulcanian and effusive deposits are zoned in F and Cl at the crystal rims, whereas those from the sub-plinian deposits and included in clinopyroxenes are not zoned. This indicates that the zoning is related to syn- or immediately pre-eruptive changes of Cl-F-H&lt;sub&gt;2&lt;/sub&gt;O during magma ascent towards the surface and can thus be modelled as diffusive reequilibration of the crystal and the melt. We obtained maximum diffusion timescales of &lt;8 hours for the unzoned apatite in sub-plinian deposits, timescales of 20&amp;#8211;22 hours for apatite in vulcanian deposits, and 600&amp;#8211;1500 hours for those in the lava flow. Thus, the time scales increase with decreasing explosivity of the eruptions, as it could be expected if magma ascent rate played the key role of eruption dynamics. However, the degassing timescales of the effusive phase are significantly longer than the eruption duration itself, which can be explained if the magma started rising in the system 1&amp;#8211;3 months prior to the onset of the eruption. The volatile-rich, fast-rising magma that fed the initial sub-plinian phase propagated through, disturbed and remobilized the shallower, more degassed batch of magma, which was erupted during the following effusive phase. Deeper, volatile-poor magma that kept moving up the open conduit, was responsible for the late vulcanian explosions.&lt;/p&gt;&lt;p&gt;Our results show that apatite is a powerful tool for probing slight changes in magma volatile chemistry and ascent rates that can vary between different phases of the same eruption and produce different eruption styles.&lt;/p&gt;&lt;p&gt;&amp;#160;&lt;/p&gt;&lt;p&gt;[1] Li and Costa, 2020, GCA [2] Li et al. 2020, EPSL&lt;/p&gt;