Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# Magma pressure discharge induces very long period seismicity

## Abstract

Volcano seismicity is one of the key parameters to understand magma dynamics of erupting volcanoes. However, the physical process at the origin of the resulting complex and broadband seismic signals remains unclear. Syn-eruptive very long period (VLP) seismic signals have been explained in terms of the sudden expansion of gas pockets rising in the liquid melt. Their origin is linked to a magma dynamics which triggers the explosive process occurring before the explosive onset. We provide evidence based on acoustic, thermal, and ground deformation data to demonstrate that VLP signals at Stromboli are generated at the top of the magma column mainly after the explosion onset. We show that VLP amplitude and duration scale with the eruptive flux which induces a decompression of 103–104 Pa involving the uppermost ~ 250 m of the feeding conduit. The seismic VLP source represents the final stage of a ~ 200 s long charge and discharge mechanism the magma column has to release excess gas accumulated at the base of a denser and degassed magma mush. The position of the VLP seismic source coincides with the centroid of the shallow mush plug and tracks elevation changes of the magma free surface.

## Introduction

Open vent volcanoes are generally characterized by a sustained and restless degassing dynamics which generates a large variety of seismic signals. Moderate explosive activity is generally explained in terms of nucleation, ascent, and bursting of a large gas slug1, which is thought to be responsible for the origin of the seismic VLP activity on many volcanoes on the globe2,3,4,5,6,7,8,9,10,11. However, the genesis of VLP seismic signals is still matter of debate.

Seismic VLP signals are commonly recorded on a large variety of magmatic systems, ranging from high-viscosity andesitic9,12 to low-viscosity basaltic volcanoes6,13,14. However, VLP signals share similar seismic characteristics2,3,4,5,6,7,8,9,10,11,15,16. Their amplitude quickly decays moving away from the craters and are thus mainly recorded at a distance of few kilometers from the vent. Frequency range is confined in the 0.01–0.2 Hz band and have wavelength tens of kilometers long. Location of the VLP source is generally shallow and confined within the first hundreds of meters below the craters. Waveform of seismic VLP signal is generally simple and with a strong rectilinear polarization vector pointing towards the craters, typical of P-waves. These similarities suggest a link to a more general dynamics which seems to be common to many active volcanoes in the world3,4,5,13,15,16.

When seismic VLPs are observed as associated with explosive activity and/or gas puffing at the surface their origin is commonly explained as the opening-closing of a crack induced by the sudden volume expansion within the fluid melt several seconds (between 3 and 60 s) before the explosion2,3,4,5,6,7,8,9,10,11,13,15,16,17. However, seismic VLPs are not always associated with visible explosive activity18,19, thus suggesting different mechanisms acting inside the conduit5. Here we use data collected at Stromboli volcano (Aeolian Islands, Italy) to shed light on the origin of syn-eruptive VLP seismic activity using a multiparametric approach.

Explosive activity at Stromboli volcano is generally moderate (VEI = 1) and persistent through time20 repeating at a mean rate of ~ 13 events/h21 (termed ‘normal activity’20,21) which provides evidence of an efficient degassing mechanism of the magma column, typical of most basaltic open-conduit volcanoes. At Stromboli, seismic VLP waveforms repeat in time in well-constrained families2,18 indicating a source mechanism typical of a non-destructive process acting in a quasi-stationary dynamics2,18,22,23. Semblance18,24 and seismic moment tensor inversion2 locate the centroid of the VLP seismic source at a shallow ~ 220–260 m depth below the craters and at ~ 160 m NW off the vents2,18 within the Sciara del Fuoco. This stability in time suggests the presence of a persistent physical constraint on the source process and location of the VLPs. At this depth, a geometrical change in the feeding conduits is believed to induce the sudden expansion of the rising gas slug and the consequent magma decompression recorded by the seismometers as VLP signal2,25. Moment tensor inversion and VLP source location are thus interpreted as evidence of an inclined dyke which is responsible for the sudden expansion of the slug while rising to the surface2. This interpretation has contributed to picturing the geometry of the shallow feeding system at Stromboli and to explaining the VLP source process as a pressure source within an opening/closing crack embedded in the volcanic edifice2,4,26. This model times the source of the VLP to seconds before the explosion (bottom-up model2) and would require a very fast ascent gas velocity (50–100 m/s)10,27 to explain the link with the explosive process as evidenced by camera and infrasound records of strombolian eruption at Stromboli, Kilauea28, Fuego4, Aso10, and Erebus3 volcanoes.

At the scale of several years, at Stromboli, the position of the VLP seismic source, as monitored by the dip of the VLP polarization seismic vector, shows large (~ 10 degrees) relative changes18,29,30,31 that are well correlated with long-term changes in the topographic elevation of summit craters30. Larger fluctuations in the VLP polarization vector occur during lateral effusive eruptions when strombolian activity is no longer visible at the summit craters. Deepening of the polarization vector has been correlated to the discharge of the magma shallow reservoir and explained by the lowering of the magma column within the feeding conduits30,31,32.

In this paper we present evidence incorporating thermal imagery, acoustic pressure, and ground deformation data to produce a new physical model that explains the seismic VLP as the final stage of recurrent charge/discharge of the excess gas accumulated in the shallow part of the feeding conduit, connecting the VLP source to the surface of the magma column and to the eruptive flux.

## Results

### The seismic VLP wavefield

We here define VLP seismic displacement as the integration of ground velocity filtered between 0.05 and 0.2 Hz. In harmony with previous experiments2,18 we derive two different families of VLP seismic signals (Figure S1). The vertical seismic component of both family shows coherent polarity for all stations and is characterized by an initial low amplitude and compressional onset (A1 in Fig. 1) lasting ~ 5–6 s that is followed by a second, larger amplitude, dilatation phase (A2 in Fig. 1) lasting ~ 4–12 s (see also Ref.2,18,24). This energetic second part of the waveform is sometimes followed by a coda 20–30 s long (A3) characterized by one to two low amplitude oscillations.

Once the VLP wavefield is rotated into the vent-radial direction, ground displacement polarities at the stations are coherent and are consistent with a roughly isotropic source producing an alternating compression-dilatation mechanism that is nearly symmetrical with respect to the crater position2,24. Although Stromboli volcano is characterized by three main active vents (Supplementary note S1) well-separated in space (~ 200 m), the VLP radiation patterns do not seem to depend on the particular exploding vent, but are quite stable18 and oriented radially with respect to the central crater area (Figure S1).

### VLP seismicity and thermal flux

Thermal radiance reflects two different phases of the explosive dynamics. While the initial increase in temperature is due to hot ash and scoria fragments ejected during the gas thrust phase (A2, Fig. 1), the waning phase of the thermal temperatures (A3 Fig. 1) is linked to plume expansion during the buoyancy of the gas-fragments cloud33,34,35. Explosive eruptions are detected as variation of thermal radiance providing estimates of the onset, duration, and thermal mass flux associated with the explosion. In harmony with previous observations made at Stromboli2,18, as well as on other volcanoes3,11,28, thermal radiance shows that seismic VLP signals are well correlated with explosive events, as they always start at the end of VLP seismic compression phase (A1 in Fig. 1a).

The link between VLP seismicity and explosive activity is better evidenced in the long-term time scale when the daily numbers of VLP30,31,35 are compared to the number of daily explosions observed by the thermal camera (ROC in Figure S2) of the permanent network (Fig. 2b). Thermal transients and VLP rates change together over days to months and highlight the fluctuations of the explosive level during normal activity (Fig. 2b).

Using high-frame rate (50 Hz) thermal images35 taken in line of sight with the exploding vent (PSF, Figure S2), we found that the onset of the thermal signal marks the transition between the compressional (A1) and the large dilatation (A2) phase of the seismic VLP (Fig. 1). In particular, the duration Δt of the thermal gas-thrust phase (from the onset, #2, to the maximum recorded temperature, #4, in Fig. 1) shows (Fig. 3) a very good correlation of R = 0.82 (Pearson correlation coefficient) with the duration of the A2 VLP dilatation phase measured on the vertical component of the ground displacement (Fig. 3b). The correlation becomes straightforward (R = 0.89) when the heat flux Φ, and the VLP amplitude of the A2 phase are compared (Fig. 3c).

Here we have considered the maximum of the heat fluxes Φ (in W) associated with explosions which have been calculated as, $$\Phi ={A}_{FOV}\sigma \epsilon \left({T}_{E}^{4}-{T}_{BKG}^{4}\right)$$ where TE is the peak of the temperature integrated within of the camera’s field of view (AFOV), TBKG is background temperature, σ is the Stefan-Boltzmann constant, and ε is emissivity (~ 0.9637).

The thermal camera therefore indicates that explosive dynamics tracks the amplitude and duration of the A2 dilatation phase (Figs. 1 and 3a) of the seismic VLP waveform. Such temporal correlation between thermal imagery data and VLP waveform is observed nearly for all explosions and it constitutes a key evidence for linking the source process to explosive activity. These observations suggest that the violent gas expansion during the gas thrust phase (Fig. 1), produces a downward momentum transferred to the volcanic edifice that is detected by seismometers as the VLP displacement field38,39 (Fig. 1d).

### VLP seismicity and acoustic emission

Acoustic waves provide a second body of evidence to associate the seismic VLP signal with explosive process occurring at the surface. Acoustic pressure (p) is explained in terms of the rate of acoustic volumetric flux (Qi)40,41,42 generated during the violent expansion of the gas-fragments mixture within the conduit, first, and into the atmosphere afterwards. Acoustic pressure thus reflects the history by which the magmatic column is reducing its internal pressure, and is proportional to the rate at which gas expands outside the conduit35,43. Here, we show that when the acoustic signal is corrected for the propagation time between the source and the station, the onset of the acoustic wave coincides with the beginning of the seismic VLP large (A2) dilatation phase (Fig. 1). Considering the topography, the conduit geometry and the scattering effects on the infrasonic amplitude p(t), the volumetric flux inside the conduit Qi can be calculated as43:

$${Q}_{i\left(t\right)}=\frac{2\pi r}{\left(1+\left|{R}_{a}\right|\right)\rho \alpha 1{0}^{\left(\frac{IL}{20}\right)}}{\int }_{0}^{t}p\left(t+\frac{r}{c}\right)dt$$
(1)

where ρ = 1.1 kg/m44 is the air density, r (~ 300 m) is the crater-station slant distance, c = 340 m/s is the sound speed in air, IL is the insertion loss of ~ 1 dB estimated for the specific source-station path, and |Ra| (= 0.96) is a parameter which takes into account the vent exit radius and frequency content of infrasonic signal. We here considered an isotropic radiation pattern (α = 1), which has proven to be valid for explosions at Stromboli43.

Our results (Fig. 4) indicate that acoustic volumetric flux ranges between 200 and 1500 m3/s and is fully comparable with gas flux directly measured by UV camera45,46. Acoustic flux correlates (R = 0.76) with the amplitude of the VLP seismic A2 dilatation phase and its duration scales with the total duration of the acoustic signal (R = 0.68). Consistent with the thermal analysis, these observations support the hypothesis that the VLP phase A2 (Fig. 1) originates during the gas thrust regime of the explosive dynamics, and its amplitude reflects the upward volumetric flux of the explosive fluid.

### VLP seismicity and ground deformation

Seismic VLP signals are at Stromboli always associated with a well resolved deformation of the ground (Fig. 2a). This process repeats in time for every explosion (Fig. 2a) in a cyclic charge/discharge process47. Considering a propagation P-wave velocity of ~ 3500 m/s2, the wavelength of the VLP seismic radiation is tens of km (~ 40 km), which compared with the distance from the source (< 500 m), indicates that we are in a near-field condition. In such a case, seismic VLP displacement reflects the source time function and can be considered to be a time dependent quasi-static volume displacement of the source2,3,16,48,49,50. Near-field conditions also make it possible to convert the low frequency band seismic signal into ground tilt using analytical solutions47,51,52.

We argue that the VLP A2 dilatation phase is reflecting the vertical force induced by the mass flux associated with the explosive ejection of the gas-fragments mixture, and we thus apply a model of quasi-static deformation to estimate the location and the associated volumetric change. We used data recorded by a network of five temporary seismic stations located within 300 m from the active vents with a ~ 70% azimuthal coverage around the craters (Fig. 5). To the best of our knowledge, these are the first measurements of the VLP wavefield at Stromboli made inside the crater area using a broadband network. This allowed us to record VLP signals in the very near-field with the best azimuthal coverage as possible, and to constrain the source position by minimizing as much as possible the effect induced by the network geometry.

Considering the open-conduit model53, the vertical VLP displacement uz can be defined as proportional to the volumetric change ΔV of an L-extended conduit:

$${u}_{z}=\frac{5}{3}G\frac{\Delta V}{\pi L}$$
(2)

where G is the Green’s function54, which depends on the position of the conduit relative to the stations. We apply a grid search procedure on a 3D space within the seismic network to find that VLP deflation is very well represented (best fit of > 0.9) by a conduit with a centroid located near to the SW crater and extending from 750 to 500 m above sea level (that is from 0 to 250 m below the craters). The epicentral location differs by about 200 m from previous semblance analysis2,18, and it is in agreement with the radial particle motion of the VLP onset (Fig. 5) and moment tensor inversion based on the same data set55. The source that best fits the observed displacements coincides with the last L = 250 m of the conduit and is associated with a volume change ΔV of ~ 33 ± 10 m3 (Fig. 5). These deflation volumes are in good agreement with previous results inferred by seismic2,56 and ground tilt modeling48.

Assuming that this volumetric change is exerting a perturbation of the shallow magma conduit, the volume change is equivalent to a decompression ΔP:

$$\Delta P=\frac{5}{3}\mu \frac{\Delta V}{\pi {a}^{2}L}$$
(3)

of 8.5 ± 3.6 × 103 Pa, where μ = 1.3 GPa is the rigidity of the host rocks47, a = 87 m is the equivalent radius of the shallow reservoir inferred by modeling the last four effusive eruptions32. Notably, this decompression is equivalent to a force of ~ 108 N which is in the same order of magnitude as the single force calculated by moment tensor inversion2.

At the explosive onset, the momentum transferred by the explosive process to the magma column can be described as the downward reaction force per unit area $$\Delta P=\rho {U}^{2}$$ produced by the eruptive jet (gas-ash-and-lava fragments) with density ρ ejected at a velocity U. Assuming that the pressure drop, calculated by quasi-static VLP displacement, represents this reaction force ΔP, for a mean exit velocity U of gas and fragments of 40 m/s (Ref.21), the plume bulk density is ~ 5 kg/m3, which is in line with previous estimates of plume (gas-ash-and-lava fragments) density at Stromboli57. This explains the good fit (R = 0.76) between the VLP vertical displacement and the infrasonic over-pressure (Fig. 4) produced by the acoustic volumetric flux Qi (Eq. 1). Ground displacement during the dilatation A2 phase of the VLP is representing the magma depressurization of the last ~ 250 m feeding conduit (i.e., shallow reservoir) in response to the gas released during Strombolian eruptions.

## Discussion

Quasi-static displacement modelling of vertical ground motion54,57, using an unprecedented seismic network deployment in very near field conditions, indicates that the source of seismic VLP is located within the crater area and it has an extension of ~ 250 m below the craters. Families of VLPs, also when associated with different explosive craters, have in the long-term26 the same back-azimuth which does not allow to discriminate between different craters (Figure S1). We argue that particle motion of VLPs at Stromboli is indicating the centroid representing the base of a shallow crystal-rich magma plug elongated in the NE direction with equivalent radius of 87 m30,32, common to all the craters, which is excited by the explosive process (Figure S1).

Thermal radiances and infrasonic waves show how the explosive onset is preceded by a small compressional VLP seismic phase (A1) lasting ~  < 5 s (Fig. 1), which is anticipating the explosion visible at the surface by several seconds (between 3 and 6 s). This time delay between the VLP (A1) and explosive onsets (t1 in Fig. 6) suggests that seismic VLP originates before the explosion is occurring at the surface and it has been inferred to reflect the ascent, expansion, and burst of a gas slug within the last portion of the magma conduit57.

According to the “bottom-up” model2, gas rising from depth44 pushes up the magma within an inclined dyke (that resides at a moment-centroid constrained depth near 220—260 m below and ~ 160 m northwest of the vents)2, generating a volumetric expansion and the compressional onset (A1 in Fig. 1). This model, however, requires unrealistically fast fluid velocities of ~ 70 m/s27 to justify gas slug propagation from the VLP location to the surface (~ 250 m)2,18 in such a short time difference of ~ 3–5 s between the VLP and the explosion onset. The timing of the source dynamics observed at Stromboli is similar to that at other volcanoes when available video and/or acoustic sensors allow to associate the explosive process with VLP activity. At Kilauea28, Erebus3,58, Fuego4 and Aso10 volcanoes, VLP seismicity mainly develops after the onset of the explosion and the short time lag between VLP and explosion onsets is difficult to reconcile with the deep position of the seismic source using a classic fluid-dynamics model.

Ground tilt at Stromboli shows that the explosive phenomena is instead associated with a ~ 200 s (in average) long inflation-deflation deformation cycle. Ground starts to inflate ~ 150 s before the explosive onset (A0 in Figs. 6 and 7b) in response to the increase of the magmastatic head due to the uplift of a “cap”, driven by the expansion of the gas which accumulates below a crystal-rich and denser magma mush59, stagnant in the shallow portion (< 250 m) of the magma column47,54. While gas is expanding upwards in the magma mush (Fig. 7b), magmastatic pressure increases and ground inflation accelerates almost ~ 10 times in the last 10–20 s47 as the gas gets closer to the surface54,60 (A1 in Figs. 6 and 7b). Considering this dense crystal-rich mush is ~ 250 m thick, the duration of the inflation (~ 150 s) indicates that gas propagates at an average velocity of ~ 1.6 m/s25.

Seismic VLP sits on top of this inflation/deflation cycle and explosion occurs at the end of the ground inflation (Fig. 7a). The small compressional onset of the VLP represents only the final stage of the acceleration of the gas expansion just before the explosion. Laboratory experiments show, in fact, that when pressure of the gas overcomes the effective tensile strength of a granular suspension, gas propagates upwards inducing the inflation of the particle pack which accelerates while approaching the free surface61,62. Inflation/deflation tilt cycles associated with seismic VLP and explosive activity have also been observed at Fuego volcano and are explained as the pressurization/depressurization of the uppermost (~ 500 m) highly viscous plug of magma4. Similarly, in this silicic magma, seismic VLPs occur at the end of the inflation and are interpreted as the explosive depressurization of the conduit induced by the brittle failure of the viscous plug4.

Once the gas-fragments cloud is released in the atmosphere, magma level suddenly drops and the ground rapidly deflates47,54. Pressure inside the magma decreases, the conduit walls contract and the shallow crystal-mush moves downward to recover the equilibrium conditions (Fig. 7b). Ground deformation (Fig. 6a) shows that magma decompression persists (~ 30 s on average) longer than the duration of the explosions (~ 4 s on average) during which pressure oscillates (A3 in Figs. 1 and 6c) for 20–30 s before the magmastatic equilibrium is restored.

Oscillation of the VLP coda with ~ 20 s period and lasting several minutes have also been observed at Kilauea13, Erebus3 and Fuego4 volcanoes as associated with explosive phenomena and explained as induced by the magma recharge force58 and repressurization of the conduit4. The origin of the coda oscillations has been related to the dimension of the conduit or to the viscoelastic properties of the magma. At Stromboli, UV camera observations show that deflation of the ground matches with a ~ 30 s long SO2 gas emission following the explosion63 (Fig. 6b) and reflecting a long-lasting post-explosive decompression of the magma column. Coda oscillations are clear in the VLPs associated with large acoustic pressure (~ 40 Pa at 100 m) and are almost absent when VLP deflation is small and acoustic pressure low (~ 4 Pa at 100 m). This evidence (Fig. 6d and Fig. S1) suggests that when eruptive dynamics is linked to a small mass flux, the induced pressure drop is too small to trigger the oscillation of the shallow magma column. Our observations are therefore consistent with a “top-down” dynamic framework28, where VLP seismicity is reflecting only the final stage (< 20 s) of the charge/discharge degassing process of the magmatic column, and it represents a small perturbation in the long lasting (~ 200 s) pressurization cycles of the conduit (Fig. 7). This process repeats in time for every single explosion in a recurrent gas charge/discharge cyclic process.

The rate of VLP seismicity is in fact well correlated to the rate of explosions visible at the surface (Fig. 2). This correlation is lost only during effusive eruptions18,30,32, when the explosive activity at the summit craters stops. At the eruption onset the dip angle of VLP polarization moves downwards of several degrees (~ 10°), tracking the drainage of the shallow magmatic system30,31,32. Nevertheless, explosive degassing is still detected using SO2 imagery, at a rate that is well above the average (up to 60% of the total flux)36, suggesting that explosive activity is still on-going deep in the conduit. At the end of effusive eruptions, explosive activity progressively resumes at the summit crater, following the progressive refilling of the magma column within the conduit30,31,64, and the explosion-VLP correlation is then restored.

Our observations support the conclusion that VLP displacement wavefield is generated in a very shallow portion of the conduit at a rate proportional to the gas flux36,61 and it can be used to track the uppermost position of the magma column. Our result opens new pathways to the interpretation of seismic VLP source process at other open-vent basaltic volcanoes, and indicates that explosive activity can also occur when magma column is far from the surface and deep inside the conduit.

## Methods

### Experiment design and data acquisition

The geophysical data presented in this work are the result of experiments specifically designed to disclose the connection between VLP (0.05–0.2 Hz) seismicity and explosive activity. We used temporary stations located at the closest realizable distances to the craters and with the best azimuthal coverage possible considering safety constraints. These results are integrated with data from the permanent stations of the network operated by the Laboratory of Experimental Geophysics (LGS) of the University of Florence (http://lgs.geo.unifi.it) since 200365. Seismic broad-band stations, thermal cameras, and acoustic stations were placed at distances ranging from < 100 m to 500 m from Stromboli’s active craters (Figure S2). All stations were synchronized using specific calibration tests and sampled at a variable rate according to the nature of the parameters but with the aim to guarantee a minimum of 0.01 s time accuracy and consistency.

We used a total of nine seismic stations consisting of Guralp CMG-40T broadband seismometers, with 30 s and 60 s eigenperiods and sensitivities of 800 V/m/s. All seismic stations were provided with differential acoustic pressure transducers with sensitivities of 25 mV/Pa, flat frequency responses between 0.01 and 100 Hz, and a full-scale range of 250 Pa. Together with the 5-element permanent infrasonic array, the acoustic wavefield was recorded with a network of 15 acoustic stations43. Ground deformation was measured using the 3 permanent tiltmeters placed in ~ 5 m depth boreholes sampled at 1 Hz with a nominal sensitivity of 1 nanoradian47,54.

Thermal images were recorded from the ROC (permanent) and PSF sites (Figure S2) with a FLIR-A20 thermal camera fitted with a 34° × 25° optical lens (9.2 mm), and 0.1 °C thermal resolution, with focal plane array of 160 × 120 pixels, and sensitive in the 7.5–13 µm spectral range. Thermal images from PSF were collected at a frame rate of 50 Hz and were synchronized using a GPS clock with an accuracy of ~ 5 ms. The PSF thermal camera was located in direct line of sight with the NE and SW craters at a slant distance of 345 m and 327 m, respectively, and thus with field of views of 218 × 153 m and 206 × 145 m35.

### VLP localization using open conduit deformation model

The source parameters associated with VLP displacement field were estimated by integrating VLP ground velocity signals recorded at five temporary broadband seismic stations (yellow dots in Figure S2) together with the permanent stations (yellow square in Figure S2) of the LGS multiparametric network54.

Following an open conduit deformation model54,57, for an ideal Poissonian solid (Poisson modulus n = 0.25), the vertical displacement uz is related to the volumetric change ΔV of a portion $$L={z}_{2}-{z}_{1}$$ of a cylindrical conduit by Eq. (2), where $${z}_{1}$$ and $${z}_{2}$$ are the top and the bottom depths, respectively. The Green’s function G in Eq. (2) is relative to the station i located at a distance x from the conduit and is defined as53,60:

$$G=\frac{-1}{2\left({z}_{2}-{z}_{1}\right)}\left[\frac{-{z}_{2}^{2}\left({z}_{2}-{z}_{1}\right)}{{R}_{2}^{3}}+\left(2\nu -1\right)\left(\frac{{z}_{2}}{{R}_{2}}-\frac{{z}_{1}}{{R}_{1}}\right)-2\nu {z}_{1}\left(\frac{1}{{R}_{2}}-\frac{1}{{R}_{1}}\right)-\left(2\nu -1\right)ln\left(\frac{{z}_{2}+{R}_{2}}{{z}_{1}+{R}_{1}}\right)\right]$$
(4)

where $${R}_{n}={\left({x}^{2}+{z}_{n}^{2}\right)}^\frac{1}{2}$$ are the associated slant distances between the top (n = 1) and the bottom (i = 2) of the conduit and the station (Figure S3).

We applied a search algorithm for the best source conduit position within a 10 m-spaced grid of 900 × 450 m centered on the summit craters and extending on the Sciara del Fuoco, and assuming that z1 can vary, with step of 10 m, within 0 and 10,000 m below the surface, and L can vary between 10 and 1000 m with steps of 10 m. The best solution is thus estimated through a least square minimization among all the possible solutions between the ratios of the Green functions Gn/Go and the measured VLP displacement ui/uo at the i station, where Go and uo are relative to the reference station. Once the best source and conduit positions as well as the conduit length L are found, the volumetric change ΔV is calculated by the slope of Eq. (2). Predicted VLP displacements differ by less than 10% respect to the measured VLP displacements, highlighting a strong linear correlation (R > 0.9) between source solution and observations (Fig. 5).

Our open conduit deformation model does not account for topographic effects on the measured VLP displacements. However, Finite Element Modeling (FEM) accounting for topography returned very comparable results, suggesting that while topography may have a large effect on displacement direction, it plays minor role (< 15%) in the associated amplitude54.

## Data availability

All data generated or analysed during this study are included in this published article (and its Supplementary Materials) or available from the corresponding author upon reasonable request.

## References

1. 1.

Blackburn, E. A., Wilson, L. & Sparks, R. S. J. Mechanisms and dynamics of strombolian activity. J. Geol. Soc. Lond. 132, 429–440 (1976).

2. 2.

Chouet, B. et al. Source mechanisms of explosions at Stromboli Volcano, Italy, determined from moment tensor inversions of very long period data. J. Geophys. Res. 108(B1), ESE 7-1-ESE 7-25 (2003).

3. 3.

Aster, R. et al. Very long period oscillations of Mount Erebus Volcano. J. Geophys. Res. 108(B11), ESE6-1-ESE6-22 (2003).

4. 4.

Lyons, J. J. & Waite, G. P. Dynamics of explosive volcanism at Fuego volcano imaged with very long period seismicity. J. Geophys. Res. 116(B09303), 1–18 (2011).

5. 5.

Maeda, Y. & Takeo, M. Very-long-period pulses at Asama volcano, central Japan, inferred from dense seismic observations. Geophys. J. Int. 185(1), 265–282 (2011).

6. 6.

Ohminato, T., Chouet, B. A., Dawson, P. B. & Kedar, S. Waveform inversion of very long period impulsive signals associated with magmatic injection beneath Kilauea Volcano, Hawaii. J. Geophys. Res. 103(B10), 23839–23862 (1998).

7. 7.

Jolly, A. D. et al. Insights into fluid transport mechanisms at White Island from analysis of coupled very long-period (VLP), long-period (LP) and high-frequency (HF) earthquakes. J. Volcanol. Geotherm. Res. 343, 75–94 (2017).

8. 8.

Jolly, A., Lokmer, I., Christenson, B. & Thun, J. Relating gas ascent to eruption triggering for the April 27, 2016, White Island (Whakaari), New Zealand eruption sequence. Earth Planets Space 70(1), 1–15 (2018).

9. 9.

Hidayat, D., Chouet, B., Voight, B., Dawson, P. & Ratdomopurbo, A. Source mechanism of very long period signals accompanying dome growth activity at Merapi volcano, Indonesia. Geophys. Res. Lett. 29(23), 33-1-33–4 (2002).

10. 10.

Ishii, K. et al. Gas flow dynamics in the conduit of Strombolian explosions inferred from seismo-acoustic observations at Aso volcano, Japan. Earth Planets Space 71(1), 1–15 (2019).

11. 11.

Waite, G. P., Nadeau, P. A. & Lyons, J. J. Variability in eruption style and associated very long period events at Fuego volcano Guatemala. J. Geophys. Res. 118(4), 1526–1533 (2013).

12. 12.

Waite, G. P., Chouet, B. A. & Dawson, P. B. Eruption dynamics at Mount St. Helens imaged from broadband seismic waveforms: Interaction of the shallow magmatic and hydrothermal systems. J. Geophys. Res. 113(B02305), 1–22 (2008).

13. 13.

Chouet, B. A., Dawson, P. B., James, M. R. & Lane, S. J. Seismic source mechanism of degassing bursts at Kilauea Volcano, Hawaii: Results from waveform inversion in the 10–50 s band. J. Geophys. Res. 115(B09311), 1–24 (2010).

14. 14.

Almendros, J., Chouet, B., Dawson, P. & Bond, T. Identifying elements of the plumbing system beneath Kilauea Volcano, Hawaii, from the source locations of very-long-period signals. Geophys. J. Int. 148(2), 303–312 (2002).

15. 15.

Nishimura, T. et al. Source process of very long period seismic events associated with the 1998 activity of Iwate Volcano, northeastern Japan. J. Geophys. Res. 105(B8), 19135–19147 (2000).

16. 16.

Kawakatsu, H. et al. Aso94: Aso seismic observation with broadband instruments. J. Volcanol. Geotherm. Res. 101(1–2), 129–154 (2000).

17. 17.

Arciniega-Ceballos, A., Chouet, B., Dawson, P. & Asch, G. Broadband seismic measurements of degassing activity associated with lava effusion at Popocatépetl Volcano, Mexico. J. Volcanol. Geotherm. Res. 170(1–2), 12–23 (2008).

18. 18.

Marchetti, E. & Ripepe, M. Stability of the seismic source during effusive and explosive activity at Stromboli Volcano. Geophys. Res. Lett. 32, 33-1-33–4 (2005).

19. 19.

Hill, D. P. et al. Very long period volcanic earthquakes beneath Mammoth Mountain, California. Geophys. Res. Lett. 29(10), 8-1-8–4 (2002).

20. 20.

Rosi, M. et al. Stromboli volcano, Aeolian Islands (Italy): Present eruptive activity and hazards. The Aeolian islands volcanoes. Geol. Soc. Lon. Mem. 37, 475–492 (2013).

21. 21.

Ripepe, M., Delle Donne, D., Harris, A. J. L., Marchetti, E. & Ulivieri, G. Dynamics of Stromboli activity: The Stromboli Volcano: An Integrated Study of the 2002–2003 Eruption. Am. Geophs. Geoph. Mon. 182, 39–48 (2008).

22. 22.

Neuberg, J., Luckett, R., Ripepe, M. & Braun, T. Highlights from a seismic broadband array on Stromboli volcano. Geophys. Res. Lett. 21(9), 749–752 (1994).

23. 23.

Kirchdörfer, M. Analysis and quasistatic FE modeling of long period impulsive events associated with explosions at Stromboli volcano (Italy). Ann. Geofis. 42(3), 379–390 (1999).

24. 24.

Chouet, B. et al. Broadband measurements of the sources of explosions at Stromboli Volcano, Italy. Geophys. Res. Lett. 26(13), 1937–1940 (1999).

25. 25.

James, M. R., Lane, S. J. & Chouet, B. Gas slug ascent through changes in conduit diameter: Laboratory insights into a volcano seismic source process in low viscosity magmas. J. Geophys. Res. 111(B05201), 1–25 (2006).

26. 26.

Maeda, Y., Kato, A. & Yamanaka, Y. Modeling the dynamics of a phreatic eruption based on a tilt observation: Barrier breakage leading to the 2014 eruption of Mount Ontake, Japan. J. Geophys. Res. 122, 1007–1024 (2017).

27. 27.

Harris, A. J. L. & Ripepe, M. Synergy of multiple geophysical approaches to unravel explosive eruption conduit and source dynamics: A case study from Stromboli. Chem. Erde 67, 1–35 (2007).

28. 28.

Patrick, M., Wilson, D., Fee, D., Orr, T. & Swanson, D. Shallow degassing events as a trigger for very-long-period seismicity at Kīlauea Volcano, Hawai ‘i. Bull. Volcanol. 73(9), 1179–1186 (2011).

29. 29.

Giudicepietro, F. et al. Changes in the VLP seismic source during the 2007 Stromboli eruption. J. Volcanol. Geotherm. Res. 182, 162–171 (2009).

30. 30.

Ripepe, M. et al. Volcano seismicity and ground deformation unveil the gravity-driven magma discharge dynamics of a volcanic eruption. Nat. Commun. 6, 6998. https://doi.org/10.1038/ncomms7998 (2015).

31. 31.

Valade, S. et al. Tracking dynamics of magma migration in open-conduit systems. Bull. Volcanol. 78, 1–12 (2016).

32. 32.

Ripepe, M. et al. Forecasting effusive dynamics and decompression rates by magmastatic model at open-vent volcanoes. Sci. Rep. 7(1), 1–9 (2017).

33. 33.

Patrick, M. R. Dynamics of Strombolian ash plumes from thermal video: Motion, morphology, and air entrainment. J. Geophys. Res. 112, B06202 (2007).

34. 34.

Marchetti, E., Ripepe, M., Harris, A. J. L. & Delle Donne, D. Tracing the differences between Vulcanian and Strombolian explosions using infrasonic and thermal radiation energy. Earth Planet. Sci. Lett. 279(3–4), 273–281 (2009).

35. 35.

Delle Donne, D. & Ripepe, M. High-frame rate thermal imagery of Strombolian explosions: Implications for explosive and infrasonic source dynamics. J. Geophys. Res. 117, 1–12 (2012).

36. 36.

Delle Donne, D. et al. Exploring the explosive effusive transition using permanent ultraviolet cameras. J. Geophys. Res. 122(6), 4377–4394 (2017).

37. 37.

Harris, A. J. & Stevenson, D. S. Thermal observations of degassing open conduits and fumaroles at Stromboli and Vulcano using remotely sensed data. J. Volcanol. Geotherm. Res. 76(3–4), 175–198 (1997).

38. 38.

Kanamori, H., Given, J. W. & Lay, T. Analysis of seismic body waves excited by the Mount St. Helens eruption of May 18, 1980. J. Geophys. Res. 89(B3), 1856–1866 (1984).

39. 39.

Brodsky, E. E., Kanamori, H. A. & Sturtevant, B. Seismically constrained mass discharge rate for the initiation of the May 18, 1980 Mount St. Helens Eruption. J. Geophys. Res. 104(B12), 29387–29400 (1999).

40. 40.

Lighthill, M. J. Waves in Fluids (Cambridge University Press, 1978).

41. 41.

Vergniolle, S. & Brandeis, G. Strombolian explosions: 1. A large bubble breaking at the surface of a lava column as a source of sound. J. Geophys. Res. 101(B9), 20433–20447 (1996).

42. 42.

Johnson, J. B. & Ripepe, M. Volcano infrasound: A review. J. Volcanol. Geotherm. Res. 206(3–4), 61–69 (2011).

43. 43.

Lacanna, G. & Ripepe, M. Modelling the acoustic flux inside the magmatic conduit by 3D-FDTD simulation. J. Geophys. Res. 125, E2019JB018849 (2020).

44. 44.

Burton, M., Allard, P., Mure, F. & La Spina, A. Magmatic gas composition reveals the source of slug driven Strombolian explosive activiy. Science 317, 227–230 (2007).

45. 45.

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(L24804), 1–5 (2006).

46. 46.

Delle, D. D. et al. Gas mass derived by infrasound and UV cameras: Implications for mass flow rate. J. Volcanol. Geotherm. Res. 325, 169–178 (2016).

47. 47.

Genco, R. & Ripepe, M. Inflation-deflation cycles revealed by tilt seismic records at Stromboli volcano. Geophys. Res. Lett. 37, L12302 (2010).

48. 48.

Wielandt, E. & Forbriger, T. Near-field seismic displacement and tilt associated with the explosive activity of Stromboli. Ann. Geophys. 42(3) (1999).

49. 49.

Legrand, D., Kaneshima, S. & Kawakatsu, H. Moment tensor analysis of near-field broadband waveforms observed at Aso Volcano, Japan. J. Volcanol. Geotherm. Res. 101(1–2), 155–169 (2000).

50. 50.

Legrand, D. et al. Interpretation of very long period tremors at Ambrym volcano, Vanuatu, as quasi static displacement field related to two distinct magmatic sources. Geophys. Res. Lett. https://doi.org/10.1029/2004GL021968 (2005).

51. 51.

Aoyama, H. & Oshima, H. Tilt change recorded by broadband seismometer prior to small phreatic explosion of Meakandake volcano, Hokkaido, Japan. Geophys. Res. Lett. 35, L06307 (2008).

52. 52.

Lyons, J. J., Waite, G. P., Ichihara, M. & Lees, J. M. Tilt prior to explosions and the effect of topography on ultra-long-period seismic records at Fuego volcano Guatemala. Geophys. Res. Lett. 39, L08305 (2012).

53. 53.

Bonaccorso, A. & Davis, P. M. Models of ground deformation from vertical volcanic conduits with application to eruptions of Mount St. Helens and Mount Etna. J. Geophys. Res. 104, 10531–11542 (1999).

54. 54.

Ripepe, M. et al. Ground deformation reveals the scale-invariant conduit dynamics driving explosive basaltic eruptions. Nat. Commun. 12, 1683 (2021).

55. 55.

Sugimura, S. et al. Moment tensor inversion of very-long-period seismicity observed in a very-near-field condition at Stromboli volcano. Geophys. Res. Abs. 21, 1 (2019).

56. 56.

Ripepe, M., Ciliberto, S. & Della Schiava, M. Time constraints for modeling source dynamics of volcanic explosions at Stromboli. J. Geophys. Res. 106(B5), 8713–8727 (2001).

57. 57.

Pistolesi, M., Delle Donne, D., Pioli, L., Rosi, M. & Ripepe, M. The 15 March 2007 explosive crisis at Stromboli volcano, Italy: Assessing physical parameters through a multidisciplinary approach. J. Geophys. Res. 116, B12 (2011).

58. 58.

Knox, H. A., Chaput, J. A., Aster, R. C. & Kyle, P. R. Multiyear shallow conduit changes observed with lava lake eruption seismograms at Erebus volcano, Antarctica. J. Geophys. Res. 123(4), 3178–3196 (2018).

59. 59.

Suckale, J., Keller, T., Cashman, K. & Persson, P. O. Flow to fracture transition in a volcanic mush plug may govern normal eruptions at Stromboli. Geophs. Res. Lett. 43, 071 (2016).

60. 60.

Nishimura, T. Ground deformation caused by magma ascent in an open conduit. J. Volcanol. Geotherm. Res. 187, 178–192 (2009).

61. 61.

Barth, A., Edmonds, M. & Woods, A. Valve-like dynamics of gas flow through a packed crystal mush and cyclic strombolian explosions. Sci. Rep. 9, 821 (2019).

62. 62.

Oppenheimer, J. et al. Analogue experiments on the rise of large bubbles through a solids-rich suspension: A “weak plug” model for Strombolian eruptions. Earth Planet. Sci. Lett. 531, 115931 (2020).

63. 63.

Pering, T. D. et al. Conduit dynamics and post explosion degassing on Stromboli: A combined UV camera and numerical modeling treatment. Geophys. Res. Lett. 43(10), 5009–5016 (2016).

64. 64.

Ripepe, M. et al. Effusive to explosive transition during the 2003 eruption of Stromboli volcano. Geology 33(5), 341–344 (2005).

65. 65.

Ripepe, M., Delle Donne, D., Lacanna, G., Marchetti, E. & Ulivieri, G. The onset of the 2007 Stromboli effusive eruption recorded by an integrated geophysical network. J. Volcanol. Geotherm. Res. 182(3–4), 131–136 (2009).

## Acknowledgements

Some results on the seismic source presented in this work find a strong support on the accurate seismological work done by Shunsuke Sugimura to whom this work is dedicated to. This study was realized by DEVNET grant of the Italian Civil Defence. D. L. has been partially supported by the PASPA-DGAPA, UNAM program, as a sabbatical year at University of Florence in 2016 and PAPIIT-DGAPA project number IN106019.

## Author information

Authors

### Contributions

All the Authors have contributed to the manuscript. In detail, M.R., D.D.D and G.L. wrote the manuscript and realized Figures. G.L. analyzed the acoustic data. D.L. contributed to define the quasi-static displacement method. S.V. carried out experiments with thermal camera. M.R., G.L. S.V. D.L. organized the 2016 seismic experiment. All Authors have reviewed the manuscript.

### Corresponding author

Correspondence to M. Ripepe.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Ripepe, M., Delle Donne, D., Legrand, D. et al. Magma pressure discharge induces very long period seismicity. Sci Rep 11, 20065 (2021). https://doi.org/10.1038/s41598-021-99513-4

• Accepted:

• Published: