Copernicus Sentinel-1 MT-InSAR, GNSS and Seismic Monitoring of Deformation Patterns and Trends at the Methana Volcano, Greece

The Methana volcano in Greece belongs to the western part of the Hellenic Volcanic Arc, where the African and Eurasian tectonic plates converge at a rate of approximately 3 cm/year. While volcanic hazard in Methana is considered low, the neotectonic basin constituting the Saronic Gulf area is seismically active and there is evidence of local geothermal activity. Monitoring is therefore crucial to characterize any activity at the volcano that could impact the local population. This study aims to detect surface deformation in the whole Methana peninsula based on a long stack of 99 Sentinel-1 C-band Synthetic Aperture Radar (SAR) images in interferometric wide swath mode acquired in March 2015–August 2019. A Multi-Temporal Interferometric SAR (MT-InSAR) processing approach is exploited using the Interferometric Point Target Analysis (IPTA) method, involving the extraction of a network of targets including both Persistent Scatterers (PS) and Distributed Scatterers (DS) to augment the monitoring capability across the varied land cover of the peninsula. Satellite geodetic data from 2006–2019 Global Positioning System (GPS) benchmark surveying are used to calibrate and validate the MT-InSAR results. Deformation monitoring records from permanent Global Navigation Satellite System (GNSS) stations, two of which were installed within the peninsula in 2004 (METH) and 2019 (MTNA), are also exploited for interpretation of the regional deformation scenario. Geological, topographic, and 2006–2019 seismological data enable better understanding of the ground deformation observed. Line-of-sight displacement velocities of the over 4700 PS and 6200 DS within the peninsula are from −18.1 to +7.5 mm/year. The MT-InSAR data suggest a complex displacement pattern across the volcano edifice, including local-scale land surface processes. In Methana town, ground stability is found on volcanoclasts and limestone for the majority of the urban area footprint while some deformation is observed in the suburban zones. At the Mavri Petra andesitic dome, time series of the exceptionally dense PS/DS network across blocks of agglomerate and cinder reveal seasonal fluctuation (5 mm amplitude) overlapping the long-term stable trend. Given the steepness of the slopes along the eastern flank of the volcano, displacement patterns may indicate mass movements. The GNSS, seismological and MT-InSAR analyses lead to a first account of deformation processes and their temporal evolution over the last years for Methana, thus providing initial information to feed into the volcano baseline hazard assessment and monitoring system.


Introduction
It is well known that some of the most scenic landscapes and islands of Greece are due to the intense volcanic activity that occurred centuries to millennia ago [1,2]. While Santorini is the most famous and a worldwide renowned touristic site, there are other active volcanic systems (i.e., Milos and Nisyros) within the south Aegean volcanic arc, which were formed from the subduction of the African tectonic plate beneath the Eurasian plate. Of this arc, the Methana peninsula is the westernmost dormant but geodynamically and hydrothermally active volcanic system [3]. The peninsula is connected with the north-eastern coast of Peloponnesus by a narrow gooseneck-shaped isthmus and stretches for approximately 44 km 2 northward into the Saronic Gulf.
Of the around 32 andesitic and dacitic lava domes scattered across the peninsula, the Methana volcano is the largest and is the subject of the present study. The volcanic hazard in Methana is considered "low" [4]; given that the last historic eruption was registered in approximately 230 BC, volcanic products derived from highly explosive eruptions were not found [3] and no alarming signs were observed in recent times. However, several aspects provide sufficient motivation for a dedicated investigation into the Methana volcano. From a hazard point of view, the neotectonic basin constituting the Saronic Gulf area is considered seismically active [5], and the active fault systems therein were considered preferential paths for present-day geothermal fluid leakage and, as such, potential sites for magma uprising [3]. Thermal springs are indeed clear manifestations of volcanic geothermal energy in Methana, and gas exhalations suggest a mixture between a dominant hydrothermal component and mantle-derived fluids [3]. Additionally, from a risk point of view, the volcano is not far from Athens (less than 50 km south-west) and any activity at the volcano would impact the local population (about 2500 inhabitants over the year, increasing during summer due to tourism). Further proof that monitoring the Methana volcano has recently been given higher priority on the agenda of Greek institutions is provided by the installations of six seismological stations by the Greek Institute of Geodynamics of the National Observatory of Athens (NOA) [6] alongside a new NOANET Global Navigation Satellite System (GNSS) station [7]. This recent development significantly updates what was stated by [4], i.e., of the whole Greek volcanic arc of Methana-Milos-Santorini (Thera)-Nisyros, Methana is among the volcanic fields not efficiently monitored.
In this context, the present study aims to investigate and characterize the ground motions observed at the surface of Methana volcano based on the results achieved by means of multi-interferogram processing of Synthetic Aperture Radar (SAR) images collected by the Copernicus Programme Sentinel-1 satellite constellation in the period from March 2015 to August 2019. Interferometric SAR (InSAR) approaches have been largely exploited to study and monitor surface deformation at different locations across the volcanic arc (e.g., [8][9][10][11][12][13][14][15][16][17]) and have proved their effectiveness to provide a spatially distributed estimation of volcanic activity due to magma chamber processes as well as shallow deformation associated with hydrothermal activity (e.g., low temperature venting). Despite this abundant literature, to the best of our knowledge, no InSAR study has previously focused on Methana.
Given the varied land cover and presence of vegetation at Methana, we intentionally selected a Multi-Temporal InSAR (MT-InSAR) processing approach providing both Persistent Scatterers (PS) and Distributed Scatterers (DS) in order to augment the monitoring capability and the number of measurement points. In this way, we aimed to characterize not only the geodynamics of the whole volcano but also the surface deformation along the steep slopes and stream gullies of the rugged terrain of the volcanic landforms, of which the susceptibility to landslides and rockfalls has been already highlighted in the literature [18]. Satellite MT-InSAR is complemented and interpreted in combination with an analysis of regional seismicity characterizing the 2006-2019 period and with geodetic data from continuous GNSS monitoring and Global Positioning System (GPS) benchmark surveying. The latter is also used to calibrate and validate the MT-InSAR results, thus tying the satellite-based estimates to the local geodetic reference used in the region and providing evidence to corroborate MT-InSAR observations.

The Study Area
Methana peninsula is located in the western Saronic Gulf south of Athens city ( Figure 1) and is one of the three volcanic centres (i.e., Aegina, Poros and Methana) representing the north-western part of the active volcanic arc of the South Aegean. This Hellenic active volcanic arc also includes the islands Milos, Santorini and Nisyros in its extension [19]. Methana and the other volcanic centres in the Saronic Gulf are mostly monogenetic, and no composite volcanic structures are present [20].
Appl. Sci. 2020, 10, x FOR PEER REVIEW 3 of 24 satellite-based estimates to the local geodetic reference used in the region and providing evidence to corroborate MT-InSAR observations.

The Study Area
Methana peninsula is located in the western Saronic Gulf south of Athens city ( Figure 1) and is one of the three volcanic centres (i.e., Aegina, Poros and Methana) representing the north-western part of the active volcanic arc of the South Aegean. This Hellenic active volcanic arc also includes the islands Milos, Santorini and Nisyros in its extension [19]. Methana and the other volcanic centres in the Saronic Gulf are mostly monogenetic, and no composite volcanic structures are present [20]. According to [24], the Saronic Gulf is the longest active area of the South Aegean with geothermal springs which are influenced by the active tectonic structure of the area as well as the existence of the magmatic chamber. Volcanic activity in Methana lasted until 1700, when a submarine eruption occurred north of Kameni Chora. Present day, active geothermal springs are reported mainly in the eastern and northern parts of the peninsula [3]. Within Methana town, in the south-east, the thermal springs are exploited at the open-air baths of the Thermal Spa and its neoclassical dwelling built in 1930.
It is also crucial to mention the presence of the Pausanias submarine volcano, located 2 km off the north-western shore of the peninsula [25]. According to [26], this submarine volcano was discovered using seismic reflection data and volcanic activity occurred during the 3rd century BC.
The Saronic Gulf constitutes a neotectonic basin which was created 4-5 million years ago [27,28]. NE-SW normal faults dominate the broader region, while N-S, E-W and NNE-SSW striking According to [24], the Saronic Gulf is the longest active area of the South Aegean with geothermal springs which are influenced by the active tectonic structure of the area as well as the existence of the magmatic chamber. Volcanic activity in Methana lasted until 1700, when a submarine eruption occurred north of Kameni Chora. Present day, active geothermal springs are reported mainly in the eastern and northern parts of the peninsula [3]. Within Methana town, in the south-east, the thermal springs are exploited at the open-air baths of the Thermal Spa and its neoclassical dwelling built in 1930.
It is also crucial to mention the presence of the Pausanias submarine volcano, located 2 km off the north-western shore of the peninsula [25]. According to [26], this submarine volcano was discovered using seismic reflection data and volcanic activity occurred during the 3rd century BC.
The Saronic Gulf constitutes a neotectonic basin which was created 4-5 million years ago [27,28]. NE-SW normal faults dominate the broader region, while N-S, E-W and NNE-SSW striking normal faults are also recognizable in the area [25] and are characterized by low seismic activity, chiefly with shallow events [29].
The geology of Methana consists of various volcanic formations with Late Pliocene-Pleistocene small andesitic to dacitic flows and domes [30,31], which overlay the older volcanic rocks, as well as the Mesozoic sedimentary basement [21,32].
A simplified geological map (Figure 1a) shows that the grey limestones are present in the northwest (Upper Triassic-Lower Jurassic), while limestones with dolomites are present in the south (Upper Jurassic-Cretaceous). These limestones constitute the basement of Methana and are located under the volcanoclastic formations (rhyolite, dacite and andesite) [33]. Additionally, the Methana peninsula has many tectonic faults, with a major fault crossing the town of Methana from west to east.
The tectonic regime of the region combined with the volcanic activity determined a rough topography of the peninsula, which shows both flat basin areas with Quaternary fill and zones dominated by high slope gradients and narrow valleys, often modelled by erosion processes (such as in the south, between Vathi and Methana) [18]. The maximum elevation is 740 m above sea level (a.s.l.) and occurs at Helona Mountain (Figure 1b).

SAR Data and Multi-Temporal InSAR Processing
A long stack of 99 SAR scenes acquired by the Copernicus Sentinel-1A and Sentinel-1B satellites in ascending mode (relative orbit 102) was used in this study. SAR images are remote sensing satellite data increasingly used for observing, mapping and monitoring Earth's surface processes [34] and are characterized by cloud-penetrating, day and night operational capabilities.
Among the Sentinel-1 SAR constellation imaging modes, we used the Interferometric Wide (IW) swath, which is based on the novel Terrain Observation with Progressive Scans (TOPS) strategy [35]. IW images are composed of three sub-swaths (i.e., IW1, IW2 and IW3), each consisting of a series of bursts, and provide a 250-km large coverage, with pixel resolutions of 5 m and 20 m (single look) in range and azimuth, respectively [36]. The Sentinel-1 data used in this study were collected from the Sentinel Hub platform and cover the time period from March 2015 to August 2019. The ascending mode acquisition geometry was selected to undertake the analysis in order to optimize visibility of the east-facing slopes, where the town of Methana (and thus the major urban settlement of the peninsula) is located (see Figure 1).
The stack of 99 Sentinel-1 SAR scenes was processed using an advanced MT-InSAR approach, an extension of the basic technique of Differential InSAR (DInSAR; e.g., [37,38]), that is based on the phase comparison of multiple SAR images gathered at different times over the study region and allows for measurements of land deformation along the line-of-sight (LOS) direction of the SAR sensor with up to millimetre precision (e.g., [39]). Advanced DInSAR methods are increasingly being used in various applications in the field of volcanology, seismology, crustal dynamics, landslides, land subsidence and geothermal energy exploitation (e.g., [40][41][42][43][44][45]).
MT-InSAR methods exploit either Persistent (PS) or Distributed (DS) Scatterers or coherent targets. Typically, PS are artificial objects that reflect the radar signal well, such as metal structures, buildings and rock outcrops, and are used in Persistent Scatterer Interferometry (PSI), e.g., the Permanent Scatterers (PS-InSAR) technique [46,47]. In urban areas, there is a prevalence of PS, and PSI methods allow for analysis of even individual structures on the ground. Methods exploiting coherent targets include algorithms such as the Small BAseline Subset (SBAS) method [48]. DS reflects lower radar energy compared to PS targets and usually spans several pixels in high-resolution SAR images [49], which exhibit similar scattering properties and can be used together for deformation estimation. The identification and monitoring of such targets is helpful especially in suburban and more rural areas, where the density of PS can be low. In such circumstances, the combination of PS and DS is important for effective displacement monitoring using SAR Single Look Complex (SLC) data. What is certain is that the processing aiming to detect both PS and DS requires large computing resource and processing time [50]. Concerning the precision of MT-InSAR measurements, error bars are extremely complex to estimate theoretically, since they depend on several factors, including the number of SAR SLC images used, the spatial density and quality (signal-to-noise ratio levels) of the measurement targets, their distance from the reference PS, the climatic conditions at the time of SAR acquisition and the satellite repeat cycle. Recent studies provide an empirical relationship between the number of SAR interferograms used and the resulting precision of the ground deformation estimates (e.g., by using 50 ERS-1/2 interferograms, a precision of 0.15 cm/year can be achieved with the SBAS method [51]).
In this study, the GAMMA SAR and Interferometry software (GAMMA Remote Sensing AG, Gümligen, Switzerland) and the Interferometric Point Target Analysis (IPTA) method [52] were used for the interferometric processing. The MT-InSAR processing workflow can be divided into two steps: the pre-processing and the main processing using GAMMA/IPTA.
During pre-processing, parameter files for the 99 SLC images were created and used to perform orbital refinement. The SLC images were co-registered based on the image used as a reference (20/09/2017), and a single-reference stack of 98 differential interferograms was generated. The Shuttle Radar Topography Mission (SRTM) 1-arcsec (i.e., 30 m/pixel) Digital Elevation Model (DEM) [53] was used during pre-processing to subtract the topography contribution from the interferograms.
Using the 99 co-registered SLC images, a multi-reference stack was also created, containing both single-pixel and multi-look differential interferometric phases. Multi-look differential interferometric phases (DS) were first generated, and afterwards, single-pixel differential interferometric phases (PS) were estimated using low spectral diversity and temporal variability of SLC intensity. Then, single-look and multi-look phases were combined into one list of both PS and DS. The multi-reference stack consists of 291 differential interferograms, pairing each scene with the 3 subsequent scenes (triple redundancy). The perpendicular baselines of those 291 pairs were quite short, with the longest baseline being 202 m. A high-accuracy Digital Surface Model (DSM) (5.5 m/pixel) provided by the Hellenic Mapping and Cadastral Organisation (Figure 1b) was also used for the IPTA in order to improve the quality of the output results. As reference point, the closest point to the GPS benchmark "Methana South" (i.e., MESO station in Figure 1a; 37 • 34 51.02562 N, 23 • 22 12.14638 E) was selected. Then, height corrections and atmospheric phases were computed using the combined point list.
The above model was refined by updating height corrections and atmospheric phases via iteration of the last step two more times. Finally, the phases were converted into deformation time series, which contained both PS and DS. Then, the deformation time series database was split into two separate ones, one for the PS and the other for the DS targets.

GNSS Continuous Monitoring Data and GPS Benchmark Surveying
A small network of two benchmark stations was established in 2006 in Methana Peninsula ( Figure 1) [54]. These two stations were part of a larger GPS network established in the north-western part of the Hellenic Volcanic Arc (including the island of Aegina and the volcanic area of Susaki) aiming to study the ground deformation of this area. The stations in Methana were located in the northern (i.e., station MENO) and southern (MESO) parts of the peninsula and were reoccupied several times up to December 2019. All measurements took place approximately within the same period of the year (usually late October to December) in an effort to minimize seasonal phenomena. The Methana benchmark network was occupied using Leica receivers (SR9500 and AX1200). For each measurement period, the two benchmark stations were occupied for a period of 48 to 60 h, with time sampling of 15 s. Such time sampling was used for local-scale solutions with nearby continuous GNSS stations for which data were available with the same time sampling. While on regional-scale solutions, the data were resampled to 30-s data sets in order to be processed together with the freely distributed 30-s GNSS data available for European stations.
In early 2015, a continuous GNSS (cGNSS) station was also established in the Lygourio (LYGO) close to Methana peninsula as part of the HxGN SmartNet (https://gr.nrtk.eu/), that, ever since, has acted as local reference station. In addition, since August 2019, cGNSS data were available from a station (MTNA) located in the western part of Methana peninsula as part of the GNSS network of the Institute of Geodynamics, National Observatory of Athens (NOA). Horizontal velocity results from two more cGNSS stations were also available from the National Technical University of Athens (NTUA), i.e., station METH installed in 2004 in the southern sector of Methana peninsula, and from [55], i.e., station 010A located to the south-east of Methana peninsula, at Poros.
The available 30-s data from the two cGNSS stations (i.e., LYGO and MTNA) that were included in this study were processed separately on a daily basis together with several other cGNSS stations located in Greece as part of a regional network. From the time series of the station's coordinates, the cGNSS station velocity was estimated in the IGb08 reference frame (see Figures A1 and A2 in Appendix A).
The GPS data were processed using BERNESE software (v.5.2, University of Bern, Bern, Switzerland) [56]. The software allows for the estimation of epoch-wise receiver coordinates in Precise Point Positioning (PPP) mode as well as in the double difference mode. For the present research, the coordinates were calculated on a static mode. To get the best results, the PPP technique was used as a first step to get a priori coordinate values that were consequently introduced in the more precise double difference method. The GPS data from the benchmark stations and the other cGNSS stations in the broader area of Methana were processed together with several other EUropean REference Frame (EUREF) and cGNSS stations from Europe and Greece (http://www.epncb.oma.be) in order to define the local reference frame and to calculate the station coordinates. All measurement periods were processed with respect to the IGb08 reference frame (http://igscb.jpl.nasa.gov/network/refframe.html). Absolute antenna phase centre corrections were used in the data analysis. Several supporting files were introduced during the processing steps, aiming to define precise coordinates for the stations. Precise orbits from CODE analysis centre (ftp://ftp.unibe.ch/aiub/CODE) were used in the analysis, and a set of Earth orientation parameters was calculated. Crustal deformation caused by changing mass due to ocean tides and the atmosphere was taken into consideration, introducing auxiliary files. For the ocean loading effect, tide loading corrections were based on the FES2004 model (http://holt.oso.chalmers.se/loading). Atmospheric tidal loading coefficients were used from a global grid based on the Ray and Ponte (2003) model [57]. For troposphere modelling, the Vienna Mapping Function (http://ggosatm.hg.tuwien.ac.at) as well as the Neill mapping function were used. Ambiguities were solved using several resolution strategies based on the baseline length and the occupation time: (i) the wide-lane and narrow-lane techniques for the medium baselines (<10-200 km), where initially the L5 linear combination is processed (wide-lane) and introduced as known (fixed) in the subsequent run (narrow-lane) where the L3 linear combination is processed; (ii) the SIGMA strategy for the very short baselines (in between the two Methana stations, where both L1 and L2 data were available and used); and finally, (iii) the quasi-ionosphere-free (QIF) strategy for the long baselines (<1000-2000 km) were the main adopted ambiguity resolution strategies.
The long occupational time of the two GPS benchmark stations in Methana, during all measurement periods, resulted into low uncertainties (at a 95% confidence level) with overall errors of about 2.4-3.1 mm and 2.5-5.1 mm for the horizontal and vertical coordinates, respectively. The processing performed on daily sessions allowed for the extraction of a set of coordinates for each 24-h long period. The final set of coordinates was the mean value of all the sessions, and the errors led to quantification of the deviation of the daily coordinate solutions from the mean value.

Seismological Data and Processing
The digital seismological data used in this study were obtained from the database of the seismological laboratory of the National and Kapodistrian University of Athens (NKUA; http://dggsl. geol.uoa.gr/en_index.html).
The optimization of epicentres for the period 2006-2019 was carried out using the HYPOINVERSE software (v.1.4, United States Geological Survey: Menlo Park, CA, USA) [58]. The final quality of the relocated epicentres is very satisfactory due to the use of the one-dimension local velocity model developed by Karakonstantis et al. [59]. In particular, the root mean square (RMS) travel-time residual for the whole study period was determined to be 0.46 s, while the means of the horizontal and vertical direction uncertainties for the events' hypocentres are 1.51 km and 5.82 km, respectively.
The seismological data were divided into two time periods. The first one coincided with the time period that the GNSS data cover, i.e., January 2006 to August 2019, while the second period covered the time span of the interferometric processing, i.e., March 2015 to August 2019 (see Section 3.1).

Evidence from the GNSS Data Analysis
The observed velocities (i.e., V E , V N and V U , along the east-west, north-south and vertical direction, respectively) and associated standard deviations (i.e., STDV) from all cGNSS and local benchmark stations are presented in Table 1. The LYGO cGNSS station exhibited SSE motion (Figure 2), calculated for a period of more than five years, that describes the anticipated regional motion of the area with respect to IGb08 [55] and is in good agreement with the two other cGNSS stations: METH and 010A. The newly established cGNSS station in the western Methana peninsula (i.e., MTNA) showed significantly different velocity due to the short period that the station was operating (i.e., from August to December 2019). The limited operational period reveals only a short portion of the station's annual trend and therefore could be biased by interannual variations. Consequently, these data have been excluded from further investigation but still briefly presented for future comparison (see Figure A1 in Appendix A). The velocity results from the two local benchmark stations in Methana peninsula (i.e., MENO and MESO) for a period of almost fourteen years (2006-2019) appeared consistent with the regional velocity field concerning the horizontal component ( Figure 2). However, the southern station (i.e., MESO) revealed significant subsidence that has to be correlated with the MT-InSAR results (see Section 4.5). Aiming to better define the local deformation in Methana peninsula, the regional velocity vector (as it was estimated in LYGO cGNSS station) was subtracted from the two local GPS benchmark stations as well as the two other cGNSS sites (i.e., METH and 010A). Although the resulting horizontal velocity vectors are very small (<1 mm/year) with significant errors (Table 2) indicating a uniform deformation of the broader Methana area, there is a pattern of differential motion in the area (Figure 3). The northern station (i.e., MENO) exhibits a westward motion, while the southern sites (i.e., MESO and METH) showed a northward motion. However, the latter is quite arbitrary due to the small horizontal velocity values and the resulting errors. The more evident differential behaviour between these two stations is in the vertical component (Table 2), with the southern station (i.e., MESO) exhibiting noticeable subsidence, as it was mentioned previously.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 8 of 24 the southern sites (i.e., MESO and METH) showed a northward motion. However, the latter is quite arbitrary due to the small horizontal velocity values and the resulting errors. The more evident differential behaviour between these two stations is in the vertical component (Table 2), with the southern station (i.e., MESO) exhibiting noticeable subsidence, as it was mentioned previously.

Observations from the Seismic Data Analysis
The analysis of the seismic data shows that the Saronikos Gulf region is mainly characterized by moderate seismic events and shallow depths, as already discussed in the literature [29].

Observations from the Seismic Data Analysis
The analysis of the seismic data shows that the Saronikos Gulf region is mainly characterized by moderate seismic events and shallow depths, as already discussed in the literature [29].

Observations from the Seismic Data Analysis
The analysis of the seismic data shows that the Saronikos Gulf region is mainly characterized by moderate seismic events and shallow depths, as already discussed in the literature [29].

PS and DS Distribution
After MT-InSAR processing, 4769 PS targets and 6234 DS targets were retrieved ( Figure 6). The spatial densities of the targets in the area were calculated as 95.14 PS/km 2

PS and DS Distribution
After MT-InSAR processing, 4769 PS targets and 6234 DS targets were retrieved ( Figure 6). The spatial densities of the targets in the area were calculated as 95.14 PS/km 2

MT-InSAR Data Precision, Calibration and Validation
The present study is based on ascending SAR dataset only (see Section 3.1). Therefore, all results and values reported hereinafter are LOS values if not explicitly stated otherwise.
The uncertainty observed in the estimated VLOS (i.e., standard deviation of the time series) is 2.8 mm/year on average for the PS dataset and 2.6 mm/year for the DS dataset. These values provide an indication of the precision of the MT-InSAR results and suggest that the ±3.0 mm/year VLOS interval can be considered the velocity range, indicating stability.
On the other hand, in order to estimate the accuracy of the MT-InSAR results, these were validated using the local velocity values of the GPS benchmark stations MESO and MENO in Methana Peninsula with respect to LYGO cGNSS station ( Table 2). A frame of reference transformation was applied to the velocity vectors in order to calculate the estimated LOS velocity of the two GPS benchmark stations and to allow for comparison with the MT-InSAR data.
Knowledge of the Euler angles (i.e., ωx, ωy and ωz; Table 3) of the Sentinel-1 LOS along the ascending relative orbit 102 allowed us to generate the rotation table for each axis: R1, R2 and R3 via Equations

MT-InSAR Data Precision, Calibration and Validation
The present study is based on ascending SAR dataset only (see Section 3.1). Therefore, all results and values reported hereinafter are LOS values if not explicitly stated otherwise.
The uncertainty observed in the estimated V LOS (i.e., standard deviation of the time series) is 2.8 mm/year on average for the PS dataset and 2.6 mm/year for the DS dataset. These values provide an indication of the precision of the MT-InSAR results and suggest that the ±3.0 mm/year V LOS interval can be considered the velocity range, indicating stability.
On the other hand, in order to estimate the accuracy of the MT-InSAR results, these were validated using the local velocity values of the GPS benchmark stations MESO and MENO in Methana Peninsula with respect to LYGO cGNSS station ( Table 2). A frame of reference transformation was applied to the velocity vectors in order to calculate the estimated LOS velocity of the two GPS benchmark stations and to allow for comparison with the MT-InSAR data.
Knowledge of the Euler angles (i.e., ω x , ω y and ω z ; Table 3) of the Sentinel-1 LOS along the ascending relative orbit 102 allowed us to generate the rotation table for each axis: R 1 , R 2 and R 3 via Equations (1)-(3), respectively: Table 3. Euler angles of the Sentinel-1 LOS along ascending mode relative orbit 102.

Euler Angle Degrees Physical
Interpretation Declination from the north (azimuth) The right-hand rule was used to calculate the Euler angles because the reference frames are right-handed. The incidence angle for the area where the two GPS benchmark stations are found (namely, ω y ) was estimated using the European Space Agency (ESA)'s SeNtinel Application Platform (SNAP) software -Sentinel-1 toolbox S1TBX (v.7.0, SkyWatch, Waterloo, Canada),while the azimuth angle of the ascending relative orbit 102 in the area (namely, ω z ) was estimated using QGIS software (v.3.14, QGIS.ORG, Switzerland).
Then, using Equation (4), the rotation matrix (5) was composed and finally exploited to transform the reference frame into the LOS geometry: The velocity vectors of GPS benchmark stations MENO, i.e., (6), and MESO, i.e., (7), were transformed to the new reference of frame using Equation (8): The LOS velocity of the two benchmark stations was finally estimated: V LOS = −0.16 mm/year for MENO and V LOS = −3.61 mm/year for MESO.
We calibrated the PS and DS datasets by tying the estimated V LOS at the closest scatterer to the MESO station (i.e., 40 m away) to the value measured at the GPS station (i.e., −3.61 mm/year). As the scatterer originally exhibited a V LOS of +0.26 mm/year, the adjustment meant a V LOS shift of −3.87 mm/year, which was applied to all PS and DS targets.
At MENO station, the calibrated V LOS from the MT-InSAR analysis equals −0.19 mm/year, thus matching extremely well with the −0.16 mm/year V LOS measured at the GPS station and providing validating evidence for the MT-InSAR dataset.
To account for the velocity shift due to calibration of the PS and DS targets to MESO station, their time series were also adjusted as follows [61]: where d i (t n ) is the adjusted position of target i at time t n , d i old (t n ) is its original position at time t n , V sh is the V LOS shift and t 0 is the first date of the time series.

Interpretation of the Observed Ground Deformation
The maximum V LOS away from the sensor (i.e., negative values) that is observed after calibrations is −8.5 mm/year from the PS dataset and −18.1 mm/year from the DS dataset ( Figure 6). PS and DS with the strongest movement away from the sensor are located on the eastern flank of the peninsula. Considering the ascending mode-looking geometry and the orientation of the LOS, this could plausibly indicate the occurrence of downslope mass movements. Compared with the topography displayed in Figure 1b, these movements concentrate on the volcanic ranges and gully erosion landforms above the urban settlements of Methana, Kipseli, Agii Teodori and Agios Georgios. Photo interpretation using Google Earth imagery, however, did not highlight changes of specific landforms for the 2012-2019 period.
The maximum V LOS towards the sensor (i.e., positive values) from the PS dataset is +7.5 mm/year, and that from the DS dataset is +3.8 mm/year ( Figure 6). Some values above +3.0 mm/year are observed for very few PS and DS scatterers, mainly located on the western flank of Methana peninsula. These scatterers are, however, characterized by a high uncertainty of their deformation rate (i.e., up to 4.0-5.0 mm/year). Their estimated velocities are therefore not considered as robust as those of other scatterers characterized by higher precision.
The northern portion of Methana town is built onto alluvium, and it seems more stable (V LOS = −1.8 mm/year) than the central zone of the town, which is moving away from the satellite at −2.5 mm/year (Figure 8a), although both zones show velocities within the range of stability. The central zone is built onto volcanoclasts (Figure 8c), and the southern zone onto limestone (−2.9 mm/year). The southern zone is moving away from the satellite with relatively low V LOS . Overall, it appears that the majority of the urban area footprint (Figure 8d) is not affected by significant deformation, while somehow, greater rates are found in the suburban zones, as revealed by the DS targets.
The deformation scenario at Vathi and Paleokastro (Figure 9a) shows that there is a slight movement away from the satellite also in this sector of the peninsula. Similar to what was observed for Methana town, PS targets in the region indicates that the urban areas have a mean value of −0.9 mm/year, while the DS targets for the suburban and more rural zones show −3.1 mm/year. Annual LOS velocities do not reveal any significant deformation patterns across the sparsely vegetated outcrops at Mavri Petra (Figure 10a), which mostly shows both PS and DS with V LOS in the range of ±3.0 mm/year and absence of any clear deformation patterns from the sole observation of the annual LOS velocities. However, evidence of seasonal deformation can be detected in the time series at several locations across the peninsula, including at Mavri Petra ( Figure 10b) and at many other sites, for instance at Vathi and Paleokastro (Figure 9b) and at Methana town (Figure 8b). However, in the latter, this happens with less prominent amplitudes. −1.8 mm/year) than the central zone of the town, which is moving away from the satellite at −2.5 mm/year (Figure 8a), although both zones show velocities within the range of stability. The central zone is built onto volcanoclasts (Figure 8c), and the southern zone onto limestone (−2.9 mm/year). The southern zone is moving away from the satellite with relatively low VLOS. Overall, it appears that the majority of the urban area footprint (Figure 8d) is not affected by significant deformation, while somehow, greater rates are found in the suburban zones, as revealed by the DS targets.  The detected seasonal fluctuations at Mavri Petra are significant, both for PS with LOS velocity close to 0 mm/year and for PS with higher LOS velocity. Those displacements have an amplitude of approximately 5 mm and periodicity of 1 year (Figure 10b). Similar amplitudes are found at Vathi and Paleokastro (Figure 9b). The peaks of the fluctuations generally occur in the period October-March, i.e., rainfall season [62]. This suggests the presence of a seasonal component in the deformation behaviour that is potentially influenced by rainfall and overlaps onto the long-term, generally linear, trend observed across the peninsula in 2015-2019. In this regard, seasonal signals are commonly observed in surface deformation records in other sites worldwide (e.g., [63]), with annual or even semi-annual periods and different amplitudes from negligible to significant, depending on their driving causes. Future ground deformation investigations in Methana could be focused on identification and removal of such periodical terms via advanced time series processing. This operation would better enhance the long-term deformation behaviour of the observed sites which, in this case, may be partly concealed by these fluctuations.
Higher LOS velocities (e.g., −5.0 to −10.0 mm/year) have been observed along steep valleys, mainly on slopes with east orientation (Figure 11a), where visibility of the ascending mode LOS eases the detection of downslope mass movements. The observed displacement could therefore be related to extremely to very slow mass movements and slope instability affecting rugged terrain and volcanic landforms that are recognized as prone to landslides and rock-falls [18]. Some seasonal fluctuations are also detected in the time series (Figure 11b).

of 24
The deformation scenario at Vathi and Paleokastro (Figure 9a) shows that there is a slight movement away from the satellite also in this sector of the peninsula. Similar to what was observed for Methana town, PS targets in the region indicates that the urban areas have a mean value of −0.9 mm/year, while the DS targets for the suburban and more rural zones show −3.1 mm/year.  Figure  8b). However, in the latter, this happens with less prominent amplitudes.
The detected seasonal fluctuations at Mavri Petra are significant, both for PS with LOS velocity close to 0 mm/year and for PS with higher LOS velocity. Those displacements have an amplitude of approximately 5 mm and periodicity of 1 year (Figure 10b). Similar amplitudes are found at Vathi and Paleokastro (Figure 9b). The peaks of the fluctuations generally occur in the period October-March, i.e., rainfall season [62]. This suggests the presence of a seasonal component in the deformation behaviour that is potentially influenced by rainfall and overlaps onto the long-term, generally linear, trend observed across the peninsula in 2015-2019. In this regard, seasonal signals are commonly observed in surface deformation records in other sites worldwide (e.g., [63]), with annual or even semi-annual periods and different amplitudes from negligible to significant, When the evidence from the MT-InSAR investigation is integrated with the outcomes of the GNSS data analysis (see Section 4.1), it can be suggested that the ground motions that are observed from 2015 to 2019 are compatible with the anticipated "low" volcanic activity of Methana. This comes out very evidently from the PS dataset and the generally low V LOS . Despite the low point density and spatial coverage of the volcanic edifice, the overall spatial distribution of the PS deformation values does not reveal a deformation field that could be reliably attributed to a typical volcanic inflation/deflation dynamic of the whole Methana, as found in the literature for other active volcanic areas. On the other hand, the DS dataset may suggest the presence of more widespread deformation patterns along the eastern flank. However, there are only few cases where the observed motions associate with specific landforms (e.g., narrow valleys, erosion gullies and superficial slides). In this regard, the photo interpretation of the slopes through Google Earth optical imagery did not highlight specific circumstances worth further investigation.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 16 of 24 depending on their driving causes. Future ground deformation investigations in Methana could be focused on identification and removal of such periodical terms via advanced time series processing. This operation would better enhance the long-term deformation behaviour of the observed sites which, in this case, may be partly concealed by these fluctuations. Higher LOS velocities (e.g., −5.0 to −10.0 mm/year) have been observed along steep valleys, mainly on slopes with east orientation (Figure 11a), where visibility of the ascending mode LOS eases the detection of downslope mass movements. The observed displacement could therefore be related to extremely to very slow mass movements and slope instability affecting rugged terrain and volcanic landforms that are recognized as prone to landslides and rock-falls [18]. Some seasonal fluctuations are also detected in the time series (Figure 11b).
When the evidence from the MT-InSAR investigation is integrated with the outcomes of the GNSS data analysis (see Section 4.1), it can be suggested that the ground motions that are observed from 2015 to 2019 are compatible with the anticipated "low" volcanic activity of Methana. This comes out very evidently from the PS dataset and the generally low VLOS. Despite the low point density and spatial coverage of the volcanic edifice, the overall spatial distribution of the PS deformation values does not reveal a deformation field that could be reliably attributed to a typical volcanic inflation/deflation dynamic of the whole Methana, as found in the literature for other active volcanic areas. On the other hand, the DS dataset may suggest the presence of more widespread deformation patterns along the eastern flank. However, there are only few cases where the observed motions associate with specific landforms (e.g., narrow valleys, erosion gullies and  Due to the lack of strong seismic events (Mw > 4) or intense seismicity in the area of Methana peninsula (see Section 4.2), the observed ground deformation could not be associated with the recorded seismicity. There are active seismic zones in the vicinity of the peninsula; however, neither the location, the depth nor the magnitude of these events could be directly associated with the observed motions. Even if a seismic contribution to triggering superficial mass movements was present, this would be plausibly limited to nearby seismic swarms. However, there is no hard evidence from ground observations to support further such a hypothesis.
Accounting for the geothermal activity and warm springs at various locations in Methana, the Due to the lack of strong seismic events (Mw > 4) or intense seismicity in the area of Methana peninsula (see Section 4.2), the observed ground deformation could not be associated with the recorded seismicity. There are active seismic zones in the vicinity of the peninsula; however, neither the location, the depth nor the magnitude of these events could be directly associated with the observed motions. Even if a seismic contribution to triggering superficial mass movements was present, this would be plausibly limited to nearby seismic swarms. However, there is no hard evidence from ground observations to support further such a hypothesis.
Accounting for the geothermal activity and warm springs at various locations in Methana, the only situation where motions were found close to a thermal spring is in the immediate vicinity of the open-air thermal baths of the Thermal Spa in the southern part of Methana town (see Figure 8a). However, the LOS deformation rate values found in the PS dataset most likely associate with rock-fall accumulation overlooking the road (Figure 12), while those in the DS dataset cannot be easily separated from the wider deformation pattern extending across the above topography (see Figure 11). The seasonality observed in the MT-InSAR time series in this area as well as at many other locations of the peninsula appear to be more plausibly due to soil moisture and groundwater level oscillations, partly controlled by precipitation, as discussed above. On the other hand, another process that could explain such seasonality could be linked with fluctuations in the geothermal reservoir located in the centre of Methana at a depth of approximately 2 km [64], though more investigation would be needed to verify the potential occurrence of such a process.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 24 process that could explain such seasonality could be linked with fluctuations in the geothermal reservoir located in the centre of Methana at a depth of approximately 2 km [64], though more investigation would be needed to verify the potential occurrence of such a process.

Conclusions
In the context of the abundant literature of satellite InSAR studies on the south Aegean volcanic arc and the published geological and geochemical research on Methana, the present study is the first that undertakes a multi-temporal InSAR analysis of the most recent ground motions in the peninsula in the period 2015-2019 and attempts to find correlations with local GPS benchmark, cGNSS and regional seismicity data collected since 2006. Methana is located in the seismically active neotectonic basin constituting the Saronic Gulf area, and recent efforts have been made by Greek

Conclusions
In the context of the abundant literature of satellite InSAR studies on the south Aegean volcanic arc and the published geological and geochemical research on Methana, the present study is the first that undertakes a multi-temporal InSAR analysis of the most recent ground motions in the peninsula in the period 2015-2019 and attempts to find correlations with local GPS benchmark, cGNSS and regional seismicity data collected since 2006. Methana is located in the seismically active neotectonic basin constituting the Saronic Gulf area, and recent efforts have been made by Greek institutions to specifically monitor the peninsula, among which is the nomination of the Permanent Scientific Committee "Greek Volcano Arc Monitoring" by the administrative board of the Earthquake Planning and Protection Organization (EPPO) in June 2020. In this context, this study is timely in starting to put together initial geodetic and satellite datasets to create a virtual monitoring network that, in the future, could be developed into a proper volcano observatory system.
Based on the integrated data analysis, the general conclusion of this research is that the ground motions observed in the Sentinel-1 MT-InSAR ascending mode data are compatible with the anticipated "low" volcanic activity of Methana. Although the volcano does not show clear signs of activity in recent years, it is a hotspot of potential concern given that it is adjacent to very high-density population areas, including Athens and the surrounding major towns. The interpreted low level of activity is clearly the initial baseline above which further investigations need to be conducted and will need to be verified once the whole archive of SAR datasets dating back to the early 1990s is analysed (e.g., ERS-1/2 and ENVISAT) and the ascending geometry of Sentinel-1 made available for this research is complemented with same InSAR analysis of the matching descending dataset acquired over the peninsula since 2014. It is evident that the use of one SAR geometry only limits the extent to which the observed LOS ground motions can be interpreted. Therefore, future research may include processing of SAR data (even at higher spatial resolution) from both ascending and descending geometries and reconstruction of the 3D deformation field as well as longer time series of geodetic data to improve our understanding of the long-term deformation behaviour of the peninsula.
Funding: Remeasurement of the GPS benchmarks MESO and MENO was financed by "HELPOS-Hellenic Plate Observing System" (MIS 5002697), which is implemented under the action "Reinforcement of the Research and Innovation Infrastructure", funded by the Operational Programme "Competitiveness, Entrepreneurship and Innovation" (NSRF 2014-2020) and co-financed by Greece and the European Union (European Regional Development Fund).