Review of Works Combining GNSS and InSAR in Europe

: The Global Navigation Satellite System (GNSS) and Synthetic Aperture Radar Interferometry (InSAR) can be combined to achieve different goals, owing to their main principles. Both enable the collection of information about ground deformation due to the differences of two consequent acquisitions. Their variable applications, even if strictly related to ground deformation and water vapor determination, have encouraged the scientiﬁc community to combine GNSS and InSAR data and their derivable products. In this work, more than 190 scientiﬁc contributions were collected spanning the whole European continent. The spatial and temporal distribution of such studies, as well as the distinction in different ﬁelds of application, were analyzed. Research in Italy, as the most represented nation, with 47 scientiﬁc contributions, has been dedicated to the spatial and temporal distribution of its studied phenomena. The state-of-the-art of the various applications of these two combined techniques can improve the knowledge of the scientiﬁc community and help in the further development of new approaches or additional applications in different ﬁelds. The demonstrated usefulness and versability of the combination of GNSS and InSAR remote sensing techniques for different purposes, as well as the availability of free data, EUREF and GMS (Ground Motion Service), and the possibility of overcoming some limitations of these techniques through their combination suggest an increasingly widespread approach. collapse regulated by lateral magma over-flow [25] of the Bàrdarbunga stratovolcano in the South-eastern sector of the Iceland island; and a conference paper combining InSAR and GNSS data [74] that focuses on ground deformation measurements of the Northern Volcanic Zone of the island. In addition, Wadge et al. [17] published an article on Mt. Etna focused on the measurement of the tropospheric water value to better consider the contribution of the atmosphere to the differential radar interferograms of ERS1/2. Atmosphere applications (10.4%), unlike subsidence, landslide activity, and tectonic and volcanic phenomena, generally involve wide areas, ranging from regional to national scales. Atmospheric applications of the combination of GNSS and SAR data encompass


Introduction
The Global Navigation Satellite System (GNSS) and Synthetic Aperture Radar Interferometry (InSAR) can be used for various applications owing to their principal functions. GPS enables the collection of three-dimensional (3D) deformation information, considering the coordinates of surveyed points. Since the 1970s, InSAR [1] has permitted the investigation of the phase difference (interferogram) between two images acquired over the same area at different times.
The GNSS term refers to any satellite constellation that provides information globally about the positioning, timing, and navigation data of elements to a GNSS receiver that then commutes this information to a precise location. The GNSS family includes several regional satellite systems, e.g., GPS (Global Positioning System), Galileo, and GLONASS (GLObal'naya NAvigatsionnaya Sputnikovaya System). GNSS performances are (i) accurate in the determination of real position, velocity and time of travel between the transmitter and receiver; (ii) intact, allowing the insertion of an alarm in case of anomaly positioning; (iii) able to work continuously without interruptions; and (iv) able to fulfil all the criteria of accuracy, integrity, and continuity. This system provides three-dimensional along the East-West, North-South, and up-down direction for ongoing deformation through the precise determination of the receiver position with millimetric accuracy in order to assess the main components. GNSS is very useful for investigating regionally low spatial frequencies, e.g., tectonic or geodynamic deformation fields, but local patterns or movements cannot be appreciated due to the large distance between adjacent GNSS stations. The e.g., 250 km with Sentinel-1 [33], the atmospheric effects combined with the reported high precision of the available InSAR approaches could be challenging, especially if the available datasets of SAR images are limited [34,35].
This review focuses on the collection and analysis of scientific contributions that combined GNSS/GPS data and InSAR approaches, i.e., both DInSAR and MTInSAR techniques, for pursuing different perspectives, such as ground deformation recognition and analysis or atmospheric contribution assessment over the entire European continent, excluding Turkey and Russia. In addition, a detailed investigation on Italy, the country most represented by scientific contributions combining GNSS and InSAR data, was conducted. The state-of-theart of the various applications of these two techniques can improve the knowledge of the scientific community and help in the further development of new approaches or additional applications in different fields.
In this review, the DInSAR and MTInSAR approaches were both considered to cover a wider range of field applications, to both rapid and slow phenomena, but analyzing the scientific contribution separately. In contrast, GNSS, GPS, cGPS (continuous GPS), and DGPS (Differential GPS) are considered unique approaches for avoiding discrimination with low statistical soundness considering that the basic approach is the same.

Data Collection
Data collection was conducted through Google Scholar's freely accessible web search engine to gather a high number of scientific contributions published that considered both InSAR and GNSS/GPS data.
To collect the scientific articles, book chapters and conference abstracts published on the Google Scholar's search engine were based on the following criteria. The first keyword adopted was the name of each European country to which the following two lists of keywords related to GNSS and SAR were connected. For the SAR group, the keywords used were "InSAR" (and related acronyms as "MT-InSAR", Multi Temporal InSAR, "DInSAR", Differential InSAR, or "A-DInSAR", Advanced DInSAR), "PS" (Persistent Scatterers), and "PSI" (Persistent Scatterers Interferometry). For the global absolute measurements, the terms "GNSS" and "GPS" (with related derived names as cGPS) were used.
The use of Google Scholar for this type of research review is considered an appropriate choice, as this search engine can be considered "essentially a superset of WoS and Scopus, with substantial extra coverage" [36]. The same author states that, in Google Scholar, approximately 95% of the available citations correspond with those contained in the Web of Science (WOS) database and approximately 92% correspond with those of Scopus in addition to a relevant number of other citations. For this reason, Google Scholar is considered an adapted search engine and a powerful tool for researching the existing literature considering the topic of the combination of GNSS and InSAR data. However, the source of the contributions highlighted by research on Google Scholar has to be validated since this system is less automatized than others, and its database can be completed with information added by the users [37].
The collected scientific contributions, including articles, reviews, book chapters and congress proceedings, were arranged in different folders based on the various searched countries. Then, the contributions were examined iteratively and catalogued in a table considering the following derivable information: -Title of the scientific contribution; -List of authors; -Type of submission (journal, book chapter, conference paper); - Year of the publication; -Country; -Area of interest (province, municipality, local toponym); -Scale of the presented work (national, regional local); -Field of application (i.e., Atmosphere, Glacial; Infrastructure, Landslide, Subsidence; Subsidence/Uplift; Tectonic; Uplift; Vulcanic, Other, and Not Specified); How the two data were used (combination, comparison, validation).

Spatial Distribution
A total of 191 scientific contributions, among which 12 were about the theorical combination of GNSS and InSAR data and 179 focused on the combined techniques, were collected and analyzed considering all European mainlands. Spatially, not all countries have been the subject of combined GNSS and InSAR studies (considering DInSAR or MTInSAR approaches). In fact, 15 European nations, in addition to Gibraltar, Andorra, San Marino, and Vatican City, are not represented by publications regarding the GNSS and InSAR data combination ( Figure 1). The higher number of publications have focused on Italy, with 47 contributions divided into 36 journal articles (e.g., [17,27,[38][39][40][41][42][43]), 4 book chapters [44][45][46], and 7 conference proceedings [47][48][49][50][51][52]. This number of contributions, which is much higher than that collected for all the other countries, can be justified when considering the high number of geohazards affecting the territory, e.g., landslide events, subsidence phenomena, geothermal processes and mining activities, and the higher use of InSAR within the Italian scientific community.
Following the classification, Netherlands is represented by 8 scientific works equally divided into international journal articles and conference papers. Switzerland and Romania provided 7 scientific contributions, 4 peer-reviewed articles, and 3 meeting proceedings, and 5 publications in international journals and 2 extended abstracts, respectively. The United Kingdom is named in 6 scientific works, of which 1 involves the Scotland territory, divided into 4 international peer-reviewed journals and 2 conference proceedings, while France only includes 3 scientific international articles and 1 conference abstract. and Malta (two scientific peer-reviewed journals) were represented by 2 scientific contributions, while Liechtenstein, Hungary, Czech Republic, and Cyprus were closely ranked, with only one scientific work.

Temporal Distribution
The first scientific contribution outlining the combined use of GNSS and InSAR data dates back to 1999 [42]. This article focused on a model used to estimate the geometry and slip distribution of the fault plane originating from the 26 September 1997, Colfiorito earthquakes (central Italy) by combining SAR interferometry and GPS measurements. In the following years, since 2000, the scientific community maintained a greater interest in data integration, also involving other European countries, such as Iceland, Spain, and Greece. Since 2009, excluding 2008, with a high number of 8 scientific contributions, a slight increment of scientific contributions was identifiable until 2012, with more than 10 articles, books or conference proceedings. The maximum number of scientific contributions was reached in 2015, with 29 contributions from 14 different countries. Next, in 2018 and 2012, a total of 20 and 17 works were published, respectively. It is worth noting that Italian case studies were published annually, ranging from the first application in 1999 to 2019, only excluding 2009, 2016, and 2020 ( Figure 2). Albania (1 article published in international journal and 1 congress abstract), Lithuania (1 publication in scientific per-reviewed journal and 1 congress abstract), and Slovakia and Malta (two scientific peer-reviewed journals) were represented by 2 scientific contributions, while Liechtenstein, Hungary, Czech Republic, and Cyprus were closely ranked, with only one scientific work.

Temporal Distribution
The first scientific contribution outlining the combined use of GNSS and InSAR data dates back to 1999 [42]. This article focused on a model used to estimate the geometry and slip distribution of the fault plane originating from the 26 September 1997, Colfiorito earthquakes (central Italy) by combining SAR interferometry and GPS measurements. In the following years, since 2000, the scientific community maintained a greater interest in data integration, also involving other European countries, such as Iceland, Spain, and Greece. The analysis of the temporal evolution in relation to the field of application of scientific contributions ( Figure 3) revealed that early studies were focused on the analysis of tectonic-induced deformations. The rise of an anticline beneath an urban area in Catania in 2001 [87] and the modeling of the dislocation of the edifices involved in the 1997 Umbria-Marche seismic sequence [88] were analyzed by integrating GNSS and InSAR data. In 2001, the first atmospheric investigation application was presented, showing primary results regarding the calibration of atmospheric effects on SAR interferograms by GPS and local atmosphere models [38]. In 2002, the joint application of GPS and DInSAR for volcano and seismic monitoring was carried out in different cases in Spain for the 1992-2002 period [55]. In 2005, the first application of the combined use of GNSS and SAR data for recognizing subsidence was published. In particular, five different measurement techniques, among which are the DGPS and MTInSAR approaches, were exploited for mapping and monitoring land subsidence in the Venice lagoon [43]. A year later, the first application to landslide studies was published dealing with the Triesenberg-Triesen landslide in Leichteinstein, where GPS monitoring data were compared with 3D models generated using an MTInSAR approach on ERS datasets [89]. Among the less represented fields of application, in 2002, the first application to glacial research was proposed, where an interferogram was generated to produce an ice surface motion map of the Gjàlp volcano in Vatnajokull, Iceland [29]. Regarding infrastructure and uplift case studies, the first key scientific contribution focused on the monitoring of the urban area of Bratislava, Slovakia [26] and on the investigation of an area close to the geothermal site of Landau, Germany [19,63].
In recent years, the combined use of GNSS and SAR data has been adopted for largescale analyses. Indeed, ground deformation maps have been produced at national scale in the United Kingdom [90] and at continental scale covering all of Europe [91]. In this case, GNSS data were implemented to identify and filter out possible residual atmospheric artefacts that may affect the quality of the employed MTInSAR data. The analysis of the temporal evolution in relation to the field of application of scientific contributions ( Figure 3) revealed that early studies were focused on the analysis of tectonicinduced deformations. The rise of an anticline beneath an urban area in Catania in 2001 [87] and the modeling of the dislocation of the edifices involved in the 1997 Umbria-Marche seismic sequence [88] were analyzed by integrating GNSS and InSAR data. In 2001, the first atmospheric investigation application was presented, showing primary results regarding the calibration of atmospheric effects on SAR interferograms by GPS and local atmosphere models [38]. In 2002, the joint application of GPS and DInSAR for volcano and seismic monitoring was carried out in different cases in Spain for the 1992-2002 period [55]. In 2005, the first application of the combined use of GNSS and SAR data for recognizing subsidence was published. In particular, five different measurement techniques, among which are the DGPS and MTInSAR approaches, were exploited for mapping and monitoring land subsidence in the Venice lagoon [43]. A year later, the first application to landslide studies was published dealing with the Triesenberg-Triesen landslide in Leichteinstein, where GPS monitoring data were compared with 3D models generated using an MTInSAR approach on ERS datasets [89]. Among the less represented fields of application, in 2002, the first application to glacial research was proposed, where an interferogram was generated to produce an ice surface motion map of the Gjàlp volcano in Vatnajokull, Iceland [29]. Regarding infrastructure and uplift case studies, the first key scientific contribution focused on the monitoring of the urban area of Bratislava, Slovakia [26] and on the investigation of an area close to the geothermal site of Landau, Germany [19,63].
In recent years, the combined use of GNSS and SAR data has been adopted for largescale analyses. Indeed, ground deformation maps have been produced at national scale in the United Kingdom [90] and at continental scale covering all of Europe [91]. In this case, GNSS data were implemented to identify and filter out possible residual atmospheric artefacts that may affect the quality of the employed MTInSAR data.

Field of Application
The joint application of SAR and GNSS data has been adopted in 9 different fields, classified as follows: (i) Tectonic; (ii) Atmosphere; (iii) Glacial; (iv) Volcanic; (v) Subsidence; (vi) Landslide; (vii) Infrastructure; (viii) Subsidence/Uplift; and ix) Uplift. In addition, another class has been added to include all nonconventional applications, as well as a Not Specified class (N.S.) for applications in which it was not possible to assign a specific field. A higher percentage of publications deal with subsidence phenomena (20.2%), followed by landslide events (19.7%) and tectonic analyses or seismic event back analyses (18.6%) ( Figure 4).
Subsidence phenomena are spatially distributed in Europe, with a higher concentration of studies (i) in the Upper Rhine Graben, at the boundary between France and Germany [62,[92][93][94]; (ii) in Italy, over the Pianura Padana Plain (central Italy), in the Emilia Romagna Region [44,47,48] and the areas of the Po Plain [20,45,49,95,96], in the Venice lagoon [41,43,97,98], in the Friuli Venezia Giulia coastal plain (North-East of Italy) [99], the Firenze-Prato-Pistoia basin (Tuscany Region, central Italy) [100], and the Sibari plain (Calabria Region, South Italy) [101]; and (iii) in Spain, the Alto Guadalentín [102,103] (South-East of Spain), in Lorca town [54] and in the Cardona salt mine [60] (North-East of Spain). In Portugal, instead, only two scientific contributions have been reported in connection to the subsiding area of Lisbon [70,104]. Eastern Europe shows a widespread distribution of land subsidence issues, which have been analyzed by combining GNSS and InSAR data, in the coastal portion of Lithuania, close to Kaszuby Lakeland [105], in the agglomeration of Warsaw, Poland [106], in the Bucharest metropolitan area in Romania [107,108], and close to the Tatra Mountains in Slovakia [109]. A specific focus on mininginduced land subsidence events and collapses was provided by research which studied the Upper Silesia Coal basin, covering a part of the territory of Poland [110] and the Czech Republic [111]. North Europe is represented only by the surface motion of peatland in Flow County, Scotland [112], and the localized ground subsidence in the urban settings of London [113] for the United Kingdom and the analysis in the prediction of long-term settlement in the Trondheim harbor, mid-Norway [76].

Field of Application
The joint application of SAR and GNSS data has been adopted in 9 different fields, classified as follows: (i) Tectonic; (ii) Atmosphere; (iii) Glacial; (iv) Volcanic; (v) Subsidence; (vi) Landslide; (vii) Infrastructure; (viii) Subsidence/Uplift; and ix) Uplift. In addition, another class has been added to include all nonconventional applications, as well as a Not Specified class (N.S.) for applications in which it was not possible to assign a specific field. A higher percentage of publications deal with subsidence phenomena (20.2%), followed by landslide events (19.7%) and tectonic analyses or seismic event back analyses (18.6%) ( Figure 4).
Subsidence phenomena are spatially distributed in Europe, with a higher concentration of studies (i) in the Upper Rhine Graben, at the boundary between France and Germany [62,[92][93][94]; (ii) in Italy, over the Pianura Padana Plain (central Italy), in the Emilia Romagna Region [44,47,48] and the areas of the Po Plain [20,45,49,95,96], in the Venice lagoon [41,43,97,98], in the Friuli Venezia Giulia coastal plain (North-East of Italy) [99], the Firenze-Prato-Pistoia basin (Tuscany Region, central Italy) [100], and the Sibari plain (Calabria Region, South Italy) [101]; and (iii) in Spain, the Alto Guadalentín [102,103] (South-East of Spain), in Lorca town [54] and in the Cardona salt mine [60] (North-East of Spain). In Portugal, instead, only two scientific contributions have been reported in connection to the subsiding area of Lisbon [70,104]. Eastern Europe shows a widespread distribution of land subsidence issues, which have been analyzed by combining GNSS and InSAR data, in the coastal portion of Lithuania, close to Kaszuby Lakeland [105], in the agglomeration of Warsaw, Poland [106], in the Bucharest metropolitan area in Romania [107,108], and close to the Tatra Mountains in Slovakia [109]. A specific focus on mining-induced land subsidence events and collapses was provided by research which studied the Upper Silesia Coal basin, covering a part of the territory of Poland [110] and the Czech Republic [111]. North Europe is represented only by the surface motion of peatland in Flow County, Scotland [112], and the localized ground subsidence in the urban settings of London [113] for the United Kingdom and the analysis in the prediction of long-term settlement in the Trondheim harbor, mid-Norway [76].
of an early-warning system for marine geohazards along the Black Sea coast [153]. Furthermore, Zoran [154] combined geospatial information, i.e., GPS and SAR images, with in situ information to assess the seismic hazard associated with the Vrancea area in Romania. Other interesting research in the field of tectonism utilized earthquake observations with remote sensing data from the Mw 6.4 Durres (Albania) earthquake that occurred in 2019 [155] and the Mw 6.5 Lorca earthquake (Spain) [56]. Studies on seismic [156] and crustal deformation [157], sometimes connected to strain rate analysis [66], or interactions with volcanic systems [68,69], connected to earthquake seismic swarms or active tectonic faults, are more common. Volcanic (10.9%) applications are mainly concentrated in (i) Greece, focused on the investigation of Santorini volcano [158,159] and Nisyros volcano [83,86]; (ii) Italy, focused on the evaluation of the deformation of Mt. Etna, Sicily Region (Southern Italy) [160,161] and the Campi Flegrei, Campania Region (Southern Italy) [40,51]; (iii) Portugal, focused on the assessment of the geodynamics and displacement mechanisms affecting the Azores islands [13,14,69]; and (iv) Spain [55], focused on the volcanic deformation of the Canary archipelago [58] and Tenerife island [24,57,162,163], inevitably depending on the presence of active volcanoes. Iceland is represented by two scientific articles, one focusing on the lateral dyke grown [72] and the gradual caldera collapse regulated by lateral magma overflow [25] of the Bàrdarbunga stratovolcano in the South-eastern sector of the Iceland island; and a conference paper combining InSAR and GNSS data [74] that focuses on ground deformation measurements of the Northern Volcanic Zone of the island. In addition, Wadge et al. [17] published an article on Mt. Etna focused on the measurement of the tropospheric water value to better consider the contribution of the atmosphere to the differential radar interferograms of ERS1/2.
Atmosphere applications (10.4%), unlike subsidence, landslide activity, and tectonic and volcanic phenomena, generally involve wide areas, ranging from regional to national scales. Atmospheric applications of the combination of GNSS and SAR data encompass A more widespread distribution of scientific contributions from landslide studies with the joint application of GNSS and SAR data can be observed in Europe. A total of 14 nations have at least one reported application. The country with the highest number of contributions is Norway, represented by the cases of the Åknes landslide [114], the Gamanjunni rock slope [78], the Osmundneset landslide in Sognog Fjordane [115], the periglacial landscape activity at Nordnesfjellet [116], and the mapping [79] and monitoring [80] of landslides in North and West Norway. Norway is followed by Switzerland, with 6 contributions, including 3 peer reviewed articles in international journals and 3 proceedings, which focused on rock glacier detection and monitoring [117,118], slope instability monitoring events in the Western Alps [119,120], and deep-seated landslide recognition by Osco [121] and Aletsch Forest [122]. The third most represented country is Italy, with 5 peer-reviewed articles and a proceeding focus on specific case studies of the Bosmatto landslides in the Northwestern Alps [39], investigations over the wider area of the Northern Apennines [123], and of the Assisi landslides [23,52], in addition to a quantitative hazard and risk assessment of slow-moving landslides in the Arno Basin in Tuscany [124]. Next, there are 4 contributions associated with Spain concerning the monitoring of the Vallcebre landslide [125], the Tena Valley [126,127], which was affected by very slow landslides, and the Portalet landslide area [53]. In Slovenia, 3 works have been dedicated to the deformation monitoring of the Potoška Planina [128][129][130], while 2 peer-reviewed publications set in Malta deal with slow-moving coastal landslides [131] and landslide susceptibility modeling [132] of the Northwestern coast. Scientific contributions are also found investigating the Triesenberg-Triesen landslide in Liechtenstein [89], the Trifon Zarezan landslide in Bulgary [133], the ground deformation affecting the Choirokoitia UNESCO World Heritage site in Cyprus [134], the La Vallette landslide in France [135], the deformation phenomena that occurred in the Tatra Mountains in Slovakia [109], the ground deformation affecting Kulcs village in Hungary [136], and the investigation of the Grande da Pipa river basin in Portugal [67].
The analysis of tectonic movements (including pre-, sin-, and post-event movements) via GNSS and SAR is mainly concentrated in the seismic areas of Europe. The scientific contributions mainly come from Italy, with 16 publications, Greece and Iceland, with 5, and Romania, with 4. Italy constitutes a large representation due to the various seismic sequences that have occurred within the last 20 years, such as the 2009 L'Aquila earthquake [46,137,138], the 2012 Emilia Romagna earthquake [139], and the 2016 Central Italy Remote Sens. 2021, 13, 1684 9 of 27 seismic swarm [21,50,[140][141][142][143], in addition to the old 1997 Colfiorito earthquake [42,88]. Furthermore, in Italy, scientific publications on the investigation of the geodynamic processes active throughout the whole Italian territory [20] or focused on the Southeastern Po Plain [144] and the geodynamic movements of the Mount Etna [87] or the strain rate of the Hyblean Plateau [145] in Sicily were developed by exploiting a combination of GNSSand InSAR-derived information. The scientific literature devoted to Greece has reported interesting research on the investigation of pre-earthquake deformation processes by examining GPS, DInSAR, MTInSAR, and ancillary data collected over Zakynthos Island [146]. In addition, the seismological analysis of the critical estimation of a future strong seismic event in the broader area of Cephalonia [82], the post-seismic investigation of a seismic sequence that occurred during the period of January-February 2014 [147] or a specific earthquake, such as the Lefkada involving the Cephalonia Transformation Fault (CTF), close to Cephalonia Island [81], and the post-seismic motion analysis of one of the most active extending grabens, the Corinthian Gulf [85], were collected. Iceland is represented by five contributions. One is the investigation of a M w (moment magnitude) 6.5 earthquake that occurred in 2000 [71,148]. The others are on the identification of extensional and interseismic deformation recorded in the Northern Vulcanic Zone [149], as well as the investigation of ground deformation close to geothermal power plants in Hengill due to the interaction of regional tectonics and volcanic deformation [150] and the post-seismic viscous relaxation in the Southern Iceland [151]. In the South Carpatians, Romania, crustal deformation was investigated by detecting small-magnitude tectonic processes [152], and GNSS data, combined with SAR interferometry, were adopted for the analysis of the abnormal seismic behaviors of Izvoarele-Galat , i and for the development of an early-warning system for marine geohazards along the Black Sea coast [153]. Furthermore, Zoran [154] combined geospatial information, i.e., GPS and SAR images, with in situ information to assess the seismic hazard associated with the Vrancea area in Romania. Other interesting research in the field of tectonism utilized earthquake observations with remote sensing data from the M w 6.4 Durres (Albania) earthquake that occurred in 2019 [155] and the M w 6.5 Lorca earthquake (Spain) [56]. Studies on seismic [156] and crustal deformation [157], sometimes connected to strain rate analysis [66], or interactions with volcanic systems [68,69], connected to earthquake seismic swarms or active tectonic faults, are more common.
Volcanic (10.9%) applications are mainly concentrated in (i) Greece, focused on the investigation of Santorini volcano [158,159] and Nisyros volcano [83,86]; (ii) Italy, focused on the evaluation of the deformation of Mt. Etna, Sicily Region (Southern Italy) [160,161] and the Campi Flegrei, Campania Region (Southern Italy) [40,51]; (iii) Portugal, focused on the assessment of the geodynamics and displacement mechanisms affecting the Azores islands [13,14,69]; and (iv) Spain [55], focused on the volcanic deformation of the Canary archipelago [58] and Tenerife island [24,57,162,163], inevitably depending on the presence of active volcanoes. Iceland is represented by two scientific articles, one focusing on the lateral dyke grown [72] and the gradual caldera collapse regulated by lateral magma overflow [25] of the Bàrdarbunga stratovolcano in the South-eastern sector of the Iceland island; and a conference paper combining InSAR and GNSS data [74] that focuses on ground deformation measurements of the Northern Volcanic Zone of the island. In addition, Wadge et al. [17] published an article on Mt. Etna focused on the measurement of the tropospheric water value to better consider the contribution of the atmosphere to the differential radar interferograms of ERS1/2.
The Subsidence/Uplift class (4.4%) includes surface deformation analyses performed at a large, regional or national scale, in which both subsidence and uplift phenomena can be identified. Case studies were identified in Netherlands, such as surface deformation analyses performed over the whole national territory [176][177][178][179] or monitoring works of the water defence structure of Waddenzee and IJsselmeer [180,181]. Two other scientific contributions focused on the Polish and Lithuanian Baltic coastal areas [105] and the surface displacement pattern affecting the permafrost areas in northern Norway [75].
The categories of Uplift, Infrastructure, and Glacial are poorly represented. In fact, three peer-reviewed articles, focused on the Tower Hamlets Council area of London [182] and on the Geothermal site Landau in Germany [19,63], which only investigated uplift phenomena resulting from the swelling of clayey soil and geothermal consequences, respectively, can be cited. The lower number of studies focused on the investigation of infrastructure can be attributed to several factors, e.g., the resolution of SAR images, the distribution of the GNSS network compared to the infrastructure dimensions, the particularity of the application and the difficulty associated with interpreting the results. An investigation of the nonlinear deformation of infrastructure via SAR data and GNSS comparison was performed in Bratislava, Slovakia [26]. Other infrastructure applications were indirect, including the monitoring of subway construction in Bucharest, Romania [183], ground deformation detection focused on a new geothermal power plant at Reykjanes, Iceland [73], observation of ground deformation associated with the Kozloduy Nuclear Power Plants (NPPs), Bulgary [184], and experimental monitoring of localized deformation on the Roman aqueducts in Rome, Italy [27]. Glacial applications are represented by specific case studies in Iceland, which focused on the 3D surface motion of the glacier surface of the Gjàlp volcano in Vatnajokull [28,29]; France, with a time-series measure of the Argentière glacier on the Mont Blanc Massif [185]; and northern Scotland (United Kingdom) through the investigation of the Glacial Isostatic Adjustment (GIA) [186].
In addition, some very particular studies of the combined use of GNSS and SAR data, categorized as Other (Figure 4), recorded the following: (i) sea level change detection along the coastline of Brest, France [187], along the North Sea and Baltic Sea Coast in Germany [188] and in the Northern Mediterranean Sea [189]; (ii) flood inundation modeling of Northamptonshire [190]; (iii) mine collapse monitoring in Rudna Mine, Poland [191]; (iv) surface and roughness and sediment texture characterization [192]; and (v) a new approach description of precise datum collection [193].

Satellite SAR Sensors
Considering all the aforementioned scientific contributions, whether they are peerreviewed articles, book chapters or conference papers, 32% of them exploit DInSARbased data, while the remaining 68% relied on one of the previously described MTInSAR approaches. In detail, more than 100 publications use PSInSAR (Persistent Scatterers Interferometry SAR) data [194], while a minor portion utilized the SPN (Stable Point Network [195]), SqueeSAR [196], SBAS (Small BAseline Subset [197]), iSBAS (Intermittent SBAS [198]), PSP-DINFSAR (Persistent Scatterer Pairs Differential InSAR [199]), CPT (Coherent Pixel Technique [200]), SPINUA (Stable Point Interferometry over Unurbanized Areas [201]), and TSIA (Two-Scale Interferometric Analysis, processing chain that performs a sequence of low-resolution (small-scale) and full-resolution (large-scale) processing [23], based on the SBAS [197] approach) algorithms. It is impossible to establish a precise number of times each algorithm was used since, in some articles, more than one approach was used for different datasets, and in many contributions, the processing algorithm was not clearly specified.
Considering the satellite wavelength of the whole database collected, approximately 150 C-band datasets were implemented, as well as 35 X-band and 15 L-band datasets. For 15 papers, the C-and X-band datasets were combined, while only 7 C-and L-band datasets were used in parallel. None of the scientific studied viewed combined X-and L-bands. Finally, applications analyzing all three bands were implemented 5 times. Publications taken into consideration were also categorized according to the satellite constellation of the SAR data. Most of the paper shows a massive application of ENVISAT data (31.1%), covering the period 2002-2010, and ERS data (27.2%), covering the period from 1992 to 2000, followed by the data of Sentinel-1 (14.2%), acquiring SAR images since April 2014, which was largely diffused over the last few years due to its short and constant revisiting time, due to its scientific contribution aim and free availability [33]. All these satellite constellations acquire in the C-band (5.6 cm wavelength), confirming its more widespread use in more common applications as a good compromise for urban and nonurbanized areas. The minor frequencies of other satellites are registered, listed from the high to low: TerraSAR-X (X-band, approximately 10%), RadarSAT (X-band, approximately 6%), Alos (L-band, approximately 6%), and COSMO-SkyMed (X-band, approximately 4.7%).
Satellite data are also used in combination with others to obtain longer timespans of analysis or to compare the performances of different bands. The combination of different constellation datasets is reported in the satellite graphs of Figure 5. It is interesting to note that, in at least one paper, all constellation satellites were combined with the other three constellations, except for RadarSAT, which did not show application when combined with the other 3 datasets. In addition, it is possible to state that the aforementioned constellation satellites are rarely combined with the other two and commonly with other datasets. Exceptions are also detectable for the COSMO-SkyMed, X-band, and Alos, L-band, whose data were more often used combined with two other constellations than standalone or only with one other dataset.

Temporal and Spatial Distribution in Italy
A more specific analysis of the temporal and spatial distribution of previous studies, as well as within the field of application, was conducted focusing on the Italian literature, since Italy is the most represented country, accounting for the 47 published scientific con-

Temporal and Spatial Distribution in Italy
A more specific analysis of the temporal and spatial distribution of previous studies, as well as within the field of application, was conducted focusing on the Italian literature, since Italy is the most represented country, accounting for the 47 published scientific contributions on the integrated use of GNSS and InSAR data.
The highest number of scientific contributions is recorded in Emilia-Romagna and Umbria, which are two regions where tectonic, seismic, subsidence, and landslide events were analyzed (Figure 6). In the Valle d'Aosta Region, in Northwestern Italy, a large landslide in a high alpine environment was studied by GNSS, InSAR, and GB-InSAR remote sensing data to reduce uncertainties. The joint combination of these techniques provided a comprehensive view of the deformation field of the landslide [39].
Landslide investigations are also identified in the Umbria Region, central Italy, for the analysis and monitoring of the slow-moving Assisi landslide [23,52]. The geodetic and interferometric products of four deep-seated landslides reactivated by the excavation of a double road tunnel that induced deformation were monitored in the Northern Apennine Despite the higher number of articles, books, and conference proceedings with case studies performed over the Italian territory, not all the regions are covered. In fact, six Italian regions do not exhibit any phenomena investigated by either GNSS or SAR data. One article [20] took into consideration the entire Italian territory using both techniques to study Earth surface deformation to provide the displacements with respect to different ground point position components.
In the Valle d'Aosta Region, in Northwestern Italy, a large landslide in a high alpine environment was studied by GNSS, InSAR, and GB-InSAR remote sensing data to reduce uncertainties. The joint combination of these techniques provided a comprehensive view of the deformation field of the landslide [39].
Landslide investigations are also identified in the Umbria Region, central Italy, for the analysis and monitoring of the slow-moving Assisi landslide [23,52]. The geodetic and interferometric products of four deep-seated landslides reactivated by the excavation of a double road tunnel that induced deformation were monitored in the Northern Apennine [123], while the same type of data was used in the Arno River basin, Tuscany Region, to update previous hazards and risk map realized by [124], Catani, et al. [202].
The highest number of applications, 15 in total, is represented by subsidence phenomena investigations. Seven scientific contributions on this theme, i.e., 2 peer-reviewed articles, 2 book chapters, and 3 conference proceedings, target the Emilia-Romagna Region, central Italy, 4 of which are specifically on the Po Plain. Only one is focused on the Tuscany Region, central Italy, more specifically in the Firenze-Prato-Pistoia Plain [100], as one on the Sibari Plain in the Calabria Region, South Italy [101], and on the plain that lies between the Tagliamento and Isonzo Rivers in the Friuli-Venezia Giulia Region, North-East of Italy [99]. The remaining 5 scientific contributions focus on the vertical deformation of the Venice lagoon [41,43,97,98] or, more generally, on the subsidence phenomenon affecting the North Adriatic Sea [45].
The high number of scientific contributions sourced from the Umbria Region, central Italy, resulted from the Assisi landslide analysis and monitoring, as well as the investigation of the earthquakes that occurred in 1997, 2009, and 2016 involving Umbria and the neighboring regions. In total, 16 contributions targeted the investigation of tectonic and seismic investigations.
In contrast to the landslide or subsidence investigations, the seismic analyses typically involved more than one region as earthquakes affect larger regions. For this reason, 4 articles and one conference proceeding focus on more than one region of central Italy. The 2009 L'Aquila earthquake, Abruzzo Region, was also represented in an article and in a book chapter analyzing the coseismic rupture [46] and slip distribution [138] concentrating the study area close to the epicenter. The 2012 Emilia-Romagna Region, instead, did not have relevant effects on the boundary regions. One of these contributions addresses the modeling of the influence of fluids and pore-pressure changes on surface displacements and on the Coulomb Failure Function (CFF) in the co-seismic and post-seismic period related to the 2012 mainshock [139]. Wang et al. [142] looked for the source parameters and triggering link between the earthquake sequences of 2009 and 2016 recorded in central Italy. In addition, 4 works combining InSAR and GNSS data were developed for determining the source parameter of the three main shocks recorded during the 2016 Central Italy earthquake sequence [50,140], for a modern observation of the spatio-temporal evolution of the seismic sequence [143] and for a novel record of near-field co-seismic fault slip measurement [21].
The spatial distribution of 5 volcanic scientific contributions, including 4 journal articles and one symposium proceeding, is strictly related to the location of volcanic cones in Italy. In fact, 3 publications focus on the analysis of nonlinear deformation [51], the evaluation of fault reactivation [203] and the investigation of magma injection [204] of the Campi Flegrei, Campania Region, and two on the ground deformation pattern of the Mt. Etna volcano in Sicily [160,161].
The Mt. Etna volcano is also the subject of two peer-reviewed articles investigating the varied tropospheric compensation expectations between the GPS and SAR interferograms [38] and the dynamic models of atmospheric movement for calculating the delays affecting radar processing [17]. The other area where the atmospheric delay analysis was undertaken was centered over the Como province in Lombardia Region, North Italy, presenting the results in an international scientific journal [18] and in an international symposium [170].
In conclusion, the Italian territory has been a relevant case study since 1999 [42], and it was continuously under investigation by the scientific community until 2019 that analyzed both GNSS and SAR data for different purposes (Figure 7).

Discussion
A substantial number of scientific contributions using remote sensing data derived from GNSS and SAR systems to investigate and analyze ground deformation in European countries, covering a time span of 31 years (from 1999 to 2020), were collected.
The spatial distribution of all the scientific contributions, considering scientific articles, book chapters or conference proceedings, mirrors the distribution of the main relevant geohazards affecting each country. Italy and Spain are well represented due to the high number of geohazards affecting their territory, especially considering the subsidence occurring in both Italy [205] and Spain [9] and the landslides affecting Italy [206]. The limited number of scientific contributions combining GNSS and InSAR data to study landslides is surprising considering the relevance of these natural hazards and the major consequences occurring throughout Europe [207]. Nonetheless, the Italian territory is also affected by frequent seismic swarms and strong earthquakes in addition to the presence of volcanoes, and both phenomena are deeply investigated by the joint use of GNSS and SAR data. In the same way, it is surprising how low the number of applications recorded in the Netherlands is considering the importance of subsidence phenomena in that country (e.g., [208]), which exert social and economic consequences and induce damage to structures and infrastructure [209,210]. Another well-represented nation is Norway, which is strongly affected by landslides, often triggered in quick clay soils [211], and characterized by the presence of permafrost, which can cause four types of landforms [212]: (i) palsas (e.g., [213]), (ii) rock glaciers (e.g., [214,215]), (iii) ice-cored moraines (e.g., [214,216]), and (iv) ice-wedge polygons (e.g., [217]).
In addition to scientific papers dedicated to the analysis of displacement, general articles employing atmospheric analyses, in which GNSS products are implemented to correct interferograms or remove the geodetic influence on InSAR displacement measure-

Discussion
A substantial number of scientific contributions using remote sensing data derived from GNSS and SAR systems to investigate and analyze ground deformation in European countries, covering a time span of 31 years (from 1999 to 2020), were collected.
The spatial distribution of all the scientific contributions, considering scientific articles, book chapters or conference proceedings, mirrors the distribution of the main relevant geohazards affecting each country. Italy and Spain are well represented due to the high number of geohazards affecting their territory, especially considering the subsidence occurring in both Italy [205] and Spain [9] and the landslides affecting Italy [206]. The limited number of scientific contributions combining GNSS and InSAR data to study landslides is surprising considering the relevance of these natural hazards and the major consequences occurring throughout Europe [207]. Nonetheless, the Italian territory is also affected by frequent seismic swarms and strong earthquakes in addition to the presence of volcanoes, and both phenomena are deeply investigated by the joint use of GNSS and SAR data. In the same way, it is surprising how low the number of applications recorded in the Netherlands is considering the importance of subsidence phenomena in that country (e.g., [208]), which exert social and economic consequences and induce damage to structures and infrastructure [209,210]. Another well-represented nation is Norway, which is strongly affected by landslides, often triggered in quick clay soils [211], and characterized by the presence of permafrost, which can cause four types of landforms [212]: (i) palsas (e.g., [213]), (ii) rock glaciers (e.g., [214,215]), (iii) ice-cored moraines (e.g., [214,216]), and (iv) ice-wedge polygons (e.g., [217]).
In addition to scientific papers dedicated to the analysis of displacement, general articles employing atmospheric analyses, in which GNSS products are implemented to correct interferograms or remove the geodetic influence on InSAR displacement measurements over wide areas, were collected and analyzed. This approach is even more relevant due to the nature of the wide areas investigated and monitored by InSAR approaches.
The combination of DInSAR and GNSS data has not always been categorical, but various trends and approaches utilized over the years can be illustrated.
The various works analyzed have shown the potential of combining both data for different purposes. For instance, GNSS data were used to improve the accuracy of the interferograms generated [11], topographically correcting the interferograms and the elevation data to measure ice surface motion maps along the Gjalp Volcano. Most of the works jointly used GNSS and DInSAR data to produce a more precise and reliable analysis of the movements detected by either combining or validating the available data. Indeed, in most of the pioneering works, authors worked to validate DInSAR data, since in most of the cases, the latter were used for the first time and their results needed to be validated. For instance, in Raucoules et al. [135], in which 3D data derived from TerraSAR-X processing over the La Vallette landslide (southern France) were validated with GNSS data; in other cases, such as in Bovenga et al. [23] and Sakkas et al. [146], GNSS data, along with geological and geotechnical data, were used to ground truth C-and X-band imagery.
Approximatively half of the scientific works are focused on single case studies or events with a local influence, while 45% involve regional areas. Temporally, more recent papers involve more large areas due to various factors: (i) the possibility to have more and improved algorithms for processing large SAR datasets; (ii) the technological progress of the last few years, allowing us to have available advanced computational resources, e.g. cloud computing by virtual machine system and parallel or distributed computing architectures used to process larger dataset over wide areas; (iii) the willingness of different radar satellite constellations, also in different bands, i.e., X-band, C-band and L-band; and (iv) the ease of obtaining free SAR images at a medium resolution, such as the Sentinel-1 constellation, by the Copernicus Access Hub of the ESA (European Space Agency), ensuring systematic worldwide covering acquisition since April 2014, or the availability of limited imagery datasets at a high resolution, i.e., X-band, by means of dedicated scientific projects in open calls. Only a few articles, 10 of the 174 examined, have as a main aim the analysis of more case studies in the same nation, within very wide areas, or throughout entire countries. In detail, four articles are noteworthy for the extension and the relevance of the work. The first two focused on the ground deformation occurring in Netherlands [176,178], while the remaining focus on the generation of atmosphere Precipitable Water Vapor (PWV) maps over large areas in Portugal [169] and the feasibility assessment of a national InSAR ground deformation map over the United Kingdom [90], which both took place in 2017. In 2020, the first continental-scale processing of Sentinel-1 imagery was described in Lanari et al. [91], which provided a calibrated ground deformation map of Europe with the velocity of deformation and cumulative displacement corrected by the GNSS network data.
The authorship, taking into account the affiliation of the first author of the scientific contributions on ground deformation or atmosphere contribution analyses over European countries, demonstrates that approximately 80% of the total refers to the same nations. For the remaining approximately 20%, 10 authorships are referred to Countries outside of Europa (e.g., 4 from USA, 5 from China, and 1 from Canada), while the other main authorships are of Italian (8), Portugal (6), and French (5) researchers. This can be justified by considering (i) collaborations or (ii) projects from different universities, institutions, or entities and (iii) the PhD abroad periods or (iv) research abroad time of researchers. Italy, with 9 contributions authored by non-Italian investigators, is the most investigated nation, not only by Italian researchers but also those from European (2) and international (5) investigators. Next, Portugal is unique in that it was only analyzed by Portugal researchers that also investigate Romania (6 contributions to 7 are from Portugal researchers), followed by Spain, with 5 contributions from European researchers.
The correlation between the temporal evaluation of European and Italian studies is interesting. For the first few years, this relation is easily comprehensible considering that the scientific contribution collected for the first three years, i.e., 1999-2000 and 2001, involved the Italian territory. It is interesting to note, comparing Figures 2 and 7, that the trend is very similar, with an incremental increase that progressed until the peak contribution in 2015, which was then followed by a reduction. The high number of scientific contributions recorded over the Italian territory severely influences the trend of the temporal distribution that occurred over Europe. In contrast, the studies published in 2020 do not follow this role; in fact, the last analyzed year exhibits a relevant number of scientific contributions of 2019 even if no publications were published involving Italy. A similar correlation can also be identified when considering the European and Italian temporal distributions and the associated field of application.
Combining the information regarding the year of publication of the scientific contributions collected and the investigated SAR satellite, it is possible to state that the later publications aimed to investigate the evolution of a phenomenon as long as possible. Obviously, this cannot be affirmed for publications involving atmospheric interactions and corrections or for those investigating the co-seismic earthquake effects of DInSAR.
The incremental increase in the number of publications reported demonstrates that the remote sensing techniques investigated in this review, i.e., geodetic and InSAR approaches, were being readily employed by the scientific community as a result of their benefits and the consequences for stakeholders and end-users as environmental managers. In recent years, this interest began to be traduced in the publication of national and regional GMS (Ground Motion Services) due to the free availability of InSAR data and, in some cases, the already GNSS-calibrated model. In fact, a higher number of applications, and thus a greater number of papers, was collected and recorded in those countries in which a Ground Motion Service (GMS) already exists or is nearly ready for presentation. In 2002, Italian law allocated funds for a Not-ordinary Plan of Environmental Remote Sensing aiming at monitoring considered at risk hydro-geomorphological areas, providing nationwide ground deformation maps calibrated by GNSS data. ERS and ENVISAT interferometric products, processed by PSInSAR and PSP-DINSFAR approaches, were published in the first phase of the project (2008-2009) and then updated in a second phase (2010-2011), while COSMO-SkyMed data were implemented in the same project in the third and final phases (2013-2015) [218,219]. Next, the first nation that made available updated GMS data over the entire territory was Norway in 2018; this publication was provided by the WebGIS portal Radasat-2 data (2010-2018) and Sentinel-1 (2015-2019), not calibrated by GNSS information [220]. In 2019, a year later, BodenBewegungsdienst Deutschland presented the Germany GMS, making the Sentinel-1 PS data available, GNSS-calibrated, processed thanks to the PSI approach [221] by the DLR German Aerospace Center research institute. In addition, Danish GMS and Dutch GMS are under development and will be published soon. More locally, Tuscany (central Italy) and Veneto Region (North-East Italy) have produced regional GMSs that are providing Sentinel-1 MTInSAR data free of charge processed by the SqueeSAR [196] algorithm. These three regions are the object of continuous monitoring, which is systematically updated and interpreted every 12 days. In addition, the application of a data-mining algorithms also highlight the PS with relevant trend changes [222,223]. Taking into account the usefulness and versatility demonstrated by this review regarding the combination of GNSS and InSAR data in ground deformation analysis, in the future, PS data will become available, and the national and regional GMS should be GNSS-calibrated. Currently, this procedure is sometimes conducted by researchers as a further step for specific investigations with good results [20,100].
Consequently, the promising results obtained by the scientific community and the free availability of data, which permit drastic cost reductions, have drawn increasing interest from the administrative managing office for the mapping and monitoring of the ground deformation issue [223][224][225].
In addition, since 2016, Europe has been working on an initiative pushing to implement the available data provided by MTInSAR products at continental scale in a project named European Ground Motion Service (EGMS) [226]. The three main products provided by the EGMS are: (i) the ground deformation maps and time-series along the LOS direction; (ii) an advanced ground deformation map in LOS direction combined with the GNSS data of the EUREF permanent network (EPN) [227]; and (iii) the two main deformation components derivable from the combination of the two orbits, i.e., the horizontal East-West and the vertical up-down deformation. All these products will be updated annually and will open several other opportunities for investigation with a consequent advancement in the knowledge in various fields, e.g., ground deformation prediction [228], mapping [224,225,229], and monitoring [230] of natural hazards, and atmospheric delay analysis.
Another use of GNSS data, which has become increasingly widespread over recent years, is aimed at the retrieval of 3D surface displacement, thus overcoming the limitations associated with DInSAR LOS displacement measurements, e.g., References [92,128]. Some limits of the geodetic and InSAR approaches can be overcome by the combination of GNSS and InSAR data, but some limitations cannot be removed. In the latest year, GNSS data were mostly used for completing SAR-derived information and filling in points were data were missing, thus improving the temporal coverage of the deformation time series, as performed in a pervious study [100], or for investigating the "real" ground deformation velocity, thereby remedying the geodetic displacement that affects the ordinary InSAR measurements. Furthermore, the areal coverage of results can be improved by joining GNSS stations and derived-PS points by PSI approaches, which can obtain deformation data over areas invisible to SAR systems, e.g., Reference [39]. It enables the obtainment of information in areas that are not covered by direct measurements or by the use of a single technique. On the other hand, GNSS stations are placed in strategic positions to create an approximatively regular network and cannot completely overpass the InSAR limits, such as the coverage of shadow areas due to the interaction between the topography and the LoS of SAR satellites. To overcome this limitation, a traditional survey with an operator by GPS must be conducted, reducing the easy and fast repeatability and considerably enlarging the time and monetary expenses.
The substantial improvement of both technical and technological approaches in recent decades, as well as the demonstrated useful and versability of the combination of GNSS and InSAR remote sensing techniques for different purposes, by more than 180 scientific contributions suggests that the incremental use of these approaches in the future will be further enhanced and widespread. The free availability of Sentinel-1 images and InSAR products by GMS will help in the increasing application in different fields and the relative spreading and dissemination.

Conclusions
A review of the published scientific contributions that focused on the European continent, excluding Turkey and Russia, combining GNSS (Global Navigation Satellite System) and InSAR (Interferometric Synthetic Aperture Radar) data, is reported here. A total of 191 scientific contributions were collected and analyzed in detail. The high number of contributions includes peer-reviewed articles and book chapters published in international journals, as well as the abstracts, conference proceedings, and extended abstracts of national and international congresses and symposiums.
An analysis of the spatial and temporal distribution of the scientific contributions was completed for Europe, while a specific focus was dedicated to Italy, since it is the most represented country, with 47 scientific contributions. Italy is also the first nation that reported a study combining GNSS and InSAR data. In fact, the first research work was recorded in 1999 and investigated the geometry and slip distribution of the fault plane originated Colfiorito earthquakes (central Italy), occurred on 26 September 1997.
Most of the joint applications in Europe, as also recognizable for Italy, focus on landslide, subsidence, and seismic/tectonic case studies. The spatial distribution mostly corresponds to the distribution of geohazards in Europe, e.g., subsidence phenomena investigated in Netherlands and volcano case studies in southern Italy, Iceland, and Spanish islands.
Considering the free access of InSAR data by Ground Motion Services (GMS), a good correspondence can be recognized. This research is interesting since it is a good demonstration that the scientific community and administrative managing offices are interested in the remotely sensed analysis of ground deformation. This can be a good starting point for the continued development of improvements in this branch of remote sensing research and in the free accessibility of data policy. Scientific advancement by the scientific community, and consequently the local administrator and the entire population, could receive a great benefit from open access data both from the perspective of research and in confidence in the products.