Abstract
Ground-Based Interferometric Synthetic Aperture Radar (GBInSAR) is an efficient technique for capturing short, subtle episodes of conduit pressurization in open vent volcanoes like Stromboli (Italy), because it can detect very shallow magma storage, which is difficult to identify using other methods. This technique allows the user to choose the optimal radar location for measuring the most significant deformation signal, provides an exceptional geometrical resolution and allows for continuous monitoring of the deformation. Here, we present and model ground displacements collected at Stromboli by GBInSAR from January 2010 to August 2014. During this period, the volcano experienced several episodes of intense volcanic activity, culminated in the effusive flank eruption of August 2014. Modelling of the deformation allowed us to estimate a source depth of 482â±â46âm a.s.l. The cumulative volume change was 4.7â±â2.6âÃâ105âm3. The strain energy of the source was evaluated 3â5 times higher than the surface energy needed to open the 6â7 August eruptive fissure. The analysis proposed here can help forecast shifts in the eruptive style and especially the onset of flank eruptions at Stromboli and at similar volcanic systems (e.g. Etna, Piton de La Fournaise, Kilauea).
Similar content being viewed by others
Introduction
Ground-Based Interferometric Synthetic Aperture Radar (GBInSAR) system offers an opportunity to image deformation due to shallow dike intrusions, shallow magma movement, opening/closing of eruptive fissures/ephemeral vents as well as landslides1,2,3,4,5,6 Since GBInSAR allows the user to choose the optimal radar location in terms of distance and incidence and azimuth angles, the deformation can be measured with exceptional geometrical resolution, while the very high sampling rate allows for continuous monitoring of the deformation1,6,7. Ground deformation at Stromboli, Italy, due to conduit processes may be difficult to detect and possibly related to strong explosions8,9,10 or syn-effusive deflations11. Both the 2002â2003 and 2007 flank eruptions completely drained the shallow storage system, as testified by the funnel-like vertical failure of the summit crater terrace associated with this eruption11,12,13,14,15,16,17,18 and to the lowering of the VLP source location19. The complete absence of strombolian activity is considered further evidence of the drainage of the shallower part of the Stromboli storage system19. According to ref. 11, the first two hours of the 2007 flank eruption drained 2.8â4.8âÃâ105âm3 of magma. Evidence of new deep magmatic input, together with the rebuilding of the summit terrace, was observed beginning in late 200918,19. During the 2010â2014 period, Stromboli experienced several episodes of intense volcanic activity2,5,6,19. Towards the end of May 2014, activity at Stromboli increased as the magma level inside the conduit rose. The volcano had several small effusions between June and August 2014 and peaked on August 7, 2014 with the opening of a fissure and an effusive eruption at â650âm a.s.l., â100âm below the NE the crater20.
GBInSAR Data
The NE portion of the summit area of Stromboli has been continuously monitored since January 20036 by a GBInSAR system located on a stable section of the flank, about 1.5âkm away from the crater terrace (Fig. 1). GBInSAR Line Of Sight (LOS) is mostly sensitive to the N-S horizontal component of displacement (average azimuth angleâ=â15°). Negative and positive values of displacement indicate, respectively, a movement toward and away from the sensor. Given the location of the system (Fig. 1), flank instability could be a possible interpretation for movements towards the sensor but not for movements away from the sensor. Since we have observed both movements away and forward the sensor and the modelling indicates that the location and depth of the source are stable, we believe that the changes in the LOS correspond to either inflation (negative values â downslope movement toward the sensor) or deflation (positive values â upslope movement away from the sensor) of the summit area1,2,3,4,5,6.
The GBInSAR system consists of a transmitting and a receiving antenna moving along a rail (3âm long in the configuration deployed at Stromboli)1. GBInSAR measures ground displacement along the LOS by computing, via cross correlation, the phase differences between the backscattered signals associated with two consecutive synthetic aperture radar (SAR) images. The ability of GBInSAR to measure volcano deformation depends on the persistence of phase coherence over time. The loss in coherence is primarily due to ground movements, e.g., lava flows or rock avalanches1,2. A coherence threshold equal to or above 0.8 is required to recognize deformation areas from a GBInSAR interferogram2. Due to the short time (11âmin) between two subsequent measurements, interferometric displacements are usually smaller than half wavelength and phase unwrapping procedures21 are not necessary. Both the range and cross-range resolutions are on average 2âmâÃâ2âm, with a precision in displacement measurements of less than 1âmm3. Displacement rates are computed by differencing the displacements obtained from two consecutive images and dividing by the time spanned. Displacement time series (Fig. 2) are acquired using an algorithm to sum, pixel by pixel, the displacements for every consecutive pair of images and then average that rate over an 8-hour interval3,22. Displacement time series of selected points (averaged over 10 pixels) are obtained from cumulative displacement maps with a precision in the displacement measurement of 0.5âmm. Daily displacement can be calculated from the time series.
Cumulative displacement maps (Fig. 1) indicate that the NE flank of Stromboli moved upslope-downslope with the major deformation localized in the upper part of the crater terrace (Fig. 1). GBInSAR displacement time series (1 January 2010â7 August 2014) reveal fluctuations in the deformation of the summit area, with three distinct periods of inflation (JulyâDecember 2011; September 2012âMay 2013; May-August 2014; orange stripes in the time series in Fig. 2), corresponding to periods of high-intensity eruptive activity, characterized by frequent and strong explosions and overflows6. Time series in the areas characterized by the largest displacement (Fig. 2b) reveal that the inflation of the three periods had different rates of deformation: 0.8âmm/day in JulyâDecember 2011, 1.2âmm/day in September 2012âMay 2013 and 2.1âmm/day in May 2014âAugust 2014. The peaks in the daily displacement rate occurred on 22 August 2011, 28 March 2013 and 2 August 2014 (Fig. 2b,c). The discrete Fourier transform of the GBInSAR time series of the daily displacement of the external rim of the crater terrace (Fig. 2b), shows that the most energetic peaks correspond to a period of 256 days (Supplementary Information). We observed that the GBInSAR displacement time series of areas affected by continuous debris deposition/erosion below the NE crater (NEC) have the same trend of the upper part of the crater terrace (Fig. 2c). Cross-correlation analysis between the two daily-averaged time series reveals that the displacement fluctuation occurred simultaneously (highest absolute-value correlation 0.69 at lag time tâ=â0). This may mean that the either the cone of the NEC is more affected by the deformation of the conduit, regardless of the erosional/depositional processes or that the swell induced by magma in the conduit is also accompanied by depositional processes on the NEC flank.
Modelling the Deformation
To determine the main parameters of the deformation source, we inverted the GBInSAR displacement employing the software dMODELS23. A number of source geometries (spherical source24, prolate spheroid25, horizontal penny-shaped source26 and tensile dislocation27), all in a flat, elastic, homogeneous, isotropic half-space are available in dMODELS23 for several geodetic techniques: leveling, tilt, GPS and InSAR. In the case of InSAR measurements, the software models the changes in range along the radar LOS28. Actual volcanic sources are not embedded cavities of simple shape but we assume that these models may reproduce the strain field created by actual storage areas and transport pathways. Given the location of the GBInSARsystem, the angle of incidence of the radar (LOS direction almost perpendicular to the volcano flank) and the fact that the system monitors only one flank of the volcano (Fig. 1), we can neglect (in first approximation) the volcano topography and model the deformation as a LOS displacement over a flat half space. The dMODELS software employs a nonlinear inversion algorithm to determine the best-fit parameters for the deformation source by searching the minimum of following the cost function23:
where, N is the number of data points, P the number of model parameters, dk are the experimental data, mk the modelling results and Ïk the data uncertainties. The non-linear inversion algorithm is a combination of local optimization (interior-point method29) and random search. This approach is more efficient for hyper-parameter optimization than trials on a grid30.
We use the empirical variogram, a measure of spatial correlation31 to determine which one the proposed source geometries best fit the deformation (see also Table 1). When two sources would fit the data with a similar precision (e.g, the deformation episode from 08 August 2012 to 08 March 2013 can be explained by either inverse faulting or deflation of a spherical source), we choose the source with the least number of parameters32.
Examples of the inversion of InSAR measurements to determine the location of the deformation sources are shown in Fig. 3. We inverted the run-up phases for each period of high-intensity eruptive activity, choosing cumulative maps in the existing dataset (two maps per day; see Table 1). To minimize the influence from sources not related to magma accumulation, we discarded maps affected by atmospheric disturbance. The models reveal a substantial stability of the deformation source over the considered time interval. The best fitting source geometry for the 130 January 2010âAugust 2014 inflation is a sphere 149â±â41âm beneath the volcano flank (Table 1; Fig. 4). The source depths can be transformed from a depth h relative (normal) to the volcano flank to a depth dasl relative to the sea level (a.s.l.) by the simple geometric transformation
where Hasl is the height of the summit crater terrace above the sea level and L is the distance between the source surface location and the summit crater terrace. Equation (2) allowed us to estimate a source depth daslâ=â482â±â46âm a.s.l. (Fig. 5).
Discussion and Conclusive Remarks
Modelling of the deformation confirmed the presence of a very shallow reservoir consistent with the persistency of magma within Stromboliâs conduits, whose existence has been proposed before from the analysis of geochemical data33,34, that broadly corresponds with the source location of syn-explosive deformation (350â600âm a.s.l.)35. The presence of a very shallow reservoir has been suggested by seismic and deformation data19 and by petrological studies of lithic ejecta and in particular from the evidence of pyrometamorphism in tephra accumulated within the crater terrace during persistent activity36. At Stromboli, zones of magma accumulation at different depth have been identified. The comparison of the source location recognized by this study with geophysical8,9,10,19,38, gas chemistry39 and petrological33,34,39,40 data argues that Stromboliâs magma plumbing configuration is a multiple-zone storage system (Fig. 6) composed by: i) a deep storage area that feed the most primitive magmas towards the surface (residence time>55 years)40; ii) an intermediate storage, mainly activated during energetic explosions and flank effusions (residence timeâ=â2â10 years)40; iii) a shallow storage that is involved in all the surface and near-surface phenomena, including explosive activity, central and flank effusions and non-eruptive dike injection (residence timeâ=â10â213 days)33,34. The comparable period of the GBInSAR displacement cycles and residence time in the shallow storage system33,34 suggests that the ground displacement in the crater terrace area is controlled by the accumulation of magma in the shallow storage system.
The modelling of GBInSAR data allows us to estimate the volume of the shallow magma accumulation and to evaluate its energy budget. The January 2010 â August 2014 unrest was characterized by the accumulation of 2.76âÃâ 104âm3 of magma. This value can be corrected to take into account both the effect of the magma compressibility (βm) and host rock (βc) stiffness41,42. Compressibility of basaltic magma falls in the range 0.4â2âÃâ10â10âPaâ1 ref. 41, while the host rock stiffness has been evaluated following the relationships:
where μ is the shear modulus, E the Youngâs modulus and Ï the Poissonâs ratio. A range of values for E and Ï can be obtained from ref. 43, considering both intact and damaged basalts (βc in the range 0.5â1Ãâ10â10âPaâ1). The corrected cumulative magma volume during January 2010âAugust 2014 is 4.7â±â2.6âÃâ105âm3 (Table 1), in agreement with previous estimates (7â±â2 105 m3) based on geochemical data33,34. These values are also of the same order of magnitude of the volume of magma drained in the starting phase of the 2002-03 and 2007 flank eruptions11,17,44, suggesting that the ephemeral vents first depleted the very shallow source. The accumulation rate fluctuated over time (Table 1), with an average value of 4.4âÃâ10â3âm3 sâ1 and a standard deviation of 2.5âÃâ10â3âm3 sâ1. A maximum accumulation rate of 1.2â±â0.7âÃâ 10â2âm3 sâ1 was estimated for the 28 May 2014-06 August 2014 inflation episode (Table 1). Based on the works of ref. 45 and ref. 46, at Stromboli 1â2âÃâ10â2âkm3 yearâ1 of magma is degassed, of which only â7âÃâ10â5âkm3 yearâ1 is erupted47. Our calculation would account for 1.5â±â0.8âÃâ10â4âkm3 yearâ1 of magma stored in the upper portion of the volcano. This leads to the conclusion that the larger portion of the magma fed at the surface falls back in the conduit due to convection39,48,49. Using the model proposed in ref. 50, we evaluated that the very shallow source recognized in this work to remain molten and open at a very shallow depth, need heat to its margins at a rate of <1âMW, much lower than heat release at the crater terrace (423â±â226âMW) evaluated by ref. 46.
The 2014 Stromboli flank eruption was characterized by the development of an eruptive fissure that moved the eruptive centre from the summit crater terrace at 750âm a.s.l. to an ephemeral vent at 650 m a.s.l. The initiation and propagation of a fracture depend primarily on a first order on the potential energy stored in a volcanic edifice when it is loaded51,52. Here the loading is primarily related to inflation of the shallow magma storage area. We calculate the strain energy of the shallow intrusions using the strain-nucleus Un52
where Vre is the volume of magma in a single deformation episode, Vc is the total magma volume, assumed here as the volume accumulated during the previous time intervals, a is the magma body radius and γ is the melt fraction (on average 50%)39,53. The strain energy in each analyzed interval is presented in Table 1. The total strain energy stored during the period January 2010âAugust 2014 was Unâ=â3.3â±â1.8âÃâ1014âJ (Table 1). The surface energy needed to open the 6â7 August eruptive fissure, considering a 170âm long (strike dimension, planar distance between the NE crater and the ephemeral vent) and 100âm tall (dip dimension, high difference between the NE crater and the ephemeral vent) fracture and an energy density of about 1.3â43âÃâ106âJâmâ2 (see ref. 51), is in the range of 109â1011âJ, much lower than the total energy stored at Stromboli during the January 2010âAugust 2014 inflation. The maximum rate and strain energy accumulated at Stromboli coincide with the inflation preceding the flank eruption of August 2014 (Fig. 7).
Ground displacement in the crater terrace area is controlled by the accumulation of magma in the shallow storage system. Fracture opening and propagation is controlled by an increase in the magma accumulation rate that allows building up significant potential energy in a short amount of time. Variations in GBInSAR time series reflect variation in strain energy stored in the shallow source and could be used to forecast the shift between summit and flank eruptions.
Additional Information
How to cite this article: Di Traglia, F. et al. Shifts in the eruptive styles at Stromboli in 2010â2014 revealed by ground-based InSAR data. Sci. Rep. 5, 13569; doi: 10.1038/srep13569 (2015).
References
Casagli, N. et al. Deformation of Stromboli Volcano (Italy) during the 2007 eruption revealed by radar interferometry, numerical modelling and structural geological field data. J. Volcanol. Geoth. Res. 182, 182â200 (2009).
Traglia, F. et al. Groundâbased InSAR reveals conduit pressurization pulses at Stromboli volcano. Terra Nova, 25, 192â198 (2013).
Nolesini, T., Di Traglia, F., Del Ventisette, C., Moretti, S. & Casagli, N. Deformations and slope instability on Stromboli volcano: Integration of GBInSAR data and analog modeling. Geomorphology, 180, 242â254 (2013).
Intrieri, E. et al. Flank instability of Stromboli volcano (Aeolian Islands, Southern Italy): integration of GB-InSAR and geomorphological observations. Geomorphology. 201, 60â69 (2013).
Di Traglia, F. et al. The ground-based InSAR monitoring system at Stromboli volcano: linking changes in displacement rate and intensity of persistent volcanic activity. B. Volcanol. 76, 1â18 (2014).
Di Traglia, F. et al. Review of ten years of volcano deformations recorded by the ground-based InSAR monitoring system at Stromboli volcano: a tool to mitigate volcano flank dynamics and intense volcanic activity. Earth-Sci. Rev. 139, 317â335 (2014).
Wadge, G. et al. Lava dome growth and mass wasting measured by a time series of groundâbased radar and seismicity observations. J. Geophys. Res.-Sol. Ea. (1978â2012). 113, 10.1029/2007JB005466 (2008).
Bonaccorso, A. et al. Stromboli 2007 eruption: Deflation modeling to infer shallowâintermediate plumbing system. Geophys. Res. Lett. 35, 10.1029/2007GL032921 (2008).
Mattia, M., Rossi, M., Guglielmino, F., Aloisi, M. & Bock, Y. The shallow plumbing system of Stromboli Island as imaged from 1 Hz instantaneous GPS positions. Geophys. Res. Lett. 31, 10.1029/2004GL021281 (2004).
Bonaccorso, A., Calvari, S., Linde, A., Sacks, S. & Boschi, E. Dynamics of the shallow plumbing system investigated from borehole strainmeters and cameras during the 15 March, 2007 Vulcanian paroxysm at Stromboli volcano. Earth Planet. Sci. Let. 357, 249â256 (2012).
Marchetti, E., Genco, R. & Ripepe, M. Ground deformation and seismicity related to the propagation and drainage of the dyke feeding system during the 2007 effusive eruption at Stromboli volcano (Italy). J. Volcanol. Geoth. Res. 182, 155â161 (2009).
Falsaperla, S., Neri, M., Pecora, E. & Spampinato, S. Multidisciplinary study of flank instability phenomena at Stromboli volcano, Italy. Geophys. Res. Lett. 33, 10.1029/2006GL025940 (2006).
Falsaperla, S., Maiolino, V., Spampinato, S., Jaquet, O. & Neri, M. Sliding episodes during the 2002â2003 Stromboli lava effusion: insights from seismic, volcanic and statistical data analysis. Geochem. Geophys. Geosyst. 9, Q04022, 10.1029/2007GC001859 (2008).
Acocella, V., Neri, M. & Scarlato, P. Understanding shallow magma emplacement at volcanoes: orthogonal feeder dikes during the 2002â2003 Stromboli (Italy) eruption. Geophys. Res. Lett. 33, 10.1029/2006GL026862 (2006).
Neri, M., Lanzafame, G. & Acocella, V. Dyke emplacement and related hazard in volcanoes with sector collapse: the 2007 Stromboli (Italy) eruption. J. Geol. Soc. 165, 883â886 (2008).
Neri, M. & Lanzafame, G. Structural features of the 2007 Stromboli eruption. J. Volcanol. Geoth. Res. 182, 137â144 (2009).
Calvari, S. et al. The 2007 Stromboli eruption: event chronology and effusion rates using thermal infrared data. J. Geophys. Res.-Sol. Ea. (1978â2012), 115, 10.1029/2009JB006478 (2010).
Calvari, S. et al. Major eruptive style changes induced by structural modifications of a shallow conduit system: the 2007â2012 Stromboli case. B. Volcanol. 76, 1â15 (2014).
Ripepe, M. et al. Volcano seismicity and ground deformation unveil the gravity-driven magma discharge dynamics of a volcanic eruption. Nat. Commun. 6, 10.1038/ncomms7998 (2015).
Rizzo, A. L. et al. The 2014 effusive eruption at Stromboli volcano (Italy): Inferences from soil CO2 flux and 3He/4He ratio in thermal waters. Geophys. Res. Lett. 42, 2235â2243 (2015).
Ghiglia, D. C. & Romero, L. A. Robust two-dimensional weighted and unweighted phase unwrapping that uses fast transforms and iterative methods. JOSA A. 11, 107â117 (1994).
Di Traglia, F., Cauchie, L., Casagli, N. & Saccorotti, G. Decrypting geophysical signals at Stromboli Volcano (Italy): Integration of seismic and GroundâBased InSAR displacement data. Geophys. Res. Lett., 41, 2753â2761 (2014).
Battaglia, M., Cervelli, P. F & Murray, J. R. dMODELS: A MATLAB software package for modeling crustal deformation near active faults and volcanic centers. J. Volcanol. Geoth. Res. 254, 1â4 (2013).
McTigue, D. F. Elastic stress and deformation near a finite spherical magma body: resolution of the point source paradox. J. Geophys. Res.-Sol. Ea. (1978â2012). 92, 12931â12940 (1987).
Yang, X. M., Davis, P. M. & Dieterich, J. H. Deformation from inflation of a dipping finite prolate spheroid in an elastic halfâspace as a model for volcanic stressing. J. Geophys. Res.-Sol. Ea. (1978â2012). 93, 4249â4257 (1988).
Fialko, Y., Khazan, Y. & Simons, M. Deformation due to a pressurized horizontal circular crack in an elastic half-space, with applications to volcano geodesy. Geophys. J. Int. 146, 181â190 (2001).
Okada, Y. Surface deformation due to shear and tensile faults in a half-space. B. Seismol. Soc. Am. 75, 1135â1154 (1985).
Tizzani, P. et al. Magma and fluid migration at Yellowstone Caldera in the last three decades inferred from InSAR, leveling and gravity measurements. J. Geophys. Res.-Sol. Ea. 120, 2627â2647 (2015).
Bonnans, J. F., Gilbert, J. C., Lemaréchal, C. & Sagastizábal, C. A. Numerical optimization: theoretical and practical aspects. Berlin, Springer-Verlag. (2006).
Bergstra, J. & Bengio, Y. Random search for hyper-parameter optimization. J. Mach. Learn. Res. 13, 281â305 (2012).
Deutsch, C. V. & Journel, A. G. Geostatistical Software Library and Userâs Guide. Oxford University Press, New York (1998).
Gordon, R. G., Stein, S., DeMets, C. & Argus, D. F. Statistical tests for closure of plate motion circuits. Geophys. Res. Lett., 14, 587â590 (1987).
Gauthier, P. J. & Condomines, M. 210Pbâ226Ra radioactive disequilibria in recent lavas and radon degassing: inferences on the magma chamber dynamics at Stromboli and Merapi volcanoes. Earth Planet. Sci. Let. 172, 111â126 (1999).
Gauthier, P. J., Le Cloarec, M. F. & Condomines, M. Degassing processes at Stromboli volcano inferred from short-lived disequilibria (210Pb-210Bi-210Po) in volcanic gases. J. Volcanol. Geoth. Res. 102, 1â19 (2000).
Chouet, B., Dawson, P. & Martini, M. Shallow-conduit dynamics at Stromboli Volcano, Italy, imaged from waveform inversions. Geol. Soc. London Spec. Publ. 307, 57â84 (2008).
Del Moro, S., Renzulli, A., Landi, P., La Felice, S. & Rosi, M. Unusual lapilli tuff ejecta erupted at Stromboli during the 15 March 2007 explosion shed light on the nature and thermal state of rocks forming the crater system of the volcano. J. Volcanol. Geoth. Res. 254, 37â52 (2013).
Apuani, T., Corazzato, C., Cancelli, A. & Tibaldi, A. Physical and mechanical properties of rock masses at Stromboli: a dataset for volcano instability evaluation. Bull. Eng. Geol. Env. 64, 419â431 (2005).
Bonaccorso, A. Evidence of a dykeâsheet intrusion at Stromboli Volcano Inferred through continuous tilt. Geophys. Res. Lett. 25, 4225â4228 (1998).
Aiuppa, A. et al. A model of degassing for Stromboli volcano. Earth Planet. Sci. Let. 295, 195â204 (2010).
Bragagni, A., Avanzinelli, R., Freymuth, H. & Francalanci, L. Recycling of crystal mush-derived melts and short magma residence times revealed by U-series disequilibria at Stromboli volcano. Earth Planet. Sci. Let. 404, 206â219 (2014).
Gudmundsson, A. Lateral magma flow, caldera collapse and a mechanism of large eruptions in Iceland. J. Volcanol. Geoth. Res. 34, 65â78 (1987).
Rivalta, E. & Segall, P. Magma compressibility and the missing source for some dike intrusions. Geophys. Res. Lett., 35, L04306, 10.1029/2007GL032521 (2008).
Heap, M. J., Vinciguerra, S. & Meredith, P. G. The evolution of elastic moduli with increasing crack damage during cyclic stressing of a basalt from Mt. Etna volcano. Tectonophysics, 471, 153â160 (2009).
Ripepe, M. et al. Effusive to explosive transition during the 2003 eruption of Stromboli volcano. Geology. 33, 341â344 (2005).
Allard, P., Carbonnelle, J., Metrich, N., Loyer, H. & Zettwoog, P. Sulphur output and magma degassing budget of Stromboli volcano. Nature. 368, 326â330 (1994)
Harris, A. J., & Stevenson, D. S. Magma budgets and steadyâstate activity of Vulcano and Stromboli. Geophys. Res. Lett., 24, 1043â1046 (1997).
Harris, A. & Ripepe, M. Synergy of multiple geophysical approaches to unravel explosive eruption conduit and source dynamicsâA case study from Stromboli. Chem. Erde-Geochem. 67, 1â35 (2007).
Burton, M. R., Mader, H. M. & Polacci, M. The role of gas percolation in quiescent degassing of persistently active basaltic volcanoes. Earth Planet. Sci. Let. 264, 46â60 (2007).
Witham, F. Conduit convection, magma mixing and melt inclusion trends at persistently degassing volcanoes. Earth Planet. Sci. Let., 301, 345â352 (2011).
Giberti, G., Jaupart, C. & Sartoris, G. Steady-state operation of Stromboli volcano, Italy: constraints on the feeding system. B. Volcanol. 54, 535â541 (1992).
Gudmundsson, A. Toughness and failure of volcanic edifices. Tectonophysics. 471, 27â35 (2009).
Gudmundsson, A. Strengths and strain energies of volcanic edifices: implications for eruptions, collapse calderas and landslides. Nat. Hazard Earth Sys. 12, 2241â2258 (2012).
Lautze, N. C. & Houghton, B. F. Linking variable explosion style and magma textures during 2002 at Stromboli volcano, Italy. B. Volcanol. 69, 445â460 (2007).
Acknowledgements
This work has been financially supported by the âPresidenza del Consiglio dei Ministri â Dipartimento della Protezione Civileâ (Presidency of the Council of Ministers - Department of Civil Protection) within the framework of the SAR.net, SAR.net2, InGrID and InGrID2015 projects; this publication, however, does not reflect the position and the official policies of the Department. Data and sensor used are designed and produced by Ellegi LLC using proprietary GBInSAR technology by LiSALab LLC, a European Commission Joint Research Centre spin-off.
Author information
Authors and Affiliations
Contributions
All the authors contributed to the manuscript. In detail: F.D.T: conception and realization of the manuscript, GBInSAR data analysis and modelling, field validation, GBInSAR maintenance. M.B: realization of the manuscript, implementation of modelling software (dMODELS), GBInSAR data modelling. T.N: realization of the manuscript, GBInSAR data analysis, field validation, GBInSAR maintenance. D.L: realization of the manuscript, GBInSAR data analysis. N.C: realization of the manuscript, financial support, responsible of the SAR.net, SAR.net2, InGrID and InGrID2015 projects.
Ethics declarations
Competing interests
The authors declare no competing financial interests.
Electronic supplementary material
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Di Traglia, F., Battaglia, M., Nolesini, T. et al. Shifts in the eruptive styles at Stromboli in 2010â2014 revealed by ground-based InSAR data. Sci Rep 5, 13569 (2015). https://doi.org/10.1038/srep13569
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/srep13569