Analysis of Thermal Anomalies in Volcanic Areas Using Multiscale and Multitemporal Monitoring : Vulcano Island Test Case

Surface temperatures derived by 208 ASTER and L8 satellite imagery were analysed to test multiscale and multitemporal capability through available sets of thermal data to support the volcanic monitoring of Vulcano Island in Italy. The analysis of thermal historical series derived by ASTER and L8 shows that two are the main thermally active areas: La Fossa crater and the mud pool of Fangaia. In this work we aimed to assess the correlation between the satellite-retrieved temperatures with those measured during the daytime ground field campaign conducted within the same time period and, in particular cases, simultaneously. Moreover, nighttime data acquired by an airborne and field campaign were processed with the same methodology applied to satellite data for a multiscale approach verification. Historical meteorological data acquired from a weather station were also considered. Statistically significant correlations were observed between nighttime acquisitions and meteorological data. Correlations were also significant for temperature measured during the airborne campaign, while differences up to 50% with daytime acquisition during the ground field campaigns were observed. The analysis of the results suggests that within nighttime data acquisition, differences between satellite-derived temperatures and ground temperature measurements are considerably reduced; therefore nighttime data acquisition is recommended to detect thermal anomalies.


Introduction
The aim of volcanic surveillance is the interpretation of the observed data in order to highlight changes of volcano activity and to define the potential risks for human health and activities.The geochemical monitoring procedures are involved in activities of volcanic observatories and include the detection and quantification of fluids that are continuously released from hydrothermal and/or magmatic sources both during quiescence and eruptive periods [1][2][3][4][5][6].
The thermal energy released by a volcanic system is a parameter of primary importance for the monitoring of volcanic activity.During inter-eruptive periods, the thermal emissions by hydrothermal fluids, the emission of steam through the fumaroles and even the diffuse heat release from the ground are a significant part of the total energy released by the volcano.This slow and restless energy release is, in the long term, comparable to the seismic release of volcano-tectonic origin [6][7][8][9][10][11].For many years, the amount of convective energy released on the surface of the Island of Vulcano has been estimated by direct samplings and measurements of the fumaroles of La Fossa cone [7,12], and the exhaling temperatures and the steam flows in some sample points have been measured.Thanks to the long-term monitoring, the time variations of fumaroles temperatures showed a combination of periodic characters and episodic anomalies, related to both geochemical changes from the source, and external perturbation, due to the variability of the atmospheric conditions [13][14][15].Some major increases of fumaroles temperature were highlighted a few months before the increase in local seismic activity [11,13] and even months before the related deformation phases which affected the northern flank of the active cone [16].The volcanic system is very dynamic and it is affected by volcano tectonic processes which could change the areal distribution and intensity of thermal anomalies.This process requires an accurate thermal mapping of the ground, which needs to be extended beyond the area where there is major evidence of fumarole release.The InfraRed (IR) optical sensors promise good results for extended thermal mapping, once detection limits are ascertained and correlated with ground measurements.
Remote sensing allows us to observe the Earth surface at a large scale, with a defined acquisition frequency, and with a favorable cost-benefit ratio; moreover, it allows us to acquire data in areas that are difficult to reach.Current satellite missions, providing imagery in the Thermal InfraRed (TIR) spectral region at high spatial resolution, give the possibility to estimate the surface temperature and highlight the main surface changes, potentially related to the buried energy sources.
To improve the systematic use of satellite data in the monitoring procedures of Volcanic Observatories, suitable integration and validation strategies are needed in current satellite missions to provide TIR data with a spectral range suitable for detecting low thermal anomalies but with ground spatial resolution which is too coarse; conversely, direct monitoring could instead validate satellite measurements [1,6,7,[10][11][12][13][14][15][16].
The improved thermal surveillance of active volcanoes enhances hazard assessment and substantially contributes to understanding the modulation of endogenic heat released through the hydrothermal processes, throughout all the volcanic phases of activity.The systematic monitoring highlights changes of volcanic phenomenology and contributes to define the related risk for human health and activities.Several volcanoes are monitored through thermal surveillance based on a thermal camera (e.g., in Italy, Mt.Etna and Solfatara di Pozzuoli near Naples, and in USA, Hawaiian volcanoes); the information obtained by this monitoring instrument is included in the local hazard plan [17][18][19][20][21].Because of the complexity of Earth Observed (EO) data, the use of medium-high spatial resolution spaceborne remote sensed data is less common.
The use of TIR satellite sensors is a consolidated technique [22] to identify volcanic activity.Different types of spatial resolution are adopted in literature: low (i.e., Meteosat Second Generation (MSG) or the Geostationary Environmental Satellite (GOES) with 4 km ground resolution) [23,24], medium (i.e., MODIS, AVHRR, Sentinel 3 with 1 km ground resolution) [25][26][27] and high (Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER), and Landsat 8 (L8) with 100 m ground resolution) [28,29].In the first case, the low spatial resolution of satellite sensors does not allow the detection of thermal anomalies in small areas of the order of a hundred meters, even if the high frequency of repeated revisiting time allows daily monitoring.In the second case, the medium spatial resolution of sensors provides TIR measurements daily and can be used for quantitative analysis including measurement of surface temperatures [30,31], effusion rates [30][31][32][33], and heat flux measurements [31,34,35], but also in this case, the spatial resolution does not allow the detection of small thermal anomalies.In the last case, the high spatial resolution sensors are the most useful for quantitative measurements and analysis [21,28] even if their low temporal resolution revisiting time of 16 days makes them ineffective for real-time monitoring due to the transient nature of many volcanic processes.
In this work, ASTER and L8 time series are processed in order to identify thermal anomalies on the Island of Vulcano, which could potentially be related to changes in volcanic activity.In the literature, thermal anomalies are extracted from data with a spatial resolution of around 1 km (e.g., ATSR, MODIS, AVHRR) [36][37][38].In this work, we use data with a spatial resolution of 90-100 m.We consider this resolution comparable with the extension of the crater area and therefore suitable to retrieve sufficient variation of the surface temperature sampled in detail, which seems comparable to the ground measurements performed along the volcanic area.

Study Area
Vulcano Island is an active volcanic system and represents the southernmost sector of the Aeolian archipelago located north of Sicily, along the southern margin of the Tyrrhenian Sea (Figure 1).The Aeolian archipelago includes seven islands and is an arc-shaped structure extending from the Island of Alicudi to Stromboli Island.Vulcano, Salina and Lipari Islands form a NNW-SSE elongated volcanic belt that cross-cuts the archipelago in its central sector along the Tindari-Letojanni strike-slip tectonic system [39][40][41][42][43][44][45][46].
Remote Sens. 2018; 10; x FOR PEER REVIEW 3 of 32 ATSR, MODIS, AVHRR) [36][37][38].In this work, we use data with a spatial resolution of 90-100 m.We consider this resolution comparable with the extension of the crater area and therefore suitable to retrieve sufficient variation of the surface temperature sampled in detail, which seems comparable to the ground measurements performed along the volcanic area.

Study Area
Vulcano Island is an active volcanic system and represents the southernmost sector of the Aeolian archipelago located north of Sicily, along the southern margin of the Tyrrhenian Sea (Figure 1).The Aeolian archipelago includes seven islands and is an arc-shaped structure extending from the Island of Alicudi to Stromboli Island.Vulcano, Salina and Lipari Islands form a NNW-SSE elongated volcanic belt that cross-cuts the archipelago in its central sector along the Tindari-Letojanni strike-slip tectonic system [39][40][41][42][43][44][45][46].Vulcano Island is actually in a quiescent stage characterized by intense and variable solphataric activity.The existing monitoring network has been set up over the last 30 years and new tools are still improving.Vulcano is composed of four main volcanic structures which become gradually younger moving toward the north and NW [48]: the south Vulcano center, the Lentia volcanic complex, Vulcanello and the La Fossa cone which is the most recent structure [49].The eruptive history of Vulcano (from 127 ka to present) has been well described by De Astis et al. [49] and its eruptive history can be summarized in eight episodes, between 127 ka and historical times (AD 1888-1890) that show a long and complex eruptive history.The eruptive episodes are followed by variable periods of quiescence associated with volcano-tectonic collapse that gave origin to the calderas of Il Piano and La Fossa [49].Each eruptive episode shows different styles ranging from lava flow to pyroclastic deposits [49,50] and also different chemical composition varying from shoshonitic to high-K calcalkaline and potassic series.The activity of La Fossa shows explosive phreatic and phreatomagmatic eruptions producing wet and dry pyroclastic surges, pumice fall deposits, and highly viscous lava flows [51].
Two caldera-type structures occur in the central sector of Vulcano and their formation can be related to events producing large volumes of ash-flow deposits and/or relatively moderate-sized explosive eruptions [41].
The structural pattern of Vulcano is dominated by the major NW-SE-to NNW-SSE-striking fault system [52,41,42,40,45], which is the superficial expression of the Tindari-Letojanni system [53,44].Vulcano Island is actually in a quiescent stage characterized by intense and variable solphataric activity.The existing monitoring network has been set up over the last 30 years and new tools are still improving.Vulcano is composed of four main volcanic structures which become gradually younger moving toward the north and NW [48]: the south Vulcano center, the Lentia volcanic complex, Vulcanello and the La Fossa cone which is the most recent structure [49].The eruptive history of Vulcano (from 127 ka to present) has been well described by De Astis et al. [49] and its eruptive history can be summarized in eight episodes, between 127 ka and historical times (AD 1888-1890) that show a long and complex eruptive history.The eruptive episodes are followed by variable periods of quiescence associated with volcano-tectonic collapse that gave origin to the calderas of Il Piano and La Fossa [49].Each eruptive episode shows different styles ranging from lava flow to pyroclastic deposits [49,50] and also different chemical composition varying from shoshonitic to high-K calcalkaline and potassic series.The activity of La Fossa shows explosive phreatic and phreatomagmatic eruptions producing wet and dry pyroclastic surges, pumice fall deposits, and highly viscous lava flows [51].
Two caldera-type structures occur in the central sector of Vulcano and their formation can be related to events producing large volumes of ash-flow deposits and/or relatively moderate-sized explosive eruptions [41].
The Piano and La Fossa calderas are the most important morpho-tectonic depressions of Vulcano Island and are bordered by NNW-SSE-, N-S-and NE-SW-trending lineaments, interpreted as a pull-apart basin with minor volcano-tectonic or volcano-related components [40,45].
Our study focuses on the area of the mud pool of "Fangaia", near Faraglione area, and La Fossa crater (Figures 1 and 2), which were formed starting from the latest eruptive episode (<8 ka) [49].The Piano and La Fossa calderas are the most important morpho-tectonic depressions of Vulcano Island and are bordered by NNW-SSE-, N-S-and NE-SW-trending lineaments, interpreted as a pull-apart basin with minor volcano-tectonic or volcano-related components [4042,4547].
Our study focuses on the area of the mud pool of "Fangaia", near Faraglione area, and La Fossa crater (Figure 1, Figure 2), which were formed starting from the latest eruptive episode (<8 ka) [4951].La Fossa is a small stratovolcano that lies in the central part of the island; it reaches an elevation of 391 m a.s.l. and it has been active since 5.5 ka [5557].La Fossa was formed by phreatomagmatic eruptions that generated dilute pyroclastic density currents [5658] and references therein.
The thermal aquifers and the widespread CO2 soil degassing in these two areas testify to the actual volcanic activity [5759, 5860--6062].
Intense fumarolic degassing persists in La Fossa crater, and significant variations observed in the fumarole vents temperature at La Fossa crater have beenwere documented .In 1926, fumarole temperatures reached more than 600 °C and the highest value, 670 °C, was measured in 1993 [6163].Until March 1995, La Fossa temperatures were higher than 500 ºC, while fumaroles located outside the crater vents ranged between 350 °C and 500 °C.After 1995, the temperature of fumaroles started to decrease while the surface of active vents increased, and many new vents opened at the edge of the active cone [1212,6264].In particular, La Fossa crater is characterized by thermal anomalies that, during the seismic activity which occurred from 2004 to 2006 [6365, 6466], Faraglione is located near Levante harbour and belongs to the Vulcanello islet.It is geologically characterized by brecciated pyroclastic material, deeply altered by the strong hydrothermal release and locally hardened (Faraglione formation) with a leucite-bearing shoshonite composition [54], ascribed to the early Vulcanello activity [49].
La Fossa is a small stratovolcano that lies in the central part of the island; it reaches an elevation of 391 m a.s.l. and it has been active since 5.5 ka [55].La Fossa was formed by phreatomagmatic eruptions that generated dilute pyroclastic density currents [56] and references therein.
The thermal aquifers and the widespread CO 2 soil degassing in these two areas testify to the actual volcanic activity [57][58][59][60].
Intense fumarolic degassing persists in La Fossa crater, and significant variations observed in the fumarole vents temperature at La Fossa crater were documented [13][14][15].In 1926, fumarole temperatures reached more than 600 • C and the highest value, 670 • C, was measured in 1993 [61].Until March 1995, La Fossa temperatures were higher than 500 ºC, while fumaroles located outside the crater vents ranged between 350 • C and 500 • C.After 1995, the temperature of fumaroles started to decrease while the surface of active vents increased, and many new vents opened at the edge of the active cone [12,62].In particular, La Fossa crater is characterized by thermal anomalies that, during the seismic activity which occurred from 2004 to 2006 [63,64], were marked by strong increases of the temperature of the fumaroles and an increase of the size of the surface area covered by the fumarolic field.
In the summer of 2010, the fumaroles temperature reached its minimum value of 260 • C.

Satellite Data
In this study, two satellite datasets were used: L8 with the on-board Thermal Infrared Sensor (TIRS), launched on February 2013 [29], and the ASTER sensor launched in December 1999 on the Terra satellite [65,66].The TIR data time series were collected excluding cloudy scenes.The time series analyzed are radiance calibrated and geometrically corrected images (level 1T) of ASTER five TIR bands covering the period from 2000 to 2017 and level 1T L8 TIRS data covering the period from 2013 to 2017.In order to have the same spatial resolution for both data sets, the 100 m ground resolution TIR L8 data were re-sampled to correspond to the ASTER TIR spatial resolution of 90 m using the nearest neighbor resampling method.
In total, 208 nighttime images (170 images of ASTER data acquired from 2000 to 2017 and 38 images of L8 data acquired from 2013 to 2017) were used (Figure 3) to perform long-term surface thermal analysis which avoids the solar reflected radiance contribution, and allows comparison with the thermal camera and surface temperature collected by meteorological stations.Four daytime ASTER and L8 collected data sets were used to compare the results of a discontinuous ground campaign aimed to analyse the fumarole fields.In this second case, the solar insulation effect has not been separated in the retrieved surface temperature.ASTER and L8 data were atmospherically corrected.Then, the processor that produces the temperature map was applied.
Remote Sens. 2018; 10; x FOR PEER REVIEW 5 of 32 In the summer of 2010, the fumaroles temperature reached its minimum value of 260 °C .

Satellite Data
In this study, two satellite datasets were used: L8 with the on-board Thermal Infrared Sensor (TIRS), launched on February 2013 [29], and the ASTER sensor launched in December 1999 on the Terra satellite [65,66].The TIR data time series were collected excluding cloudy scenes.The time series analyzed are radiance calibrated and geometrically corrected images (level 1T) of ASTER five TIR bands covering the period from 2000 to 2017 and level 1T L8 TIRS data covering the period from 2013 to 2017.In order to have the same spatial resolution for both data sets, the 100 m ground resolution TIR L8 data were re-sampled to correspond to the ASTER TIR spatial resolution of 90 m using the nearest neighbor resampling method.
In total, 208 nighttime images (170 images of ASTER data acquired from 2000 to 2017 and 38 images of L8 data acquired from 2013 to 2017) were used (Figure 3) to perform long-term surface thermal analysis which avoids the solar reflected radiance contribution, and allows comparison with the thermal camera and surface temperature collected by meteorological stations.Four daytime ASTER and L8 collected data sets were used to compare the results of a discontinuous ground campaign aimed to analyse the fumarole fields.In this second case, the solar insulation effect has not been separated in the retrieved surface temperature.ASTER and L8 data were atmospherically corrected.Then, the processor that produces the temperature map was applied.Since the L8 TIRS data are delivered in quantized and calibrated scaled Digital Numbers (DN) representing all the multispectral image, the B10 data were automatically extracted and converted to Top Of Atmosphere (TOA) radiance by using the following equation: where: LTOA is the Radiance at the top of the Atmosphere W•m-2 •sr−1 •µm−1; Since the L8 TIRS data are delivered in quantized and calibrated scaled Digital Numbers (DN) representing all the multispectral image, the B 10 data were automatically extracted and converted to Top Of Atmosphere (TOA) radiance by using the following equation: where: L TOA is the Radiance at the top of the Atmosphere W•m −2 •sr −1 •µm −1 ; M L is the multiplicative rescaling factor (3.342 × 10 −4 ) for L8 B 10 metadata, from the MTL file furnished with the image data; Q cal is the quantized and calibrated standard digital pixel number (DN) value; A L is the additive rescaling factor (0.1) for L8 B 10 metadata, from the MTL file.
Starting from the B 10 TOA spectral radiance (L TOA ), the influence of the atmosphere has been removed using the "CIRILLO" module [69] and the radiance emitted by the surface at specific wavelength (L λ ) has been obtained.CIRILLO is homemade software for the correction of spaceborne images acquired in the VIS-TIR spectral range.The atmospheric terms needed by CIRILLO are computed by the known radiative transfer models Modtran4.0[70,71] and 6S [72].Moreover, CIRILLO takes into account the illumination changes due to relief computing and corrects them starting from digital terrain models and solar position.
Using these data, the spectral radiance (L λ , W•m −2 •sr −1 •µm −1 ) emitted by surface at temperature T(K) could be calculated by rearranging the following equation from [73]: L TOA the radiance received by B 10 of the sensor with brightness temperature T 10 ; L λ the radiance from the surface at B 10 ; L u the upwelling atmospheric radiance components (obtained using Modtran [70]); L d the downwelling atmospheric radiance components (obtained using Modtran [70]); ε the surface emissivity for B 10 extracted by the ASTER 05 emissivity data [74]; τ the atmospheric transmittance for B 10 when the view zenith angle is θ (obtained using Modtran [70]).TIRS is treated as nadir viewing since the view angle is no more than 7.5 • .
Atmospheric profiles in terms of pressure, temperature, geopotential height, relative humidity needed for obtaining the atmospheric components in (2) were downloaded from the University of Wyoming (WYO) website [75].
Since the focus of the investigation is the estimation of a surface temperature map, in order to invert equation (2), we have assumed the emissivity value provided by the ASTER05 reference map.The surface temperature has been obtained by applying the formula reported in [76]: where T is the surface temperature ( • K); L λ is the surface spectral radiance; K 1 is the Band-specific thermal conversion constant provided in the MTL file; K 2 is the Band-specific thermal conversion constant provided in the MTL file; Once we obtained the temperature, the conversion to degrees Celsius was done in order to better compare with the measurements collected during the field campaigns.

Processing of ASTER Data
ASTER covers a wide spectral region with 14 bands, distributed along three sensors; the Thermal Infrared (TIR) sensor with 5 bands in the spectral range from 8.125 to 11.65 µm at spatial resolution of 90 m has been used in this work.
The standard TIR ASTER L1T products were used to estimate the surface temperature from the ASTER data using the Temperature Emissivity Separation (TES) algorithm [77].The atmospheric correction was performed using "CIRILLO" atmospheric correction tools [69] as described in the previous section.

Direct Measurements from Continuous Monitoring Stations
Vulcano is constantly monitored by a ground network of instruments for continuous acquisition of fumaroles present on top of the crater.Sensors position equipped with thermometers are located on the main active fumaroles, to track changes from the high temperature of fumaroles (HTF), shown in Figure 4 (modified from [13]).

Remote Sens. 2018; 10; x FOR PEER REVIEW 7 of 32
Vulcano is constantly monitored by a ground network of instruments for continuous acquisition of fumaroles present on top of the crater.Sensors position equipped with thermometers are located on the main active fumaroles, to track changes from the high temperature of fumaroles (HTF), shown in Figure 4 (modified from [13]).The monitoring network components are thermal probes, inserted within the fumarole vents at a depth of 50 cm, consisting of thermocouples (K type).The buried part of the thermo-couple wire is protected by a cylindrical case, made by stainless steel, while a head made of teflon protects a part of the electronic components of the system from the direct exposure to acidic fluids [13][14][15].
Many other temperature measurements were occasionally carried out on the La Fossa crater with a portable low-cost radiometer assembled into the Pyrogips system [78,79].This instrument supplies the advantage of in situ methods (i.e.accuracy of direct measurement by a thermocouple) producing a detailed grid of surface temperatures that can be very useful for establishing a network of sampling points or installing a new site for geochemical monitoring.The string of data returns the GPS position, surface temperature (derived from the pyrometer and thermocouple), ambient temperature, barometric pressure and air moisture.The surface temperature acquired during four field campaigns (2012, two in 2014 and 2017) was used to compare them with the temperatures estimated by satellite data acquired, considering dates and time as close as possible to the ground survey.In Table 1 the dates for each campaign and available satellite data are reported.Since the surveys were conducted during the morning, daytime satellite acquisitions (about 09:30 UTC) were processed by using the same methodology described in § 2.1.and its subsections..The monitoring network components are thermal probes, inserted within the fumarole vents at a depth of 50 cm, consisting of thermocouples (K type).The buried part of the thermo-couple wire is protected by a cylindrical case, made by stainless steel, while a head made of teflon protects a part of the electronic components of the system from the direct exposure to acidic fluids [13][14][15].
Many other temperature measurements were occasionally carried out on the La Fossa crater with a portable low-cost radiometer assembled into the Pyrogips system [78,79].This instrument supplies the advantage of in situ methods (i.e., accuracy of direct measurement by a thermocouple) producing a detailed grid of surface temperatures that can be very useful for establishing a network of sampling points or installing a new site for geochemical monitoring.The string of data returns the GPS position, surface temperature (derived from the pyrometer and thermocouple), ambient temperature, barometric pressure and air moisture.The surface temperature acquired during four field campaigns (2012, two in 2014 and 2017) was used to compare them with the temperatures estimated by satellite data acquired, considering dates and time as close as possible to the ground survey.In Table 1 the dates for each campaign and available satellite data are reported.Since the surveys were conducted during the morning, daytime satellite acquisitions (about 09:30 UTC) were processed by using the same methodology described in Section 2.1 and its subsections.However, the cost to benefit ratio of the in situ measurements cannot be suggested for the geochemical surveillance of Vulcano since it is less convenient in comparison to the growing technical facilities offered by thermal remote sensing both proximal (by drone or by airplane) and by satellite [7].

Proximal Thermocamera Imaging with Hand-Held IR Camera
The infrared thermography is based on radiation heat transfer which is an energy transport mechanism that occurs under the form of electromagnetic waves.During this transport, radiant energy may partially be absorbed and reflected by a body.Consequently, the intensity of radiation received at the optical sensors depends on several parameters such as emissivity, viewing angle, atmospheric effects instrumental noise and, particularly for volcanic targets, it depends on the thermal heterogeneity of the target [18,62,[80][81][82][83][84][85].
Hand-held IR camera technology has considerably improved over the last decade, resulting in a proliferation of commercially available instruments, which allow volcanic observatories to achieve indirect thermal monitoring.
The camera features are reported in Table 2.
Table 2. Hand-held IR camera features.

IR Camera Features
Temperature range 20 The camera was calibrated before and after the survey held on 16 July 2014, with a variable temperature target.This initial pre-processing of the TIR imagery required the conversion of at-sensor-radiance values to surface temperature using a first-order correction for atmospheric absorption, surface emissivity and ambient temperature effects using ResearchIR Software.This correction includes the distance from the target to the camera, the atmosphere temperature and the humidity content, all measured in situ.An emissivity value of 0.95 was used for the first-order correction that is an appropriate value for this kind of volcanic rock, i.e., basalt with a smooth surface ranged from 0.971 to 0.975 [84].
Although the cameras are intrinsically temperature-calibrated, numerous extrinsic field conditions can compromise the correct temperature measurements [84][85][86][87].The IR camera usually measures an apparent temperature for several reasons lower than the actual surface temperature.If there is much condensed water in the plume, the target surface may be entirely or partially blocked by fume, attenuating the radiant signal of the surface significantly [85].The pixel footprint size also plays a part in the temperature accuracy-and effective sensitivity-of the cameras [84] because an individual pixel temperature is a composite of all the temperatures within the pixel footprint.If the hot target is small and the majority of the pixel area is cooler, this "mixed pixel" issue results in the pixel-integrated temperature being significantly lower than that of the hot target.

Meteorological Data
Data from the weather stations located near the study areas were compared to the different time series of ground temperatures to verify if the trends of ground temperatures observed and both direct and indirect measurements are related to external changes (Section 2.2.1).
Prior to interpreting the observed temperature ground variations, due to the complexity of volcanic and hydrothermal system, it is useful to discriminate the different influences of the external components as they affect primarily the air temperature, but they also influence the ground temperatures.
In previous studies, the rainfall, barometric pressures, and normal air temperatures have been obtained by the Servizio Informativo Agrometeorologico Siciliano (SIAS, Regione Siciliana), that is the regional network of weather stations [14,15].The two weather stations in Figure 5 were considered, because they are located both at a safe distance from the hydrothermal influence and the local volcanic activity, while their locations reflect the normal weather conditions of the region.Data from the weather stations located near the study areas were compared to the different time series of ground temperatures to verify if the trends of ground temperatures observed and both direct and indirect measurements are related to external changes ( §2.2.1).
Prior to interpreting the observed temperature ground variations, due to the complexity of volcanic and hydrothermal system, it is useful to discriminate the different influences of the external components as they affect primarily the air temperature, but they also influence the ground temperatures.
In previous studies, the rainfall, barometric pressures, and normal air temperatures have been obtained by the Servizio Informativo Agrometeorologico Siciliano (SIAS, Regione Siciliana), that is the regional network of weather stations [14,15].The two weather stations in Figure 5 were considered, because they are located both at a safe distance from the hydrothermal influence and the local volcanic activity, while their locations reflect the normal weather conditions of the region.The historical weather data were acquired from a weather station near Vulcano harbor, under public access provided by the Oplao web site [88].The web site offers a data archive of weather information from July 2008 until today.These temperatures were used to compare the trend of surface temperature estimated by satellite data with the local climate, as validation measurements and for representative of the soil-atmosphere interface temperature.When the satellite data were available during the 2008-2017 period, the temperature, reported at 21:00 UTC was compared daily with the results of the satellite-processed data, by a paired dataset.
This comparison is aimed to verify if the background temperature trends, derived by satellite image processing, and validation measurements provided by the weather station near Vulcano harbor, have similar behavior, fostering the capability to extend local measurements to a wider area.

Airborne Data
In order to introduce an intermediate observational scale and instrumental resolution between the direct measurements and satellite data, multispectral airborne data were acquired and elaborated.The measurement campaign took place in the period 10-12 September 2010 [89].This campaign was characterized by measurements made by the Airborne Thematic Mapper (ATM) multispectral plane sensor, and contemporaneous on-site validation measurements made by ground instrumentation on both marine and terrestrial surfaces.The airborne remote sensing mission was carried out by the crew of a Coast Guard aircraft, an ATR 42 aircraft, with a SENSYTECH AADS 1268 ATM-E scanning system installed on board.This sensor consists of a passive scanner for aerial remote sensing that records the reflected and emitted terrestrial radiation, divided into 11 spectral The historical weather data were acquired from a weather station near Vulcano harbor, under public access provided by the Oplao web site [88].The web site offers a data archive of weather information from July 2008 until today.These temperatures were used to compare the trend of surface temperature estimated by satellite data with the local climate, as validation measurements and for representative of the soil-atmosphere interface temperature.When the satellite data were available during the 2008-2017 period, the temperature, reported at 21:00 UTC was compared daily with the results of the satellite-processed data, by a paired dataset.
This comparison is aimed to verify if the background temperature trends, derived by satellite image processing, and validation measurements provided by the weather station near Vulcano harbor, have similar behavior, fostering the capability to extend local measurements to a wider area.

Airborne Data
In order to introduce an intermediate observational scale and instrumental resolution between the direct measurements and satellite data, multispectral airborne data were acquired and elaborated.The measurement campaign took place in the period 10-12 September 2010 [89].This campaign was characterized by measurements made by the Airborne Thematic Mapper (ATM) multispectral plane sensor, and contemporaneous on-site validation measurements made by ground instrumentation on both marine and terrestrial surfaces.The airborne remote sensing mission was carried out by the crew of a Coast Guard aircraft, an ATR 42 aircraft, with a SENSYTECH AADS 1268 ATM-E scanning system installed on board.This sensor consists of a passive scanner for aerial remote sensing that records the reflected and emitted terrestrial radiation, divided into 11 spectral bands ranging from Visible to Thermal Infrared.The flight on Vulcano Island was carried out according the acquisition characteristics reported in Table 3.

Statistical Indices
The analysis of the mean of surface temperature retrieved by satellite data in an area of about 0,13 Sq.km for the La Fossa crater and about 0.1 Sq.km near the "Fangaia" mud pool, and the temperature measured by the meteorological station located in the Vulcano harbor offers a very interesting cross comparison and allows us to validate the methodology used for the satellite data surface temperature.In order to emphasize the agreement among the analysed datasets, we choose to apply the following statistical indices: Pearson correlation (Equation ( 4)), the Normalized Standard Error (NSE, Equation ( 5)) and the Normalized Bias (NB, Equation ( 6)).These indices give quantitative measurements of the relationship between the two measurements and demonstrate a high degree of correlation.
where N is the total number of processed data points, Si and Mi are the satellite surface temperature and meteorological station temperature of the i-th measurements, and S and M are the average of the mean of Si and Mi, respectively.

Surface Temperature Retrieved by Satellite Data and Weather Station Time Series
Here we focus on the time series obtained by analysing 170 nighttime images of ASTER data acquired from 2000 to 2017 and 38 nighttime images of L8 data acquired from 2013 to 2017.The meteorological station time series were collected considering the same day and time of satellite acquisition in order to compare the information.Two areas of Vulcano Island were considered: La Fossa crater and the mud pool "Fangaia" areas.These sites have been selected because they have an apparent temperature higher than the surrounding area.In Figures 6 and 7, the scatter plots of the surface temperature retrieved by the satellite and those registered by the meteo station for the two test areas and the three correlation indices are displayed.Concerning the ASTER time series, we can reasonably state that the strong agreement between the meteorological observation and satellite data starting from 2008 can validate the methodology for the evaluation of temperature by using satellites.Moreover, even if no meteorological data are available before July 2008, the agreement showed in Figure 6 and Figure 7 allows us to extend the validity of the ASTER time series for the period ranging from 2000 to 2008.As L8 acquisition started in 2013, it is not possible to extend this analysis backward.
During a field campaign on 3 September 2015, surface temperature was measured on the top of the La Fossa crater by using a thermometer installed in a multisensor instrument.The measured value [90] is coherent with the Oplao web site value and with the temperature estimated by L8 data acquired only two days later (Table 4).This comparison (even if related to only one day) can be Concerning the ASTER time series, we can reasonably state that the strong agreement between the meteorological observation and satellite data starting from 2008 can validate the methodology for the evaluation of temperature by using satellites.Moreover, even if no meteorological data are available before July 2008, the agreement showed in Figures 6 and 7 allows us to extend the validity of the ASTER time series for the period ranging from 2000 to 2008.As L8 acquisition started in 2013, it is not possible to extend this analysis backward.
During a field campaign on 3 September 2015, surface temperature was measured on the top of the La Fossa crater by using a thermometer installed in a multisensor instrument.The measured value [90] is coherent with the Oplao web site value and with the temperature estimated by L8 data acquired only two days later (Table 4).This comparison (even if related to only one day) can be considered as further confirmation of the methodology used to estimate the temperature from satellites.The temperature time series, plotted in Figure 8, are the result of the long-term monitoring of 4 fumarole vents, which in the last 30 years have shown temperature variations ranging on the whole, from 670 • C to 120 • C. Since September 1991, 24 temperatures a day have been recorded and validated for the sites selected for volcanic surveillance.The output temperatures of the main fumaroles located on the top of La Fossa cone have not been stationary, a few changes in trends and many impulsive variations affected all the monitored sites [13].The temperature time series, plotted in Figure 8, are the result of the long-term monitoring of 4 fumarole vents, which in the last 30 years have shown temperature variations ranging on the whole, from 670 °C to 120 °C.Since September 1991, 24 temperatures a day have been recorded and validated for the sites selected for volcanic surveillance.The output temperatures of the main fumaroles located on the top of La Fossa cone have not been stationary, a few changes in trends and many impulsive variations affected all the monitored sites [13].The linear trends of emission temperatures recorded by the monitoring network are highlighted in Figure 8, and their best-fitting equations have been reported in the legend.The monitored fumaroles located on the rim zone (HTF and New HTF in Figure 4) have shown two positive trends for more than ten years, until the end of 2013 (Figure 8).At least from the year 2004, and 4 years later, the fumaroles located on the inner slope have shown the same time variations as the fumarole on the rim zone (compare the time variations of HTF and Inner Slope in Figure 8).Then, the trend of temperatures on the inner slope (Inner Slope in Figure 4) dramatically diverged from the trend of the rim zone, during the year 2007.The beginning of 2007 was characterized by strong decreases recorded in all the monitored fumaroles but, while the upper rim recovered the positive trending variations of temperature, the inner slope continued to show a decreasing trend.Afterwards, until The linear trends of emission temperatures recorded by the monitoring network are highlighted in Figure 8, and their best-fitting equations have been reported in the legend.The monitored fumaroles located on the rim zone (HTF and New HTF in Figure 4) have shown two positive trends for more than ten years, until the end of 2013 (Figure 8).At least from the year 2004, and 4 years later, the fumaroles located on the inner slope have shown the same time variations as the fumarole on the rim zone (compare the time variations of HTF and Inner Slope in Figure 8).Then, the trend of temperatures on the inner slope (Inner Slope in Figure 4) dramatically diverged from the trend of the rim zone, during the year 2007.The beginning of 2007 was characterized by strong decreases recorded in all the monitored fumaroles but, while the upper rim recovered the positive trending variations of temperature, the inner slope continued to show a decreasing trend.Afterwards, until the end of 2017, the main features shown by the fumarole temperature of the inner slope have been the strongest perturbations of external origin, due to the seasonal changes of weather conditions.
On the active cone of La Fossa [91], the experimental acquisition of ground temperatures was performed at 4 different sites during three years, from 2004 to 2007.The measurements of ground temperature, made by contact probes, showed that three of these sites were characterized by anomalous profiles of temperature, suggesting a heating source located below.These stations highlight a process of steam heating within the ground and a thermal source, and they allow us to evaluate the heat flux from the soil.Under dry conditions, that is in the absence of steam or water within the monitored profile, the conductive heat flux is the only form of heat transfer, thus the temperature of the ground is the only detectable variable in order to evaluate the heat flux coming from the soil heated by a steam source located at the bottom of the monitoring profile.

Proximal Yhermocamera Imaging with a Hand-Held IR Camera
The geothermal flux and volcanic gas activity generate positive anomalies mainly at the surface, located along the rim and in the inner flank of La Fossa cone.A simplified, steady-state thermal budget model was used to map the thermal anomalies of this sector of Vulcano.The ground measurements show well-localized and thermally active zones.The thermal infrared images present a visualization of the heat flux associated with these structures.
Generally, two physical characteristics of the exposed surfaces are the main causes of the space variability of temperatures: the thermal capacity of the exposed surfaces and the shallow permeability of the ground.
The surface of the cone of La Fossa on the top is not homogeneous.It is covered by later La Fossa Tuffs, pyroclastic breccias, including large bombs.Sand tuffs and breccias show interrelation of ashes, some with accretionary lapilli, although subordinate, authigenic-juvenile material erupted in the form of thin pumice layers and some larger scoriae bombs.Moreover there is a diffuse interbedding of argillitic alteration residual and fumaroles deposits.The enlargement of the exhaling area and/or an increase of thermal anomaly surrounding the main fumarole vents (due to steam heating from the bottom source) can highlight significant increases of thermal release even if the maximum temperature of fumarole fluids falls.Responding to thermodynamic basic principles, the exchange of energy drives each component towards the thermal equilibrium.
During the survey on 16 July 2014, soil temperatures ranged from 20 • C to 114 • C on the fumarole field at a distance of 300 at the bottom to 400 m on the top of the rim. Figure 9 shows the thermographic map of the main diffusive and thermal manifestations after optical and atmospheric corrections required; a mixed area with diffuse heat flux (yellow) and an area with high temperature vents (red) are widespread in a region with low temperature (blue).
Remote Sens. 2018; 10; x FOR PEER REVIEW 13 of 32 the end of 2017, the main features shown by the fumarole temperature of the inner slope have been the strongest perturbations of external origin, due to the seasonal changes of weather conditions.On the active cone of La Fossa [91], the experimental acquisition of ground temperatures was performed at 4 different sites during three years, from 2004 to 2007.The measurements of ground temperature, made by contact probes, showed that three of these sites were characterized by anomalous profiles of temperature, suggesting a heating source located below.These stations highlight a process of steam heating within the ground and a thermal source, and they allow us to evaluate the heat flux from the soil.Under dry conditions, that is in the absence of steam or water within the monitored profile, the conductive heat flux is the only form of heat transfer, thus the temperature of the ground is the only detectable variable in order to evaluate the heat flux coming from the soil heated by a steam source located at the bottom of the monitoring profile.

Proximal Yhermocamera Imaging with a Hand-Held IR Camera
The geothermal flux and volcanic gas activity generate positive anomalies mainly at the surface, located along the rim and in the inner flank of La Fossa cone.A simplified, steady-state thermal budget model was used to map the thermal anomalies of this sector of Vulcano.The ground measurements show well-localized and thermally active zones.The thermal infrared images present a visualization of the heat flux associated with these structures.
Generally, two physical characteristics of the exposed surfaces are the main causes of the space variability of temperatures: the thermal capacity of the exposed surfaces and the shallow permeability of the ground.
The surface of the cone of La Fossa on the top is not homogeneous.It is covered by later La Fossa Tuffs, pyroclastic breccias, including large bombs.Sand tuffs and breccias show interrelation of ashes, some with accretionary lapilli, although subordinate, authigenic-juvenile material erupted in the form of thin pumice layers and some larger scoriae bombs.Moreover there is a diffuse interbedding of argillitic alteration residual and fumaroles deposits.The enlargement of the exhaling area and/or an increase of thermal anomaly surrounding the main fumarole vents (due to steam heating from the bottom source) can highlight significant increases of thermal release even if the maximum temperature of fumarole fluids falls.Responding to thermodynamic basic principles, the exchange of energy drives each component towards the thermal equilibrium.
During the survey on 16 July 2014, soil temperatures ranged from 20°C to 114°C on the fumarole field at a distance of 300 at the bottom to 400 meters on the top of the rim. Figure 9 shows the thermographic map of the main diffusive and thermal manifestations after optical and atmospheric corrections required; a mixed area with diffuse heat flux (yellow) and an area with high temperature vents (red) are widespread in a region with low temperature (blue).The statistical analysis applied to sequences of the thermal images analysed aims to define the extension of thermal anomalies, contributing to the evaluation of surface flux.The Normal quantile-quantile plot (QQ plot, Figure 10) is a graphical tool to verify if a set of data came from The statistical analysis applied to sequences of the thermal images analysed aims to define the extension of thermal anomalies, contributing to the evaluation of surface flux.The Normal quantile-quantile plot (QQ plot, Figure 10) is a graphical tool to verify if a set of data came from some theoretical distribution, such as a normal or exponential one.The data values are ordered and the cumulative distribution values are calculated as (i-0.5)/nfor the i-th ordered value out of n total values.The scatterplot is created by plotting the real set of quantiles against the normal distributed population of the same quantiles; a line that is roughly straight corresponds to a normal distribution of the real dataset.Any diversion from the straight line indicates a diversion of the dataset from the normal distribution.The probability plots reveal different populations of data.According to these populations the space variability observed on the surface temperatures has been inferred from the main components of the heat transport (radiative, convective conductive) which participated in the heat exchange processes according to the different rate of flux reaching the ground surface.
The temperature distribution on a probability plot, as shown in Figure 10, was the visual method used for discriminating the different groups, and also for choosing the scale of values for the temperature thematic maps.We interpreted the lowest class of ground temperature (<20 The temperature distribution on a probability plot, as shown in Figure 10, was the visual method used for discriminating the different groups, and also for choosing the scale of values for the temperature thematic maps.We interpreted the lowest class of ground temperature (<20 °C) as falling in the water saturated layer; the second class of values (included in the range from 20 to 27 °C) was interpreted as the surface condition dominated by conductive heat transfer, absence of liquid water and very low steam content; the third class (temperatures in the range from 28 to 48 °C) was for the steam heated soil; finally, the hottest class of values included the low temperature fumaroles ( T > 48), where the convective transfer dominates.
Figure 10.qq-plot of the infrared temperature measurements (La Fossa crater rim).The red line identifies the thermal quantile groups of dataset and the range temperature; the black line is the reference line to assess the Normal Gaussian distribution assumption).

Airborne Campaign
The airborne survey carrying the multi-spectral infrared and visible imaging spectrometer (ATM) was carried out with the primary objective of obtaining high-spatial-resolution thermal images of the Vulcano Island.The predawn airborne thermal data of channel 10.5 µ m were examined in order to retrieve sea and land surface temperatures, with particular attention to fumarole areas emitting at higher temperatures (about 50 °C).Thermal data were corrected for atmospheric and reflection effects by using equation (2).Parameters of the equation were obtained using the Modtran model atmosphere [70,92].Using a sea-land mask we successively identify the thermal anomalies within the area (Figure 11).The striping effect showed in Figure 11 is an artifact due to the generation of the image and it does not influence the quality and results of the applied procedure.Sea temperature retrieved from airborne data was compared with data obtained from contemporaneous ground field measurements using a portable thermo radiometer and handheld thermal camera.

Airborne Campaign
The airborne survey carrying the multi-spectral infrared and visible imaging spectrometer (ATM) was carried out with the primary objective of obtaining high-spatial-resolution thermal images of the Vulcano Island.The predawn airborne thermal data of channel 10.5 µm were examined in order to retrieve sea and land surface temperatures, with particular attention to fumarole areas emitting at higher temperatures (about 50 • C).Thermal data were corrected for atmospheric and reflection effects by using equation (2).Parameters of the equation were obtained using the Modtran model atmosphere [70,92].Using a sea-land mask we successively identify the thermal anomalies within the area (Figure 11).The striping effect showed in Figure 11 is an artifact due to the generation of the image and it does not influence the quality and results of the applied procedure.Sea temperature retrieved from airborne data was compared with data obtained from contemporaneous ground field measurements using a portable thermo radiometer and handheld thermal camera.Gas emissions are observed both in the crater and in the beach area (Fumarole beach) at the northern foot of the La Fossa crater; vents also occur at different sites off the island.The presence of submarine gaseous emissions along the coast of "Fangaia" area is evidenced.
In Figure 12 the temperature profile derived from the map thermal distribution of Figure 11 is shown.The water temperature profile shows a regular decrease in temperature starting from the open sea to the coast.This is due to the slow night sea cooling down, since the water has higher specific heat capacity than land.Near the coast, temperatures increase suddenly and sharply due to the Gas emissions are observed both in the crater and in the beach area (Fumarole beach) at the northern foot of the La Fossa crater; vents also occur at different sites off the island.The presence of submarine gaseous emissions along the coast of "Fangaia" area is evidenced.
In Figure 12  Gas emissions are observed both in the crater and in the beach area (Fumarole beach) at the northern foot of the La Fossa crater; vents also occur at different sites off the island.The presence of submarine gaseous emissions along the coast of "Fangaia" area is evidenced.
In Figure 12 the temperature profile derived from the map thermal distribution of Figure 11 is shown.The water temperature profile shows a regular decrease in temperature starting from the open sea to the coast.This is due to the slow night sea cooling down, since the water has higher specific heat capacity than land.Near the coast, temperatures increase suddenly and sharply due to the The water temperature profile shows a regular decrease in temperature starting from the open sea to the coast.This is due to the slow night sea cooling down, since the water has higher specific heat capacity than land.Near the coast, temperatures increase suddenly and sharply due to the presence of submarine gas emissions.The average temperature rises from 20 • C in the sea to about 28 • C in the land.
In the case of the submarine fumarole area, it was possible to compare temperatures derived from remote sensing data and contemporaneous field measurements.Figure 13 shows the water surface temperature measured from thermal camera.A temperature of 21.9 • C was measured which is in good agreement with the ATM-derived temperature of 21. presence of submarine gas emissions.The average temperature rises from 20 °C in the sea to about 28 °C in the land.
In the case of the submarine fumarole area, it was possible to compare temperatures derived from remote sensing data and contemporaneous field measurements.Figure 13 shows the water surface temperature measured from thermal camera.A temperature of 21.9 °C was measured which is in good agreement with the ATM-derived temperature of 21.6 °C.

Discussion
In the next sections the comparison between multiscale and multitemporal will be described.

Satellite Data Compared with Ground Data
We have generated subsets of collected ground temperature considering those included in a pixel of 90x90 m (corresponding to L8 and ASTER pixel size).For each of these subsets the average ground temperature has been calculated and compared with surface temperature estimated by satellites (Table 5, Table 6, Table 7).Taking into account the satellite spatial resolution, 15 satellite pixels drape the ground distribution of collected points.
In Figure 14 -Figure 17, the draping of the 15 identified satellite pixels is overlapped on the path followed during the field campaign reported in Table 1.

Discussion
In the next sections the comparison between multiscale and multitemporal will be described.

Satellite Data Compared with Ground Data
We have generated subsets of collected ground temperature considering those included in a pixel of 90 × 90 m (corresponding to L8 and ASTER pixel size).For each of these subsets the average ground temperature has been calculated and compared with surface temperature estimated by satellites (Tables 5-7).Taking into account the satellite spatial resolution, 15 satellite pixels drape the ground distribution of collected points.
In Figures 14-17, the draping of the 15 identified satellite pixels is overlapped on the path followed during the field campaign reported in Table 1.
We have generated subsets of collected ground temperature considering those included in a pixel of 90x90 m (corresponding to L8 and ASTER pixel size).For each of these the average ground temperature has been calculated and compared with surface temperature estimated by satellites (Table 5, Table 6, Table 7).Taking into account the satellite spatial resolution, 15 satellite pixels drape the ground distribution of collected points.
In Figure 14 -Figure 17, the draping of the 15 identified satellite pixels is overlapped on the path followed during the field campaign reported in Table 1.A first comment can be made about the spatial resolution of the measurements: "point" for ground data and "area" for data retrieved by satellite.Furthermore, it is important to remark that the average ground estimated temperature is higher than the temperature estimated by using satellite data.The ground measurements were collected near the hottest point where fumaroles were present and, as shown in Figures 14 and 18, their contribution is spread at satellite pixel scale.Similarly, the draping of the identified satellite pixels, overlapping the path followed during the field campaign held on the 15/09/2014 survey, is reported in Figure 19.Similarly, the draping of the identified satellite pixels, overlapping the path followed during the field campaign held on the 15/09/2014 survey, is reported in Figure 19.The Comparison between ground measurements and satellite estimation temperatures in the mud pool "Fangaia", is reported in Table 8.Similarly, the draping of the identified satellite pixels, overlapping the path followed during the field campaign held on the 15/09/2014 survey, is reported in Figure 19.Table 8.Comparison between ground measurements and satellite estimation temperatures in the mud pool "Fangaia".During the field campaign on 16 July 2014, a nighttime IR map was produced by means of the Thermal Flir camera 620T (Figure 20).This map has been compared with the surface temperature map obtained by processing satellite acquisitions the day before (ASTER acquired on 15 July at 20:45 UTC).As the resampling of the satellite data results in a resolution of 90 × 90 m, a regularly-spaced 90 × 90 m grid was created for the IR map to compare the temperatures of satellite images and IR ground image.For the comparison, we have considered as representative satellite pixels the ASTER and L8 pixels, since they have the hottest temperature compared to the other pixels: these hottest pixels are 2 for ASTER and 2 for L8 and they are fully overlapping.As reported in Table 9, the measured areas have a good correlation with the ground measurements.map obtained by processing satellite acquisitions the day before (ASTER acquired on 15 July at 20:45 UTC).As the resampling of the satellite data results in a resolution of 90x90 m, a regularly-spaced 90x90 m grid was created for the IR map to compare the temperatures of satellite images and IR ground image.For the comparison, we have considered as representative satellite pixels the ASTER and L8 pixels, since they have the hottest temperature compared to the other pixels: these hottest pixels are 2 for ASTER and 2 for L8 and they are fully overlapping.As reported in Table 9, the measured areas have a good correlation with the ground measurements.In order to constrain the hottest pixels in the La Fossa crater area, an analysis based on satellite data time series shows that the northern sector of the La Fossa crater inner slope is characterized by higher temperature (Figure 20) with respect to a wide area as represented in Figure 21.
In order to constrain the hottest pixels in the La Fossa crater area, an analysis based on satellite data time series shows that the northern sector of the La Fossa crater inner slope is characterized by higher temperature (Figure 20) with respect to a wide area as represented in Figure 21.By analysing ASTER time series, we pointed out that not all the pixels sensed by satellites are characterized by a significant thermal signature.A restricted number of pixels, highlighted by satellites, is systematically warmer than nearby.We determine that pixels numbered 7, 8, 12 and 13 are those for which the temperature is higher than the surrounding pixels, i.e. 192 count over 208 (Table 10 and Figure 22).By analysing ASTER time series, we pointed out that not all the pixels sensed by satellites are characterized by a significant thermal signature.A restricted number of pixels, highlighted by satellites, is systematically warmer than nearby.We determine that pixels numbered 7, 8, 12 and 13 are those for which the temperature is higher than the surrounding pixels, i.e., 192 count over 208 (Table 10 and Figure 22).Similarly, the geographical positions of the hottest pixels derived by time series satellite data (Table 11) for the "Fangaia" mud pool were detected (Figure 23, Figure 24).Similarly, the geographical positions of the hottest pixels derived by time series satellite data (Table 11) for the "Fangaia" mud pool were detected (Figures 23 and 24).

Satellite Data Compared with Airborne Data
In this section the discussion is based on the comparison between the retrieved ASTER thermal data series and the thermal data obtained by the ATM.The night-time satellite passes considered are those more temporally closed to the ATM airborne campaign.This results in ASTER acquisitions on 6 September and 22 September 2010.The comparison has been performed co-locating ATM and ASTER processed data for ground temperature and dividing the "Fangaia" site by the La Fossa crater

Satellite Data Compared with Airborne Data
In this section the discussion is based on the comparison between the retrieved ASTER thermal data series and the thermal data obtained by the ATM.The night-time satellite passes considered are those more temporally closed to the ATM airborne campaign.This results in ASTER acquisitions on 6 September and 22 September 2010.The comparison has been performed co-locating the ATM and ASTER processed data for ground temperature and dividing the "Fangaia" site by the La Fossa crater areas.Two different comparisons were performed.In a first comparison the mean of the temperature values retrieved for the two areas has been calculated considering the whole thermal anomalies derived by ATM and ASTER ground temperatures.Figure 25 shows the ATM thermal data and the areas containing the thermal anomalies considered to calculate the temperature means.In Table 12, the comparison of the mean temperatures retrieved is shown.The error in temperature values has been assigned considering the error of the processing algorithm.
areas.Two different comparisons were performed.In a first comparison the mean of the temperature values retrieved for the two areas has been calculated considering the whole thermal anomalies derived by ATM and ASTER ground temperatures.Figure 25 shows the ATM thermal data and the areas containing the thermal anomalies considered to calculate the temperature means.In Table 12, the comparison of the mean temperatures retrieved is shown.The error in temperature values has been assigned considering the error of the processing algorithm.The comparison highlights, as expected, that ATM has higher temperature values than ASTER.In the "Fangaia" area there is quite good agreement between the ASTER acquired on 6 September and the ATM-retrieved temperature (Figure 26).In the La Fossa Crater area, the ATM-derived temperature is higher for both ASTER acquisitions.In order to better understand such difference, a second comparison has been performed.The comparison of ASTER acquired on 6 September and ATM thermal anomalies has been performed pixel by pixel for the La Fossa Crater (Figure 26).Many thermal anomalies derived from ATM cannot be discriminated at an ASTER ground sampling resolution of 90 m, which is too low to capture small fumarolic areas, which instead ATM airborne data with 3 m of ground resolution can identify.In order to make a closer comparison, the ATM mean temperature has been recalculated pixel by pixel and only the ASTER and ATM common pixels are considered.Referring to Figure 26, we have obtained an average value of 25.8 ± 1.4 • C from 6 September ASTER data and an average temperature of 27.6 ± 2.4 • C from ATM data.
temperature is higher for both ASTER acquisitions.In order to better understand such difference, a second comparison has been performed.The comparison of ASTER acquired on 6 September and ATM thermal anomalies has been performed pixel by pixel for the La Fossa Crater (Figure 26).Many thermal anomalies derived from ATM cannot be discriminated at an ASTER ground sampling resolution of 90 m, which is too low to capture small fumarolic areas, which instead ATM airborne data with 3 m of ground resolution can identify.In order to make a closer comparison, the ATM mean temperature has been recalculated pixel by pixel and only the ASTER and ATM common pixels are considered.Referring to Figure 26, we have obtained an average value of 25.8±1.4°C from 6 September ASTER data and an average temperature of 27.6±2.4°C from ATM data.

Conclusions
This study is based on the comparison among satellite data with measurements collected in different campaigns (airborne and ground).By analysing 2000-2017 spaceborne remotely sensed time series data of surface temperature covering the whole island, we have detected hot pixels and used them as focal points.Assuming that all pixels on the island are affected by the same seasonal conditions, the hottest pixels identified will be warmer independently of seasonal and daily variations.During the daytime, the reflected solar radiation is difficult to remove from thermal channels.Moreover, the error caused by the solar correction term can affect the retrieval accuracy, especially under low reflectance in thermal infrared regions.Therefore, in our analysis, the reflected solar radiation was neglected.Considering that thermal anomalies may occur on the ground surface with low intensity, an accurate surface thermal mapping should rely only on nighttime remotely sensed images rather than daytime measurements.

Conclusions
This study is based on the comparison among satellite data with measurements collected in different campaigns (airborne and ground).By analysing 2000-2017 spaceborne remotely sensed time series data of surface temperature covering the whole island, we have detected hot pixels and used them as focal points.Assuming that all pixels on the island are affected by the same seasonal conditions, the hottest pixels identified will be warmer independently of seasonal and daily variations.During the daytime, the reflected solar radiation is difficult to remove from thermal channels.Moreover, the error caused by the solar correction term can affect the retrieval accuracy, especially under low reflectance in thermal infrared regions.Therefore, in our analysis, the reflected solar radiation was neglected.Considering that thermal anomalies may occur on the ground surface with low intensity, an accurate surface thermal mapping should rely only on nighttime remotely sensed images rather than daytime measurements.
The multiscale approach has been offered by measurement of the surface temperature data collected during field campaigns and remotely sensed data from airborne and spaceborne platforms.The multitemporal characteristics have been approached by comparing long-term time series (satellite data) and discontinuous measurements (ground and airborne campaigns).The potential of these new applications could improve the results of direct monitoring systematically performed for volcanic surveillance since 1990 by a network of stations for ground temperature measurements.Surface temperature estimated by remote sensing data has been validated using independent measurements provided by the Oplao website.
In order to validate the accuracy of satellite surface temperature, we compared ASTER/L8 data, close in time, to the airborne and ground campaigns, paying attention to the hottest pixels highlighted in the above-mentioned procedure.The statistical analysis through Pearson correlation index, NSE and NB has shown a strong agreement among the data with a Pearson correlation index varying from 0.8956 to 0.9276.At the same time, the NSE and NB indices show very high correlation considering the values ranging from 0.1117 to 0.1377 and from −0.0157 to 0.0531, respectively.The reliability of the surface temperature retrieval methodology has allowed to compare these values with those acquired by field and airborne measurements revealing differences of about 5% with nighttime acquisition and up to 50% with daytime acquisition.
After detecting the temperature changes affecting the fumaroles areas, a further outcome of this study would be to extend the comparison of data obtained through direct and indirect acquisition to the areas affected by a diffused thermal release.With this constraint the thermal monitoring of heat flux of La Fossa cone, emitted from both the fumaroles and the diffuse outgassing, could provide important support for the detection of the early warning phenomena that precede a volcanic crisis.Moreover, the geophysical and geochemical signals for the last 29 years have been located within the La Fossa caldera, the summit area of the La Fossa Cone and the village located at its foot.Indeed, during the monitored period, the interval of time with anomalous energy output always correlated with the temperature change in the natural volatile emanations, and they could all be updated by the continuous monitoring of the summit fumaroles.
In order to highlight the main changes in the thermal release from a solphataric area, it is suggested to implement a systematic procedure based on the integration and comparison of the surface temperatures retrieved using remote sensing with direct method measurements.
In this work we have analysed the possibility of observing and estimating thermal anomalies on Vulcano island using a multiscale and multitemporal approach.By considering the reliability of the proposed approach, the authors believe that it can be also used for monitoring volcanic areas where ground networks are not present or because the volcanoes lie in remote areas.

Figure 1 .
Figure 1.Location of the Island of Vulcano in the Aeolian archipelago, southern Tyrrhenian Sea and zoomed image of Vulcano based on a Lidar-derived DEM-shaded relief image [47].

Figure 1 .
Figure 1.Location of the Island of Vulcano in the Aeolian archipelago, southern Tyrrhenian Sea and zoomed image of Vulcano based on a Lidar-derived DEM-shaded relief image [47].

Figure 2 .
Figure 2. The two sites selected for the analysis: La Fossa crater (a), "Fangaia" mud pool (b) and Vulcano weather station (c).In the box the two areas viewed by L8 on 8 July 2017.Faraglione is located near Levante harbour and belongs to the Vulcanello islet.It is geologically characterized by brecciated pyroclastic material, deeply altered by the strong hydrothermal release and locally hardened (Faraglione formation) with a leucite-bearing shoshonite composition [5456], ascribed to the early Vulcanello activity [4951].La Fossa is a small stratovolcano that lies in the central part of the island; it reaches an elevation of 391 m a.s.l. and it has been active since 5.5 ka[5557].La Fossa was formed by phreatomagmatic eruptions that generated dilute pyroclastic density currents [5658] and references therein.The thermal aquifers and the widespread CO2 soil degassing in these two areas testify to the actual volcanic activity [5759, 5860--6062].Intense fumarolic degassing persists in La Fossa crater, and significant variations observed in the fumarole vents temperature at La Fossa crater have beenwere documented.In 1926, fumarole temperatures reached more than 600 °C and the highest value, 670 °C, was measured in 1993[6163].Until March 1995, La Fossa temperatures were higher than 500 ºC, while fumaroles located outside the crater vents ranged between 350 °C and 500 °C.After 1995, the temperature of fumaroles started to decrease while the surface of active vents increased, and many new vents opened at the edge of the active cone [1212, 6264].In particular, La Fossa crater is characterized by thermal anomalies that, during the seismic activity which occurred from 2004 to 2006 [6365, 6466],

Figure 2 .
Figure 2. The two sites selected for the analysis: La Fossa crater (a), "Fangaia" mud pool (b) and Vulcano weather station (c).In the box the two areas viewed by L8 on 8 July 2017.

Figure 3 .
Figure 3. Histogram with the number of satellite acquisitions per year.

Figure 3 .
Figure 3. Histogram with the number of satellite acquisitions per year.2.1.1.Processing of L8 Data L8 carries two sensor payloads: the OLI (Operational Land Imager) and the TIRS (Thermal Infrared Sensor).In this work, band from 10.6 to 11.2 µm (B 10 ) acquired by TIRS sensor is used; B 11 known calibration problems in the dual thermal bands of L8 [67,68] were taken into account and therefore are not considered.Since the L8 TIRS data are delivered in quantized and calibrated scaled Digital Numbers (DN) representing all the multispectral image, the B 10 data were automatically extracted and converted to Top Of Atmosphere (TOA) radiance by using the following equation:

Figure 4 .
Figure 4. Location on the top of La Fossa Cone of the sensors for the thermal monitoring included in the INGV network (on a Lidar-derived DEM-shaded relief image [47]).The interval of acquisition is one hour.

Figure 4 .
Figure 4. Location on the top of La Fossa Cone of the sensors for the thermal monitoring included in the INGV network (on a Lidar-derived DEM-shaded relief image [47]).The interval of acquisition is one hour.

Figure 6 .
Figure 6.Scatter plots of surface temperature retrieved by ASTER (left) and L8 (right) data and temperature measured by the meteo station on La Fossa crater.

Figure 6 .
Figure 6.Scatter plots of surface temperature retrieved by ASTER (left) and L8 (right) data and temperature measured by the meteo station on La Fossa crater.

Figure 7 .
Figure 7. Scatter plots of surface temperature retrieved by ASTER (left) and L8 (right) data and temperature measured by the meteo station in "Fangaia" mud pool.

Figure 7 .
Figure 7. Scatter plots of surface temperature retrieved by ASTER (left) and L8 (right) data and temperature measured by the meteo station in "Fangaia" mud pool.

Figure 8 .
Figure 8.Time series of temperature selected from the continuous monitoring network of INGV.

Figure 8 .
Figure 8.Time series of temperature selected from the continuous monitoring network of INGV.

Figure 9 .
Figure 9. Thermographic map of the La Fossa crater rim acquired on16/07/2014 with a hand-held IR camera.The red boxes mark the craters showing high thermal anomalies.

Figure 9 .
Figure 9. Thermographic map of the La Fossa crater rim acquired on16/07/2014 with a hand-held IR camera.The red boxes mark the craters showing high thermal anomalies.
• C) as falling in the water saturated layer; the second class of values (included in the range from 20 to 27 • C) was interpreted as the surface condition dominated by conductive heat transfer, absence of liquid water and very low steam content; the third class (temperatures in the range from 28 to 48 • C) was for the steam heated soil; finally, the hottest class of values included the low temperature fumaroles (T > 48), where the convective transfer dominates.Remote Sens. 2018; 10; x FOR PEER REVIEW 14 of 32 some theoretical distribution, such as a normal or exponential one.The data values are ordered and the cumulative distribution values are calculated as (i-0.5)/nfor the i-th ordered value out of n total values.The scatterplot is created by plotting the real set of quantiles against the normal distributed population of the same quantiles; a line that is roughly straight corresponds to a normal distribution of the real dataset.Any diversion from the straight line indicates a diversion of the dataset from the normal distribution.The probability plots reveal different populations of data.According to these populations the space variability observed on the surface temperatures has been inferred from the main components of the heat transport (radiative, convective conductive) which participated in the heat exchange processes according to the different rate of flux reaching the ground surface.

Figure 10 .
Figure10.qq-plot of the infrared temperature measurements (La Fossa crater rim).The red line identifies the thermal quantile groups of dataset and the range temperature; the black line is the reference line to assess the Normal Gaussian distribution assumption).

Figure 11 .
Figure 11.Surface temperature retrieved by ATM from the 12 September 2010 predawn flight.Higher temperatures relate to fumaroles distributed within the La Fossa crater and along the coast of the Volcano bay.The red line corresponds to the temperature profile of Figure 12 and the white dot to the measurement location of Figure 13.

Figure 12 .
Figure 12.Sea water temperature profile along the red line shown in Figure 11.The red line intersects the temperature value estimated at that location of Figure 13.

Figure 11 .
Figure 11.Surface temperature retrieved by ATM from the 12 September 2010 predawn flight.Higher temperatures relate to fumaroles distributed within the La Fossa crater and along the coast of the Volcano bay.The red line corresponds to the temperature profile of Figure 12 and the white dot to the measurement location of Figure 13.

32 Figure 11 .
Figure 11.Surface temperature retrieved by ATM from the 12 September 2010 predawn flight.Higher temperatures relate to fumaroles distributed within the La Fossa crater and along the coast of the Volcano bay.The red line corresponds to the temperature profile of Figure 12 and the white dot to the measurement location of Figure 13.

Figure 12 .
Figure 12.Sea water temperature profile along the red line shown in Figure 11.The red line intersects the temperature value estimated at that location of Figure 13.

Figure 12 .
Figure 12.Sea water temperature profile along the red line shown in Figure 11.The red line intersects the temperature value estimated at that location of Figure 13.

Figure 13 .
Figure 13.Sea water measured from the thermo camera.

Figure 13 .
Figure 13.Sea water measured from the thermo camera.

Figure 20 .
Figure 20.IR map acquired on La Fossa crater during the 16/07/2014 with a hand-held IR camera with satellite pixel.Red dots indicate the number of pixels as described in Figure 14-Figure 17.

Figure 20 .
Figure 20.IR map acquired on La Fossa crater during the 16/07/2014 with a hand-held IR camera with satellite pixel.Red dots indicate the number of pixels as described in Figures 14-17.

Figure 21 .
Figure 21.Red points represent the geographical position of the hottest pixels derived by time series satellite data for La Fossa crater.

Figure 21 .
Figure 21.Red points represent the geographical position of the hottest pixels derived by time series satellite data for La Fossa crater.

Figure 22 .
Figure 22.Histogram reports the numbering of the satellite pixel along the x-axis and the occurrence of the hottest pixels along the y-axis for La Fossa crater.

Figure 22 .
Figure 22.Histogram reports the numbering of the satellite pixel along the x-axis and the occurrence of the hottest pixels along the y-axis for La Fossa crater.

Figure 23 .
Figure 23.Red points represent the geographical position of the hottest pixels derived by time series satellite data for the "Fangaia"mud pool.

Figure 23 .
Figure 23.Red points represent the geographical position of the hottest pixels derived by time series satellite data for the "Fangaia"mud pool.

Figure 23 .
Figure 23.Red points represent the geographical position of the hottest pixels derived by time series satellite data for the "Fangaia"mud pool.

Figure 24 .
Figure 24.Histogram reports the numbering of the satellite pixel along the x-axis and the occurrence of the hottest pixels along the y-axis for "Fangaia" mud pool.

Figure 24 .
Figure 24.Histogram reports the numbering of the satellite pixel along the x-axis and the occurrence of the hottest pixels along the y-axis for "Fangaia" mud pool.

Figure 25 .
Figure 25.Thermal anomalies detected by ATM Thermal data (Red dots).The two green boxes are the area considered for the temperature meaning.Brightness temperature of ATM Thermal channel 11 over the Vulcano island is represented in grey scale with shaded relief represented in green.

Figure 25 .
Figure 25.Thermal anomalies detected by ATM Thermal data (Red dots).The two green boxes are the area considered for the temperature meaning.Brightness temperature of ATM Thermal channel 11 over the Vulcano island is represented in grey scale with shaded relief represented in green.

Figure 26 .
Figure 26.ASTER acquired on 6 September and ATM-retrieved temperature comparison of the "Fangaia" area (left) and La Fossa Crater (right).

Figure 26 .
Figure 26.ASTER acquired on 6 September and ATM-retrieved temperature comparison of the "Fangaia" area (left) and La Fossa Crater (right).

Table 1 .
Ground and Satellite data.

Table 1 .
Ground and Satellite data.

Table 3 .
Main characteristics of the predawn ATM flight.

Table 4 .
Comparison of temperature derived from the web site, field campaign and L8 data.

Table 4 .
Remote Sens. 2018; 10; x FOR PEER REVIEW 12 of 32 considered as further confirmation of the methodology used to estimate the temperature from satellites.Comparison of temperature derived from the web site, field campaign and L8 data.

Table 8 .
Comparison between ground measurements and satellite estimation temperatures in the mud pool " Fangaia".

Table 9 .
ASTER data of 15 July 2014 at 20:45 UTC, temperature measured with a hand-held IR camera and rescaled on ASTER pixel size, percentage of ground measured points falling in the 90x90 m ASTER pixel.

Table 9 .
ASTER data of 15 July 2014 at 20:45 UTC, temperature measured with a hand-held IR camera and rescaled on ASTER pixel size, percentage of ground measured points falling in the 90 × 90 m ASTER pixel.N • Pixel ASTER T ( • C) Thermal Camera T ( • C)

Table 10 .
Number of satellite (ASTER and L8) images for which the points are the hottest for the La Fossa area.

Table 10 .
Number of satellite (ASTER and L8) images for which the points are the hottest for the La Fossa area.

Table 11 .
Number of satellite (ASTER and L8) images for which the points are the hottest for "Fangaia" area.

Table 11 .
Number of satellite (ASTER and L8) images for which the points are the hottest for "Fangaia" area.

Table 12 .
Comparison of the mean temperature derived from ATM and ASTER considering the areas of La Fossa crater and "Fangaia".

Table 12 .
Comparison of the mean temperature derived from ATM and ASTER considering the areas of La Fossa crater and "Fangaia".Site ASTER Mean T ( • C)