Stromboli Research Articles

Pre-eruptive Conditions of the 3 March 2015 Lava Fountain of Villarrica Volcano (Southern Andes)

Villarrica or Rukapillan (35.9°S; 2,847 m a.s.l.) is one of the most active volcanoes in South America and is the highest-risk volcano in Chile. It has an open conduit with a persistent lava lake. On the 3 March 2015, Strombolian activity rapidly progressed into a 1.5-km-high lava fountain, erupting at least ∼ 2.4 × 106 m3 of tephra. Soon after, the activity returned to mild Strombolian “background” explosions, which lasted until early 2017. Understanding the pre-eruptive conditions of such paroxysmal events is fundamental for volcanic hazard assessment. We present major and trace element geochemistry for glass and crystalline phases of basaltic andesite paroxysm pyroclasts (52–56 wt.% SiO2), and for the subsequent Strombolian “background” activity through February 2017 (54–56 wt.% SiO2). The lava fountain source magma was initially stored in a deeper and hotter region (9.4–16.3 km; ca. 1140 °C) and was then resident in a shallow (≤ 0.8 km) storage zone pre-eruption. During storage, crystallising phases comprised plagioclase (An66–86), olivine (Fo75–78) and augite (En46–47). Equilibrium crystallisation occurred during upper-crustal magmatic ascent. During storage in the shallower region, magma reached H2O saturation, promoting volatile exsolution and over-pressurization, which triggered the eruption. In contrast, subsequent “background” explosions involving basaltic-andesite were sourced from a depth of ≤ 5.3 km (ca. 1110 °C). Pre-eruptive conditions for the 2015 lava fountain contrast with historical twentieth-century eruptions at Villarrica, which were likely driven by magma that underwent a longer period of mixing to feed both effusive and explosive activity. The rapid transition to lava-fountaining activity in 2015 represents a challenging condition in terms of volcano monitoring and eruption forecasting. However, our petrological study of the pyroclastic materials that erupted in 2015 offers significant insights into eruptive processes involving this type of eruption. This aids in deciphering the mechanisms behind sudden eruptions at open conduit systems.

Silent Very Long Period Seismic Events (VLPs) at Stromboli Volcano, Italy

AbstractAt Stromboli Volcano, Italy, very long period (VLP) seismic signals and Strombolian eruptions have been attributed to the unsteady flow of gas slugs through the shallow plumbing system followed by explosive slug bursting at a free surface. In data from a 2018 seismo‐acoustic deployment, ∼92% of events in two main VLP multiplets do not coincide in time with impulsive infrasonic signals (the expected signal of explosive slug bursting); we term these “silent VLPs.” The lack of infrasonically detected explosions relative to repeating VLPs does not support the commonly invoked “gas slug” model. We propose that VLPs may be generated when gas bubbles move into a weak semi‐solid plug in the uppermost portion of the conduit. The plug then acts as a mechanical filter in which pathways vary and guide or trap ascending gas slugs, allowing for passive (silent) gas release and explosive escape mechanisms decoupled in time from VLPs.

Deformation at the open-vent Erebus volcano, Antarctica, from more than 20 years of GNSS observations

Erebus volcano on Ross Island in Antarctica is an iconic open-vent volcano that has hosted a convecting anorthoclase-rich phonolite lava lake for at least 50 years. Recent magnetotelluric observations have imaged a conduit system originating at least 60km below Erebus. This terminates in a seismically-imaged shallow magma reservoir about 500m NW of the crater with an upper surface at about 500m depth. A narrow and inclined terminal conduit system connects the shallow reservoir to the lava lake. Bomb-ejecting Strombolian eruptions from the lava lake and sporadically active adjacent vents are common. Larger eruptive activity with locally substantially elevated hazard has also been observed, including exceptionally energetic Strombolian activity in 1984 and an Inner Crater phreatic explosion in 1993. Despite sustained degassing and frequent eruptions, geochemical data show the composition of the lava has remained stable for the last 17ka, consistent with the long-lived transcrustal magmatic system underlying the lake.Global Navigation Satellite System (GNSS) data collection on Ross Island began in the late 1990s with campaign observations. In the early 2000s, additional benchmarks were added closer to the summit of Erebus and continuous GNSS (cGNSS) sites were co-located with a seismic network. Today, seven cGNSS sites operate on the summit plateau and flanks of the volcano, and a network of eight campaign benchmarks has been surveyed episodically.We present the first comprehensive geodetic data analysis and modeling results integrating these more than two decades of data collected at Erebus. We resolve long-term subsidence of Ross Island, which a simple viscoelastic model links to the long-term growth of Erebus over the last 20ka. The data also show multi-year cycles of inflation and deflation in the summit area, consistent with activity in the shallow summit magmatic system. These small amplitude (several mm/yr) transient dynamics suggest multi-year pulses of pressurization and depressurization, or geometric changes within the shallow magmatic system that we can reproduce with analytical source models. The most recent inflation pulse lasted from November 2020 until March 2022 when several stations moved radially away from the shallow magmatic system and upwards at 10–15mm/yr. Based on prior cycles, this deformation may result in increased eruptive activity suggesting that continued and enhanced surveillance of Erebus is warranted. These observations contrast with the long-term general stability of the lava lake, but reflect Inner Crater dynamics, which can include changes in lava lake elevation and associated topographic changes of over 20m on multi-year time scales. Our results emphasize the value of long-term and campaign-aided high-accuracy GNSS observations at open-vent volcanoes. This is especially true for volcanoes like Erebus which are remote and may only be accessible for a few months a year, and that deform at amplitudes and periods that may be difficult to resolve with satellite-based radar.

Automatic Detection of the Explosive Activity of the Mt. Etna Volcano through Doppler Radar Monitoring

Improving the capability to detect volcanic explosive activity could be strategic for the task of a volcano observatory to inform civil protection authorities and air traffic controllers. The detection of explosive volcanic activity can be done in real time and also under bad visibility conditions by using the radar remote sensing technique. Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo (INGV-OE) installed an S-band pulse Doppler radar in a shelter located at about 3 km south of the active volcanic vents in order to enhance the Etna volcano’s surveillance. Here, we describe the realisation of a system that exploits such device, aimed at continuously monitoring the explosive eruptive activity at the Mt. Etna summit craters through an automatic processing flow of the radar data. We analysed the signals recorded during 23 eruptive episodes that occurred at the Etna South-East Crater during the second half of 2021; these episodes were characterised by an opening Strombolian activity and the subsequent evolution into a lava fountain. To identify the onset of both volcanic phenomena, empirical thresholds of radar time series were extracted with the help of thermal and visible images acquired by the INGV-OE cameras’ network. The resulting monitoring tool automatically operates 24/7 for volcanic surveillance, providing real-time data to the INGV-OE control room.

Long range infrasound monitoring of Yasur volcano

The atmospheric injection of gas and material produced by an explosive volcanic eruption determines a rapid compression of the atmosphere, which subsequently propagates as longitudinal elastic waves (sound). The size of the source, generally greater than tens of meters, and its duration, longer than a few seconds, result into an emitted signal that is particularly rich in low frequencies (f < 20 Hz), thus determining an efficient infrasound radiation. Thanks to the low spectral content and the reduced attenuation in the atmosphere, infrasound is capable of propagating for very large distances.In this study we show how the infrasonic monitoring at regional distances (> 100 km) is efficient in recording and characterizing volcanic events. For the purpose of our study, infrasound signal radiated from Yasur volcano (Tanna Island, Vanuatu) was studied for a period of twelve years (January 2008 – December 2019). Signals from Yasur were registered at a source-to-receiver distance of 400 km by the IS22 infrasound array, located in New Caledonia, part of the Comprehensive Nuclear-Test-Ban Treaty (CTBT) International Monitoring System (IMS). The predominantly explosive Strombolian activity of this volcano makes it a perfect subject to be studied by infrasound technology. Detections of volcanic infrasound are modulated according to the seasonal variation of stratospheric winds and corrected for attenuation accounting for real atmospheric specification between the source and the receiver to retrieve the pressure at the source. Next, they are used to evaluate long-term (yearly) and short term (hourly) variations of activity over the period of analysis. Results are compared with thermal anomalies recorded by the MODerate resolution Imaging Spectroradiometer (MODIS) installed on NASA's Terra and Aqua satellites.We show that even at regional distances (400 km) it is possible to follow the long term (yearly) fluctuations of ordinary explosive activity during periods of optimal propagation of infrasonic waves in the atmosphere. In addition, we show that the time resolution retrieved from the signal analysis allows to follow variations of volcanic activity at hourly time scale, thus representing a valuable source of information, in particular in areas where local geophysical observations are missing.

The recent volcanism of Flores Island (Azores): Stratigraphy and eruptive history of Funda Volcanic System

Monogenetic volcanoes occur in many different geotectonic settings and are usually small and short-lived. They can experience a variety of eruptive styles, even during the same eruption. In monogenetic volcanic fields, volcanism usually migrates to different locations over time, making volcanic hazard assessment very challenging. The eruptive history of a volcanic region, including the size, style, and location of previous eruptions, provides valuable information to help predict the behaviour of future volcanic events and their associated hazards. Here, we reconstruct for the first time the eruptive history of the Funda Volcanic System (FVS), one of the most recent (∼3 ka) monogenetic eruptive centres of Flores Island (Azores), based on a detailed tephrostratigraphic work coupled with geochemical analysis of glass shards and radiocarbon dating. We identified at least three volcanic events at FVS spaced by time intervals of ∼100 yr. The first event (3430 cal yr BP) was a small Strombolian eruption, the second event (3330 cal yr BP) started as a violent Strombolian eruption and may have ended as phreatomagmatic, and the third event (3250 cal yr BP) was exclusively phreatomagmatic. Our results demonstrate that volcanism at the FVS was more prolonged and recurrent than previously reported. Moreover, we show that the FVS experienced different eruptive styles in a short timeframe, ranging from small basaltic eruptions to violent explosive phreatomagmatic events. Such diversity of eruptive styles results in different volcanic products, which have different hazard implications. Our new results contribute to the knowledge of the recent volcanic activity of Flores Island, and we anticipate them to be of paramount importance for future volcanic hazard assessment studies.

Spatio-temporal changes in degassing behavior at Stromboli volcano derived from two co-exposed SO2 camera stations

Improving volcanic gas monitoring techniques is central to better understanding open-vent, persistently degassing volcanoes. SO2 cameras are increasingly used in volcanic gas studies, but observations are commonly limited to one single camera alone viewing the volcanic plume from a specific viewing direction. Here, we report on high frequency (0.5 Hz) systematic measurements of the SO2 flux at Stromboli, covering a 1-year long observation period (June 2017-June 2018), obtained from two permanent SO2 cameras using the same automated algorithm, but imaging the plume from two different viewing directions. Our aim is to experimentally validate the robustness of automatic SO2 camera for volcano monitoring and to demonstrate the advantage of using two co-exposed SO2 camera stations to better capturing degassing dynamics at open-vent volcanoes. The SO2 flux time-series derived from the two SO2 camera stations exhibit good match, demonstrating the robustness of the automatic SO2 camera method. Our high-temporal resolution SO2 records resolve individual Strombolian explosions as transient, repetitive gas bursts produced by the sudden release of over pressurized gas pockets and scoriae. Calculations show that explosive degassing activity accounts for ∼10% of the total SO2 emission budget (dominated by passive degassing) during mild regular open-vent activity. We show that the temporal variations of the explosive SO2 flux go in tandem with changes in total SO2 flux and VLP seismicity, implicating some commonality in the source processes controlling passive degassing and explosive activity. We exploited the spatial resolution of SO2 camera to discriminate degassing at two distinct regions of the crater area, and to minimize biases due by the station position respect to the target plume. We find that the SO2 fluxes from southwest-central (SWCC) and northeast (NEC) crater areas oscillate coherently but those from the NEC are more sensitive to the changes in the volcanic intensity. We interpret this as due to preferential gas/magma channeling into the structurally weaker north-eastern portion of the crater terrace in response to increasing supply rate of buoyant, bubble-rich magma in the shallow plumbing system.

Energy threshold changes in volcanic activity at Mt. Etna (Italy) inferred from volcanic tremor

From the 2010s on, pattern classification has proven an effective method for flagging alerts of volcano unrest before eruptive activity at Mt. Etna, Italy. The analysis has been applied online to volcanic tremor data, and has supported the surveillance activity of the volcano that provides timely information to Civil Protection and other authorities. However, after declaring an alert, no one knows how long the volcano unrest will last and if a climactic eruptive activity will actually begin. These are critical aspects when considering the effects of a prolonged state of alert. An example of longstanding unrest is related to the Christmas Eve eruption in 2018, which was heralded by several months of almost continuous Strombolian activity. Here, we discuss the usage of thresholds to detect conditions leading to paroxysmal activity, and the challenges associated with defining such thresholds, leveraging a dataset of 52 episodes of lava fountains occurring in 2021. We were able to identify conservative settings regarding the thresholds, allowing for an early warning of impending paroxysm in almost all cases (circa 85% for the first 4 months in 2021, and over 90% for the whole year). The chosen thresholds also proved useful to predict that a paroxysmal activity was about to end. Such information provides reliable numbers for volcanologists for their assessments, based on visual information, which may not be available in bad weather or cloudy conditions.

Inferences on the magmatic plumbing system at Stromboli volcano (Italy) from trace element geochemistry of matrix glasses and minerals in different types of explosive eruptions

Stromboli (Italy) is a basaltic volcano characterized by persistent, mild strombolian activity, occasionally interrupted by lava effusion and more violent explosive events, named major explosions and paroxysms depending on their intensity and magnitude. The normal activity is fed by a shallow and degassed highly porphyritic (HP) shoshonitic basalt carrying about 50 vol.% crystals settled in a shoshonitic glassy matrix (K2O > 3.8 wt.%). The more energetic explosions erupt a deep, volatile-rich, low-porphyritic (LP) magma with < 10 vol.% crystals in a shoshonitic basaltic glassy matrix (K2O < 2.4 wt.%). Products with intermediate glass composition are also found in the more violent explosive events. In this study, we present a new data set of major and trace element contents in matrix glasses and minerals performed in products from different types of explosive activity that occurred at Stromboli between 1998 and 2020. This large data set is used to put constraints on the evolution and architecture of the intermediate plumbing system, where the transformation from LP to HP occurs. Results indicate that, compared to paroxysms, the glassy matrices of the LP pumices from major explosions are richer in incompatible trace elements (and K2O wt.%) due to < 15 wt.% fractionation of clinopyroxene and olivine. This points to a chemical zoning of the deep reservoir and suggests that major explosions are fed by magmas residing in its upper part. Among the major explosions, the homogeneous intermediate glasses in the products from the 19 July 2020 event originate from the interplay of mixing and crystal fractionation processes. The crystallization of euhedral microphenocrysts of An-rich plagioclase suggests that batches of magma can pond and crystallize for few days (< 11) at the base of the intermediate zone of the plumbing system, at pressure coinciding with the entering of plagioclase into the system (< 100 MPa). As a relevant point for understanding the pre- and syn-eruptive magma dynamics, data indicate a positive correlation between the magnitude of the explosions and the depth of the supply magma.

Crater-rim collapse at Stromboli volcano: understanding the mechanisms leading from the failure of hot rocks to the development of glowing avalanches

The failures of volcanic crater-rims frequently lead to the development of avalanches of glowing rocks, a hybrid phenomenon between rock avalanches and pyroclastic density currents, reaching considerable distances from the eruptive centres and therefore be a serious threat for inhabited areas. The collapse conditions have been here taken in to account by means of slope stability analysis using a Limit Equilibrium Method analysis on the crater-rim of Stromboli volcano (Italy). A Stromboli, crater-rim collapses occurred frequently (at least seven events in the last two decades) and were always associated with high-level of magma within the conduits, testified by the increased eruptive activity and ground-deformation. The more frequent/intense eruptive activity produced a greater accumulation of volcaniclastic material, whereas the high level of magma increased magmastatic thrust on the deposits. Volcaniclastic material was modelled by combining the nonlinear failure envelopes as the Generalized Hoek and Brown criterion, with the addition of the failure's envelope of the rockfill-like material described by the Barton-Kjaernsli criterion, taking into consideration the presence of discontinuities within proximal, partially welded, volcaniclastic masses. In addition to the lithological and morphological characteristics of the crater terrace rim and the magmatic thrust, the effects of the explosions in terms of seismic ground acceleration and disturbance factor (D) of the volcaniclastic material were also considered here. While the ground acceleration compatible with the explosive activity of Stromboli has little influence on the stability of the crater terrace rims, the increase in D increases the proneness for failure.

Integration of Remote Sensing and Offshore Geophysical Data for Monitoring the Short-Term Morphological Evolution of an Active Volcanic Flank: A Case Study from Stromboli Island

The Sciara del Fuoco (SdF) collapse scar at Stromboli is an active volcanic area affected by rapid morphological changes due to explosive/effusive eruptions and mass-wasting processes. The aim of this paper is to demonstrate the importance of an integrated analysis of multi-temporal remote sensing (photogrammetry, COSMO-SkyMed Synthetic Aperture Radar amplitude image) and marine geophysical data (multibeam and side scan sonar data) to characterize the main morphological, textural, and volumetric changes that occurred along the SdF slope in the 2020–2021 period. The analysis showed the marked erosive potential of the 19 May 2021 pyroclastic density current generated by a crater rim collapse, which mobilized a minimum volume of 44,000 m3 in the upper Sciara del Fuoco slope and eroded 350,000–400,000 m3 of material just considering the shallow-water setting. The analysis allowed us also to constrain the main factors controlling the emplacement of different lava flows and overflows during the monitored period. Despite the morphological continuity between the subaerial and submarine slope, textural variations in the SdF primarily depend on different processes and characteristics of the subaerial slope, the coastal area, the nearshore, and “deeper” marine areas.

Forecasting and communicating the dispersion and fallout of ash during volcanic eruptions: lessons from the September 20, 2020 eruptive pulse at Sangay volcano, Ecuador

Volcanic ash is a hazard to human health and activities. Accurate and timely forecasts coupled with effective communication about the dispersion and fallout of volcanic ash during explosive events are essential to reduce impacts on local communities and limit economic losses. In this study, we present the first detailed description of an eruptive pulse at Sangay volcano and its eruption source parameters. The eruptive pulse on September 20, 2020, started at about 9:20 UTC and lasted between 90 and 100 min, producing an eruptive column that rapidly separated into (1) a higher (15.2 km above sea level, asl), gas-rich cloud moving east-southeast and (2) a lower (12.2 km asl), ash-rich cloud moving west and causing ash fallout up to 280 km from the volcano. Field data collected immediately after the event allow estimating the volume of bulk tephra to be between 1.5 and 5.0 × 106 m3, corresponding to a volcanic explosivity index of 2. The eruptive pulse, identified as violent Strombolian, emitted andesitic ash that was more mafic than products ejected by Sangay volcano in recent decades. Component analysis and glass chemistry of juvenile particles support the hypothesis that this event excavated deeper into the upper conduit compared to typical Strombolian activity at Sangay volcano, while grain-size analysis allows reconstruction of the total grain-size distribution of the fallout deposit. The discrepancies between the ash fallout simulations performed with the Ash3D online tool and the actual deposit are mainly the result of inaccurate pre- and syn-eruptive configurations and highlight the importance of additional processes such as aggregation. Communication products issued during the event included (1) several standard short reports, (2) volcano observatory notices for aviation, (3) social media posts, and (4) a special report providing the results of the ash fallout simulation. Although communication was effective with the authorities and the connected population, an effort must be made to reach the most vulnerable isolated communities for future events. This study shows that forecasting ash dispersion and fallout during volcanic eruptions can guide early warnings and trigger humanitarian actions, and should become a standard in volcano observatories worldwide.

Multi-parametric study of an eruptive phase comprising unrest, major explosions, crater failure, pyroclastic density currents and lava flows: Stromboli volcano, 1 December 2020–30 June 2021

Open conduit volcanoes like Stromboli can display elusive changes in activity before major eruptive events. Starting on December 2020, Stromboli volcano displayed an increasing eruptive activity, that on 19 May 2021 led to a crater-rim collapse, with pyroclastic density currents (PDCs) that spread along the barren NW flank, entered the sea and ran across it for more than 1 km. This episode was followed by lava flow output from the crater rim lasting a few hours, followed by another phase of lava flow in June 2021. These episodes are potentially very dangerous on island volcanoes since a landslide of hot material that turns into a pyroclastic density current and spreads on the sea surface can threaten mariners and coastal communities, as happened at Stromboli on 3 July and 28 August 2019. In addition, on entering the sea, if their volume is large enough, landslides may trigger tsunamis, as occurred at Stromboli on 30 December 2002. In this paper, we present an integration of multidisciplinary monitoring data, including thermal and visible camera images, ground deformation data gathered from GNSS, tilt, strainmeter and GBInSAR, seismicity, SO2 plume and CO2 ground fluxes and thermal data from the ground and satellite imagery, together with petrological analyses of the erupted products compared with samples from previous similar events. We aim at characterizing the preparatory phase of the volcano that began on December 2020 and led to the May–June 2021 eruptive activity, distinguishing this small intrusion of magma from the much greater 2019 eruptive phase, which was fed by gas-rich magma responsible for the paroxysmal explosive and effusive phases of July–August 2019. These complex eruption scenarios have important implications for hazard assessment and the lessons learned at Stromboli volcano may prove useful for other open conduit active basaltic volcanoes.

A GIS-based multi-hazard assessment at the San Pedro volcano, Central Andes, northern Chile

Recent advances in the modeling of volcanic phenomena have allowed scientists to better understand the stochastic behavior of volcanic systems. Eruptions can produce various types of volcanic phenomena of different sizes. The size of a given volcanic phenomenon dominates its spatial distribution and is commonly represented by volume/mass parameters in the models that reproduce their behavior. Multi-hazard assessments depend on first-order parameters to forecast hazards at a given geographic location. However, few multi-hazard assessments consider the size of the eruption (e.g., tephra fallout) to co-parameterize the size of the accompanying phenomena (e.g., mass flows) in a given eruptive scenario. Furthermore, few studies simulate multi-phenomenon eruptive scenarios with semi-continuous variations in their size, something that allows a better quantification of the aleatoric variability of the system. Here, we present a multi-hazard assessment of the San Pedro volcano, a high-threat volcano from northern Chile, that produced two large-size Plinian eruptions (VEI 5 and 6) in the last 16 ka, and ten Strombolian eruptions (VEI 2) between 1870 and 2021 CE, with the latest occurring on 2 December 1960 CE. We use intra-scenarios (i.e., subdivisions of eruptive scenarios) to explore the size variability of explosive volcanic phenomena. The size of intra-scenarios is extrapolated from the largest-size deposits of each type of phenomenon from the geologic record of the San Pedro volcano. We simulate explosive intra-scenarios for tephra fallout, concentrated PDCs, and lahars, and effusive scenarios for blocky lava flows. On the local scale, mass flows are likely (66–100%) to affect transport and energy infrastructure within a 14 km radius of the volcano. On the regional scale, large-size eruptions (VEI 5) in the rainy season are about as likely as not (33–66%) to accumulate 1 cm of tephra on energy, transport, and mining infrastructure over a 50 km radius, and these same eruptions are unlikely (10–33%) to accumulate 1 cm of tephra on the city of Calama. This work shows how multi-phenomenon intra-scenarios can be applied to better quantify the aleatoric variability of the type and size of volcanic phenomena in hazard assessments.

Explosive Paroxysmal Events at Etna Volcano of Different Magnitude and Intensity Explored through a Multidisciplinary Monitoring System

Between 13 December 2020 and 21 February 2022, Etna volcano produced a sequence of 66 paroxysmal explosive eruptions, with Strombolian activity at the summit craters climaxing in lava fountains and eruption columns extending several kilometers above the craters, accompanied by minor and short-lasting lava flows from the crater rim. We selected three of these episodes that occurred within a short space of time, between 13 December 2020 and 12 March 2021, of different magnitude (i.e., erupted volume) and intensity (i.e., mass eruption rate or instantaneous eruption rate), and analyzed them from a multidisciplinary perspective. The aim was to gain insights into those parameters that mostly reveal the eruptive process for hazard assessment purposes. The multidisciplinary data consist of calibrated visible images, thermal images, seismic and infrasound data, ground deformation detected from the strainmeters, as well as satellite SEVIRI images. From these data, we obtained the timing of each paroxysmal event, the erupted volume in terms of tephra and lava flows, and the corresponding deflation of the source region, together with the development of the lava fountains and eruption columns with time. The results enabled determining that the smallest episode was that of 13 December 2020, which comprised three distinctive pulses but did not produce an eruptive column detectable from either monitoring cameras or satellites. The 28 February 2021 episode was remarkable for the short amount of time required to reach the climax, and was the most intense, whereas the 12 March 2021 event showed the longest duration but with an intensity between that of the previous two. Our results show that these three paroxysmal events display a typical trend, with the first event also being the smallest in terms of both erupted volume and intensity, the second being the most intense, and the third the one of greatest magnitude but less intense than the second. This is coherent with the end of the first paroxysmal phase on 1 April 2021, which was followed by 48 days of eruptive pause before starting again. In this context, the end of the paroxysmal phase was anticipated by a more effusive episode, thus heralding a temporary decline in the gas content within the feeding magma batch.

Shallow magma dynamics at open-vent volcanoes tracked by coupled thermal and SO2 observations

Open-vent volcanic activity is typically sustained by ascent and degassing of shallow magma, in which the rate of magma supply to the upper feeding system largely exceeds the rate of magma eruption. Such unbalance between supplied (input) and erupted (output) magma rates is thought to result from steady, degassing-driven, convective magma overturning in a shallow conduit/feeding dyke. Here, we characterize shallow magma circulation at Stromboli volcano by combining independent observations of heat (Volcanic Radiative Power; via satellite images) and gas (SO2, via UV camera) output in a temporal interval (from August 1, 2018 to April 30, 2020) encompassing the summer 2019 effusive eruption and two paroxysmal explosions (on July 3 and August 28, 2019). We show that, during the phase of ordinary strombolian explosive activity that preceded the 2019 effusive eruption, the average magma input rate (0.1-0.2 m3/s) exceeds the magma eruption rate (0.001-0.01 m3/s) by ∼2 orders of magnitude. Conversely, magma input and output rates converge to an average of ∼0.4 m3/s during the summer 2019 summit effusion, implying an overall suppression of magma recycling back into the feeding system, and hence of excess degassing. We find that, during the effusive eruption, the peak in SO2 emissions lags behind the thermal emission peak by ∼27 days, suggesting that magma output, feeding the lava flow field, initially dominates over magma input in the conduit. We propose that this conduit mass unloading, produced by this initial phase of the effusive eruption, leads to an overall decompression (of up to 30 Pa/s) of the shallow plumbing system, ultimately causing ascent of less-dense, volatile-rich magma batch(es) from depth, enhanced explosive activity, and elevated SO2 fluxes culminating into a paroxysmal explosion on August 28. Our results demonstrate that combined analysis of thermal and SO2 flux time-series paves the way to improved understanding of shallow magmatic system dynamics at open-vent volcanoes, and of the transition from explosive to effusive activity regimes.

Physical and morphological characterization of the 19 May 2021 ash cloud deposit at Stromboli (Italy)

We report on the ash cloud related to the gravitational collapse of a portion of the Stromboli volcano crater rim that occurred on 19 May 2021. The collapse produced a pyroclastic density current (PDC) that spread along the northwest flank of the volcano and propagated in the sea for about 1 km. The PDC was associated with a convective ash cloud that rapidly dispersed eastward and deposited a thin layer (< 1 mm) of very fine pinkish ash over the village of Stromboli. The deposit was sampled shortly after the emplacement (within a few hours) and prior to any significant reworking or re-sedimentation. We present a comprehensive description of the deposit including dispersal, sedimentological characteristics and textural and geochemical features. We also compare the 19 May 2021 deposit with fine-ash deposits connected to other PDCs and landslides previously occurring at Stromboli and with the distal ash of a paroxysmal explosive eruption of Mt. Etna volcano. Results indicate that the distributions of the mass on the ground and of the grain size are not correlated with the distance from the source. Also, the componentry reflects a preponderance of remobilized material ingested by the PDC. Therefore, the great amount of fine ash can be ascribed to clasts comminution processes, although the prevalence of dense crystalline components records an overall equiaxial shape, revealing a paucity of elongated clast with complex morphology. Furthermore, the outcomes of this work aim to create a collection of data of a co-PDC ash cloud that may prove useful for comparison with other deposits worldwide.

A Petrologic Insight into Transitioning Eruption Styles from the Devil’s Rock Region, Ambae, Vanuatu

AbstractAmbae Island is the largest volcano in the New Hebrides Arc with recent eruptive activity occurring primarily at the summit and along the island’s rift zone. The Devil’s Rock area forms a prominent outcrop on the SW coast. Eruptive deposits here are derived from both strombolian and phreatomagmatic eruptions that contain a similar olivine- and clinopyroxene-rich juvenile basaltic component. This study focuses on a particular transition from strombolian to phreatomagmatic activity to understand if the change in eruption style is a function of magmatic processes or properties (e.g. different composition, ascent rate, degassing history) or if it is driven purely by external factors (e.g. magma–water interaction and/or vent migration). Melts from the strombolian to phreatomagmatic phase record the same melt compositions and volatile contents, suggesting the same magma batch is involved throughout the eruption. More broadly, similarities in H2O, CO2 and S concentrations between olivine- and pyroxene-hosted melt inclusions from Devil’s Rock melt inclusions and those erupted during the 2017–2018 summit eruptions may indicate that a long-term shared magmatic reservoir exists beneath Ambae. Physical characteristics of juvenile tephra including groundmass crystallinity and porosity are combined with melt inclusion compositions to better understand the degassing and crystallisation history and melt evolution of this volcanic system across the transitioning eruptive sequence. Groundmass crystallisation is variable and negatively correlated with connected porosity of erupted scoria reflecting mixing of materials at the vent and inclusion of dense clasts from conduit margins. A direct comparison of crystallinities between strombolian and phreatomagmatic phases reveals higher crystallinity in the strombolian deposits, which is reflective of post-fragmentation crystallisation of clasts. This is particularly evident in the proximal strombolian materials. Qualitative crystallisation textures of melt inclusions are used in a similar fashion to groundmass crystallinities to assess the relative timing of cooling. These trends mirror those of the groundmass and suggest longer cooling times and more effective degassing for samples of the transitional materials. Based on our analysis of deposits at Devil’s Rock, the transition from a strombolian to a phreatomagmatic eruption style was likely driven by groundwater or seawater incursion into the shallow conduit, close to modern-day sea level. Overall, these results suggest a dynamic system where different magmatic cooling histories for strombolian versus phreatomagmatic eruptive phases are reflected in changing groundmass crystallinity. This highlights the propensity for transitions in eruption style over seemingly short time intervals and significantly enhancing eruption explosivity.

Performance of a Rotational Sensor to Decipher Volcano Seismic Signals on Etna, Italy

AbstractVolcano‐seismic signals such as long‐period events and tremor are important indicators for volcanic activity and unrest. However, their wavefield is complex and characterization and location using traditional seismological instrumentation is often difficult. In 2019 we recorded the full seismic wavefield using a newly developed 3C rotational sensor co‐located with a 3C traditional seismometer on Etna, Italy. We compare the performance of the rotational sensor, the seismometer and the Istituto Nazionale di Geofisica e Vulcanologia‐Osservatorio Etneo (INGV‐OE) seismic network with respect to the analysis of complex volcano‐seismic signals. We create event catalogs for volcano‐tectonic (VT) and long‐period (LP) events combining a STA/LTA algorithm and cross‐correlations. The event detection based on the rotational sensor is as reliable as the seismometer‐based detection. The LP events are dominated by SH‐type waves. Derived SH phase velocities range from 500 to 1,000 m/s for LP events and 300–400 m/s for volcanic tremor. SH‐waves compose the tremor during weak volcanic activity and SH‐ and SV‐waves during sustained strombolian activity. We derive back azimuths using (a) horizontal rotational components and (b) vertical rotation rate and transverse acceleration. The estimated back azimuths are consistent with the INGV‐OE event location for (a) VT events with an epicentral distance larger than 3 km and some closer events, (b) LP events and tremor in the main crater area. Measuring the full wavefield we can reliably analyze the back azimuths, phase velocities and wavefield composition for VT, LP events and tremor in regions that are difficult to access such as volcanoes.

Inception and Evolution of La Corona Lava Tube System (Lanzarote, Canary Islands, Spain)

AbstractGrowing interest in studying large terrestrial lava tubes is motivated in part by their analogy with their extra‐terrestrial counterparts. However, on Earth, the formation of such structures is still poorly understood. Here, the lava tube system of La Corona (Lanzarote, Canary Islands, Spain) is studied to identify how pre‐existing stratigraphy can govern a lava tube's evolution. Combining terrestrial laser scanner technology with field observations and geochemical analyses of the pre‐existing lava enabled us to reconstruct the three‐dimensional geometry of the lava tube system, the paleo‐surface trough which it developed, and the volcanic series into which it carved its path. We show that a pyroclastic layer played a key role in the development of the lava tube. The layer—Derived from late Quaternary Strombolian activity—Is traceable along almost the full length of the tube path and defines the paleo‐topography. The excavation process mostly happens because of the mechanical strength of the substrate, that controls the widening of the growing lava tube. Other influential parameters controlling erosion include slope variations of the paleo‐surface (i.e., knickpoints), and the lava physical properties. Since weak layers such as regolith are a common feature of extra‐terrestrial lava flows, the processes seen at La Corona to the may be highly relevant to the development of planetary lava tube systems.