Monitoring and Mapping of Shallow Landslides in a Tropical Environment Using Persistent Scatterer Interferometry: A Case Study from the Western Ghats, India

Persistent Scatterer Interferometry (PSI) techniques are now well established and accepted for monitoring ground displacements. The presence of shallow-seated landslides, ubiquitous phenomena in the tropics, offers an opportunity to monitor and map these hazards using PSI at the regional scale. Thus, the Western Ghats of India, experiencing a tropical climate and in a topographically complex region of the world, provides an ideal study site to test the efficacy of landslide detection with PSI. The biggest challenge in using the PSI technique in tropical regions is the additional noise in data due to vegetation. In this study, we filtered these noises by utilizing the 95-percentile of the highest coherence data, which also reduced the redundancy of the PSI points. The study examined 12 landslides that occurred within one of the three temporal categories grouped as Group 1, Group 2, and Group 3, categorized in relation to PSI monitoring periods, which was also further classified into east- and west-facing landslides. The Synthetic Aperture Radar (SAR) data is in descending mode, and, therefore, the east-facing landslides are characterized by positive deformation velocity values, whereas the west-facing landslides have negative deformation values. Further, the landslide-prone areas, delineated using the conventional factor of safety (FS), were refined and mapped using PSI velocity values. The combination of PSI with the conventional FS approach helped to identify exclusive zones prone to landslides. The main aim of such an attempt is to identify critical areas in the unstable category in the map prepared using FS and prioritizing the mitigation measures, and to develop a road map for any developmental activities. The approach also helps to increase confidence in the susceptibility mapping and reduce false alarms.


Introduction
Landslides are natural hazards that pose severe social and economic consequences because they can cause widespread damage to life and property [1]. Landslide triggering mechanisms fall into three critical categories [2]: geological (e.g., seismic activity from other hazards like earthquakes or volcanoes), geomorphological (e.g., changes in slope features through surficial processes), and hydrological (e.g., precipitation, groundwater, or freeze/thaw) [3]. The latter, hydrological triggering mechanisms, are the most prevalent cause of landslides around the world, especially in tropical regions and areas that receive heavy precipitation [4][5][6][7]. Oversaturated, non-cohesive soils and variegated debris may

Study Area
The Western Ghats, a mountain range, and prominent escarpment paralleling the southwestern coast of India ( Figure 1a) is a hotspot for hydrologically-triggered landslides [29]. About 8% (roughly 1400 km 2 ) of the Western Ghats of Kerala is classified as a critical zone for mass movements [11,30].
The present study is concentrated in the Ernakulam, Idukki, and Kottayam districts of central Kerala (shown in Figure 1b) and lies between 9 • 40 N to 10 • 20 N latitude and 76 • 40 E to 77 • 10 E longitude covering an area of 2221 km 2 . The average annual rainfall is between 250 cm and 425 cm; the temperature varies between 21 • C to 29 • C. Landslides here are shallow-seated and relatively small, typically less than 2 m in depth [11] (Figure 1c). Rugged hills and deep valleys characterize this section of the Western Ghats. Holocene sediments and soils rest unconformably atop Precambrian bedrock composed of gneiss (e.g., biotite, hornblende-biotite, and more felsic granite gneiss), charnockite, granulite, and quartzite [31]. The stratigraphic column of the state is shown in Table 1. The contact between sediments and bedrock usually forms the landslide rupture surface. Table 1. Stratigraphic column for the state of Kerala (Source: Nandakumar and Balakrishnan [32]).

Age
Formation Lithology  Granitic gneiss, Schists etc. Charnockites and associated rocks Khondalite suite of rocks and its associates

Data and Methodology
Twenty-two descending single look complex (SLC) synthetic aperture radar (SAR) images acquired between March 2016 and June 2017 from the Sentinel-1A satellite ( Table  2), courtesy of the European Space Agency (ESA), were downloaded from the Alaska Satellite Facility at the University of Alaska Fairbanks [33]. Sentinel-1A acquires SAR images at an approximate wavelength of 5.6 cm (C-band) over an area with dimensions 250 by 180 km at a revisiting cycle of 12 days. The incidence angle for the Interferometric Wide Swath (IWS) mode ranges from 29.1 degrees to 46 degrees, with a central incidence angle

Data and Methodology
Twenty-two descending single look complex (SLC) synthetic aperture radar (SAR) images acquired between March 2016 and June 2017 from the Sentinel-1A satellite ( Table 2), courtesy of the European Space Agency (ESA), were downloaded from the Alaska Satellite Facility at the University of Alaska Fairbanks [33]. Sentinel-1A acquires SAR images at an approximate wavelength of 5.6 cm (C-band) over an area with dimensions 250 by 180 km at a revisiting cycle of 12 days. The incidence angle for the Interferometric Wide Swath (IWS) mode ranges from 29.1 degrees to 46 degrees, with a central incidence angle of approximately 37.55 degrees. Sentinel-1A images were processed using the PSI technique to measure ground deformation rates with an accuracy up to 1 mm/yr [34] following the default approach provided with the ENVI SARscape software [19,20,35]. The methodology behind landslide mapping and monitoring is shown in Figure 2. of approximately 37.55 degrees. Sentinel-1A images were processed using the PS nique to measure ground deformation rates with an accuracy up to 1 mm/yr [34] foll the default approach provided with the ENVI SARscape software [19,20,35]. The m ology behind landslide mapping and monitoring is shown in Figure 2.

Selection of Landslide Location
Twelve different landslide locations were collated from various media so mostly newspaper articles, and verified through field studies. The details of these slides are shown in Table 3. All these landslides were confined to roadcuts. These slides were categorized into three groups based on activity periods that align wi data acquisition: Landslides that occurred before March 2016, i.e., representing lan that occurred before the SAR data collection, are referred to as 'Group 1'; between 2016 and June 2017, i.e., landslides that occurred during the SAR data collection period, are referred to as 'Group 2', and after June 2017, i.e., landslides after the SA collection time-period, are referred to as 'Group 3' (Figure 3). Individual landslide studied against the PSI deformation values. The dates of data were selected prior 2018 anomalous rainfall that occurred in Kerala which resulted in one of the biggest and 4726 landslides, so that more significant movements associated with this anom rainfall should not manifest in our results, thereby, any anomaly can be avoided.

Selection of Landslide Location
Twelve different landslide locations were collated from various media sources, mostly newspaper articles, and verified through field studies. The details of these landslides are shown in Table 3. All these landslides were confined to roadcuts. These landslides were categorized into three groups based on activity periods that align with PSI data acquisition: Landslides that occurred before March 2016, i.e., representing landslide that occurred before the SAR data collection, are referred to as 'Group 1'; between March 2016 and June 2017, i.e., landslides that occurred during the SAR data collection time-period, are referred to as 'Group 2', and after June 2017, i.e., landslides after the SAR data collection time-period, are referred to as 'Group 3' (Figure 3). Individual landslides were studied against the PSI deformation values. The dates of data were selected prior to the 2018 anomalous rainfall that occurred in Kerala which resulted in one of the biggest floods and 4726 landslides, so that more significant movements associated with this anomalous rainfall should not manifest in our results, thereby, any anomaly can be avoided.

Monitoring of Individual Landslides
These three categories of landslides were correlated with the temporal deformation of PSI points. These three groups were further subdivided into two based on the eastfacing and west-facing slope. This classification is performed to account for the Line of Sight (LoS) of Sentinel-1A satellite that plays a vital role in the deformation direction measured using PSI. These Sentinel-1A images were acquired in descending mode. Thus, landslides can show up in PSI data in two ways: 1. Negative velocity values on west-facing slopes, indicating downward or westward (downslope) movement, and 2. Positive velocity values on east-facing slopes, indicating downward or eastward (downslope) movement.
This method will show good results for landslides with E-W oriented slopes (±45°) [36]. The east and west-facing slopes were determined using the 30 m spatial resolution Shuttle Radar Topographic Mission (SRTM) digital elevation model (DEM) downloaded from the United States Geological Survey (USGS) website [37].

Filtering of PSI Points
The PSI technique can produce a large number of points depending on the choice of a coherence threshold. Although all the points convey the same amount of information,

Monitoring of Individual Landslides
These three categories of landslides were correlated with the temporal deformation of PSI points. These three groups were further subdivided into two based on the east-facing and west-facing slope. This classification is performed to account for the Line of Sight (LoS) of Sentinel-1A satellite that plays a vital role in the deformation direction measured using PSI. These Sentinel-1A images were acquired in descending mode. Thus, landslides can show up in PSI data in two ways: 1.
Negative velocity values on west-facing slopes, indicating downward or westward (downslope) movement, and 2.
Positive velocity values on east-facing slopes, indicating downward or eastward (downslope) movement.
This method will show good results for landslides with E-W oriented slopes (±45 • ) [36]. The east and west-facing slopes were determined using the 30 m spatial resolution Shuttle Radar Topographic Mission (SRTM) digital elevation model (DEM) downloaded from the United States Geological Survey (USGS) website [37].

Filtering of PSI Points
The PSI technique can produce a large number of points depending on the choice of a coherence threshold. Although all the points convey the same amount of information, they are affected differently by random noise and the strength of the signal being measured. To identify and isolate only the points that convey the maximum amount of information and avoid potential noise and artifacts (especially the influence of vegetation in the tropical region), we undertook a selection process looking for points with a high signal to noise ratio and signal quality utilizing the coherence value. For the filtering, we select only Persistent Scatterer (PS) points that have higher coherence than the 95-percentile.

Creation of Slope Stability Map
Geologic, hydrologic, and geomorphologic data, including the normalized difference vegetation index (NDVI), slope angle, groundwater level, and soil strength parameters, were used to create a landslide susceptibility map. NDVI was generated using the Near-Infrared (NIR) and Red bands of Landsat 8 [36]. The formula for NDVI generation is: and was reclassified based on Holben [38]. The soil map was taken from the National Bureau of Soil Survey data [39], which used the Unified Soil Classification System (USCS). The slope map was prepared using SRTM data [40]. The slope map was classified based on the Bureau of Indian Standard [41]. All thematic maps used for the creation of landslide susceptibility map are shown in Figure 4. The spatial extent of the landslide susceptibility map was then reduced by a factor of safety (FS) analysis to identify and target landslideprone slopes. An FS map was prepared based on input factors used in the Geographic Information System Tool for Infinite Slope Stability Analysis (GIS-TISSA) [42] and Probabilistic Infinite Slope Analysis model (PISA-m) [43]. The equation for calculating the FS is: where c r is root cohesive strength (pressure); c s is soil cohesive strength (pressure); q t is tree surcharge (pressure); γ m is saturated unit weight of the soil (force/volume); γ sat is moist unit weight of the soil (force/volume); D is soil thickness; H w is pore pressure coefficient (0 ≤ H w ≤ 1); β is the angle of topographical slope; ϕ is the angle of internal friction; In this study, we assumed a partially saturated soil condition (H w = 0.5) for the slope stability analysis.

Refinement of Slope Stability Map
A slope aspect map of the study area was prepared, which was further edited as an east and west aspect map, by taking all pixels with values 0-180 as east and 180-360 as west. PSI points with positive values were extracted from eastern slopes, and PSI points with negative values were extracted from west-facing slopes. All PSI points belonging to a slope of 5 degrees or less were avoided as this area was deemed flat. A refinement matrix was created to refine the conventional FS map prepared using GIS-TISSA by taking into consideration the PS LoS velocities (V LoS ).
Overlaying of these two derivatives, i.e., one prepared using GIS-TISSA and other using V LoS resulted in the combinations of different slope conditions. This final landslideprone zone map includes three hazard classifications: unstable, critical, and stable.

Refinement of Slope Stability Map
A slope aspect map of the study area was prepared, which was further edited as an east and west aspect map, by taking all pixels with values 0-180 as east and 180-360 as west. PSI points with positive values were extracted from eastern slopes, and PSI points with negative values were extracted from west-facing slopes. All PSI points belonging to a slope of 5 degrees or less were avoided as this area was deemed flat. A refinement matrix was created to refine the conventional FS map prepared using GIS-TISSA by taking into consideration the PS LoS velocities (VLoS).

Monitoring Shallow Debris Flows
PSI results using the Sentinel-1A SLC SAR image are shown in Figure 5a. The small triangular cut-over portion is due to that fact that one Sentinel-1 image is composed of 27 smaller bursts. When such data is used for delineating the area of interest (AOI), AOI may overlap with multiple bursts. Thus, it is a pre-processing artifact. Figure 5a shows the V LoS (both positive and negative values). The PSI analysis provided 153,876 PS points, and each point yields average velocity values (mm/yr), incremental displacement values (mm) and total displacement values (mm). The velocity varies between −24.38 mm/yr to +28.22 mm/yr. 95-percentile coherence PS points were only selected for further analysis to avoid noise. The resulting PSI points were 9206 (Figure 5b). And these are the PSI points that were used in this study. The velocity of these points varies between −13.56 mm/yr to +11.28 mm/yr. The 12 landslide locations, collected from different sources and field verification (shown in Figure 1) were placed into three temporal groups: eight landslides in "Group 1", two landslides in "Group 2", and two landslides in "Group 3" and further classified into east and west-facing (Table 3).

Monitoring Shallow Debris Flows
PSI results using the Sentinel-1A SLC SAR image are shown in Figure 5a. The small triangular cut-over portion is due to that fact that one Sentinel-1 image is composed of 27 smaller bursts. When such data is used for delineating the area of interest (AOI), AOI may overlap with multiple bursts. Thus, it is a pre-processing artifact. Figure 5a shows the VLoS (both positive and negative values). The PSI analysis provided 153,876 PS points, and each point yields average velocity values (mm/yr), incremental displacement values (mm) and total displacement values (mm). The velocity varies between −24.38 mm/yr to +28.22 mm/yr. 95-percentile coherence PS points were only selected for further analysis to avoid noise. The resulting PSI points were 9206 (Figure 5b). And these are the PSI points that were used in this study. The velocity of these points varies between −13.56 mm/yr to +11.28 mm/yr. The 12 landslide locations, collected from different sources and field verification (shown in Figure 1) were placed into three temporal groups: eight landslides in "Group 1", two landslides in "Group 2", and two landslides in "Group 3" and further classified into east and west-facing (Table 3). PSI monitoring of landslides was done for these 12 landslides. Sentinel-1A on its descending mode characterizes the west-facing slope with negative velocity values, indicating down-slope or lateral movement to the west, and positive velocity values on eastfacing slopes, indicating down-slope or lateral movement to the east [44][45][46]. Buffering was carried out as the landslides were demarcated based on local reports and hence may not represent landslide initiation points. Buffering was performed based on the run-out length of landslides, which is 101 m. The east-facing landslide having west-face slope in the buffer was removed and vice-versa, for getting more accurate results. PSI monitoring of landslides was done for these 12 landslides. Sentinel-1A on its descending mode characterizes the west-facing slope with negative velocity values, indicating down-slope or lateral movement to the west, and positive velocity values on east-facing slopes, indicating down-slope or lateral movement to the east [44][45][46]. Buffering was carried out as the landslides were demarcated based on local reports and hence may not represent landslide initiation points. Buffering was performed based on the run-out length of landslides, which is 101 m. The east-facing landslide having west-face slope in the buffer was removed and vice-versa, for getting more accurate results.
Out of the three different temporal periods, the Group 2 category shows the maximum deformation value (St Mary's Orthodox Church with −35.67 mm) followed by Group 3 (Kamakshi village landslide with a deformation value of −28.03 mm) and Group 1 categories (NH-49 Salaf masjid landslide with −21.07 mm deformation). These three typical cases are shown in Figure 6. The rest of the temporal deformation graphs for other landslides are shown in Supplementary Material. All these landslides are west-facing. The maximum deformation recorded for the east-facing landslides is +10.11 mm/yr for Munnar-Gudaral Road landslide that belongs to Group 2. This was followed by +8.66 mm deformation for the Adimali-Rajakkad Road landslide that belongs to Group 3 category landslide. Group 1 has a maximum deformation of +8.53 mm, both for Painavu-Kattapana landslide and Thadiyampadu landslide. These deformations indicate that the maximum deformation is for the landslides that occurred during the time of PSI monitoring (Group-2), whereas the landslides that are going to happen in the future have smaller deformation compared to Group-2. It is also noted that sites that have experienced landslides in the past (Group-1) are also showing signs of deformation, which is indicative of reactivation and/or consolidation/settlement of landslides.
cases are shown in Figure 6. The rest of the temporal deformation graphs for other landslides are shown in Supplementary Material. All these landslides are west-facing. The maximum deformation recorded for the east-facing landslides is +10.11 mm/yr for Munnar-Gudaral Road landslide that belongs to Group 2. This was followed by +8.66 mm deformation for the Adimali-Rajakkad Road landslide that belongs to Group 3 category landslide. Group 1 has a maximum deformation of +8.53 mm, both for Painavu-Kattapana landslide and Thadiyampadu landslide. These deformations indicate that the maximum deformation is for the landslides that occurred during the time of PSI monitoring (Group-2), whereas the landslides that are going to happen in the future have smaller deformation compared to Group-2. It is also noted that sites that have experienced landslides in the past (Group-1) are also showing signs of deformation, which is indicative of reactivation and/or consolidation/settlement of landslides.

Mapping and Refining the Conventional FS Map
The conventional FS map prepared using GIS-TISSA (Figure 7a) was refined using PSI velocity values (Figure 7b,c). Since PSI measures the slant range displacement and velocities of ground points, landslide movements appear differently in the eastern and

Mapping and Refining the Conventional FS Map
The conventional FS map prepared using GIS-TISSA (Figure 7a) was refined using PSI velocity values (Figure 7b,c). Since PSI measures the slant range displacement and velocities of ground points, landslide movements appear differently in the eastern and western slopes over any satellite orbital geometry (ascending or descending). Since all the images downloaded were of the descending mode, the landslides are observed as negative velocity points in the western slope and positive velocity points in the eastern slope. Thus, the PS points were divided on the basis of slope aspect as east and west. All points falling in the flat aspect were avoided as this area is physiographically plain and free from landslides. The rest of the points, after the 95-percentile filtering, were extracted.
Upon extracting, the positive PS points in the eastern slope comprise 3852 points, and the negative PS points in the western slope amounted to 5023. It is to be noted that these numbers were generated after reducing the points on flat ground. The PS points in the east slope range from 0 to +11.28 mm/yr, whereas in the west slope, it spans from 0 to −23.03 mm/yr. This was divided into the equal class interval. The values between 0 and ±8 mm/yr are considered as low V LOS ; velocity values between ±8 and ±16 mm/yr as medium V LOS and the rest as high V LOS . Thus, the eastern slopes do not have a high V LOS class. These points were used to refine the corresponding GIS-TISSA pixels in the eastern and western slopes using a conversion matrix (Figure 8a). Such a matrix was used by Ciampalini et al. [47]. Pixels in the unstable class were not changed. In the critical class, pixels coinciding with high-velocity PS points were moved to unstable class. Pixels coinciding with medium velocity were unchanged, while those that coincided with low velocity were moved to stable class. In the stable class, pixels coinciding with high-velocity PS points were moved to the critical class. At the same time, those coinciding with medium and low velocity was retained in the same class. After refining, it was found that there was a marked decrease in the percentage of pixels in the critical class (Figure 8b). Most pixels in the critical classes were deemed as stable, thereby reducing the number of false positives.
Geomatics 2021, 1, FOR PEER REVIEW 10 western slopes over any satellite orbital geometry (ascending or descending). Since all the images downloaded were of the descending mode, the landslides are observed as negative velocity points in the western slope and positive velocity points in the eastern slope. Thus, the PS points were divided on the basis of slope aspect as east and west. All points falling in the flat aspect were avoided as this area is physiographically plain and free from landslides. The rest of the points, after the 95-percentile filtering, were extracted. Upon extracting, the positive PS points in the eastern slope comprise 3852 points, and the negative PS points in the western slope amounted to 5023. It is to be noted that these numbers were generated after reducing the points on flat ground. The PS points in the east slope range from 0 to +11.28 mm/yr, whereas in the west slope, it spans from 0 to −23.03 mm/yr. This was divided into the equal class interval. The values between 0 and ±8 mm/yr are considered as low VLOS; velocity values between ±8 and ±16 mm/yr as medium VLOS and the rest as high VLOS. Thus, the eastern slopes do not have a high VLOS class. These points were used to refine the corresponding GIS-TISSA pixels in the eastern and western slopes using a conversion matrix (Figure 8a). Such a matrix was used by Ciampalini et al. [47]. Pixels in the unstable class were not changed. In the critical class, pixels coinciding with high-velocity PS points were moved to unstable class. Pixels coinciding with medium velocity were unchanged, while those that coincided with low velocity were moved to stable class. In the stable class, pixels coinciding with high-velocity PS points were moved to the critical class. At the same time, those coinciding with medium and low velocity was retained in the same class. After refining, it was found that there was a marked decrease in the percentage of pixels in the critical class ( Figure 8b). Most pixels in the critical classes were deemed as stable, thereby reducing the number of false positives.

Discussion
The purpose of this study is to show the efficiency of PSI in monitoring and mapping shallow-seated landslides in tropical areas. Deformation caused by landslides, both local and regional scale, as well as both slow-or extremely slow-moving and fast-moving landslides can be studied by radar and can be evaluated using PSI [48]. In this study, 12 landslides, separated into three temporal categories based on PSI monitoring as Group 1, Group 2, and Group 3, which were further divided into east and west-facing landslides, were studied in detail to understand how PSI data can be utilized for studying shallowseated landslides. The analyzed PSI data spanned from March 2016 through July 2017.
Landslides of these three periods also show characteristic deformation rates. This is because the monitoring period covers all the sequence of development of shallow-seated landslide viz. the monsoon season, the subsequent increase in pore-water pressure, and the final failure. The west-facing slopes showed the largest deformation values (e.g., St. Mary's Orthodox Church landslide with a value of −35.67 mm), whereas the maximum for the east-facing landslides is +10.11 mm/yr for the Munnar-Gudaral Road landslide. Both belong to the Group 2 landslides. The second largest deformation values demonstrated category is Group 3 landslides. These deformations indicate that the maximum deformation is for the landslides that occurred during the time of PSI monitoring (Group 2), whereas the landslides that are going to happen in the future (Group 3) have smaller deformation compared to Group 2. It is also noted that sites that have experienced landslides in the past (Group 1) are also showing signs of deformation, which is indicative of reactivation and/or settlement/compaction of landslides. The West-facing Kamakshi village landslide with −28.03 mm and the east-facing Adimali-Rajakad Road landslide with +8.66 mm deformation belongs to Group 3. Lastly, in the category of Group 1 landslides, none of the developmental stages of the shallow-seated landslide were under surveillance, and hence the group was characterized by the least deformation value. NH-49 Salaf Masjid landslide with −21.07 mm deformation and +8.53 mm for Painavu landslide are in the

Discussion
The purpose of this study is to show the efficiency of PSI in monitoring and mapping shallow-seated landslides in tropical areas. Deformation caused by landslides, both local and regional scale, as well as both slow-or extremely slow-moving and fast-moving landslides can be studied by radar and can be evaluated using PSI [48]. In this study, 12 landslides, separated into three temporal categories based on PSI monitoring as Group 1, Group 2, and Group 3, which were further divided into east and west-facing landslides, were studied in detail to understand how PSI data can be utilized for studying shallowseated landslides. The analyzed PSI data spanned from March 2016 through July 2017.
Landslides of these three periods also show characteristic deformation rates. This is because the monitoring period covers all the sequence of development of shallow-seated landslide viz. the monsoon season, the subsequent increase in pore-water pressure, and the final failure. The west-facing slopes showed the largest deformation values (e.g., St. Mary's Orthodox Church landslide with a value of −35.67 mm), whereas the maximum for the eastfacing landslides is +10.11 mm/yr for the Munnar-Gudaral Road landslide. Both belong to the Group 2 landslides. The second largest deformation values demonstrated category is Group 3 landslides. These deformations indicate that the maximum deformation is for the landslides that occurred during the time of PSI monitoring (Group 2), whereas the landslides that are going to happen in the future (Group 3) have smaller deformation compared to Group 2. It is also noted that sites that have experienced landslides in the past (Group 1) are also showing signs of deformation, which is indicative of reactivation and/or settlement/compaction of landslides. The West-facing Kamakshi village landslide with −28.03 mm and the east-facing Adimali-Rajakad Road landslide with +8.66 mm deformation belongs to Group 3. Lastly, in the category of Group 1 landslides, none of the developmental stages of the shallow-seated landslide were under surveillance, and hence the group was characterized by the least deformation value. NH-49 Salaf Masjid landslide with −21.07 mm deformation and +8.53 mm for Painavu landslide are in the Group 1 landslides. Bianchini et al. [48] and Cigna et al. [49] also used PSI data for describing the state of activity of landslides in Calabria Region, Italy. Similarly, Mateos et al. [50] used PSI as a measure to assess the damage caused by landslides in Granada, Spain.
An attempt to improve the existing landslide susceptibility map using PSI data was also attempted in this study. Accuracy of conventional landslide susceptibility maps needs improvement as it tends to have a large area of false positives and negatives. The method deployed to improve the accuracy, i.e., melding of GIS-TISSA derived landslide susceptibility map and PSI data, resulted in the refinement of the landslide susceptibility map. The ultimate result ended with a reduction in the critical category from 17 to 2 percentage points. This will help in efficient planning of disaster risk reduction activities.
Combining the PSI information with the slope stability model has significant practical implications. Traditionally, slope stability analysis is performed using physics-based [7] or statistical models [51]. Moreover, the models developed using these approaches are highly conservative, leading to several false positives. Because of the large false positives, the community often loses its confidence in these models and get desensitized to any early warning or risk predictions based on these models. Therefore, any opportunity to refine the slope stability model will enhance community confidence in risk estimations based on these models. The proposed approach that we have demonstrated in this study, combining PSI and the slope stability modeling, is a valuable step to achieve that goal. The combination of PSI with the stability model helped to refine the stability model output. The refined output (by eliminating false positives) also has significant economic implications as well. If the agency decides to monitor or mitigate the landslide risk at the critical location, then a refined output based on PSI and stability model would require fewer locations monitored or stabilized, which will provide cost savings.

Conclusions
This study is an attempt to monitor and map shallow-landslides set in a tropical set-up. The study resulted in the following conclusions: i.
The three types of temporal landslides pertaining to different periods of PSI analysis resulted in providing three distinct PSI velocity values. These refer to landslides that occurred before the analyzed time period, during PSI analysis, and after the analysis. ii.
Traditional slope stability analysis methods often lead to significant false positives. The proposed approach combining PSI with traditional slope stability shows promise in reducing these false positives. The PSI derived velocity was used to refine the existing landslide susceptibility map. Such refinement will help to identify areas that require utmost priority in terms of monitoring or mitigation. iii.
Developing countries like India will only have limited financial provisions for adopting management practices. The freely available radar data and the provisions of identifying areas requiring urgent management are able to be deciphered through studies like the one mentioned here. iv.
Such an approach could also help the engineering community gain better community confidence as the slope stability model is more targeted and less susceptible to false positives. Institutional Review Board Statement: Not applicable.

Informed Consent Statement: Not applicable.
Data Availability Statement: Data will be shared based on request.