The 2014 Effusive Eruption at Stromboli : New Insights from In Situ and Remote-Sensing Measurements

In situ and remote-sensing measurements have been used to characterize the run-up phase and the phenomena that occurred during the August–November 2014 flank eruption at Stromboli. Data comprise videos recorded by the visible and infrared camera network, ground displacement recorded by the permanent-sited Ku-band, Ground-Based Interferometric Synthetic Aperture Radar (GBInSAR) device, seismic signals (band 0.02–10 Hz), and high-resolution Digital Elevation Models (DEMs) reconstructed based on Light Detection and Ranging (LiDAR) data and tri-stereo PLEIADES-1 imagery. This work highlights the importance of considering data from in situ sensors and remote-sensing platforms in monitoring active volcanoes. Comparison of data from live-cams, tremor amplitude, localization of Very-Long-Period (VLP) source and amplitude of explosion quakes, and ground displacements recorded by GBInSAR in the crater terrace provide information about the eruptive activity, nowcasting the shift in eruptive style of explosive to effusive. At the same time, the landslide activity during the run-up and onset phases could be forecasted and tracked using the integration of data from the GBInSAR and the seismic landslide index. Finally, the use of airborne and space-borne DEMs permitted the detection of topographic changes induced by the eruptive activity, allowing for the estimation of a total volume of 3.07 ± 0.37 × 106 m3 of the 2014 lava flow field emplaced on the steep Sciara del Fuoco slope.


Introduction
Stromboli volcano (Italy; Figure 1), a stratovolcano located at the easternmost end of the Aeolian Archipelago, experienced a flank eruption from August-November 2014 [1][2][3][4][5][6][7].In this paper, in situ and remote-sensing measurements at Stromboli between May and November 2014 are presented.Data comprise videos recorded by the visible and infrared cameras network [8,9] (Figure 1a), ground displacement recorded by the permanent sited Ku-band, Ground-Based Interferometric Synthetic Aperture Radar (GBInSAR) device [10,11], ground displacements in the seismic band (0.02-10 Hz) [12,13] (Figure 1a), and high-resolution Digital Elevation Models (DEMs) reconstructed based on tri-stereo PLEIADES-1 imagery and Light Detection and Ranging (LiDAR) data, respectively.Data allowed us to characterize the precursors and the phenomena that occurred during the eruption, as well as an estimation of the erupted volume.
The last flank eruption started at Stromboli on 7 August 2014, preceded by 2 months of increased Strombolian activity and several lava overflows from the craters expanding along the SdF [2,26].Overflows were often accompanied by landslides along the SdF [7], described as rock-falls and/or gravel slides, evolving down slope to gravel flows (Figure 2a).The onset of the 2014 flank eruption (6-7 August) involved the breaching of the summit cone with emplacement of a landslide along the SdF (Figure 2b), the opening of an eruptive fissure on the NE flank of the cone [24], and the effusion of lava from the crater rim at first and from the eruptive fissure later, feeding the 2014 lava flow field [4,5,7].The eruption was characterized by the lava effusion in the SdF from a fissure at 650 m above sea level (a.s.l.) that lasted until 13 November 2014 [4].
The 2014 lava flow field was similar to others erupted from high-elevation vents (as described for the 2002-2003 lava flow field by [27]), comprising: (i) a series of tumuli and lava flows around the effusive vent at ∼650 m a.s.l.(proximal shield); (ii) a medial zone fed by small flows, and characterized by frequent lava crumbling down slope and producing a debris field; (iii) a basal toe composed of the debris flow field emplaced above the stacked lava delta [28].The erupted volume has been estimated to be 7.4 × 10 6 m 3 by [4] and ∼5.5 × 10 6 m 3 by [5].

The Camera Monitoring Network (Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Etneo -INGV-OE)
The camera monitoring network comprises thermal infrared and visible cameras located at Il Pizzo, at ~918 m elevation and ~250 m from the craters, plus visible and thermal infrared cameras (SQV400 and SQT400, respectively) at 400 m elevation on the N flank of the SdF and ~800 m from the craters (Figure 1a).The SQV400 and SQT400 cameras allow a view from NE of the North-East Crater NEC area (Figure 1a) and of the upper eastern sector of the SdF.An additional thermal camera installed in July 2014, is located at the lower eastern end of the SdF (~190 m a.s.l., Figure 1b) and is focused on the lower portion of the slope between ~400 m and the N coast line.To obtain a description of the eruptive activity, the total number of explosive events that occurred during each day of cloud-free observation was manually counted and plotted as an integer versus time.On average, between 24 July 2014 and 10 August 2014, ~30% of the days were affected by clouds and/or by system failure.In such cases data are lacking.10 broadband seismic stations.In this work, we used data from 3 seismic stations (Figure 1b).STR1 is located on the flank of the volcano at an elevation of about 560 m and has the longest and most complete record of seismological parameters.STRA is located at about 840 m, it is the closest to the summit craters, therefore is particularly suitable to study the dynamics of explosive activity.Lastly, STR8, located at about 570 m, is the closest to the SdF, hence it is the best suited to study seismic signals related to landslides in this area.The seismological parameters considered here are: the amplitude of volcanic tremor, the amplitude of explosion quakes and the inclination of the seismic polarization in the Very-Long-Period (VLP) band (0.05-0.5 Hz).We also performed the analysis of seismic signals, using neural networks, which allows detecting signals related to landslides occurring along the SdF [30].Landslides signals have shown to be a reliable precursor of effusive activity at Stromboli.Their occurrence rate increased markedly during the onset of the previous effusive eruption in 2007 [12].To quantify the intensity of landslide activity, Esposito, et al. [31] defined a landslide index following

Seismic Network (Istituto Nazionale di Geofisica e Vulcanologia, Osservatorio Vesuviano -INGV-OV)
During the 2014 eruption the seismic broadband network of Stromboli [29] was operating with 10 broadband seismic stations.In this work, we used data from 3 seismic stations (Figure 1b).STR1 is located on the flank of the volcano at an elevation of about 560 m and has the longest and most complete record of seismological parameters.STRA is located at about 840 m, it is the closest to the summit craters, therefore is particularly suitable to study the dynamics of explosive activity.Lastly, STR8, located at about 570 m, is the closest to the SdF, hence it is the best suited to study seismic signals related to landslides in this area.
The seismological parameters considered here are: the amplitude of volcanic tremor, the amplitude of explosion quakes and the inclination of the seismic polarization in the Very-Long-Period (VLP) band (0.05-0.5 Hz).We also performed the analysis of seismic signals, using neural networks, which allows detecting signals related to landslides occurring along the SdF [30].Landslides signals have shown to be a reliable precursor of effusive activity at Stromboli.Their occurrence rate increased markedly during the onset of the previous effusive eruption in 2007 [12].To quantify the intensity of landslide activity, Esposito, et al. [31] defined a landslide index following a fuzzy logic scheme.This parameter varies smoothly between 0 and 1, with 0 indicating a low probability of ongoing landslide activity and, conversely, 1 indicating a high probability of occurrence.

GBInSAR
InSAR technique has been successfully used to catch slope deformation at active volcanoes [31][32][33], though Synthetic Aperture Radar (SAR) satellites are no longer able to follow the evolution of very fast slope deformation.Conversely, given their repeat time (minutes), GBInSAR can correctly measure very fast displacements, thus the InSAR technique is ideal for monitoring, surveillance, and early-warning applications [24,34].At Stromboli, a GBInSAR apparatus is working since 2003, located in a stable area N of the SdF (Figure 1b), and consisting of a transmitting and a receiving antenna moving along a 3 m long rail [10].The GBInSAR, working in Ku-band, has also the great advantage of penetrating the dust clouds (abundant especially during the collapse events, see Figure 2), and working in every light and atmospheric condition.The GBInSAR measures ground displacement along the line of sight (LOS; Figure 1a,b), by computing via cross correlation the phase differences between the backscattered signals associated with two SAR images.Range and cross-range resolution are on average 2 × 2 m, with a measured displacement precision lower than 1 mm [10,11].Due to the short-elapsed time (11 min) between two subsequent measurements, the interferometric displacements are usually smaller than half the wavelength (8.6 mm for the Ku-band), therefore unwrapping procedures are generally not necessary [10].Moreover, unwrapping is a time-consuming process, avoiding the operational use of the GBInSAR during crisis management.For the Stromboli GBInSAR images, a coherence mask (threshold equal to 0.8) is set to reduce the noisy areas of the interferogram [11].Displacement maps and time-series are obtained "stacking" the interferogram phase with a displacement measurement precision of 0.5 mm, obtained using 1-h averaged SAR images [35].

Topographic Data and Co-Registration
Topographic change detection of Stromboli island due to the 2014 eruption was performed by comparing pre-and post-eruptive DEMs, the former generated from airborne LiDAR data and the post-eruption DEM from PLEIADES-1 satellite data.
The pre-2014 LiDAR-DEM was obtained elaborating the 3D point cloud acquired during an airborne survey carried on May 2012 by "BLOM Compagnia Generale Ripreseaeree S.P.A." (www.blomasa.com).The data were acquired using the Leica ADS80 sensor which has instrumental vertical and horizontal accuracy of ±10/20 cm and ±25 cm, respectively (see [36] for 2012 DEM features).The acquired point cloud has a mean point density of 8 pt/m 2 .
The post-2014 DEM was derived from the very high-resolution (VHR) tri-stereo optical imagery from the PLEIADES-1 satellites.PLEIADES-1 constellation is composed by two satellites, PLEIADES-1A (PHR1A) and PLEIADES-1B (PHR1B).These satellites can sense three or more synchronous images of the same area, with angles variable between ~6• e ~28 • .The stereoscopic triplet is composed of three nearly simultaneously acquired images, one backward looking, one forward looking, plus a third near-nadir image [37][38][39][40][41]. Tri-stereo images are 100% cloud free and were acquired by PHR1A on 27 May 2017, with a total areal coverage of 58 km 2 .This acquisition mode permits to obtain a DEM using the photogrammetric processing of three images (tri-stereo mode), with a pixel xy resolution of 1 m × 1 m.The dataset comprised also VHR optical imagery (0.5 m × 0.5 m resolution for Panchromatic + 2 m × 2 m Multispectral data).To assess the accuracy of the heights and their horizontal position in the Pleiades-1 DEM, six Ground-Control Points (GCPs) were collected on the map database (Cartographic XY standard deviation: 0.15 m).Nine tie points were automatically collected in the images.A block adjustment including all the satellite scenes was performed.The block adjustment was validated when the following accuracy was achieved: (i) pixel xy bias smaller than 0.3 pixels; (ii) pixel xy standard deviation smaller than 0.3 pixels; (iii) pixel xy maximum smaller than 2 pixels.Topographic change detection using multi-temporal DEMs was performed by differencing two DEMs of the same area derived from data taken at different time.This calculation is often affected by errors depending on mismatches between the two DEMs, which lead to artefact ∆z [42].This error can be detected and reduced by measuring and minimizing the DEM differences in areas where the two DEMs are supposed to be equal, i.e., those areas not reasonably affected by relevant natural changes.
The 2012 and 2017 topographic data were here co-registered by minimizing the root mean square (RMS) error between two DEMs.The methods used iteratively vary the three angles of rotation, the translation, and the magnification or reduction factor of one DEM using a custom-made algorithm based on the MINUIT minimization library [43], as described by [44][45][46].MINUIT is a tool to find the minimum value of multi-parameter functions and can be freely downloaded (http://www.cern.ch/minuit).The LiDAR-PLEIADES DEMs co-registration was performed by calculating the minimization parameters for the whole volcano with the exclusion of the investigated areas (i.e., the SdF).The RMS error in this case is given by the MINUIT output.The 2012 LiDAR was used as reference for co-registering the PLEIADES DEM.It is worth emphasizing that this RMS error is the root mean square residuals (in elevation) between a 3D point cloud and the reference DEM, rather than a true absolute error.Co-registration decreased the RMS error from the initial 3.37 m to the final 1.27 m.
The differences between the co-registered DEMs were used to detect and outline the extent of the areas affected by topographic changes in the considered time spam, and to calculate the volume and thickness variations inside them.Moreover, the perimeters of observed changes were double-checked using the ortho-rectified PLEIADES-1 images.Field surveys also strengthened data interpretation.
Volume (V) emplaced or lost between two acquisitions was calculated from the DEM difference according to the equation: (see [42] and references therein), where x is the grid step and z i is the height variation within the grid cell i.These values were then summed up for all the cells in the selected areas in which we want to calculate volume changes.
The residual volume errors were calculated as: where A is the investigated area and σ ∆z is the co-registration residual RMS error between the two DEMs.Following [42], this formula assumes that the errors were completely correlated, and thus the volume error calculated as above is overestimated.

Results
Explosive activity at Stromboli from the crater terrace was at an average of ~10 explosions h −1 during the first four months of 2014 (INGV Reports, available at www.ct.ingv.it).Between end of June-early July, the explosion number peaked at ~25 explosions h −1 (INGV Reports, available at www.ct.ingv.it),when also several lava flows poured out from the South-West crater (SWC, Figure 1a), and from two of the three vents located in the NEC area (NEC1 and NEC-hornito, Figure 1a, Table 1, Supplementary Materials), spreading along the SdF.Between 22 June and 5 August 2014, several powerful explosions and small landslides along the SdF occurred (Supplementary Materials).
The GBInSAR recorded an increase in the LOS displacement rate directed towards the sensor in the crater terrace area, indicating inflation of the summit plumbing system, starting from 30 May 2014 (Figure 3B,C).The following months (June-early August 2014) were characterized by slightly fluctuating displacement rate, with values always above the threshold that distinguishes the anomalous from the persistent activity (>0.05 mm/h directed towards the sensor; as defined by [11]).A small deflation of the crater terrace occurred around 24 July 2014, after the output of several small overflows, then a more rapid inflation of the crater terrace, and a week later also of the NEC-talus, was observed (Figure 3C).One of the most striking features in the seismological parameters during the 2014 eruption was the progressive, but rapid increase in both volcanic tremor and explosion-quake amplitude, starting approximately on 2 August 2014 (Figure 3D,E).Both quantities showed a steady increment, followed by a sudden drop in the amplitudes occurring about on 6 August 2014 at 12:00 (Figure 3D,E).This event was associated with the onset of the effusive activity on 6 August 2014, clearly evidenced by the landslide index, approaching the value of 1 in the same period (Figure 3A).However, a detailed analysis of the VLP polarization (Figure 3F), which is a good indicator of changes in the level of the magmatic column within the conduits, revealed that the most significant drop in the magma column level occurred only at about 04:00 of 7 August 2014, just before the effusive vent opened at 650 m a.s.l.(V1 and V2; see Table 1, Figure 4, and Supplementary Materials for further details).
The displacement rate suddenly increased on 5 August 2014 for both the external flank of the crater terrace and the NEC-talus, reaching their maximum on 6 August 2014 (Figure 3B,C).An increase in the explosion rate and intensity was recorded on 6 August 2014 morning, especially from the NEC.Several overflows occurred since 08:50 (flows F1 to F3 in Figure 4), and the rim of the crater terrace started to be clearly unstable since 11:00, when three arcuate fractures on the NE crater rim formed between NEC1 and NEC-hornito (Table 1, Supplementary Materials).
The GBInSAR recorded an increase in the displacement rate (since 10:00 of 6 August 2014) in both the external flank of the crater terrace and the NEC-talus.The maximum displacement was recorded mainly in the NEC-hornito area (Figure 5a).Rock-falls started to occur at 12:22, and soon after a lava flow erupted from the saddle between the NEC-hornito and the SWC (Flow F2A in Figure 4a), expanding N along the SdF (12:32).Between 13:00 and 14:21, the GBInSAR recorded a large loss in the SAR coherence in the NEC-hornito, consistent with the occurrence of fast-moving material (landslides and/or overflows) from the NEC-hornito (Figure 5b).Three landslides (starting at 14:05, 14:06, and 14:07, Figure 6), described as rock-falls/gravel slides evolving in rapid gravel flows, were directly observed (Figure 2b).The collapse scar observed on a photo taken from helicopter (Figure 2a) had an estimated size of ~40 m width, 30 m high and 20 m depth, with a total collapsed volume estimated in ~24,000 m 3 .Effusive vents moved northward from the NEC-hornito towards the NEC1.Afterwards, the displacement towards the sensor (up to 236 mm h −1 ) shifted from the NEC-hornito towards the NEC-talus (Figure 5c).Between 01:29 and 01:32 of 7 August 2014 (Figure 5c), the GBInSAR detected the formation of a fracture upslope of the NEC-talus.At the same time, the GBInSAR recorded an inversion in the LOS displacement rate (movements directed away from the sensor) in the crater terrace area, coherent with the deflation of the shallow conduit (Figure 5d).terrace started to be clearly unstable since 11:00, when three arcuate fractures on the NE crater rim formed between NEC1 and NEC-hornito (Table 1, Supplementary materials).
The GBInSAR recorded an increase in the displacement rate (since 10:00 of 6 August 2014) in both the external flank of the crater terrace and the NEC-talus.The maximum displacement was recorded mainly in the NEC-hornito area (Figure 5a).Rock-falls started to occur at 12:22, and soon after a lava flow erupted from the saddle between the NEC-hornito and the SWC (Flow F2A in Figure 4a), expanding N along the SdF (12:32).Between 13:00 and 14:21, the GBInSAR recorded a large loss in the SAR coherence in the NEC-hornito, consistent with the occurrence of fast-moving material (landslides and/or overflows) from the NEC-hornito (Figure 5b).Three landslides (starting at 14:05, 14:06, and 14:07, Figure 6), described as rock-falls/gravel slides evolving in rapid gravel flows, were directly observed (Figure 2b).The collapse scar observed on a photo taken from helicopter (Figure 2a) had an estimated size of ~ 40 m width, 30 m high and 20 m depth, with a total collapsed volume estimated in ~ 24,000 m 3 .Effusive vents moved northward from the NEC-hornito towards the NEC1.Afterwards, the displacement towards the sensor (up to 236 mm h −1 ) shifted from the NEC-hornito towards the NEC-talus (Figure 5c).Between 01:29 and 01:32 of 7 August 2014 (Figure 5c), the GBInSAR detected the formation of a fracture upslope of the NEC-talus.At the same time, the GBInSAR recorded an inversion in the LOS displacement rate (movements directed away from the sensor) in the crater terrace area, coherent with the deflation of the shallow conduit (Figure 5d).At 05:04 V1 opened at ~650 m a.s.l.(Figure 4) and at 05:08 the NE flank of NEC1 started collapsing, as evidenced by the associated high-frequency seismic signal (Figure 7).At 05:01, the live cam detected a downslope curved fracture opened on the flank of the NEC cone, followed at 05:16 by the opening of the effusive vent V2 at ~650 m a.s.l.(Table 1 and Figure 4).The onset of the increased landslide activity occurred at 05:08, and its amplitude reached a peak at ~05:24, i.e., during the opening of V1 and V2 (Table 1 and Figure 7).With the shift of the active vent from V1 to V2, a lava flow started spreading along the Pianoro (Supplementary Material) and the flank eruption started.
We analyzed in detail 24 hours of high-frequency seismic signals related to landslide activity along the SdF starting from 10:00 on 6 August 2014 (Figure 7B).Between 10:00 and 18:00 on 6 August 2014, the repeated lava overflows triggered a moderate landslide activity along SdF (Fig. 7A).Before the opening of vent V1, at 04:01 on 6 August, we observed a progressive increase in the landslides, which dropped just after this event.However, the clear onset of the increased landslide activity occurred at 05:08 and the amplitude reached a peak at around 05:24, a few minutes after the opening of vent V2 at 05:16.After a pause of about 8 minutes, the landslide activity resumed at 05:32 and At 05:04 V1 opened at ~650 m a.s.l.(Figure 4) and at 05:08 the NE flank of NEC1 started collapsing, as evidenced by the associated high-frequency seismic signal (Figure 7).At 05:01, the live cam detected a downslope curved fracture opened on the flank of the NEC cone, followed at 05:16 by the opening of the effusive vent V2 at ~650 m a.s.l.(Table 1 and Figure 4).The onset of the increased landslide activity occurred at 05:08, and its amplitude reached a peak at ~05:24, i.e., during the opening of V1 and V2 (Table 1 and Figure 7).With the shift of the active vent from V1 to V2, a lava flow started spreading along the Pianoro (Supplementary Materials) and the flank eruption started.
We analyzed in detail 24 h of high-frequency seismic signals related to landslide activity along the SdF starting from 10:00 on 6 August 2014 (Figure 7B).Between 10:00 and 18:00 on 6 August 2014, the repeated lava overflows triggered a moderate landslide activity along SdF (Figure 7A).Before the opening of vent V1, at 04:01 on 6 August, we observed a progressive increase in the landslides, which dropped just after this event.However, the clear onset of the increased landslide activity occurred at 05:08 and the amplitude reached a peak at around 05:24, a few minutes after the opening of vent V2 at 05:16.After a pause of about 8 min, the landslide activity resumed at 05:32 and lasted for about 1 h, before reaching a stationary level (slightly higher than the value observed before 05:08) and possibly related to the ongoing lava flow (Figure 7C).
The GBInSAR recorded a large loss in coherence at the NEC-talus and the SdF area, consistent with the onset of the effusion (Figure 5D).At 05:29 the flow had travelled the Pianoro flat surface and at 05:30 started to expand along the SdF steep slope (Supplementary Materials).Here it formed at least 5 lava flow branches (Supplementary Materials), reaching the coast with a hot avalanche 30 min later (Table 1, Supplementary Materials).The landslide activity resumed at 05:32, and the collapse lasted for about 1 h, coherent with several landslides' pulses (Figure 8).By 06:24 several lava branches reached the sea, forming a lava fan ~500 m wide (Table 1, Supplementary Materials).Effusive and explosive activity declined between 7 and 8 August 2014, with only one lava flow active along the coast and deeper magma level within the feeder conduit, evidenced by more collimated jets produced during the explosive activity at the summit vents (Table 1, Supplementary Materials).
Remote Sens. 2018, 10, x FOR PEER REVIEW 11 of 21 lasted for about 1 hour, before reaching a stationary level (slightly higher than the value observed before 05:08) and possibly related to the ongoing lava flow (Figure 7C).The GBInSAR recorded a large loss in coherence at the NEC-talus and the SdF area, consistent with the onset of the effusion (Figure 5D).At 05:29 the flow had travelled the Pianoro flat surface and at 05:30 started to expand along the SdF steep slope (Supplementary Material).Here it formed at least 5 lava flow branches (Supplementary Material), reaching the coast with a hot avalanche 30 minutes later (Table 1, Supplementary materials).The landslide activity resumed at 05:32, and the collapse lasted for about 1 hour, coherent with several landslides' pulses (Figure 8).By 06:24 several lava branches reached the sea, forming a lava fan ~500 m wide (Table 1, Supplementary materials).Effusive and explosive activity declined between 7 and 8 August 2014, with only one lava flow active along the coast and deeper magma level within the feeder conduit, evidenced by more collimated jets produced during the explosive activity at the summit vents (Table 1, Supplementary materials).Loss in SAR coherence persisted until 8 August 2014 on the whole scene illuminated by the GBInSAR device.In the crater terrace, the coherence increased again since 9 August 2014 when no explosive activity was observed from the monitoring cameras, and the GBInSAR recorded displacement away from the sensor until late September indicating deflation of the summit.
The images recorded by the thermal camera located at Il Pizzo (SPT; see Figure 1B for location) showed significant morphology changes of the summit craters since 9 August 2014, showing that the apex of cones and hornitos within NEC1, NEC-hornito and SWC had been removed by the several collapses and hot "gravel flows" of the previous days.The thermal images displayed also a much lower temperature of the whole crater area, indicating a decline in the explosive activity.Likewise, the rate of effusion decreased significantly from the previous day, with only one wider lava flow starting from the eruptive fissure and replacing the 4-5 branches of the previous day.This flow was reaching the coast forming a narrow channel apparently cooling down, with overflows sometimes spilling over the channel levees.Loss in SAR coherence persisted until 8 August 2014 on the whole scene illuminated by the GBInSAR device.In the crater terrace, the coherence increased again since 9 August 2014 when no explosive activity was observed from the monitoring cameras, and the GBInSAR recorded displacement away from the sensor until late September indicating deflation of the summit.
The images recorded by the thermal camera located at Il Pizzo (SPT; see Figure 1B for location) showed significant morphology changes of the summit craters since 9 August 2014, showing that the apex of cones and hornitos within NEC1, NEC-hornito and SWC had been removed by the several collapses and hot "gravel flows" of the previous days.The thermal images displayed also a much lower temperature of the whole crater area, indicating a decline in the explosive activity.Likewise, the rate of effusion decreased significantly from the previous day, with only one wider lava flow starting from the eruptive fissure and replacing the 4-5 branches of the previous day.This flow was reaching the coast forming a narrow channel apparently cooling down, with overflows sometimes spilling over the channel levees.Until 15 October 2014 the GBInSAR recorded displacement towards the sensor, indicating inflation of the crater area.Afterwards, very low displacement rate (<0.05 mm/h) oscillating between away and towards the sensor, have been recorded.In the SdF area, after the onset of the 2014 flank eruption, the radar recorded a very low coherence for three days; this was related to the initial fast-moving lava flow.Since 10 August 2014, low coherence zones were related only to a small part of the monitored scenario, corresponding to active lava tongues, while the rest of the monitored SdF showed variable displacement rates (1-10 mm/h).
The seismic signals related to minor landslides along SdF continued to be recorded after the development of the lava flow on SdF.

Remote Sens. 2018, 10, x FOR PEER REVIEW 13 of 21
Until 15 October 2014 the GBInSAR recorded displacement towards the sensor, indicating inflation of the crater area.Afterwards, very low displacement rate (<0.05 mm/h) oscillating between away and towards the sensor, have been recorded.In the SdF area, after the onset of the 2014 flank eruption, the radar recorded a very low coherence for three days; this was related to the initial fastmoving lava flow.Since 10 August 2014, low coherence zones were related only to a small part of the monitored scenario, corresponding to active lava tongues, while the rest of the monitored SdF showed variable displacement rates (1-10 mm/h).
The seismic signals related to minor landslides along SdF continued to be recorded after the development of the lava flow on SdF.In the period between the acquisitions of the two DEMs (May 2012-May 2017), large topographic changes occurred.The areas interested in volume gain and loss are shown in Figure 9, which includes the crater terrace and the SdF depression.The largest variations were related to the emplacement of the 2014 lava flow field in the NE part of the SdF (Figure 9).The total emplaced volume is 3.07 ± 0.37 × 10 6 m 3 , calculated inside the black dashed line.The average thickness of lava was 10.45 ± 1.27 m, with the maximum value reaching 35 m (Figure 9, profile 1).In the period between the acquisitions of the two DEMs (May 2012-May 2017), large topographic changes occurred.The areas interested in volume gain and loss are shown in Figure 9, which includes the crater terrace and the SdF depression.The largest variations were related to the emplacement of the 2014 lava flow field in the NE part of the SdF (Figure 9).The total emplaced volume is In the central part of the SdF, between 0 and ~518 m a.s.l., an increased volume of 1.70 ± 0.30 × 10 6 m 3 has been estimated, with an average thickness of 7.11 ± 1.27 m.This volume is related to the accumulation of volcaniclastic material derived from the erosion of the SdF, from the overflowinduced landslides, and from the NEC-talus collapse [7].

Discussion
The run-up phase of the 2014 flank eruption was characterized by very vigorous explosive activity and episodic lava overflows from the crater terrace.This activity was preceded by an increase in the displacement rate towards the sensor of the crater terrace at anomalous level since the end of May 2014, and by frequent landslides along the SdF.Ground deformation of the crater terrace strongly increased on 26 July 2014 and was followed by increase of the volcanic tremor and explosionquake amplitude, starting approximately on 2 August 2014.The explosive activity peaked on 5-6 August 2014, anticipating by ~1 day the increase in the displacement rate towards the sensor (inflation of the crater terrace).
Different authors [1,3,6] suggested that the 2014 eruption was triggered by a gradual recharge of the shallow plumbing system and upper conduits.The SO2 flux in the persistent Stromboli plume has been measured for two months before the onset of the eruption, detecting a strong increase above normal activity in the SO2 pulses (puffing and explosions) before the effusion onset.This is consistent with the increase in gas bubble supply and magma transport rate feeding the uppermost storage system at Stromboli [47].Ground displacement is consistent with the presence of a very shallow reservoir below the crater terrace [2].Long-term (month to weeks) increase in the ground deformation recorded by the GBInSAR in the crater terrace area, as well as the long-term rise of the VLP seismic source towards the surface [5,6], are associated with the accumulation of magma in the uppermost storage system at Stromboli [2].Contrarily, the detected increase in short-term (days to hours) ground deformation and amplitude of volcanic tremor are associated with the impulsive bulk In the central part of the SdF, between 0 and ~518 m a.s.l., an increased volume of 1.70 ± 0.30 × 10 6 m 3 has been estimated, with an average thickness of 7.11 ± 1.27 m.This volume is related to the accumulation of volcaniclastic material derived from the erosion of the SdF, from the overflow-induced landslides, and from the NEC-talus collapse [7].

Discussion
The run-up phase of the 2014 flank eruption was characterized by very vigorous explosive activity and episodic lava overflows from the crater terrace.This activity was preceded by an increase in the displacement rate towards the sensor of the crater terrace at anomalous level since the end of May 2014, and by frequent landslides along the SdF.Ground deformation of the crater terrace strongly increased on 26 July 2014 and was followed by increase of the volcanic tremor and explosion-quake amplitude, starting approximately on 2 August 2014.The explosive activity peaked on 5-6 August 2014, anticipating by ~1 day the increase in the displacement rate towards the sensor (inflation of the crater terrace).
Different authors [1,3,6] suggested that the 2014 eruption was triggered by a gradual recharge of the shallow plumbing system and upper conduits.The SO 2 flux in the persistent Stromboli plume has been measured for two months before the onset of the eruption, detecting a strong increase above normal activity in the SO 2 pulses (puffing and explosions) before the effusion onset.This is consistent with the increase in gas bubble supply and magma transport rate feeding the uppermost storage system at Stromboli [47].Ground displacement is consistent with the presence of a very shallow reservoir below the crater terrace [2].Long-term (month to weeks) increase in the ground deformation recorded by the GBInSAR in the crater terrace area, as well as the long-term rise of the VLP seismic source towards the surface [5,6], are associated with the accumulation of magma in the uppermost storage system at Stromboli [2].Contrarily, the detected increase in short-term (days to hours) ground deformation and amplitude of volcanic tremor are associated with the impulsive bulk degassing from the shallow plumbing system [35,48], whereas the increase in the explosion number and intensity and amplitude of explosion quakes are indicative of the increase in the release rate and overpressure intensity of individual explosions.
One of the most interesting aspects of the data presented here is the relationship between the instability events and the effusive activity at Stromboli volcano (Figure 10).Other authors have already shown the critical state of Stromboli during the very early stages of effusive eruptions [2,5,11,22,49].During the run-up to the 2014 flank eruption, gravel flows consisting of mixed breccia and loose volcaniclastic deposits were generated mainly in two ways: (i) by the crumbling of the overflows, or (ii) by the collapse of some portions of the external part of the crater terrace and apex of the hornitos (Figure 2).The second type of collapse is a recurring phenomenon at Stromboli, occurring when there is a large accumulation of magma below the crater terrace (causing overflows, increase in the frequency of explosions, spattering activity, inflation of the whole summit).Furthermore, weakness of the crater wall (e.g., induced by the strong explosive activity on 5 August 2014), and mechanical erosion caused by vent opening and by lava fingering, may have triggered the NEC-hornito collapse on 6 August 2014 [48].Here we have documented for the first time that, after the NEC-hornito collapse, eruptive vents migrated from south to north (Figure 10a).Overflows and mass-wasting in the NEC-hornito area preceded the formation of fracture below the NEC1 (Figure 10b) by ~11 h, and the onset of the flank eruption and NEC1 collapse by ~15 h, suggesting that 11-15 h before the onset of the effusive eruption, the crater terrace was full of magma.Then, the location of effusive vents displayed a clockwise rotation along the rim of the crater terrace, from the NEC-hornito towards NEC2 and NEC1 (Figure 2), and finally magma migrated towards the NEC-talus, inducing fracture opening and landslides occurrence.This process suggested a gradual pressurization and/or expansion of the feeder dike along a SW-NE direction, which is also the main tectonic trend of the island [49,50].Lava started to outflow from a first vent (V1; Figure 10c) before the NEC-talus collapse and shifted at vent V2 simultaneously with the beginning of the NEC-talus landslide (Figure 10d).There could be different explanations for this vent shift: 1.
a structural barrier generated by the presence of a buried structure [51,52] or different stratified material (Figure 2b), which in turn has caused the magma-filled fracture to re-orientate [53,54]; 2. magma to flow from V1 at V2 in response to the tensile stress occurring in the talus produced by the downwards displacement observed by GBInSAR.
It is important to emphasize that this shift has prevented the dike from intruding within the SdF, as happened in the early stages of past eruptions (e.g., the 2007 eruption) [8].Instead during the first days of the 2002-03 flank eruption, the intrusion of magma propagating laterally from the conduit towards a lower altitude was considered the cause that destabilized the SdF slope, with the triggering of tsunamigenic landslides on 30 December 2002 [22][23][24][25].
The flank eruption onset produced the drainage of a superficial magma reservoir [5,49], inducing ground deflation [55], a deepening of the VLP seismic source [5], and the cessation of summit explosive activity [5].
The volume of the drained magma can have critical implications on the eruption style and on its transition from effusive to explosive [56].A previous author [4] calculated a total bulk volume of 7.4 × 10 6 m 3 using data derived by the new satellite Technology Experiment Carrier-1 (TET-1) by the German Aerospace Center (DLR), whereas others [5] evaluated a total bulk volume of 5.5 × 10 6 m 3 by thermal images from satellites using the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on-board the Terra and Aqua satellites.This last value was also confirmed by [57].In MODIS approach, the error is usually estimated at 50% [58].The flank eruption onset produced the drainage of a superficial magma reservoir [5,49], inducing ground deflation [55], a deepening of the VLP seismic source [5], and the cessation of summit explosive activity [5].
The volume of the drained magma can have critical implications on the eruption style and on its transition from effusive to explosive [56].A previous author [4] calculated a total bulk volume of 7.4 × 10 6 m 3 using data derived by the new satellite Technology Experiment Carrier-1 (TET-1) by the German Aerospace Center (DLR), whereas others [5] evaluated a total bulk volume of 5.5 × 10 6 m 3 by thermal images from satellites using the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor on-board the Terra and Aqua satellites.This last value was also confirmed by [57].In MODIS approach, the error is usually estimated at 50% [58].
To improve the volume estimation, we calculated the topographic change by comparing the 2012 and 2017 DEMs.This approach considers the volume variations due to erosion and accumulation occurred between the two acquisitions, therefore also the topographic changes occurred before or after the eruption and the volume of the lava flow emplaced below sea level during the first days of the eruption.However, in this area there was no evidence of significant erosive phenomena [7], and moreover the evidence of large accumulation (i.e., more than 1-2 × 10 6 m 3 ) of lava or a debris flow under the sea is still missing [28].
We estimated that during the 2014 flank eruption, a volume of 3.07 ± 0.37 × 10 6 m 3 of lava was emplaced on the steep SdF slope, i.e., most of the volume emplaced in the subaerial part of the volcano.This value is close to the previous estimations of 5.5 × 10 6 m 3 [5,57], which however comprises also the minor portion emplaced below sea level.To improve the volume estimation, we calculated the topographic change by comparing the 2012 and 2017 DEMs.This approach considers the volume variations due to erosion and accumulation occurred between the two acquisitions, therefore also the topographic changes occurred before or after the eruption and the volume of the lava flow emplaced below sea level during the first days of the eruption.However, in this area there was no evidence of significant erosive phenomena [7], and moreover the evidence of large accumulation (i.e., more than 1-2 × 10 6 m 3 ) of lava or a debris flow under the sea is still missing [28].
We estimated that during the 2014 flank eruption, a volume of 3.07 ± 0.37 × 10 6 m 3 of lava was emplaced on the steep SdF slope, i.e., most of the volume emplaced in the subaerial part of the volcano.This value is close to the previous estimations of 5.5 × 10 6 m 3 [5,57], which however comprises also the minor portion emplaced below sea level.
Some authors [57] proposed that the critical factor for triggering the most violent explosive events (paroxysms, [59][60][61]) is the decompression rate, suggesting that rates >10 Pa s −1 can potentially favor the ascent of gas-rich magma batches responsible for the paroxysmal explosions that occurred at Stromboli during the 2002-03 and 2007 flank eruptions.
A critical aspect is that the paroxysmal explosions of both 5 April 2003 and 15 March 2007 occurred while effusive activity was still going on, and when the lava effusion rate (related to the magma decompression rate) was decreasing with respect to the initial phases of the eruptions.Moreover, the 2002-03 and 2014 flank eruptions were characterized by the same mean effusion rate [57] but, differently from the 2002-03, the 2014 eruption did not produce any paroxysm during the effusive activity.
A different approach to understand the transition from effusive to explosive style can be obtained by estimating the cumulated effusive volume.In fact, a previous author [56] proposed that the erupted volume during flank eruptions at Stromboli is very important for hazard estimation, because it is considered a key factor for the triggering of paroxysmal explosive activity.Calvari et al. [56] proposed that there is an erupted magma bulk volume limit of 6.5 ± 1 × 10 6 m 3 .Once that this volume is poured out, for example by drainage due to the opening of a flank fissure, the decompression induced on the shallow storage can trigger paroxysmal activity, by bringing at shallow levels the gas-rich magma that is in the lower conduit [62][63][64].Calvari et al. [56] point out that the rate of effusion, and consequently the decompression rate, is less crucial to trigger powerful explosive activity.This interesting aspect has also been observed trough laboratory experiments allowing direct observations of explosive expansion of a slowly decompressed magma analog [65].The analog experiments highlighted that under these conditions an explosive eruption would only occur when sufficient magma was spilled from the conduit, in agreement with the considerations of [56].The inferred volume of 2014 emitted lava (3.07 ± 0.37 × 10 6 m 3 ), even adding 1-2 × 10 6 m 3 of underwater volume not measurable with our DEMs approach, is below the limit of 6.5 ± 1 × 10 6 m 3 suggested by [56] and, therefore, consistent with the absence of the paroxysmal explosion.This consideration also explains the different behavior of the 2014 eruption (with absence of paroxysmal activity), respect to the 2002-03 eruption (with paroxysmal explosion).In fact, although the two eruptions had the same mean effusion rate, the 2014 emitted volume did not reach the critical limit (6.5 × 10 6 m 3 ) that was instead overcome by the 2002-03 eruption.

Conclusions
Integration of in situ and remote-sensing measurements can strongly support the management of eruptive crises at active volcanoes.In this paper, data from the in situ (visible and infrared) live-cam and seismic monitoring network, measuring the ground motion in the seismic band (0.02-10 Hz), have been integrated with ground displacement recorded by a GBInSAR device.Furthermore, to constrain the erupted lava volume, DEMs from LiDAR and PLEIADES-1 tri-stereo were compared, also supported by the high-resolution PLEIADES-1 optical images, supporting the areal mapping of the 2014 lava flow field.The main findings are: • live-cam and explosion-quake data revealed that the explosive activity peaked between 5 and 6 August 2014, whereas the GBInSAR device recorded a drastic increase in the displacement rate since the morning of 6 August, consistent with a strong inflation of the crater terrace; • ground displacement started to show evidence of sliding in the crater terrace after the 6 August 2014 evening, as also evidenced by seismic signals; • the onset of the 2014 flank eruption involved the breaching of the summit cone with emplacement of a landslide along the SdF (anticipated by the GBInSAR measurements, observed by the live cam, and recorded by the seismic data), the opening of an eruptive fissure on the NE flank of the cone (observed by the live cam), and the effusion of lava from the crater rim at first and from the eruptive fissure later, feeding the 2014 lava flow field; • the eruption was characterized by the lava effusion along the SdF from a fissure at 650 m above sea level (a.s.l.) that lasted until 13 November 2014, with a total volume of 3.07 ± 0.37 × 10 6 m 3 of lava emplaced on the steep SdF slope; • this volume is below the limit of 6.5 ± 1 × 10 6 m 3 expected for triggering a paroxysmal explosion.
This work highlights the importance of considering data from in situ sensors and remote-sensing platforms in monitoring active volcanoes.Comparison of data from live-cams, seismic monitoring, and ground displacements recorded by GBInSAR devices provide information about the eruptive activity, nowcasting the shift in eruptive style from effusive to explosive.At the same time, the landslide activity during the run-up and onset phases could be forecasted and tracked using the integration of

Figure 1 .
Figure 1.(a) Geographic setting of the Island of Stromboli, located at the easternmost end of the Aeolian Archipelago, and names of the summit vents within the crater terrace.NEC: North-East Crater area; CC: Central Crater area; SWC: South-West Crater area (b) 3D-view of the investigated sectors observed from north, highlighting the location of the measurement stations used in this work.STR1, STR8, and STRA: broadband seismic stations.SQT400 and SQV400: thermal and visible live cam, respectively.In both (a) and (b), the line of sight (LOS) displacement map produced by the GBInSAR apparatus between 30 May 2014 and 6 August 2014 is also shown, highlighting the strong inflation of the crater terrace and the stability of the Sciara del Fuoco before the 2014 flank eruption.The measure-areas of the time-series are indicated by white circles.All Digital data were collected in the Projected coordinate system: WGS 1984 UTM zone 33 Projection: Transverse Mercatore.Maps were generated using ESRI ArcGIS CAMPUS (Università degli Studi di Firenze Licence; http://www.siaf.unifi.it/vp-1275-arcgis-licenza-campus.html).

Figure 1 .
Figure 1.(a) Geographic setting of the Island of Stromboli, located at the easternmost end of the Aeolian Archipelago, and names of the summit vents within the crater terrace.NEC: North-East Crater area; CC: Central Crater area; SWC: South-West Crater area (b) 3D-view of the investigated sectors observed from north, highlighting the location of the measurement stations used in this work.STR1, STR8, and STRA: broadband seismic stations.SQT400 and SQV400: thermal and visible live cam, respectively.In both (a) and (b), the line of sight (LOS) displacement map produced by the GBInSAR apparatus between 30 May 2014 and 6 August 2014 is also shown, highlighting the strong inflation of the crater terrace and the stability of the Sciara del Fuoco before the 2014 flank eruption.The measure-areas of the time-series are indicated by white circles.All Digital data were collected in the Projected coordinate system: WGS 1984 UTM zone 33 Projection: Transverse Mercatore.Maps were generated using ESRI ArcGIS CAMPUS (Università degli Studi di Firenze Licence; http://www.siaf.unifi.it/vp-1275-arcgis-licenza-campus.html).

Figure 2 .
Figure 2. (a) Gravel flow associated with the 7 July 2014 overflow.The landslide is composed by mixed breccia from the overflow front and the (mainly fine-grained) debris eroded along the Sciara del Fuoco (photo taken by the Università degli Studi di Firenze personnel during field survey).(b) Landslides associated with the collapse of the NEC-hornito, occurred on 6 August 2014 (photo taken by UNIFI personnel during boat survey).

Figure 2 .
Figure 2. (a) Gravel flow associated with the 7 July 2014 overflow.The landslide is composed by mixed breccia from the overflow front and the (mainly fine-grained) debris eroded along the Sciara del Fuoco (photo taken by the Università degli Studi di Firenze personnel during field survey).(b) Landslides associated with the collapse of the NEC-hornito, occurred on 6 August 2014 (photo taken by UNIFI personnel during boat survey).
-13:00 Incandescent blocks accumulation along the coast 06 Aug 2014 at 13:08 Overflow from NEC-hornito reached the coast 06 Aug 2014 14:05-14:08 Three hot "gravel flows" from the NEC-hornito along the SdF, reached the coast and went on spreading along the sea surface for several tens of meters 06 Aug 2014 at 14:50 Overflow from NEC-hornito reached the coast 06 Aug 2014 at 15:46 Overflow from NEC-hornito reached the coast 06 Aug 2014 at 16:02 Hot "gravel flows" from the NEC-hornito along the SdF 06 Aug 2014 at 16:08 Overflow from NEC-hornito reached the coast 07 Aug 2014 ~02:30 Hot "gravel flows" from the NEC-hornito along the SdF 07 Aug 2014 ~03:00 -NEC-hornito lava decreased -Explosive activity from NEC increased -Increased size of the hot avalanche deposit 07 Aug 2014 at 03:40 Hot "gravel flows" from the NEC onto the Pianoro flat area 07 Aug 2014 04:01 NE flank of NEC1 started to collapse Lava flow from the NEC towards the Pianoro flat area Hot "gravel flows" from the NEC-talus onto the Pianoro flat area and along the SdF 07 Aug 2014 05:01 Downslope curved fracture opened on the flank of the cone 07 Aug 2014 05:04-05:16 V1 opened at 05:04 and V2 at 05:16, both at ~650 m a.s.l.Landslides onto the Pianoro flat area and along the SdF Lava flow onto the Pianoro flat area 07 Aug 2014 05:30 Lava flow along the SdF 07 Aug 2014 ~06:02 Hot "gravel flows" from the NEC-talus reached the coast 07 Aug 2014 06:24 Lava flow from V1 reached the coast Remote Sens. 2018, 10, x FOR PEER REVIEW 8 of 21

Figure 3 .
Figure 3. Relevant geophysical time-series from May to August 2014.(A) Landslide index derived from seismic data; GBInSAR cumulative displacement data for (B) the crater terrace and (C) for the NEC-talus.(D) Volcanic tremor amplitude at station STR1.(E) Explosion-quake amplitudes at station STR1.(F) VLP polarization inclination at station STRA.

Figure 3 .
Figure 3. Relevant geophysical time-series from May to August 2014.(A) Landslide index derived from seismic data; GBInSAR cumulative displacement data for (B) the crater terrace and (C) for the NEC-talus.(D) Volcanic tremor amplitude at station STR1.(E) Explosion-quake amplitudes at station STR1.(F) VLP polarization inclination at station STRA.

Figure 4 .
Figure 4. (a) The external flank of the crater terrace after the onset of the 2014 flank eruption (photo taken by INGV personnel on 11 August 2014 from helicopter survey).Lava overflows emitted during the run-up phase (F1 to F3) and lava flow produced during and after the onset (F4A and F4B) are also reported.(b) Detail of the NEC-talus area highlighting the scar produced during the onset of the eruption (7 August 2014), the first lava flow after crater breaching (F4A, erupted from vent V1), the upper portion of the 2014 lava flow field (F4B), and the vent of the 2014 flank eruption (V2) (photo taken by UNIFI personnel on 8 October 2014 from field survey).It is interesting to note the internal stratification in the NEC-talus remnants (consisting of volcaniclastic levels) that are continuous in all the inner walls of the landslide scar.

Figure 4 .
Figure 4. (a) The external flank of the crater terrace after the onset of the 2014 flank eruption (photo taken by INGV personnel on 11 August 2014 from helicopter survey).Lava overflows emitted during the run-up phase (F1 to F3) and lava flow produced during and after the onset (F4A and F4B) are also reported.(b) Detail of the NEC-talus area highlighting the scar produced during the onset of the eruption (7 August 2014), the first lava flow after crater breaching (F4A, erupted from vent V1), the upper portion of the 2014 lava flow field (F4B), and the vent of the 2014 flank eruption (V2) (photo taken by UNIFI personnel on 8 October 2014 from field survey).It is interesting to note the internal stratification in the NEC-talus remnants (consisting of volcaniclastic levels) that are continuous in all the inner walls of the landslide scar.

Figure 5 .
Figure 5. GBInSAR interferograms formed using a pair of images recorded between (a) the 12:09 and 12:20 on 6 August 2014, evidencing the instability of the NEC-hornito; (b) the 13:59 and 14:10 on 6 August 2014, evidencing the instability of the NEC-talus; (c) the 01:29 and 01:40 on 7 August 2014, evidencing the concomitant instability of the NEC-talus and the deflation of the crater terrace, coherent with the sliding of the NEC-talus induced by the NE-ward magma migration from the crater terrace; (d) the 05:26 and 05:37 on 7 August 2014, characterized by loss in SAR coherence due to fast movement both in the crater area (sliding processes) and along the SdF (lava flow emplacement).Interferograms in (a), (b), and (c) are affected by phase wrapping.

Figure 5 .
Figure 5. GBInSAR interferograms formed using a pair of images recorded between (a) the 12:09 and 12:20 on 6 August 2014, evidencing the instability of the NEC-hornito; (b) the 13:59 and 14:10 on 6 August 2014, evidencing the instability of the NEC-talus; (c) the 01:29 and 01:40 on 7 August 2014, evidencing the concomitant instability of the NEC-talus and the deflation of the crater terrace, coherent with the sliding of the NEC-talus induced by the NE-ward magma migration from the crater terrace; (d) the 05:26 and 05:37 on 7 August 2014, characterized by loss in SAR coherence due to fast movement both in the crater area (sliding processes) and along the SdF (lava flow emplacement).Interferograms in (a-c) are affected by phase wrapping.

Figure 6 .
Figure 6.Thermal images recorded by the SCT camera located at Labronzo, ~100 m a.s.l.(see Figure 1b for location), showing the hot "gravel flows" that reached the sea on 6 August 2014, at 14:05 (A-D), 14:06 (E-H) and 14:07 (I-L), spreading on the sea surface.Temperature is in °C.The velocities of the three landslides of hot blocks in the frontal region were 9.8 m s −1 , 5.9 m s −1 and 9.8 m s −1 , respectively, based on web-cameras measurements.

Figure 6 .
Figure 6.Thermal images recorded by the SCT camera located at Labronzo, ~100 m a.s.l.(see Figure 1b for location), showing the hot "gravel flows" that reached the sea on 6 August 2014, at 14:05 (A-D), 14:06 (E-H) and 14:07 (I-L), spreading on the sea surface.Temperature is in • C. The velocities of the three landslides of hot blocks in the frontal region were 9.8 m s −1 , 5.9 m s −1 and 9.8 m s −1 , respectively, based on web-cameras measurements.

Figure 7 .
Figure 7. Plot of the high-frequency (>10 Hz) seismic amplitudes averaged over 10 s windows (station STR8, EW seismic component) for: (A) details referred to the time interval 10:00-18:00 of 6 August 2014, highlighting the signals associated with landslides that occurred soon after the onset of the overflows, whose onset is marked by orange vertical lines; (B) the time interval 10:00 on 6 August 2014-10:00 on 7 August 2014; (C) details referred to the time interval 03:00-10:00 of 7 August 2014, highlighting the signals associated with landslides that occurred at the opening V1 and V2 vents (red vertical lines).

Figure 7 .
Figure 7. Plot of the high-frequency (>10 Hz) seismic amplitudes averaged over 10 s windows (station STR8, EW seismic component) for: details referred to the time interval 10:00-18:00 of 6 August 2014, highlighting the signals associated with landslides that occurred soon after the onset of the overflows, whose onset is marked by orange vertical lines; (B) the time interval 10:00 on 6 August 2014-10:00 on 7 August 2014; (C) details referred to the time interval 03:00-10:00 of 7 August 2014, highlighting the signals associated with landslides that occurred at the opening V1 and V2 vents (red vertical lines).

Figure 8 .
Figure 8. Lava flows (white) and hot avalanches (yellow, green, and purple dotted outline) spreading along the lower slope of Sciara del Fuoco and to the sea, as observed by the thermal images recorded by the SCT camera located at Labronzo, ~100 m a.s.l. on the east side of the Sciara on 7 August 2014.See Figure 1b for camera location.(A-C): lava flow spreading to the coast, flanked by a hot avalanche (yellow dotted outline) at 5:59; (D-F): another hot avalanche spreading at 6:00 (yellow dotted line); (G-I): a much bigger hot avalanche spreading at 6:01 (red dotted outline) forming a conspicuous dust cloud).

Figure 8 .
Figure 8. Lava flows (white) and hot avalanches (yellow, green, and purple dotted outline) spreading along the lower slope of Sciara del Fuoco and to the sea, as observed by the thermal images recorded by the SCT camera located at Labronzo, ~100 m a.s.l. on the east side of the Sciara on 7 August 2014.See Figure 1b for camera location.(A-C): lava flow spreading to the coast, flanked by a hot avalanche (yellow dotted outline) at 5:59; (D-F): another hot avalanche spreading at 6:00 (yellow dotted line); (G-I): a much bigger hot avalanche spreading at 6:01 (red dotted outline) forming a conspicuous dust cloud).

Figure 9 .
Figure 9. Stromboli topographic changes detected comparing 2012 and 2017 DEMs.Dashed black line encloses the area covered by the 2014 lava flow.Dashed yellow lines, labelled with P1, P2 and P3, show the position of the profiles along the flow represented in the insets.

Figure 9 .
Figure 9. Stromboli topographic changes detected comparing 2012 and 2017 DEMs.Dashed black line encloses the area covered by the 2014 lava flow.Dashed yellow lines, labelled with P1, P2 and P3, show the position of the profiles along the flow represented in the insets.

Figure 10 .
Figure 10.Schematic sketches representing the onset of the 2014 flank eruption as reconstructed on the basis of the new data presented in this paper.(a) The run-up phase of the 2014 flank eruption was characterized by high magma level beneath the crater terrace, which peaked on 6 August 2014 morning, with the occurrence of a large overflow from the NEC-hornito, triggering its collapse.Then the activity migrated from the NEC-hornito towards the NEC1.(b) The onset of the eruption coincided with the lateral dike propagation from the crater terrace along a NE-SW direction, triggering the instability of the NEC-talus, as evidenced by the increase in displacement rate at NECtalus (until the loss of coherence in the GBInSAR data, compare with Figure 5) and landslides number (Figure 7).(c) Then the shift in the effusive vent from V1 to V2 occurred, (d) closely related to the NEC-talus collapse.Magmatic and mass-wasting phenomena are represented in red and in black, respectively.

Figure 10 .
Figure 10.Schematic sketches representing the onset of the 2014 flank eruption as reconstructed on the basis of the new data presented in this paper.(a) The run-up phase of the 2014 flank eruption was characterized by high magma level beneath the crater terrace, which peaked on 6 August 2014 morning, with the occurrence of a large overflow from the NEC-hornito, triggering its collapse.Then the activity migrated from the NEC-hornito towards the NEC1.(b) The onset of the eruption coincided with the lateral dike propagation from the crater terrace along a NE-SW direction, triggering the instability of the NEC-talus, as evidenced by the increase in displacement rate at NEC-talus (until the loss of coherence in the GBInSAR data, compare with Figure 5) and landslides number (Figure 7).(c) Then the shift in the effusive vent from V1 to V2 occurred, (d) closely related to the NEC-talus collapse.Magmatic and mass-wasting phenomena are represented in red and in black, respectively.

Table 1 .
Summary of the eruptive activity during the run-up and onset phases of the 2014 flank eruption, as derived from the camera monitoring network.