Rapid provision of maps and volcanological parameters: quantification of the 2021 Etna volcano lava flows through the integration of multiple remote sensing techniques

At active volcanoes recurring eruptive events, erosive processes and collapses modify the edifice morphology and impact monitoring and hazard mitigation. At Etna volcano (Italy) between February and October 2021, 57 paroxysmal events occurred from the South-East Crater (SEC), which is currently its most active summit crater. Strombolian activity and high lava fountains (up to 4 km) fed lava flows towards the east, south and south-west, and caused fallout of ballistics (greater than 1 m in diameter) within 1–2 km from the SEC. The impacted area does not include permanent infrastructure, but it is visited by thousands of tourists. Hence, we rapidly mapped each lava flow before deposits became covered by the next event, for hazard mitigation. The high frequency of the SEC paroxysms necessitated integration of data from three remote sensing platforms with different spatial resolutions. Satellite (Sentinel-2 MultiSpectral Instrument, PlanetScope, Skysat and Landsat-8 Operational Land Imager) and drone images (visible and thermal) were processed and integrated to extract digital surface models and orthomosaics. Thermal images acquired by a permanent network of cameras of the Istituto Nazionale di Geofisica e Vulcanologia were orthorectified using the latest available digital surface model. This multi-sensor analysis allowed compilation of a geodatabase reporting the main geometrical parameters for each lava flow. A posteriori analysis allowed quantification of bulk volumes for the lava flows and the SEC changes and of the dense rock equivalent volume of erupted magma. The analysis of drone-derived digital surface models enabled assessment of the ballistics’ distribution. The developed methodology enabled rapidly and accurate characterisation of frequently occurring effusive events for near real-time risk assessment and hazard communication.


Introduction
Frequent eruptive events at an active volcano, coupled with erosive processes, instability episodes and subsequent structural failure, continuously modify the edifice morphology.These phenomena have major implications for the development of monitoring and hazard mitigation strategies (McGuire 1996).Near real-time monitoring and forecasting is a necessary step in mitigating the impact of an eruption and of lava flows.Sharing the produced maps with local civil protection in the first few hours following the start of an eruption helps to build effective mitigation strategies (Chevrel et al. 2022).Topographic changes at an active volcano necessitate updating topographic datasets, re-running lava flow emplacement models and quantification of volume variations (gains and losses) to support hazard assessment 58 Page 2 of 25 and risk mitigation (Felpeto et al. 2007;Cappello et al. 2012;Turner et al. 2017;Neal et al. 2019;Chevrel et al. 2021;Barsotti et al. 2022;Pedersen et al. 2022).The topographic datasets that comprise digital surface models (DSMs) and orthophotos can be updated by processing data acquired by different sensors (thermal and visible cameras, and light detection and ranging, i.e.LIDAR) and platforms (satellites, aeroplanes, helicopters, unoccupied aerial systems, i.e.UASs, and from the ground).Satellite data, and in particular Pleiades optical data, have been previously processed to produce DSMs of volcanic areas and to quantify the eruptive deposits by comparing multi-temporal datasets (Pieri and Abrams 2004;Bagnardi et al. 2016;Ganci et al. 2018;Ganci et al. 2020;Pinel et al. 2020;Proietti et al. 2020;Gouhier et al. 2022;Ganci et al. 2022).In addition to satellite data that allow DSM extraction, satellite shortwave and thermal infrared images have been widely used for the monitoring of active volcanoes since the 1960s, providing repeated lowcost and synoptic observations with a known revisit time (e.g.Francis 1979;Francis and Rothery 1987;Oppenheimer et al. 1993;Wooster and Rothery 1997;Lachlan-Cope et al. 2001;Patrick et al. 2005;Gray et al. 2019;Ganci et al. 2020).Thermal data can be provided at low-spatial (1 to 4 km) but high-temporal (30 min to 3 day) resolution by the advanced very-high-resolution radiometer, the along-track scanning radiometer and the geostationary operational environmental satellite or at high-spatial (30 to 120 m) but lowtemporal resolution (16 days) by the Landsat 7 enhanced thematic mapper plus (Harris et al. 1997(Harris et al. , 1999(Harris et al. , 2000(Harris et al. and 2001)).Data from aerial vehicles have been traditionally used for volcanological monitoring and mapping, also applying the structure from motion processing method (Baldi et al. 2000;Harris et al. 2007;Marsella et al. 2009;Marsella et al. 2014;De Beni et al. 2015;Dvigalo et al. 2016;Neri et al. 2017;Thiele et al. 2017).UASs have been used extensively to derive topographic data for mapping volcanic areas, supporting lava hazard analyses, estimating lava flow and dome volumes and performing high-resolution morphometric analysis (e.g.Turner et al. 2017;Dietterich et al. 2017;Favalli et al. 2018;De Beni et al. 2019;Wakeford et al. 2019;James et al. 2020;Lormand et al. 2020;Andaru et al. 2021;Civico et al. 2021;Carr et al. 2022;Civico et al. 2022;Román et al. 2022).LIDAR data have been acquired on active volcanoes for measuring lava flows and domes as well as crater rim changes (e.g.Favalli et al. 2009;Marsella et al. 2009;Jozkow et al. 2016;Alsadik and Remondino 2020;Dietterich et al. 2021;Gomez et al. 2022).Satellite data allow production of DSMs from sub-metre (WorldView and Pleiades) to metre (Ikonos and Quickbird) sampling distances, relatively high temporal frequency (hours to days, if using multiple platforms), metre-level accuracy and swath width up to two tens of kilometres (Proietti et al. 2020).Surveys with aircraft, helicopters and UAS can be used instead of, or in addition to, satellite data as they can provide very highresolution (cm-pixel) data for mapping and DSM generation (Turner et al. 2017).The advantages of UAS surveys over other remote sensing techniques are as follows: (1) low cost of survey; (2) low-altitude flight operations that result in a centimetre-scale spatial data; (3) automated mission planning that allows acquisition of consistent and repeat datasets for the same area; and (4) increased safety for the pilot compared to in situ measurements (Turner et al. 2017).Moreover, remotely piloted systems represent a platform capable of delivering rapid, reactive and repeatable data acquisition for any location and benefit from a high degree of flexibility in terms of flight planning (Turner et al. 2017).
Advanced remote observation methods now play an important role in the understanding, tracking and monitoring of volcanic processes for both remote volcanoes where field surveys are challenging and for well-monitored sites (Pyle et al. 2013).Volcanic activity generally affects large areas (for example on Etna between 0.03 and 40 km 2 , Coltelli et al. 2017;Branca et al. 2013), and may develop rapidly (seconds to hours).Therefore, it is desirable to integrate data acquired by different in-orbit and airborne remote sensing platforms that have different spatial and temporal resolutions, and images acquired by ground-based sensors (Honda and Nagai 2002;Torres et al. 2004;James et al. 2006James et al. , 2007James et al. , 2012James et al. and 2014;;Ernst et al. 2008;Diefenbach et al. 2012;Ganci et al. 2013 and2020 ;Coltelli et al. 2017;Carr et al. 2019;Di Traglia et al. 2020;Calvari et al. 2021;Caballero et al. 2022;Musacchio et al. 2023).However, critical gaps in the data acquisition capabilities within each of all remote platforms persist due to the rapid, unpredictable and dangerous nature of many volcanic processes, and poor meteorological conditions.For these reasons, a multi-sensor and multi-platform approach is valuable when there is a need to collect data with rapid recurrence, such as to map multiple eruptions over 24 h.On Etna, such behaviour has been well documented since the 1970s, particularly since 2011 and into 2021 (e.g.Tanguy and Patanè 1984;Chester et al. 1985;Behncke et al. 2006;Behncke et al. 2014;De Beni et al. 2015;Ganci et al. 2019;Proietti et al. 2020;Ganci et al. 2022).
Here, we focus on Etna's 2021 eruptive activity which, due to the short intervals between events, required rapid mapping and quantification of each lava flow field (before deposits became covered by the next event), evaluation of the SEC cone growth and estimates of proximal and distal fallout deposit volumes.Our analysis is based on combining data with a range of spatial and temporal resolutions from different multiple remote sensing techniques and platforms (satellite, ground and UAS).

Etna's summit activity
Etna is a composite basaltic stratovolcano, located along the Ionian coast of eastern Sicily (Italy), and characterised by quasi-persistent activity at its four summit craters: Voragine (VOR), North-East Crater (NEC), Bocca Nuova (BN) and South-East Crater (SEC, Fig. 1).The summit craters gave rise to gas emissions, Strombolian to lava fountain activity, and lava overflows.Paroxysms are among the most powerful explosive events to occur at Etna's summit craters (Bombrun et al. 2016;Corsaro et al. 2017).They typically last no longer than a few hours (Barsotti et al. 2010) and are associated with violent Strombolian activity or pulsating to sustained lava fountains, and produce eruptive columns of ash, lapilli and steam, rising up to 14 km above sea level (Corsaro et al. 2017).In most cases, these events also produce fast-moving lava flows that extend a few kilometres from the summit and may continue after the end of the explosive phase (e.g.Alparone et al. 2003;Behncke 2009;Behncke et al. 2014;De Beni et al. 2015;Ganci et al. 2018;Proietti et al. 2020;Corsaro and Miraglia 2022).Summit lava flows are not life-threatening due to their slow propagation, but danger can result from associated phenomena, e.g. the interaction between lava flow and snow or water to create localised phreatic explosions (Chester et al. 1985;Blong 1984).Lava fountains and sustained Strombolian activity produce fallout of ballistic juvenile and lithic bombs in areas less than 5 km from the vent.This fallout represents a danger to tourists visiting the summit of the volcano and sportspeople (both in winter for skiing and in the other seasons for trekking), as well as guides and researchers (Andronico et al. 2015 and2021 ;Osman et al. 2019).
Plumes produced during paroxysms are commonly dispersed by the prevailing westerly winds to deposit tephra over the eastern flanks of the volcano.Lapilli-sized particles and even blocks and bombs can be transported downwind up to several kilometres, while ash-sized particles may be transported over hundreds of kilometres.In a few of Etna's paroxysms, for example on 23 November 2013, ash and tephra fallout was observed up to Calabria and Puglia regions in southern Italy, and as far as 400 km from the volcano (Andronico et al. 2015;Donnadieu et al. 2016).Smallvolume pyroclastic density currents (PDCs, < 2 × 10 6 m 3 ) may also be produced during or after a paroxysm.PDCs not directly associated with the fountain can be due to (1) collapse of pyroclastic fountains; (2) the failure and/or sliding of pyroclastic deposits emplaced on steep slopes during intense lava fountaining; (3) the gravitational collapse of lava flow fronts; (4) the lava interaction/mixing with of wet rocks or snow (Behncke 2009;Norini et al. 2009;Andronico et al. 2018).
Approximately 30 paroxysmal events occurred at Etna's summit craters between 1900 and 1970.More than 180 summit paroxysms have occurred since the 1970s, suggesting they are occurring at an increased frequency (Behncke and Neri 2003).A particularly high number of events (66) were observed between January and August 2000, and most of them occurred at the SEC (Alparone et al. 2003;Behncke et al. 2006).NEC generated its last paroxysms in 1995-1996and 1998, and BN in 1999(Allard et al. 2006;Behncke et al. 2003).A lower number of paroxysms have occurred at VOR, but these are more violent and briefer than those at the SEC, NEC and BN.The VOR paroxysms of 22 July 1998 and 4 September 1999 were among Etna's most powerful eruptive events of the past century (Corsaro et al. 2017).The first event generated a sub-Plinian eruption column that ascended 10 km, resulting in widespread ash fall, mainly to the south (Aloisi et al. 2002).The latter event formed a 5-cm-thick tephra layer in Giarre, ~10 km from the vent, but it was not well documented (Andronico et al. 2015).

Etna's explosive activity 2010-2020
In 2010, Etna exhibited mild degassing activity, followed in 2011 by Strombolian activity from both the BN and the SEC.Between 2011 and 2013, discontinuous and vigorous Strombolian activity at the BN built a small cone on the crater floor and fed a new lava flow (Proietti et al. 2020).At the VOR, two brief periods of eruptive activity in early 2013 and early 2015 consisted of moderate Strombolian explosions and ash emissions (Corsaro et al. 2017).In May 2016, seven paroxysmal events at VOR formed ash plumes up to 14 km high.The related deposits filled the crater and the adjoining BN and overflowed on the western flank from BN lower rim (Proietti et al. 2020).Between September 2019 and mid-May 2020, discontinuous Strombolian and effusive activity on the VOR crater floor caused the growth of an intra-crater cone, and fed lava flows from VOR into BN.Between December 2020 and April 2021, explosive activity occurred again at the VOR, and was associated with sporadic lava effusion (Pioli et al. 2022).
At the SEC, fifty paroxysms occurred between 2011 and 2014 and dramatically changed its morphology with the formation of a new scoria cone, initially called New South-East Crater (NSEC).The overlapping of all the lava flows, emplaced during each paroxysm, produced composite lava flow fields along the east, south and west flanks.After 2013, the eruptive activity led to the coalescence of the NSEC with the SEC (Vicari et al. 2011;Andronico et al. 2014;Behncke et al. 2014;De Beni et al. 2015;Falsaperla and Neri 2015;Acocella et al. 2016;Andronico et al. 2018;Proietti et al. 2020).Between 2017 and 2020, discontinuous vigorous Strombolian and effusive activity occurred at the SEC.In 2020, four paroxysmal events at the SEC produced lava flows that extended 2.8 km to the SE, into the Valle del Bove, and to the SW as well as 1.6 km to the NE, into the Valle del Leone (De Beni et al. 2021;Calvari and Nunnari 2022).This activity culminated in the SEC explosive-effusive events of 2021.

The 2021 SEC activity
In 2021, the SEC gave rise to 58 paroxysmal events in just over 9 months, from 18 January to 23 October 2021.These events fed lava flows that rapidly propagated to the S-SW and into the Valle del Bove and the Valle del Leone, formed eruptive columns up to 11 km in height and caused fallout of ballistics greater than 1 m in diameter within 1-2 km from the SEC, and moderate to light tephra fallout at increasing distances from the proximal to distal areas.The proximal area is that including the original volcanic landform (i.e. the source vent and the SEC cone) and the 'distal area' represents the area covered by a fallout deposit that continuously mantles the topography (Nèmeth and Martin 2007).The first lava fountain occurred on 18 January and was preceded, on 17 and 20 January, by effusive activity at the SEC that fed lava flows towards the Valle del Bove (Andronico et al. 2021;Marchese et al. 2021;Calvari and Nunnari 2022).About a month later, a series of paroxysmal events began at the SEC, with 17 lava fountains occurring between 16 February and 1 April 2021, each lasting between 0.8 and 10.5 h, with an average duration of 2.9 h.The time interval between two successive events ranged from just over 1 day to just under 8 days (Table 1).On 19 May 2021, paroxysms resumed after a month and a half of quiescence.Until 23 October, the SEC produced 40 paroxysms that lasted between 0.5 and 4.9 h, with an average duration of 1.6 h.The Table 1 Progressive number of the paroxysmal event (PE) of the first cycle; event date (day and month, the year is 2021 for all events); start time of lava flow (LFl); start, end time and duration of lava foun-tain (LFo); time interval between two lava fountains (LFos); dates are shown, along with start and end times for events that occurred on two or more days; -denotes missing information  These events are described in the weekly 'Bollettino Etna' activity reports (https:// www.ct.ingv.it/ index.php/ monit oragg io-e-sorve glian za/ prodo tti-del-monit oragg io/ bolle ttinisetti manali-multi disci plina ri? filter% 5Bsea rch% 5D= Etna).
Lastly, images acquired by the INGV-OE permanent network of surveillance cameras were orthorectified and analysed.
To produce a DSM of Etna that quantifies the pre-paroxysm topography with a spatial resolution of 1 m and a vertical accuracy of 1 m, we processed a triplet of Pleiades satellite visible images acquired on 22 August 2020 (image spatial resolution of 0.5 m, acquired on demand) using the MicMac open-source library (Rupnik et al. 2017).DSM accuracy was evaluated by measuring the standard deviation of the difference between the 2020 DSM and a 2015 Pleiades DSM in areas where no topographic changes occurred (Ganci et al. 2019).
We analysed orthorectified images acquired by Skysat and PlanetScope (spatial resolution of 0.5 and 3.7 m, respectively; acquired on demand and every day, respectively) by setting thresholds on the near-infrared and red (~0.8 μm and ~0.6 μm, respectively) bands (Ganci et al. 2019).
Active or recently erupted lava flows were identified in orthorectified images acquired by Sentinel-2 MSI (spatial resolution of 10 m and 20 m, for visible and thermal images respectively; temporal resolution of 2-3 days) and Landsat 8 OLI sensors (spatial resolution of 30 m; temporal resolution of 7-9 days).The flows were identified by thresholding in the short-wave infrared (SWIR) bands (Bands 7 and 12 for Landsat 8 and Sentinel-2, respectively) centred at ~2 μm (Ganci et al. 2018;Cappello et al. 2019).
If satellite or UAS images were not available or adequate for satisfactory mapping, we used images acquired by Monte Cagliato and Schiena dell'Asino thermal cameras (EMCT and ESR, frame rate of 2 frame/s; Fig. 1).These two sensors belong to a permanent network of ground-based surveillance cameras (visible and thermal) that is installed on Etna's flanks and is managed by the Istituto Nazionale di Geofisica e Vulcanologia-Osservatorio Etneo (INGV-OE https:// www.ct.ingv.it/ index.php/ monit oragg io-e-sorve glian za/ segna li-intempo-reale/ video-sorve glian za-vulca nica-etna).At the end of the paroxysms, we selected the images that effectively showed the lava fields.We orthorectified these images by re-projecting them onto the 2020 DSM using an algorithm that considers the position and orientation of each camera (Ganci et al. 2013;Calvari et al. 2021), to obtain georeferenced thermal images (spatial resolution of 5 m).For image orthorectification, it is of primary importance to use an upto-date DSM of the volcano since morphological changes will increase error in the image georeferencing and orthorectification (Ganci et al. 2013).
In 2021, we carried out nine UAS surveys with the aim of monitoring eruptive activity and acquiring data useful for photogrammetry processing (Table 3).Six surveys were performed over the SEC crater, using a DJI Mavic Enterprise.The high-resolution visible (12 Mpx) and thermal images (160 × 120 px) collected enabled monitoring the morphostructural evolution of the SEC summit portion, which was deeply fractured and potentially unstable.Two surveys were performed over the Valle del Bove, where most of the lava flows were emplaced, using a DJI Mavic Enterprise and a DJI Mavic 2 Pro to collect visible images (20 Mpx).One UAS survey was performed in the area of Cisternazza, a pit crater that formed during the 1792 eruption (Fig. 1), using a DJI Phantom 4 Pro to collect high-resolution visible images (20 Mpx).A previous survey in the Cisternazza area had been carried out in 2020 using the same UAS (Table 3).
The UAS surveys only partially covered the areas to be monitored; images collected above the Valle del Bove covered only a part of the 2021 lava flows, and the SEC surveys generally covered only its south and east sides.These partial coverages reflect the difficulty of UAS surveying in hostile environments such as the summit area of Etna and the Valle del Bove.Due to the need to operate from a safe location, it was only possible for the UAS to take off from the 2002-2003 scoria cone (Fig. 1).From this takeoff point, the north side of the SEC is not visible and UAS flight regulations did not allow operation beyond visual line of sight.Not all flyovers yielded useful images for photogrammetric analysis due to the presence of gas and clouds that obscured the ground.Suitable visible images (Table 3) were processed through a structure-from-motion photogrammetric software (Agisoft Metashape v1.6.6, https:// www.agiso ft.com.) and we extracted DSMs and orthomosaics with resolutions between 4 and 55 cm/px.During data processing, the SEC and the Valle del Bove datasets were georeferenced using only the aircraft position measured by the on-board GNSS.This is because we were unable to position the ground control points (GCPs) around the SEC cone and inside the Valle del Bove due to volcanic hazard and the extremely rough terrain.In contrast, we georeferenced and processed the data acquired on the Cisternazza area using 6 GCPs, whose positions were measured through a Glonass GNSS with centimetric to sub-centimetric precision.

Rapid lava flow mapping
We rapidly mapped almost all the lava flows in ESRI Arc-GIS software.For seven events (PEs from 22 to 24, from 26 to 28, and 36 in Table 2), we did not have satellite, UAS or surveillance camera images available owing to poor weather conditions.Therefore, we did not map the corresponding lava flows.In three more cases, we had only one reference dataset for two paroxysmal events that were very close in time (PEs 30 and 31; PEs 42 and 43; PEs 44 and 45 in Table 2).Because we were unable to distinguish and separately map the individual lava flows, we produced three maps, each covering the cumulative deposits from two paroxysms.For each mapped lava flow, we measured geometrical parameters of maximum length, minimum front elevation and covered area.We recorded these measurements into a geodatabase GIS along with the analysed datasets (Tables 4 and 5), the start time of lava fountains (LFo) and lava flows (LFl), as well as the end time of lava fountains (Tables 1 and 2).These times are determined by INGV-OE control room staff from surveillance camera images, and are reported in the "Bollettino Etna" activity reports.The end of both lava effusion and of the front advance is not exactly identifiable from image analysis due to the gradual nature of the cooling process.After the emplacement of each lava flow, an initial bulk volume estimate was recorded in the geodatabase, based on a planimetric approach, i.e. multiplying the measured area by a plausible average thickness of 2 m, estimated from ground surveys (Stevens et al. 1999).We overlaid all the mapped lava flows to produce a 2021 resurfacing map for Etna (Proietti et al. 2011), which indicates how many times each cell (of the 2020 DSM) was covered by lava in 2021.We also sequentially superimposed the map of each lava flow, starting from that of 16 February, to follow the temporal evolution of the 2021 resurfacing map.

A posteriori cross-section analysis of the lava flows
Because high-resolution DSMs were not available for each event, we performed (in ESRI ArcGIS) an a posteriori analyses to evaluate the bulk volume of each of the 2021 lava flows by applying the cross-section approach (De Beni et al. 2021).This method is based on two assumptions: (1) the thickness of a simple lava field remains constant on a surface with a constant slope; (2) comparable average thicknesses are associated with comparable slopes.Each of the 2021 eruptive events lasted from a few hours to less than 3 days, thus formed simple lava flows (sensu Walker 1971Walker , 1973) ) whose emplacement and thickening are essentially governed by slope changes.We evaluated the slope from the 2020 DSM and defined seven polygons that envelope areas with homogeneous slope values (Fig. 3a).Then we used these polygons to cut each simple lava flow into parts.We drew eleven cross-sections along the lava flows emplaced on 23 and 31 March 2021 in Valle del Bove and towards S-SE, respectively (Fig. 3a and b).The cross-sections were located in two slope areas (0-25° and 26-34°) and were plotted where the flow morphology and their thicknesses are unambiguously measurable on the 27 March and the 16 September 2021 DSMs, respectively.The first DSM was acquired towards the end of the first eruptive cycle, shortly after the emplacement of the 23 March lava field.The S-SE branch of the March 31 lava flow was not covered by successive events, so it could be measured on the latter DSM.Interpolation of the preeruption surface at the eleven cross-sections allowed us to estimate the average thickness to be assigned in each slope area and to the flow portions inside them (Fig. 3c-g).We calculated the bulk volume of each flow through the planimetric approach (Stevens et al. 1999).We estimate an error of 20% (Table 6) for these volumes, which is the magnitude commonly associated with this approach (Calvari et al. 1994).For each lava flow that we did not map (owing to poor weather, # in Table 5), we assigned a volume equal to that measured for an event with the same flow direction and lava fountain duration.For each pair of paroxysms very close in time with only one reference dataset (*, **, *** in Table 5), a total volume (V tot ) was calculated by multiplying the measured area (which is cumulative for both events) by doubling the average single-flow thickness, as evaluated from the cross-section analysis (described above).This approach might overestimate lava volume, but we assume to less than the estimated error of 20% (Table 6).The total volume was then distributed between the two events (1 and 2) proportionally to their durations (D 1 and D 2 ), e.g.:

A posteriori topographic analysis of the lava flows
We also estimated bulk lava flow volumes through a posteriori analysis of the available DSMs in ESRI ArcGIS.We applied the topographic approach that calculates the elevation differences between datasets acquired on different dates before, during or after an eruption (Stevens et al. 1999).In particular, we compared the 27 March 2021 DSM and the 2014 Etna's DSM (pixel size of 1.5 m, De Beni et al. 2015) in the Valle del Bove.The 2014 DSM is less up-to-date than the 2020 one; however, after its acquisition, only one lava flow was emplaced on the analysed zone (06-08 December 2015, Proietti et al. 2020).We preferred to use this reference surface because, inside the Valle del Bove, it is more accurate and less noisy than the 2020 one and therefore it allows for a better estimation of the lava flow thicknesses.The topographic comparison of the two DSMs allowed measuring the thickness accumulated during the paroxysmal events from 1 to 16 (Fig. 4a and Table 1).The measured value was divided by the number of overlapped lava flows (Fig. 4b) to evaluate an average thickness of 2.0 m for each flow unit.This value is congruent with the average thickness estimated by ground observation and measurements, which were used for the preliminary volume evaluation.It is also congruent with the average thicknesses assigned to slope areas D and E based on the cross-section analysis (Fig. 3g).We analysed a cross-section cutting the Monti Centenari scoria cone (formed during the 1852 eruption), we identified the overlap area of the different lava flows and we distinguished the different contributions to lava thickness (Fig. 3e).By comparing the 2021 and 2014 DSMs, we also measured the lava thickness accumulated within the Monti Centenari scoria cone (formed during the 1852 eruption) during the 9 and 12 March 2021 events (Fig. 3f).

A posteriori topographic analysis of the SEC cone
We measured the bulk volume changes of the SEC cone and its elevation through topographic comparison (Stevens et al. 1999) in ESRI ArcGIS, of the 16 September 2021 and the August 2020 DSMs.In this case, we preferred the 2020 DSM over that from 2014, because it better represents the pre-2021 topographic surface.The 2021 DSM does not include the NW portion of the SEC, which corresponds to 35% of the cone basal area (as measured with ESRI ArcGIS).However, this side of the SEC underwent minor morphology changes with respect to the south and east ones because most tephra are generally dispersed towards the SE-ESE by the prevailing westerly wind.The missing elevation changes were estimated by applying the Natural Neighbor interpolation tool of ArcGIS (https:// deskt op.arcgis.com/ en/ arcmap/ 10.6/ manage-data/ terra ins/ hidden-natur al-neigh bor-inter polat ion.htm) to the values measured in adjacent areas of the SEC, where no lava flows were emplaced.The interpolated volume corresponds to 11% of the total volume changes.To improve the quantification of the SEC volume changes, we subtracted the value measured for the portion of lava flows emplaced on the cone to that obtained from the DSM comparison.
In the topographic approach, the maximum volume error is linearly dependent on the standard deviation of the elevation change (σ Δz ), calculated from regions where no change occurred (De Beni et al. 2019).However, the September 2021 DSM does not include unchanged areas for measuring σ Δz .This DSM is limited to the SEC cone, which was changed significantly by the paroxysmal events, including its southern base where ballistic and fallout products were deposited.Volume changes of the SEC cone, between 2011 and 2016, were previously measured through the topographic approach, with an associated error of 4% (Proietti et al. 2020).Therefore, we assume that the SEC volume error ranges between 4 and 20%, as previously found for topographic and planimetric analyses (Table 6).It should also be considered that in the natural neighbor interpolation of a smooth surface, such as that of the NW flank of the SEC, the volume error can be as high as 40% (Table 6; De Beni et al. 2019).Considering that the contribution to SEC volume change measured from data interpolation is 11% of the total change, this gives an error of 4.4%, which should be added to that from topographic analysis, resulting in a total error of about 8%.
The error on volumes evaluated by comparing two point clouds, without GCPs or over an area where no natural change occurred, has been previously estimated with the Cloud Compare software through the Multiscale Model to Model Cloud Comparison and its successive precision map modification (M3C2 and M3C2-PM, Lague et al. 2013;James et al. 2017).We cannot apply the precision map calculation (M3C2-PM) because we do not have precision estimates for the 2020 Pleiades DSM (generated by the MicMac software).We therefore applied the plugin M3C2, considering diameter values between 5 and 20 m for both the normal scale and projection scale and we restricted the differencing to the vertical direction.The obtained distance uncertainties were rasterised to the same spatial resolution of the two DSMs (1 m).However, the obtained confidence intervals are smaller than the 1% of the volume changes evaluated from DSMs comparison.
Ultimately, considering that we calculated the total volume of the SEC cone by combining multiple methodologies with different uncertainties, we believe that an error of 20% (already considered for lava volumes, Table 6) is a conservative estimate for our analysis.

A posteriori analysis of the fallout deposits and ballistics
We measured the thickness of the proximal fallout deposits at the SEC south base by comparing the 16 September 2021 UAS-derived DSM and the 22 August 2020 Pleiades DSM, in ESRI ArcGIS.This deposit accumulated from 53 events (of the 57 paroxysms occurring between the two surveys) with plume dispersion towards the WSW-ESE.The analysis was limited to an area (~2.7 × 10 4 m 2 ) over which the lava flows had not spread.The thickness of the distal fallout deposits from 27 events (with plume dispersion towards the S-ESE) was determined using DSMs from the two UAS surveys (05 August 2020 and 16 September 2021) above and around the Cisternazza pit crater (Fig. 1 and Table 3), which is about 2.5 km SSE from the SEC.This comparison covered the central flat area inside the Cisternazza (~0.3 × 10 4 m 2 ) and a ring (~4.0 × 10 4 m 2 ) around the pit crater.To estimate the quality of the co-registration between the two datasets, we identified 16 easily recognisable morphological elements, such as overhanging and steep rocks, for which wind would have prevented any fallout deposit accumulating.We measured their 3D coordinates (X, Y and Z) in both DSMs and subtracted the 2020 values from the 2021 values.We obtained mean differences around zero (0.04 m, −0.06 m and −0.02 m for X, Y and Z, respectively), with associated standard deviations of 0.09 m, 0.10 m and 0.14 m.We thus considered 0.14 m as the error associated with thickness measurements in the Cisternazza area.
Finally, we performed a visual analysis of the orthomosaics extracted from the 16 September 2021 UAS survey for an area at the south base of the SEC, between ~3050 and 2900 m of elevation.This analysis enabled us to identify ballistics larger than 20 cm and to record their location in a point-type shapefile.Their spatial distribution was then assessed using the Kernel Density tool (in ESRI ArcGIS), which calculates the local neighbourhood feature density.

Quantification of the volume of magma erupted
We measured the bulk volumes of all the lava fields and the changes of the SEC cone due to the 2021 paroxysmal events by applying the cross-section and topographic analyses.We converted bulk volumes to DRE volumes using vesicularity values from the literature of 20% and 50% for the lava flows and the cone, respectively (Behncke et al. 2006, Table 6).A range of vesicularity values has been previously assumed for Etna effusive and explosive products (Andronico et al. 2018) but 20% and 50% values have been used previously for the 2011-2016 SEC eruptive activity, which was similar to 2021 (Behncke et al. 2014;De Beni et al. 2015;Corsaro et al. 2017;Proietti et al. 2020).We also estimated the total DRE volume of distal tephra by considering the percentage of lava flows (73%), proximal (23%) and distal (4%) pyroclasts previously found for the 25-26 October 2013 event (Andronico et al. 2018, Table 6).The authors proposed that a similar distribution of deposits, showing a predominant output of lava flows and coarse tephra deposited around the cone, with respect to the distal tephra, was observed during the paroxysmal phase of most of the Etna lava fountains (Andronico et al. 2018).We finally combined the DRE volumes measured for the 2021 composite lava flow fields and for the SEC changes with the estimate of the distal fallout deposits, to quantify the DRE volume of the erupted magma.

Results
For the first cycle (from 16 February to 1 April) and in the Valle del Bove, we find that the longest lava flow reaching the lowest front elevation (4.3 km and 1700 m a.s.l., respectively) was fed by an effusive vent on 24 February (PE 6, Fig. 2, Table 1 and Table 4).For the same cycle but on the SW flank, the longest lava flow reached the lowest front elevation (3.5 km and 2600 m a.s.l., respectively) on 31 March (PE 17, Table 1 and Table 4).For the second cycle (from 19 May to 23 October), the longest lava flow in the Valle del Bove (2.8 km and 1950 m a.s.l.) was fed by an effusive vent on 09 August (PE 54, Table 2 and Table 5) and for the SW flank (2.9 km and 2350 m a.s.l.) was on 08 July (PE 50, Table 2 and Table 5).
The ranges of the bulk volume of each paroxysm are (1.2-4.8)× 10 6 m 3 and (0.1-2.1 × 10 6 m 3 ), for the first and the second eruptive cycles, respectively.The corresponding ranges of the DRE volumes are (1.0-3.8)× 10 6 m 3 and (0.1-1.6 × 10 6 m 3 ), respectively (Table 6, Table 7, and Fig. 5).For the first and the second eruptive cycles, the paroxysmal events that erupted the minimum lava volumes occurred on 19 March and 22 May, respectively.Those that erupted the maximum lava volumes occurred on 17 February and 9 August, respectively (Fig. 5).The bulk volume of the composite lava flow field from the first cycle (17 paroxysms) is (42 ± 8) × 10 6 m 3 and corresponds to an average volume per event of (2.5 ± 0.5) × 10 6 m 3 (Table 7).For the second cycle (40 paroxysms), it is (29 ± 6) × 10 6 m 3 and corresponds to an average volume per event of (0.7 ± 0.1) × 10 6 m 3 (Table 8).These measurements give a total bulk volume for 2021 of (71 ± 14) × 10 6 m 3 , corresponding to an average volume per event of (1.2 ± 0.3) × 10 6 m 3 , and a DRE lava volume of (57 ± 11) × 10 6 m 3 (Table 9).The elevation changes to the SEC scoria cone, derived from comparison of the 16 September 2021 UAS and the 2020 Pleiades DSMs, resulted in a maximum increase of 135 m.The top of the SEC, located on the crest of its south-east rim, reached on 16 September 2021 an altitude of 3348 ± 1 m a.s.l.(Fig. 6).The measured changes of the SEC result in   a bulk volume of (27 ± 5) × 10 6 m 3 (Table 9), corresponding to an average volume per event of (0.5 ± 0.1) × 10 6 m 3 , and a DRE cone volume of (13 ± 3) × 10 6 m 3 (Table 9).The total DRE volume of magma erupted in 2021 is (73 ± 15) × 10 6 m 3 and it include the lava fields, the SEC change and the estimated distal fallout (3 ± 1) × 10 6 m 3 .The combination of the lava flow maps produced for each paroxysm provides the 2021 Etna resurfacing map.This map counts how many times each pixel has been covered during all the 2021 paroxysmal events and provides a quantitative description of the spatio-temporal evolution of the most frequently covered areas (Fig. 8a; an animation of its temporal evolution is available as online resource 1).
The proximal fallout deposits accumulated at the SEC south base have average and maximum thicknesses of 19 and 26 m (Fig. 6a) and become thinner further from the cone.This material was deposited during 53 paroxysms with plume dispersion towards the WSW-ESE, so the corresponding average deposited thickness per event is 0.4 m.The distal fallout from the 27 events with plume dispersion towards the S-ESE was derived from comparison of the 2021 and 2020 DSMs around the Cisternazza pit crater (Fig. 7).Its thickness has a maximum measured value of 3 m in the base of the Cisternazza (average value and standard deviation of 0.80 m and 0.18 m, respectively).This corresponds to an average deposit thickness per event of up to 0.11 m.The fallout is thinnest in the area around the pit crater, where the maximum measured value is 0.50 m (average value and standard deviation of 0.15 m and 0.07 m, respectively).This corresponds to an average deposit thickness per event up to 2 cm.
From the analyses of the high-resolution orthomosaic (16 September 2021 UAS survey), we observe that the SEC south base was highly impacted by ballistics.In particular, ballistics between 600 and 900 m from the vent have a discontinuous distribution and a maximum size of 4 m.The larger ones rolled down the slope for at least 750 m (Fig. 6b, c).Moreover, ballistics smaller than 1.5 m are dispersed discontinuously from 900 up to 1200 m from the vent and intercept the tourist track.We present the density map of ballistics drawn for the area at the SEC south base, through a kernel density analysis of the position of ballistics with a linear dimension greater than 0.2 m (Fig. 8b).

Discussion
This work presents the results of the timely mapping of the SEC cone and of the lava flows emplaced on Etna volcano, from 57 paroxysmal events in 2021.These events occurred with a minimum interval of 4 h and mapping was based on the synergistic integration of multi-sensor and multiplatform data, i.e. visible, thermal, near-infrared and SWIR images, acquired from satellites, UAS surveys and the INGV-OE ground-based surveillance cameras.This ensured a constant data flow at multiple spatial and temporal resolutions (from 1.9 cm/px, for UAS orthophotos, to 30 m/px for Landsat-8 orthoimages; from less than 1 day to less than 6 days after an event has occurred), allowing detailed spatial and temporal coverage and understanding of Etna's 2021 eruptive events.The lava flow maps and relative measurements were promptly released to the Italian Civil Protection Department through the "Bollettini and Comunicati" (standardised documents issued by the INGV-OE, to communicate the alert level and briefly describe volcanic activity).We quantify the bulk volume of each lava flow, of the composite flow fields (for the two eruptive cycles and the whole 2021 activity) and of the SEC cone.From these, we derive the average volume per events (for both the flows and the SEC), and the DRE volume of erupted magma.We also measure the average thickness per event of fallout deposits (both proximal and distal).Finally, we evaluate the lava flow resurfacing map and the density map of ballistics at the SEC south base.
The higher average volume per event measured for the first cycle with respect to the second one (2.5 × 10 6 m 3 versus 0.7 × 10 6 m 3 ) is mainly due to the February paroxysms that erupted almost double the lava volume of the subsequent events.In March 2021, erupted volumes decreased but they were still generally greater than those of the second cycle (Fig. 5).The total DRE volume of magma erupted during the 57 SEC paroxysms (16 February-23 October 2021) is (73 ± 15) × 10 6 m 3 and corresponds to an average output rate of about 3.4 m 3 /s.This value is higher than the corresponding values measured between 22 August 2020 and 29 June 2022 (2.88; Ganci et al. 2023) and much higher than Etna's long-term steady state output (0.8-0.9 m 3 /s; Wadge and Guest 1981;Bonaccorso and Calvari 2013).However, the 2021 higher output rate can be explained by the low eruptive activity observed between 26 May 2016 (end of the 2011-2016 eruptive period, Proietti et al. 2020) and 17 February 2021 (before the beginning of the 2021 first eruptive cycle).In this period, Etna produced only short-lived lava flows and five paroxysms.The erupted total DRE volume is less than 30 × 10 6 m 3 (Cappello et al. 2019;De Beni et al. 2021;"Bollettino Etna" (n.d.) which is much lower that the value estimable for a steady-state accumulation over a period of about 4.7 years (Bonaccorso and Calvari 2013).
The high-resolution DSM extracted from the 16 September 2021 UAS survey (Table 3) enabled us to measure that the SEC cone summit (3348 m a.s.l.) is the new peak of Etna volcano (Fig. 6), now higher than the NEC (3321 ± 1 m a.s.l., in the 2020 DSM).
Our resurfacing map shows that, at the end of the 2021 eruptive activity, the most affected area was the southwest base of the SEC, where up to 35 flows overlap.This zone is less than 1 km from the guides' information point, which is frequently visited by tourists, that is to Fig. 8a.This result highlights the need to provide near real-time monitoring and forecasting data to civil protection authorities.Such information is pivotal for managing access for volcanologists, stakeholders, and tourists and mitigating the hazard related to summit lava flows.In fact, summit lava flows are slow enough to allow people to move away from their fronts but, in some cases, they can be hazardous to any people near them.For example, the interaction between lava flows and snow or water can create localised phreatic explosions (Chester et al. 1985;Blong 1984).The resurfacing map could support the assessment of the areas most exposed to inundation by Etna's summit lava flows and can be used as input and checking data of lava flow numerical models (Del Negro et al. 2013).
The density map of ballistics (Fig. 8b) should only be interpreted qualitatively because it underestimates the actual ballistics number density.Many ballistics have been mantled subsequently by ash and lapilli deposition.In addition, on impact, many ballistics break into parts too small to be identified in the orthophoto.Nevertheless, the density map highlights the risk associated with ballistics especially near the tourist track that was intercepted at about 2900 m.Our analysis shows that the SEC south base, which is close to the tourist footpaths leading to the 2002-2003 scoria cones and the summit craters, was highly impacted.These results could support studies aimed at modelling the trajectories of ballistics during lava fountains, with the aim of identifying the affected area in near real-time and quantifying the hazard (Alatorre-Ibargüengoitia et al. 2012;Fitzgerald et al. 2017 and2014 ;Spanu et al. 2016;Costa et al. 2023).
The measured maximum thickness of the distal fallout deposits around the Cisternazza pit crater (0.50 m) has the same order of magnitude as its (0.14 m) error.This small accumulation could be due to tephra remobilisation by strong winds.Greater thicknesses, ranging between 0.50 and 3.00 m, have accumulated on the lowest portion of the Cisternazza pit crater.However, these are likely influenced by the funnel-like crater geometry, with steep walls up to 60 m deep, acting as a natural collector for both fallout and reworked deposits, as well as the debris material eventually detached from the boundary walls.
Our volume estimates align with those of previous work on 2021 Etna's eruptive activity.The bulk volume (71 × 10 6 m 3 ) we obtained for the composite lava field formed during the 2021 paroxysms is in agreement with its recent estimates from Pleiades DSMs (Ganci et al. 2022(Ganci et al. , 2023)).This value is also comparable to that measured for long-lived (420 days) lateral eruptions, i.e. during 2008-2009(Ganci et al. 2012)).We also note that the volume estimates for each of the 2021 lava flows are well within the range of uncertainties of the corresponding values estimated from thermal satellite data (Ganci et al. 2023).The bulk lava volume we measured for the 12 March event (2.3 ± 0.5 × 10 6 m 3 ) is compatible with that (2.4 × 10 6 m 3 ) measured previously by analysing the strain changes, the images acquired from the INGV-OE surveillance cameras, and satellite data (Calvari et al. 2021).The total volume erupted by the Etna lava fountains of February-March 2021 was previously estimated as 52.5 × 10 6 m 3 from strain changes measured by high-precision borehole dilatometers (Bonaccorso et al. 2021).This result is compatible with our DRE magma volume estimate for the first cycle (43 ± 9 × 10 6 m 3 ), by adding the value measured for lava flows and that extrapolated for the SEC.The latter is obtained by redistributing the overall 2021 volume of the SEC based on the ratio between the lava volumes measured for the first cycle and all the 2021 events through: where the subscript SEC or lava indicates that the volume (V) refers to the SEC or lava flows, respectively, whereas the superscript I cycle or tot indicates that the volume (V) is related to the events of the first cycle or all the 2021 events, respectively.
The 2021 activity (average volumes per event of 1.2 × 10 6 m 3 and 0.5 × 10 6 m 3 for the lava flows and the SEC, respectively) can be compared with the 50 paroxysms occurring in 2011-2014.The 2011-2014 bulk volumes for lava flows and the New South-East Crater (NSEC) were 84 × 10 6 m 3 and 50 × 10 6 m 3 , respectively, with corresponding average volumes per event of 1.7 × 10 6 m 3 and 1.0 × 10 6 m 3 (De Beni et al. 2015).Although the lava flow volume per event estimates are similar, estimates of average SEC volume increase per event in 2021 are approximately half of those for 2011-2014.This difference could be explained by both the changing conditions of cone growth and the timing between 2011 and 2014 compared with 2021.In 2011, the NSEC started to grow from a small pit crater (about 5 × 9 m) on the lower east flank of the "old" SEC.The proximal pyroclastic deposits accumulated into a scoria cone of comparable size to the old one (Behncke et al. 2014;De Beni et al. 2015).During 2011-2014, the time intervals between consecutive events were longer than those in 2021, allowing the NSEC to be more stable than the SEC.The only lateral collapse occurred on 11 February 2014 (Andronico et al. 2018), when the NSEC cone had equalled in size the old SEC.From 13 December 2020, there were several lateral collapses of the SEC, due to the formation of fractures, combined with a frequent explosive activity, which weakened its east and south summit areas.The lateral collapses were associated with pyroclastic flows (on 16 February, 28 February, 24 March and 23 October).The 2021 proximal pyroclastic deposits accumulated on a yet-constructed cone whose flanks, having an average slope of 30% (measured from the UAS DSMs), had already reached the repose angle for dry pyroclastic deposits (32-33%, Wood 1980).Therefore, the cone could no longer accumulate further considerable volumes of pyroclastic products, which flowed down and accumulated at its base, where the slope decreases.This hypothesis is supported by the proximal deposit thickness measured at the lower SE flank of the SEC, ranging between 10 and 26 m (Fig. 6a).It is also reasonable to hypothesise that a large part of the deposits flowed into the Valle del Bove.
Our analyses confirm that the integration of satellite data with those acquired from other platforms (UAS and ground cameras) makes it possible to monitor eruptive events whose frequency is greater than the satellite's revisit time.The analysis of thermal images, acquired from one of the three remote platforms, in addition to visible images, makes it possible to distinguish active lava flows from recent cooling ones, especially when the eruption frequency is high (e.g. up to three events in the same day).UAS can provide valuable data for volcano monitoring.However, it must always be assessed whether the level of volcanic activity and weather conditions permit a flight.Weather conditions can change very quickly, going from sunny to overcast or foggy in a matter of minutes.Such abrupt changes could compromise flight performance and survey results and, in the worst case, cause a loss of control of the drone.Thus, automated return-to-base is a valuable capability for a UAS.Low temperatures (e.g.−10 °C, or lower) and moderate to very-high-speed wind (more than 20 knots, 37 km/h) can reduce battery capacity by up to 50% and require surveys to be split into multiple flights.A strong knowledge of the volcano and expertise in piloting in extreme conditions are required to maximise survey success.Finally, even if fixed ground-based cameras allow observation of only some areas of the volcano's flanks, their real-time image acquisition supports volcano surveillance and mapping when satellite and UAS data are not available.

Conclusions
The integration of multi-sensor and multi-platform remote sensing data has proven effective in monitoring and quantifying high-frequency eruptive events, such as the 2021 Etna paroxysms (with up to three events in one day).Rapid results from our analysis, which were promptly released to the Italian Civil Protection Department, are the lava flow maps for most of the eruptive events and the measurement of the related geometric parameters.We also released three datasets that could support further quantification of the hazard from summit eruptions, namely the 2021 Etna's resurfacing map, its temporal evolution and the ballistics density map.Further outcomes are the quantification of the bulk volumes for the lava field, and the change of the SEC cone, and the corresponding DRE values.The comparability of the 2021 volume with that of the 2008-2009 eruption, the high output rate of the 2021 paroxysms and the low volume erupted in the previous period (26 May 2016-17 February 2021) lead us to conclude that the occurrence of a repeated series of many low-volume events is sufficient to restore the volcano's equilibrium state, analogously to a long-lasting lateral eruption.Therefore, by the end of 2021, Etna volcano has reached its equilibrium state.Thereafter, there have been only two paroxysmal events, in February 2022, two mild and purely effusive eruptions in Valle del Bove and Valle del Leone (11 May-13 June 2022 and 27 November 2022-7 February 2023) and two paroxysmal events in May and August 2023.To date, 2 months after the last event, the volcano only shows mild degassing.The methodology implemented during the 2021 paroxysm monitoring is also applicable to a potential lateral eruption of Etna.In this case, mapping and quantification of erupted products would be more urgent, as such flows could arrive near inhabited areas, posing a greater risk than summit eruptions.Our strategy can be exported to other volcanoes and extended to different environmental monitoring.

Fig. 1
Fig. 1 Composite lava flow field from all the eruptive events of 2021, superimposed on shaded relief extracted from a Pleiades triplet of 22 August 2020; contour lines are shown every 200 m, elevation is ellipsoidal.The figure also shows the summit craters, the two thermal cameras (EMCT and ESR), and the location of effusive vents (white stars which are not to scale).The inset shows the location of Etna with respect to Sicily.SC stands for scoria cones

Fig. 2 a
Fig. 2 a Lava flows of the first eruptive cycle; b lava flows of the second eruptive cycle.In both figures the background image is the Pleiades shaded relief of 22 August 2020; contour lines are shown

Fig. 3 a
Fig. 3 a Slope map superimposed on shaded relief.Fuchsia and grey polygons delimit the lava flows of 23 and 31 March 2021, respectively; white polygons delimit the areas of homogeneous slope, black lines identify the cross-sections measured on the DSM of 27 March 2021, the blue square indicated the area enlarged in panel b; b limits of the flows for which we estimated the average thicknesses by crosssection, the colours of the polygons correspond to slope areas (see the legend in panel a); the red square indicated the area enlarged in

Fig. 4 a
Fig. 4 a Cumulative lava thickness of the 2021 lava flows (LFls) up to the 27 March 2021 plus the December 2015 flow; b corresponding resurfacing map that counts the number of accumulated flows

Fig. 5
Fig. 5 DRE volume (V) and corresponding error bar (black lines, calculated as plus or minus 20% of the corresponding value) of each lava flow of the first (a) and second cycle (b); the different ranges of the x and y axes of the two figures were chosen to facilitate the discrimination of individual events

Fig. 6 a
Fig. 6 a Elevation change of the SEC cone and thickness of the proximal deposits, the red lines are the limits of the 2021 composite flow field, the white arrow indicate the top of the SEC, and white dashed rectangle delimits the area showed in panel b; b orthomosaic (16 Sep-

Fig. 7 Fig. 8 a
Fig. 7 Thicknesses of the distal deposits measured inside and around the Cisternazza pit crater; the upper left inset locates the Cisternazza surveys (filled red polygon) with respect to the Etna summit area.Background images overlap the 16 September 2021 shaded relief and that of 22 August 2020

Table 3
Summary table of the UAS surveys; UAS model; survey target and goal indicate the area of interest and the focus of the overflight; the resolution of the derived DSMs and orthomosaics are given DSMs through the Global Mapper software (v22.0,https:// www.bluem arble geo.com/ global-mapper/) and managed the DSMs with the ESRI ArcGIS software.
We performed real-time correction of the GNSS data using a reference station of the INGV permanent GNSS network (RING http:// ring.gm.ingv.it/).Before analysing the elevation changes of the SEC and the Valle del Bove, we co-registered the different pairs of point clouds by manual and automatic alignment through the Cloud Compare (v2.12.4) open-source software (https:// www.danie lgm.net/ cc/https:// www.danie lgm.net/ cc/).In this way, we improved the correspondence between two datasets acquired at different dates (pre and post-eruption).We then exported the point clouds into

Table 4
Progressive number of the paroxysmal event (PE) of the first cycle; datasets analysed for mapping the lava flows (acquisition date: day and month, the year is 2021 for all datasets); directions, maxi-

Table 5
Progressive number of the paroxysmal event (PE) of the second cycle; datasets analysed for mapping the lava flows (acquisition date: day and month, the year is 2021 for all datasets); directions, maximum length (Max L) and minimum (Min) front elevation (ellip-soidal values) of the lava flows (LFls); *, ** and *** dataset of cumulative maps, each referring to two paroxysms very close in time; # not mapped lava flow owing to poor wheatear; -missing information

Table 6
Parameters used

Table 7
Statistics (minimum, maximum and arithmetic mean) of the average thickness (Av thick), area (A), bulk and DRE volumes (V bulk and V DRE ).These values (including the arithmetic mean of the average thickness) are calculated by considering all the lava flows, and the composite lava field, of the first cycle

Table 8
Statistics (minimum, maximum and arithmetic mean) of the average thickness (Av thick), area (A), bulk and DRE volumes (V bulk and V DRE ).These values (including the arithmetic mean of the average thickness) are calculated by considering all the lava flows, and the composite lava field, of the second cycle