Advanced SAR Interferometric Analysis to Support Geomorphological Interpretation of Slow-Moving Coastal Landslides ( Malta , Mediterranean Sea )

An advanced SAR interferometric analysis has been combined with a methodology for the automatic classification of radar reflectors phase histories to interpret slope-failure kinematics and trend of displacements of slow-moving landslides. To accomplish this goal, the large dataset of radar images, acquired in more than 20 years by the two European Space Agency (ESA) missions ERS-1/2 and ENVISAT, was exploited. The analysis was performed over the northern sector of Island of Malta (central Mediterranean Sea), where extensive landslides occur. The study was assisted by field surveys and with the analysis of existing thematic maps and landslide inventories. The outcomes allowed definition of a model capable of describing the geomorphological evolution of slow-moving landslides, providing a key for interpreting such phenomena that, due to their slowness, are usually scarcely investigated.


Introduction
The use of interferometric methods to infer earth surface deformation from satellite-borne synthetic aperture radar has become a regular tool in geoscience since the early 1990s [1].Generally, sporadic/intermittent movements can be measured, e.g., displacements higher than 1 cm occurring during major earthquakes, deformations in volcanic areas and other kinds of episodic motion [2].The satellite orbit repetition determines a temporal bracketing as two or more acquisitions are required.Non-episodic motions can be surveyed by combining long series of radar acquisitions thus, it is possible to measure, for example, the rate of subsidence due to mining [3], to ground-water extraction [4] or to slow-moving landsliding [5][6][7].The algorithms that are nowadays used to monitor and measure long-term evolution of ground deformation descend from two types of approach: the Short BAseline Subsets (SBAS) proposed by [8], and the PSI (Persistent Scatterer Interferometry) developed by [9].
A technique derived from the latter approach, named Interferometric Point Target Analysis (IPTA) [10], was applied in the overall framework of landslide hazard assessment and management.IPTA was combined with the free-ware PSTime, implemented by [11], to achieve a better comprehension of the geomorphological evolution of extremely slow mass movements.This combined analysis was applied to the north-western portion of the Island of Malta (central Mediterranean Sea), where extensive landslide processes occur [12,13].Almost 20 years of ERS and ENVISAT acquisitions, stored in the European Space Agency (ESA) archives, were exploited.Analyzing the phase of isolated coherent radar reflectors through time, the evolution of the landslides was tracked and the rate of deformation was measured down to 1-2 mm/year [14].The interferometric phase histories were then automatically classified according to kinematic classes and trend of movements.The investigation was accompanied by field surveys and with the analysis of existing thematic maps and landslide inventories.The outcomes allowed to define a model capable to describe the evolution of coastal slopes affected by slow-moving landslides, providing a key for interpreting phenomena that are usually scarcely investigated.

Geological and Geomorphological Setting
Malta is the main island of the Maltese archipelago, which is located 90 km south of Sicily and 290 km north of Africa.The archipelago rises on the submerged Malta Plateau, through which it is connected with the Hyblean region of Sicily, and constitutes the northern shoulder of the Pantelleria Rift.The Maltese islands are characterized by a Late Oligocene-Miocene stratigraphic sequence made by three carbonatic and one clayey formations.From the oldest to the youngest they are: Lower Coralline Limestone Formation (Fm.), Globigerina Limestone Fm., Blue Clay Fm. and Upper Coralline Limestone Fm.The landscape is shaped by marine action, karst and limited fluvial processes as well as relevant and widespread coastal mass movements.
The study area includes the north-western coast of Malta and is limited southerly by Il-Qarraba peninsula and northerly by Paradise Bay, extending for about 10 km 2 (Figure 1).The structural setting of this portion of the island is characterized by an ENE-WSW oriented horst and graben system.From a morphological viewpoint, the grabens correspond to valleys and lowlands, whereas horsts make up plateaus and low-altitude promontories.Along the shoreline, plunging cliffs and rocky coasts alternate with inlets hosting pocket beaches, appreciated by tourists and locals for their pristine beauties [15].From a geological viewpoint, horsts consist of Upper Coralline Limestone Fm., which overlies lowlands or gentle slopes composed by Blue Clay Fm. [16][17][18][19].The superposition of brittle limestone plateaus on plastic-ductile clayey terrains favors slope-failure processes.The anthropic print is less significant than in the southern and eastern coasts, also thanks to the institution of the Il-Majjistral Nature and History Park that allowed protecting the central sector of the study area [20].

Type and Distribution of Coastal Landslides
Most of the mass movements in the study area are block slides often related to rock spreading and characterized by slow or extremely slow velocities.Rock spreads are witnessed by the presence of several tension cracks located along the edges of the karst plateaus (Figure 2), which can exceed 200 m in horizontal persistence such as in wide sectors of Marfa Ridge plateau (aerial photo 3981 in Figure 2), in the northern portion of Anchor Bay, at Ras Il-Wahx promontory and Il-Qarraba peninsula (aerial photo 4076 in Figure 2).These tension cracks are generated by the different mechanical properties of limestone and the underlying clayey terrains, and favor rainfall infiltration, which feeds the underground karst aquifer.The combination of structural discontinuities and tension cracks, generated and enlarged by the action of rock spreads, isolates limestone pillars or large blocks, which are detached, slowly lowered or toppled from the plateau cliffs.These blocks are abundant and form wide accumulations (named as rdum by the locals), which slowly move towards the sea [12,13,21].

Type and Distribution of Coastal Landslides
Most of the mass movements in the study area are block slides often related to rock spreading and characterized by slow or extremely slow velocities.Rock spreads are witnessed by the presence of several tension cracks located along the edges of the karst plateaus (Figure 2), which can exceed 200 m in horizontal persistence such as in wide sectors of Marfa Ridge plateau (aerial photo 3981 in Figure 2), in the northern portion of Anchor Bay, at Ras Il-Wahx promontory and Il-Qarraba peninsula (aerial photo 4076 in Figure 2).These tension cracks are generated by the different mechanical properties of limestone and the underlying clayey terrains, and favor rainfall infiltration, which feeds the underground karst aquifer.The combination of structural discontinuities and tension cracks, generated and enlarged by the action of rock spreads, isolates limestone pillars or large blocks, which are detached, slowly lowered or toppled from the plateau cliffs.These blocks are abundant and form wide accumulations (named as rdum by the locals), which slowly move towards the sea [12,13,21].

Materials and Methods
Interferometric Point Target Analysis (IPTA) investigates the temporal and spatial characteristics of the phase backscattered from radar reflectors on the ground to accurately map surface deformation histories, terrain heights and relative atmospheric path delays.Usually, radar reflectors can be found in built environments, since walls and roofs can consistently backscatter the radar signal through time.The displaced limestone blocks over the clayey slopes of the north-western coast of Malta can also be considered suitable reflectors, as it can be noticed by the large amounts of point target candidates detected in the study area.
The processed data include 50 ERS images acquired between 1992 and 2000 and 33 ENVISAT acquisitions covering the period 2004-2009.Details of the processed datasets have been resumed in Tables 1 and 2.

Materials and Methods
Interferometric Point Target Analysis (IPTA) investigates the temporal and spatial characteristics of the phase backscattered from radar reflectors on the ground to accurately map surface deformation histories, terrain heights and relative atmospheric path delays.Usually, radar reflectors can be found in built environments, since walls and roofs can consistently backscatter the radar signal through time.The displaced limestone blocks over the clayey slopes of the north-western coast of Malta can also be considered suitable reflectors, as it can be noticed by the large amounts of point target candidates detected in the study area.
The processed data include 50 ERS images acquired between 1992 and 2000 and 33 ENVISAT acquisitions covering the period 2004-2009.Details of the processed datasets have been resumed in Tables 1 and 2.
After the co-registration to a single reference image, denoted as Master in the above mentioned tables, the corresponding interferograms were computed.The reference images were chosen among those: (i) acquired near the temporal average of the stack of the processed images, to minimize the temporal decorrelation of interferograms; (ii) provided with an orbit as close as possible to the center of the orbital tube spanned by the stack of the processed images to minimize the spatial decorrelation of interferograms (Figure 3).After the co-registration to a single reference image, denoted as Master in the above mentioned tables, the corresponding interferograms were computed.The reference images were chosen among those: (i) acquired near the temporal average of the stack of the processed images, to minimize the temporal decorrelation of interferograms; (ii) provided with an orbit as close as possible to the center of the orbital tube spanned by the stack of the processed images to minimize the spatial decorrelation of interferograms (Figure 3).These candidate scenes were then processed with selected images of the stack to generate interferograms characterized by small perpendicular baselines and short time intervals.With a priori knowledge of deformation signal locations and a good coherence, a very flat differential interferogram was interpreted as a clear indication of low atmospheric distortions for both acquisitions and hence in the candidate scene.The interferograms were then subtracted to a synthetic phase, derived from the 3-Arc-Seconds Shuttle Radar Topography Mission (SRTM) Digital Elevation Model (DEM), obtaining 49 ERS and 32 ENVISAT differential interferograms.A Persistent Scatterer (PS) candidates list was compiled based on two detection algorithms implemented in IPTA processing chain.The first one analyzes the spectral properties of individual Single Look Complex (SLC) images and identifies point targets from the spectral behavior and a dominating backscattering.The second exploits large SLC data stacks to detect pixels showing low intensity variability through time.The spatial differences (between pairs of point target candidates) of the differential interferometric phases were analyzed by means of a bi-dimensional regression that takes into account the perpendicular baseline and the time elapsed between two acquisitions.Differences between the measured and modeled phases are related to uncompensated topography, orbital errors, atmospheric path delays, non-linear deformations and noise.Standard deviation of the modeled phase represents a quality index that is used to reject point target candidates that do not meet IPTA requirements.The computation consisted of an iteration process.At the end of every regression, the residuals were inspected, interpreted and filtered out accordingly.In the final phase model, the term related to the displacement was isolated and decomposed in turn as follows: Displacement phase " plinear deformation phaseq `pnon-linear deformation phaseq `pnoiseq The freeware PSTime [11] has been recently used for the automatic classification of PS time series.The algorithm, based on a conditional sequence of statistical tests, automatically groups the PS into classes of predefined trends (such as uncorrelated, linear, quadratic, bilinear and discontinuous trend).PSTime was applied to classify the deformation phase histories obtained by the IPTA analysis, in order to investigate the geomorphological evolution of the coastal landslides during the last two decades.Nevertheless, field surveys, previous studies and the analysis of the deformation recorded by a GNSS (Global Navigation Satellite System) monitoring network deployed since 2005 at Anchor Bay and Ghajn Tuffieha Bay [12,13,21] indicated that the kinematic aspects of these landslide do not fit properly the PSTime trend classes.The parameters of the statistical test of PSTime were hence redefined to force the time series in a new customized version of target trends.The deformation recorded at the GNSS benchmarks locations were projected along the sensors LOS (Line of Sight) according to the following formula: where: d LOS is the deformation along the radar LOS, dE is the deformation measured by the GNSS benchmark in the east coordinate, dN is the deformation measured by the GNSS benchmark in the north coordinate, dH is the deformation measured by the GNSS benchmark in the elevation, φ is the look angle of the radar, θ is the azimuth direction.Deformations along the LOS of the GNSS benchmarks were used as a model to define six new classes (Figure 4a1-f1) that identified rock blocks as: a.

Stable b.
Deforming linearly downhill c.
Accelerating downhill e.
With discontinuous rate and trend The good agreement between the new modeled trend classes and the time series actually measured at some representative PS locations can be seen in Figure 4a2-f2.The radar reflectors classified as uncorrelated by the PSTime algorithm and those deforming linearly with a rate included in the interval (1, −1) mm/year were considered stable points in the model.

Data Analysis and Results Interpretation
Over the northern part of Malta, 2763 ERS and 6386 ENVISAT persistent scatterers (Figure 5) were identified.
Datasets were processed using the same reference point in order to avoid biased results.Most of the point targets are stable, while deformations are recorded along several sectors of the coast, where slope-failure processes are widespread.PSI processing allowed recognizing active landslides and measure their movement rates, which resulted to be extremely slow, ranging from 1 to 7 mm/year.Furthermore, it is observable that the rate of activity has changed through years.Some sectors show higher rate of deformations in the ERS (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001) acquisitions than in the ENVISAT (2004ENVISAT ( -2009) ) and vice-versa.ENVISAT dataset (Figure 5) appears more noisy then ERS because of the lower number of interferograms processed.Statistical parameters for the assessment of measurement quality are reported in Table 3.The good agreement between the new modeled trend classes and the time series actually measured at some representative PS locations can be seen in Figure 4a2-f2.The radar reflectors classified as uncorrelated by the PSTime algorithm and those deforming linearly with a rate included in the interval (1, ´1) mm/year were considered stable points in the model.

Data Analysis and Results Interpretation
Over the northern part of Malta, 2763 ERS and 6386 ENVISAT persistent scatterers (Figure 5) were identified.
Datasets were processed using the same reference point in order to avoid biased results.Most of the point targets are stable, while deformations are recorded along several sectors of the coast, where slope-failure processes are widespread.PSI processing allowed recognizing active landslides and measure their movement rates, which resulted to be extremely slow, ranging from 1 to 7 mm/year.Furthermore, it is observable that the rate of activity has changed through years.Some sectors show higher rate of deformations in the ERS (1992)(1993)(1994)(1995)(1996)(1997)(1998)(1999)(2000)(2001) acquisitions than in the ENVISAT (2004-2009) and vice-versa.ENVISAT dataset (Figure 5) appears more noisy then ERS because of the lower number of interferograms processed.Statistical parameters for the assessment of measurement quality are reported in Table 3.In order to evaluate the reliability of the analysis, the deformations recorded by homologous PS (the targets that can be identified in the two datasets) were plotted versus time.As it can be noticed from Figure 6, the results appear to be consistent with the expected kinematic behavior of the investigated processes.
Four landslide areas were chosen as test sites for interferometric analysis, based on the presence of a statistically relevant number of targets detected in the two datasets.The chosen areas are the following: (i) western sector of Marfa Ridge; (ii) stretch between Rdum id-Delli and Rdum Majjiesa; (iii) Ras Il-Wahx headland; (iv) peninsula of Il-Qarraba.

Marfa Ridge Landslides
Marfa Ridge consists of a karst plateau which displays the maximum altitude of the northern coast of Malta (125 m a.s.l.).It is crossed by an E-W oriented system of faults and Globigerina Limestone outcrops at sea level.Towards the sea, the peninsula is characterized by block slides affecting clayey slopes capped by the Upper Coralline Limestone plateau (Figure 7).Rock falls and rock topples occur at the edges of the plateau and are favored by rock spreading cf.[22].The histogram in Figure 8 highlights that the majority of ERS PS shows a downhill deformation, while several ENVISAT PS are affected by uplift.Moreover, it can be noticed that the site shows different deformation trends.Thereby, areas characterized by linear subsidence in ERS images and appearing stable in ENVISAT ones are alternated with areas that, in contrast, seem to be stable in ERS dataset and show uplift in ENVISAT acquisitions.A difference in the rate of deformation, can also be noticed from Figure 8, between the northern and the southern sectors, being the former slightly faster.In summary the two dataset allowed us to frame two different stages of an extremely slow deformation process.

Marfa Ridge Landslides
Marfa Ridge consists of a karst plateau which displays the maximum altitude of the northern coast of Malta (125 m a.s.l.).It is crossed by an E-W oriented system of faults and Globigerina Limestone outcrops at sea level.Towards the sea, the peninsula is characterized by block slides affecting clayey slopes capped by the Upper Coralline Limestone plateau (Figure 7).

Marfa Ridge Landslides
Marfa Ridge consists of a karst plateau which displays the maximum altitude of the northern coast of Malta (125 m a.s.l.).It is crossed by an E-W oriented system of faults and Globigerina Limestone outcrops at sea level.Towards the sea, the peninsula is characterized by block slides affecting clayey slopes capped by the Upper Coralline Limestone plateau (Figure 7).Rock falls and rock topples occur at the edges of the plateau and are favored by rock spreading cf.[22].The histogram in Figure 8 highlights that the majority of ERS PS shows a downhill deformation, while several ENVISAT PS are affected by uplift.Moreover, it can be noticed that the site shows different deformation trends.Thereby, areas characterized by linear subsidence in ERS images and appearing stable in ENVISAT ones are alternated with areas that, in contrast, seem to be stable in ERS dataset and show uplift in ENVISAT acquisitions.A difference in the rate of deformation, can also be noticed from Figure 8, between the northern and the southern sectors, being the former slightly faster.In summary the two dataset allowed us to frame two different stages of an extremely slow deformation process.Rock falls and rock topples occur at the edges of the plateau and are favored by rock spreading cf.[22].The histogram in Figure 8 highlights that the majority of ERS PS shows a downhill deformation, while several ENVISAT PS are affected by uplift.Moreover, it can be noticed that the site shows different deformation trends.Thereby, areas characterized by linear subsidence in ERS images and appearing stable in ENVISAT ones are alternated with areas that, in contrast, seem to be stable in ERS dataset and show uplift in ENVISAT acquisitions.A difference in the rate of deformation, can also be noticed from Figure 8, between the northern and the southern sectors, being the former slightly faster.In summary the two dataset allowed us to frame two different stages of an extremely slow deformation process.

Rudm Id-Delli-Rdum Majjiesa Landslides
This site shows different stages of evolution of rock spreading and block sliding: the northernmost part (Rdum id-Delli) is characterized by a limestone plateau from which large portions were isolated and lowered.The latter exhibit a developed network of persistent fractures.At the sea level smaller blocks are abundant.In the southernmost sector (Rdum Majjiesa), rock fall deposits are frequent at the bottom of the plateau edge and are remobilized by block slides (Figure 9).PS at Rdum id-Delli block slide show a linear downhill trend through time, except three targets that record a discontinuous rate and trend of deformation, consisting of subsidence followed by

Rudm Id-Delli-Rdum Majjiesa Landslides
This site shows different stages of evolution of rock spreading and block sliding: the northernmost part (Rdum id-Delli) is characterized by a limestone plateau from which large portions were isolated and lowered.The latter exhibit a developed network of persistent fractures.At the sea level smaller blocks are abundant.In the southernmost sector (Rdum Majjiesa), rock fall deposits are frequent at the bottom of the plateau edge and are remobilized by block slides (Figure 9).

Rudm Id-Delli-Rdum Majjiesa Landslides
This site shows different stages of evolution of rock spreading and block sliding: the northernmost part (Rdum id-Delli) is characterized by a limestone plateau from which large portions were isolated and lowered.The latter exhibit a developed network of persistent fractures.At the sea level smaller blocks are abundant.In the southernmost sector (Rdum Majjiesa), rock fall deposits are frequent at the bottom of the plateau edge and are remobilized by block slides (Figure 9).PS at Rdum id-Delli block slide show a linear downhill trend through time, except three targets that record a discontinuous rate and trend of deformation, consisting of subsidence followed by PS at Rdum id-Delli block slide show a linear downhill trend through time, except three targets that record a discontinuous rate and trend of deformation, consisting of subsidence followed by uplift.Slight differences in the rate between ERS and ENVISAT dataset can be noticed in the histograms (Figure 10).This stretch of coast shows several accelerating downhill trends.
Remote Sens. 2016, 8, 443 12 of 18 uplift.Slight differences in the rate between ERS and ENVISAT dataset can be noticed in the histograms (Figure 10).This stretch of coast shows several accelerating downhill trends.

Ras Il-Wahx Landslides
The western sector of the 100 m-high Ras Il-Wahx promontory shows spectacular landslides (Figure 11).Rock spreading affects a deeply jointed area extending inland up to 70 m from the cliff edge.The headland is surrounded by block slide accumulations, which consist of tens of piled-up blocks sometimes of relevant volume.These are partially sunk in the clayey terrains and rotated backward.

Ras Il-Wahx Landslides
The western sector of the 100 m-high Ras Il-Wahx promontory shows spectacular landslides (Figure 11).Rock spreading affects a deeply jointed area extending inland up to 70 m from the cliff edge.The headland is surrounded by block slide accumulations, which consist of tens of piled-up blocks sometimes of relevant volume.These are partially sunk in the clayey terrains and rotated backward.
According to interferometric analysis this site shows the highest deformation rates of the entire north-western coast.Moreover, this is the only sector of the coast where the rate of deformation has increased through time (see histogram in Figure 12), despite a small cluster of blocks in the northern part of the promontory that reached stability (ENVISAT dataset) after an evolution characterized by discontinuous deformation rate and trend (ERS dataset).

Ras Il-Wahx Landslides
The western sector of the 100 m-high Ras Il-Wahx promontory shows spectacular landslides (Figure 11).Rock spreading affects a deeply jointed area extending inland up to 70 m from the cliff edge.The headland is surrounded by block slide accumulations, which consist of tens of piled-up blocks sometimes of relevant volume.These are partially sunk in the clayey terrains and rotated backward.According to interferometric analysis this site shows the highest deformation rates of the entire north-western coast.Moreover, this is the only sector of the coast where the rate of deformation has increased through time (see histogram in Figure 12), despite a small cluster of blocks in the northern part of the promontory that reached stability (ENVISAT dataset) after an evolution characterized by discontinuous deformation rate and trend (ERS dataset).According to interferometric analysis this site shows the highest deformation rates of the entire north-western coast.Moreover, this is the only sector of the coast where the rate of deformation has increased through time (see histogram in Figure 12), despite a small cluster of blocks in the northern part of the promontory that reached stability (ENVISAT dataset) after an evolution characterized by discontinuous deformation rate and trend (ERS dataset).

Il-Qarraba Landslides
A peculiar landform of Malta is the head-shaped peninsula of Il-Qarraba, which is connected to inland through a narrow isthmus of clayey and marly terrains [17].The latter are overlain in the central part of the peninsula by a flat limestone slab, characterized by 10 m-high cliffs.The western and northern portions of the limestone slab are affected by rock spreads evolving into block slides characterized by the presence of huge blocks.The limestone blocks move towards the sea and are dominant in the western and northern sector of the peninsula (Figure 14).Both the ERS and ENVISAT datasets show the same trend of deformation (Figure 15): the blocks are either stable or linearly subsiding and no uplifting is recorded, except for few PS.

Il-Qarraba Landslides
A peculiar landform of Malta is the head-shaped peninsula of Il-Qarraba, which is connected to inland through a narrow isthmus of clayey and marly terrains [17].The latter are overlain in the central part of the peninsula by a flat limestone slab, characterized by 10 m-high cliffs.The western and northern portions of the limestone slab are affected by rock spreads evolving into block slides characterized by the presence of huge blocks.The limestone blocks move towards the sea and are dominant in the western and northern sector of the peninsula (Figure 14).Both the ERS and ENVISAT datasets show the same trend of deformation (Figure 15): the blocks are either stable or linearly subsiding and no uplifting is recorded, except for few PS.

Il-Qarraba Landslides
A peculiar landform of Malta is the head-shaped peninsula of Il-Qarraba, which is connected to inland through a narrow isthmus of clayey and marly terrains [17].The latter are overlain in the central part of the peninsula by a flat limestone slab, characterized by 10 m-high cliffs.The western and northern portions of the limestone slab are affected by rock spreads evolving into block slides characterized by the presence of huge blocks.The limestone blocks move towards the sea and are dominant in the western and northern sector of the peninsula (Figure 14).Both the ERS and ENVISAT datasets show the same trend of deformation (Figure 15): the blocks are either stable or linearly subsiding and no uplifting is recorded, except for few PS.

Discussion
The analysis of the results cannot be separated from a brief discussion on the constraints of the technique employed to measure deformations.Considering that there are no substantial differences in terms of wavelength and acquisition geometry between the synthetic radar mounted on board of ERS and ENVISAT satellites and that the images are acquired from the same track, it can be assumed that the sensitivity to deformations is the same in the two datasets.The fact that block slide movement direction is mainly E-W oriented makes interferometric measurement very effective [23].Taking into account that Malta was framed from a descending orbit, the expected downhill displacement of the limestone blocks resulted in an increment of the LOS distance.In contrast the shortening of the LOS distance, that in the previous chapters was simply defined as uplift, can be interpreted as a back-tilting of the limestone blocks towards the sensor favored by the swelling of the clayey deposits.The evolution of the landslides which characterize the north-western coast of Malta is in agreement with the model proposed by [21,22,24], according to which the deformation process is composed by the following three phases: (i) rock spreading affecting the intensely jointed limestone cap; (ii) progressive detachment of limestone blocks of large size from the plateau; (iii) tilting and sliding of the blocks above the underlying clays which deform visco-plastically.Within this framework, rock falls and earth flows/slides represent shallow evidence of a deeper movement involving the entire slope from the plateau to the coastline.Lateral spreading also favors the development of bulges in the clays at the foot of the cliffs because of the overburden [25].Moreover, the results here obtained properly fit with the model proposed by [21], that applied a multidisciplinary approach by means of geomorphological and geophysical surveys and in situ monitoring activities at Anchor Bay site.
IPTA analysis coupled with PSTime algorithm allowed to detect and discriminate six different phase variation patterns that can be associated to the above mentioned evolutionary model (Figure 13).
These phase trends can be interpreted from a geomorphological viewpoint as follows: 1.
Stable area.No significant deformation trend can be detected by the interferometric analysis.
The majority of stable PS can be found in correspondence of wide block slide deposits or over the less jointed portions of plateaus.It is likely that the stability of the latter can be related to the limited thickness of the limestone slabs or to the presence of extensive network of joints that do not favor the presence of a karst aquifer.These types of settings occur at Il-Qarraba for which the interferometric results confirm the above mentioned hypothesis.

2.
Downhill-moving area.A significant increment in LOS distance between the target and the sensor can be detected.The associated deformation trend can be interpolated linearly and can be constant through time or discontinuous (i.e., accelerating or decelerating).The PS showing this behavior were interpreted as blocks detached from the plateau that began to move downward towards the sea.The deformation process may be ascribable to the presence of a slip surface.This kind of movement is more common in the upper part of steep slopes as in the area of Rdum id-Delli and Rdum Majjiesa.Again the interferometric results appear consistent with this geomorphological interpretation.

3.
Discontinuous-moving area.The seaward translation of the blocks on the surface is irregular.When the steepness of the slope decreases or the slip surface progressively becomes sub-horizontal, the blocks might come to a temporary rest.The consequent uplift measured in terms of reduction of the LOS is likely to be associated to a backward tilting of the blocks towards the plateau favored by the bulging and expansion of the underlying clay.This style of deformation is evident at Marfa Ridge where blocks that slide towards the sea, according to the ERS dataset, subsequently showed a back-tilting in the interferometric analysis performed using ENVISAT images.On the other hand, according to PS measurements on ERS dataset, at Ras Il-Whax blocks first show a discontinuous deformation trend and then appear to be stable or to re-start their descent along the slope.

4.
Uplift-moving area.It is the natural evolution of the discontinuous trend, as a matter of fact it is frequently observed around cluster of blocks characterized by this type of deformation as can be clearly noted in Figures 7 and 11.This behavior is also favored by the bulging and the expansion of the underlying clay.
It is likely that the main factors controlling the style and rate of the deformations are the weight of the limestone plateaus; the relative ratios of limestone plateau and clayey terrains thickness; and the undermining at the toe of the cliff.Where the limestone caps have considerable thickness and the clayey terrains outcrops at sea level, the displacement rates are higher: perfect examples are Ras Il-Wahx site and the coast comprised between Rdum id-Delli and Rdum Majjiesa where 20 m-thick of limestone cap can be found, which is almost the double of the thickness recorded at Il-Qarraba.Another proof of this hypothesis is the case of Marfa Ridge: on the northernmost sector of the ridge, the blocks show higher rates of movement compared to the south.Again, a considerable thickness of the limestone plateau (17 m) can be found along with the outcropping of Blue Clays at sea level.This setting is not preserved in the southernmost part, where the Globigerina limestone appears along the coastline.The process of undermining is the direct consequence of marine water action and the presence of an aquifer inside the limestone slab that favors the occurrence of springs at the contact between limestone and clay.This process can be found along the stretch of coast between Rdum id-Delli and Ras Il-Wahx and where relatively high rate of displacements have been recorded from the interferometric analysis.

Conclusions
The analysis of long-term evolution of widespread slow-moving landslides is an expensive and time-consuming task to achieve using traditional monitoring techniques.Advanced SAR interferometric analyzes represent nowadays the only suitable methodology, even by virtue of the fact of being able to exploit nearly 20 years of data.Interferometric methods have been often used to measure surface displacements caused by mass movements, but in this study IPTA processing and a modified version of PSTime algorithm were coupled in order to provide more relevant information.These outcomes assisted in the geomorphological interpretation of slow-moving landslides, which are seldom investigated with such detailed accuracy.From the interferometric phase, information allowing to define a model of the Maltese coastal landslides, congruent to field observations and GNSS measurements, were extracted.The mechanism of evolution was inferred and reliable measurements of the rate and trend of deformation were provided.The importance of the results achieved is evident assuming that all the information collected from the interferometric analysis is crucial for an effective landslide hazard assessment [26].The combined methodology used in this study shows the potential for being extended to other areas affected by slow-moving landslides and characterized by similar geological settings.Further developments of this research will include the collection and analysis of the new ESA mission Sentinel-1 up-to-date datasets.

Figure 1 .
Figure 1.Geological and geomorphological setting of the study area.

Figure 1 .
Figure 1.Geological and geomorphological setting of the study area.

Figure 2 .
Figure 2. Distributions of the main cracks affecting the plateaus.(A) Detail of the northernmost part of the study area; (B) detail of the southernmost part of the study area.

Figure 2 .
Figure 2. Distributions of the main cracks affecting the plateaus.(A) Detail of the northernmost part of the study area; (B) detail of the southernmost part of the study area.

Figure 3 .
Figure 3. Spatial and temporal baselines of the interferograms.The red spot define the acquisition chosen as reference; red crosses denote noisy interferograms not used in the final processing.

Figure 3 .
Figure 3. Spatial and temporal baselines of the interferograms.The red spot define the acquisition chosen as reference; red crosses denote noisy interferograms not used in the final processing.

Figure 5 .
Figure 5. Results of the IPTA analysis performed on the ERS dataset (1999-2000) and ENVISAT dataset (2004-2009).Below the GNSS benchmark location at the site of Anchor Bay and Il-Qarraba.

Figure 7 .
Figure 7. Marfa Ridge site seen from the north.

Figure 7 .
Figure 7. Marfa Ridge site seen from the north.

Figure 7 .
Figure 7. Marfa Ridge site seen from the north.

Figure 8 .
Figure 8. Marfa Ridge landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6a.

Figure 8 .
Figure 8. Marfa Ridge landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6a.

Figure 8 .
Figure 8. Marfa Ridge landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6a.

Figure 10 .
Figure 10.Rdum id-Delli-Rdum Majjiesa landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.Blue bars are ERS PS, red are ENVISAT.The white box on the maps locate the PS whose time series are reported in Figure 6b.

Figure 10 .
Figure 10.Rdum id-Delli-Rdum Majjiesa landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.Blue bars are ERS PS, red are ENVISAT.The white box on the maps locate the PS whose time series are reported in Figure 6b.

Figure 11 .
Figure 11.Ras Il-Wahx site seen from the south.

Figure 11 .
Figure 11.Ras Il-Wahx site seen from the south.

Figure 12 .
Figure 12.Ras Il-Wahx landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6c.The white line defines the cross section of Figure 13.

Figure 12 .
Figure 12.Ras Il-Wahx landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6c.The white line defines the cross section of Figure 13.

Figure 12 .
Figure 12.Ras Il-Wahx landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6c.The white line defines the cross section of Figure 13.

Figure 15 .
Figure 15.Il-Qarraba landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6d.

Figure 15 .
Figure 15.Il-Qarraba landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6d.

Figure 15 .
Figure 15.Il-Qarraba landslide map and distribution of the PS, color coded according to the trend classes defined in Figure 4: (A) ERS satellite PS deformation trend; (B) ENVISAT satellite PS deformation trend; (C) Histogram showing the relative frequency of occurrence in each class of movement rate.The white box on the map locates the PS whose time series are reported in Figure 6d.

Table 1 .
ERS dataset.Track 494, frame 2889 descending.Images denoted in the NOTES column as noisy are those that have been removed from the final iteration due to the strong presence of atmospheric disturbances.Image denoted with MASTER is the reference images from which all the interferograms have been generated.

Table 1 .
ERS dataset.Track 494, frame 2889 descending.Images denoted in the NOTES column as noisy are those that have been removed from the final iteration due to the strong presence of atmospheric disturbances.Image denoted with MASTER is the reference images from which all the interferograms have been generated.

Table 2 .
ENVISAT dataset.Track 494, frame 2889 descending, swath I2.Images denoted in the NOTES column as noisy are those which have been removed from the final iteration due to the strong presence of atmospheric disturbances.Image denoted with MASTER is the reference images from which all the interferograms have been generated.

Table 3 .
Statistical parameters of the interferometric analysis.