Anatomy of a Paroxysmal Lava Fountain at Etna Volcano: The Case of the 12 March 2021, Episode

On 13 December 2020, Etna volcano entered a new eruptive phase, giving rise to a number of paroxysmal episodes involving increased Strombolian activity from the summit craters, lava fountains feeding several-km high eruptive columns and ash plumes, as well as lava flows. As of 2 August 2021, 57 such episodes have occurred in 2021, all of them from the New Southeast Crater (NSEC). Each paroxysmal episode lasted a few hours and was sometimes preceded (but more often followed) by lava flow output from the crater rim lasting a few hours. In this paper, we use remote sensing data from the ground and satellite, integrated with ground deformation data recorded by a high precision borehole strainmeter to characterize the 12 March 2021 eruptive episode, which was one of the most powerful (and best recorded) among that occurred since 13 December 2020. We describe the formation and growth of the lava fountains, and the way they feed the eruptive column and the ash plume, using data gathered from the INGV visible and thermal camera monitoring network, compared with satellite images. We show the growth of the lava flow field associated with the explosive phase obtained from a fixed thermal monitoring camera. We estimate the erupted volume of pyroclasts from the heights of the lava fountains measured by the cameras, and the erupted lava flow volume from the satellite-derived radiant heat flux. We compare all erupted volumes (pyroclasts plus lava flows) with the total erupted volume inferred from the volcano deflation recorded by the borehole strainmeter, obtaining a total erupted volume of ~3 × 106 m3 of magma constrained by the strainmeter. This volume comprises ~1.6 × 106 m3 of pyroclasts erupted during the lava fountain and 2.4 × 106 m3 of lava flow, with ~30% of the erupted pyroclasts being remobilized as rootless lava to feed the lava flows. The episode lasted 130 min and resulted in an eruption rate of ~385 m3 s−1 and caused the formation of an ash plume rising from the margins of the lava fountain that rose up to 12.6 km a.s.l. in ~1 h. The maximum elevation of the ash plume was well constrained by an empirical formula that can be used for prompt hazard assessment.

A model for gravity changes induced by lava fountaining at Mt Etna

<div> <p><span>Lava fountains represent a common eruptive phenomenon at basaltic volcanoes, which consist of jets of fluid lava ejected into the atmosphere from active vents or fissures. They are driven by rapid formation and expansion of gas bubbles during magma ascent. The dynamics of lava fountains is thought to be controlled by the gas accumulation in the foam layer at the top of a shallow magmatic reservoir, which eventually collapses triggering the lava fountaining. Gravity measurements taken from a location close to summit of Mt. Etna during the 2011 lava fountain episodes showed a pre-fountaining decrease of the gravity signal. The interplay between gas accumulation in the foam layer and its subsequent exsolution in the conduit has been interpreted as the mechanism producing the gravity decrease and eventually leading to the foam collapse and onset of the lava fountaining activity. Gravity measurements have proved helpful in recording the earliest phases anticipating the lava fountain episodes and inferring the amount of gas involved. However, more accurate estimates of the accumulating and ascending gas volume and total magma mass require considering the possible effect of non-spherical magma chamber geometries and magma compressibility. </span></p> </div><div> <p><span>Under task 4.4 of the H2020 NEWTON-g project, we are accomplishing a detailed study aimed to simulate the gravity signal produced in the stage prior to a lava fountain episode, through a magma chamber - conduit model. We use a prolate ellipsoidal chamber matching the inferred shape of the shallow chamber active at Mt. Etna during the lava fountain episodes, and calculate the surface gravity changes induced by inflow of new magma into the chamber-conduit system. We use a two-phase magma with fixed amount of gas mass fraction and account for magma compressibility. We find that a realistic chamber shape and magma compressibility play a key role and must be considered to produce realistic gravity changes simulations. We combine our physical model with empirical distributions of recurrence time and eruption size of the past lava fountains at Mt. Etna to stochastically simulate realistic time series of gravity changes. The final goal of this study is to develop a prediction model for the amount of magma and duration of lava fountains at Mt. Etna.</span></p> </div>

Multisensor Characterization of the Incandescent Jet Region of Lava Fountain-Fed Tephra Plumes

Explosive basaltic eruptions eject a great amount of pyroclastic material into the atmosphere, forming columns rising to several kilometers above the eruptive vent and causing significant disruption to both proximal and distal communities. Here, we analyze data, collected by an X-band polarimetric weather radar and an L-band Doppler fixed-pointing radar, as well as by a thermal infrared (TIR) camera, in relation to lava fountain-fed tephra plumes at the Etna volcano in Italy. We clearly identify a jet, mainly composed of lapilli and bombs mixed with hot gas in the first portion of these volcanic plumes and here called the incandescent jet region (IJR). At Etna and due to the TIR camera configuration, the IJR typically corresponds to the region that saturates thermal images. We find that the IJR is correlated to a unique signature in polarimetric radar data as it represents a zone with a relatively high reflectivity and a low copolar correlation coefficient. Analyzing five recent Etna eruptions occurring in 2013 and 2015, we propose a jet region radar retrieval algorithm (JR3A), based on a decision-tree combining polarimetric X-band observables with L-band radar constraints, aiming at the IJR height detection during the explosive eruptions. The height of the IJR does not exactly correspond to the height of the lava fountain due to a different altitude, potentially reached by lapilli and blocks detected by the X-band weather radar. Nonetheless, it can be used as a proxy of the lava fountain height in order to obtain a first approximation of the exit velocity of the mixture and, therefore, of the mass eruption rate. The comparisons between the JR3A estimates of IJR heights with the corresponding values recovered from TIR imagery, show a fairly good agreement with differences of less than 20% in clear air conditions, whereas the difference between JR3A estimates of IJR height values and those derived from L-band radar data only are greater than 40%. The advantage of using an X-band polarimetric weather radar in an early warning system is that it provides information in all weather conditions. As a matter of fact, we show that JR3A retrievals can also be obtained in cloudy conditions when the TIR camera data cannot be processed.

Volcanic Jet Noise from the Kilauea Fissure 8 Lava Fountain

<p>Seismic and acoustic signals are important for remote real time and post-eruption analysis of volcanic eruptions. To properly interpret these signals it is critical to connect their characteristics with eruption parameters. In this study, we present an analysis of the infrasound emissions by the sustained lava fountain at Fissure 8 during the 2018 eruption of Kilauea Volcano, Hawaii. This eruption was one of the largest and most destructive events in Hawaii’s historic times. Large (35.5 km2) lava flows covered much of the Lower East Rift Zone (LERZ) and destroyed property and infrastructure. This activity was dominated by high lava effusion rates at Fissure 8 and lava fountains up to 80 m tall. The energetic output of gas and lava produced sustained, broadband acoustic waves which were recorded by a four-element infrasound array deployed 0.6 km northwest of the fountain. The spectrum of the infrasound is similar to that of man-made jets and is termed volcanic jet noise. We compare the spectrum of the recorded infrasound signal with models developed for man-made jets such as rockets and jet engines. These models predict different spectral shapes for fine scale turbulence (FST), produced by incoherent movement of the gases, and large scale turbulence (LST), produced by coherent instability waves. The dominance of one or the other turbulent noise source is highly directional. We compare the infrasonic signals with observations of fountain properties, such as pyroclast velocity and height, to help understand the jet noise signals and determine quantitative fountain properties from the infrasound. The results of this work will contribute to the understanding of the physics of lava fountain sound generation, its dependence on eruption parameters, and ultimately provide a tool for rapid assessment of eruption style and dynamics.</p>

Unravelling the links between seismo-acoustic signals and eruptive parameters: Etna lava fountain case study

Deriving eruption source parameters from geophysical data is critical for volcano hazard mitigation, yet remains a challenging task in most volcanoes worldwide. In this work, we explored the temporal relationship between geophysical signals and eruptive parameters measured during six explosive episodes from the New South-East Crater of Mt. Etna (Italy). The quadratic reduced seismic velocity and pressure were calculated to track the temporal variation of volcanic elastic radiation, and the lava fountain height was estimated by thermal camera image processing. The temporal relationships between these geophysical and eruptive time series were studied. In particular, the first considered lava fountain exhibited a “clockwise hysteresis” pattern: higher seismic amplitude with respect to the fountain height during the waxing phase as compared to the waning phase. We also calculated the regression parameters for both linear and power laws, linking seismo-acoustic and eruptive time series. For the linear regressions, we found fairly constant values of the scaling factors in five out of six eruptive episodes, which can be considered as a promising step to derive eruption source parameters from geophysical data in real-time. Regarding power law regressions, a clear relationship was observed between the exponents determined for the power law linking quadratic reduced velocity and lava fountain height, and the time interval duration from the previous eruption. These results suggest that the condition of the uppermost part of the plumbing system (e.g. viscosity of residing magma and conduit conditions) play a key role in the seismic energy generation during the eruptions.

Changes in SO2 Flux Regime at Mt. Etna Captured by Automatically Processed Ultraviolet Camera Data

We used a one-year long SO2 flux record, which was obtained using a novel algorithm for real-time automatic processing of ultraviolet (UV) camera data, to characterize changes in degassing dynamics at the Mt. Etna volcano in 2016. These SO2 flux records, when combined with independent thermal and seismic evidence, allowed for capturing switches in activity from paroxysmal explosive eruptions to quiescent degassing. We found SO2 fluxes 1.5–2 times higher than the 2016 average (1588 tons/day) during the Etna’s May 16–25 eruptive paroxysmal activity, and mild but detectable SO2 flux increases more than one month before its onset. The SO2 flux typically peaked during a lava fountain. Here, the average SO2 degassing rate was ~158 kg/s, the peak emission was ~260 kg/s, and the total released SO2 mass was ~1700 tons (in 3 h on 18 May, 2016). Comparison between our data and prior (2014–2015) results revealed systematic SO2 emission patterns prior to, during, and after an Etna’s paroxysmal phases, which allows us to tentatively identify thresholds between pre-eruptive, syn-eruptive, and post-eruptive degassing regimes.