The Use of TERRA-ASTER Satellite for Landslide Detection

: In August 2005, numerous shallow landslides occurred in the region of Vorarlberg (Austria), particularly induced by unfavourable event-related weather conditions. Two scenes of TERRA-ASTER sensor were used for the identification of the vegetation change induced by the landslides. The focus of this study is the establishment of a reliable method, comparable to aerial-photo visual interpretation standards, able to identify accurately landslides by processing a series of medium-resolution remote sensing optical data, before and after a catastrophic event. A very intuitive workflow for a semi-automatic image classification for the detection of landslide-induced change on the image data is proposed. The accuracy and validation assessment was carried out by means of a landslide (aerial-photos derived) inventory. By taking into account the central area of investigation, the landslide detection method, which adopted an innovative double classification workflow (a first supervised followed by an unsupervised algorithm), delivered a very high producer accuracy (81.5%) coupled to a more-than-acceptable user accuracy (68.9%) and kappa coefficient (72.9%).


Introduction
Landslides are one of the widest spread types of geohazard all around the world and are mostly caused by preparatory triggering factors such as earthquakes [1], human activities [2] and extreme rain episodes [3]. Such gravitational mass movements occur nearly everywhere in Austria; in fact, in 2014, it accounted for 306 events (67 rockfalls and 239 landslides) which left five deaths and 10 injured people together with extensive damage to public and private properties. In the context of disaster management, a rapid and representative overview of the areas hit by hazardous events [4] must integrate surficial mapping, which enables a detailed classification and delimitation of the phenomena, with aerospace-borne image analysis [5,6] with its ability to provide regional scale landslide detection and accessibility to remote areas [7]. Since the launch of LANDSAT in the early seventies, optical image analysis applied to the recognition of landslide scars and deposition areas is been used as a catastrophe assessment tool [8] and proved to be useful for the update of landslide inventory archives [9,10]. Although those early works established a series of standard methods for mass movements' detection [11], the limitations of the minimumrecognizable pixel size lay on their spatial resolution. Besides the limitation of optical imageries such as clouds, shadows caused by either meteorological conditions or elevation exposition, the use of optical remote sensing analysis and image classification is very indicative of very rapid to moderate mass movements such as flow-like landslides, rock planar slides, debris slide, and rotational slides [12]. In recent years, several very high resolution (VHR) satellites have been launched and, due to their fine spatial capacity to resolve pixels, more and more scientists have become familiar with object based image analysis (OBIA) techniques applied to landslide recognition [13,14]. Since April 2015, SENTINEL-2, a new generation of high resolution satellites and high temporal resolution, is in orbit and thanks to an agreement between ESA and the Copernicus Consortium the access to this data is been secured open and free to everyone. TERRA ASTER images have provided, since 2000, with their 15 m spatial resolution (for the VNIR bands), a very suitable sensor for landslide scar and trail detection [15]. Most of the approaches used to address mass movement recognition were change-detection algorithms [1,16] and pixel-based supervised classification [8]. Several studies highlighted the limits of pixel-by-pixel against OBIA methods [17,18] and some addressed the attention as well to the use of unsupervised classification as a viable way to overcome landslide fragmentation detected by supervised analysis [19]. Another limitation of the classical supervised classification relies on the impossibility to directly remove spectral analogues classified pixels such as bare land, river beds, build-up areas and agriculture. A way to face such constraints is to use in the postclassification phase updated ancillary data, such as stream network [4] or digital elevation model derivatives [20,21], to improve the method's mapping accuracy. A new type of paradigm, connected to the availability of new sensors, now implies to develop a mixed approach, which will integrate the quickness of pixel-based algorithms with the textural and structural information exploitation capability offered by OBIA classifiers. The "super-pixel" approach proposed by Wu et al. [22] represents such an example.
The main motivation behind this work is the assessment of the quality of remote sensing optical data (with medium resolution) beside the aerial photos (with very high resolution) as a complementary tool to promptly identify, on a regional scale, spontaneous and reactivated landslides driven by extreme weather conditions. The aim of this paper is to define an intuitive workflow for a semi-automated pixel-based image classification combined with a GIS multi-parameter data handling for the detection of landslide-induced by a heavy rainfall event. The change detection is aiming at emphasizing the spectral response in vegetation change in the pre-and post-event images due to the abrupt removal of vegetation induced by the manifestation of landslide activities. The method used a pairs of Terra-ASTER images and was applied to a study case located in the region of Vorarlberg, Austria. The main difference between our research and previous studies is that a double classification scheme was applied to a bi-temporal classical pixel change detection approach for landslide recognition. This classification scheme is composed of two steps: 1. supervised classification via maximum likelihood classifier with a set of three indexes of spectral change. 2. unsupervised classification with another set of three indexes of spectral change.

Study Area
The study area of interest is located in the eastern part of the Austrian federate state of Vorarlberg ( Figure 1). The eastern part of the Voralberg region included in this study is characterised by the presence of the two main rivers of the region the Ill and the Bregenzer Ach both flowing from south-east toward north-west. The overall elevation, which varies between 580 and 2900 m, and the main slope angles, varying between 20° and 40°, plays an important role on the proneness to landslide. The day of the catastrophic event the local weather station registered 97 mm of rain which represents half of the amount (based on 30 years averaged measurements) of precipitation occurring during a typical month of August. The formation of a strong low-pressure system originated over the Adriatic Sea brought a large amount of moisture over the Balkans, Austria and southern Germany [23]. This so called 5-b weather condition led to long-lasting phases of heavy rain (4-5 days), resulting in the alpine floods events of 2005 for the northern Pre-and Central Alps [24].
The disaster occurred around the 19-23 August 2005, affected almost the entire Vorarlberg province and caused a series of damages which were calculated summing up to 100 million Euros [25]. The cause of property loss was related not only to landslides but also to erosion and flooding. In fact, several large rivers exceeded the 100-year floodreturning period (HQ100). Based on a previous work [26], a landslide inventory from multi-temporal digital aerial-photo interpretation (2005)(2006) is being used for the delimitation of the area in exam and for the validation of the classification method. A major limitation of the aerial-photo derived interpretation, represented by the fact that the aerial-photo 2005 coverage was not extended to the whole region but limited to the valley floor, was overcome by using the available ASTER data. The northern part with the molasses (pink), the penninic and helvetic zone (dark and light green); • The central part with the northern calcareous alps (blue); • The southern part with the central alps and penninc metamorphic zone (purple).
In the northern area, soft sedimentary and clastic deposits are the main lithologies prone to spontaneous and quick landslides. The central area is more subject to rock fall and rockslide (black and white points in the landslide inventory in Figure 2) due to the presence of limestone and dolomites, whereas darker and harder rocks of methamorphic origin can be found in the southern part of the study area. More in the south, the landscape is dominated by high relief cliffs, narrow and U shaped valleys, and periglacial phenomena such as rock glaciers and solifluction are common. On the other hand, heading toward the north of the study area, the landscape was reshaped by fluvial activity, the relief is flatter and the incidence of rock fall and rockslide is less dominant. By visualizing the local land cover map ( Figure 2) it is clear that the dominant cover types are haerbaceous and forest in the north, whereas in the south the reduction in vegetated areas is replaced by the presence (mostly where the relief is high) by bare rock and shrubs. As already shown in Figure 2, there are six classes of types of landslide in the reference inventory and the main classes earth flow, debris slide and earth flow /debris slide comprise the majority of them (with 94% of the total). The visually interpreted inventory was created by manually digitizing, at the scar locations, points, reunited in a geodatabase, which stores the type of process and activity class of a series of landslides visually identified in a multi-temporal set of aerial photos (1930, 1950; 1970, 1990, 2001, 2005, 2006, 2009, 2012). By means of the database, it is possible to infer that 86% of the landslides were newly triggered phenomena that mostly affected loose and earth-like sediments. Further polyline features such as scars, scarps, and crowns together with polygons delimiting single landslides are not part of the inventory.

Challenges and Chances in Remote Sensing Landslide Recognition
A series of demands are related to the use of remote sensing applied to landslide detection in mountainous and urbanised environments. Landslides are commonly associated with prolonged heavy rain episodes, which lead to the water saturation of the shallowest layer of soil and regolith and an increased groundwater level within the slope [27]. However, under such mentioned circumstances, the interaction of those slope instabilities with rivers (which act as undercutting agents) and torrents (which transport the gravitational solid mass originated on his bank along pre-existing channels) must be taken into account when satellites images are used for landslide identification. One of the challenges faced by change detection classification is the not-straightforward ease of separating processes such as earthflow / debris slide from debris flow and debris flood due to their similar spectra signatures. The main difference between the two sub-groups already mentioned is that, although earthflow and debris slide do not share the same slope failure mechanism, the initiation of debris flows and debris floods, unlike the first group, is normally due to strong interactions with water along a steep channel [12].
In Figure 3 is shown the area of Warth-Teschenberg on the river Lech where two earth flows (in yellow) were activated in 2005; on the other hand, two debris slides (in white) on the river flank were already visible on the aerial photo from 2001. In this case, the distinction between the debris flood action and the potential reactivation of previously originated landslides on the river flank represents a challenging task. The augmented temporal resolution in terms of revisiting time of the newcomer's satellites plays an important role, because the closer an image is to the date of the event, the better the chance to identify a phenomenon fresh and untouched by remedial measures. The large scene extension of satellite data represents a better chance in remote sensing to cope with the problem of aerial-photo gaps ( Figure 4) which are in fact remarkably common in case of a catastrophic event imaged by airborne acquisition. An important challenge is the evaluation of the ability of remote sensing method to detect the complete phenomenological expression of the landslide (scar, transport and accumulation zone). A typical high-resolution optical satellite offers a 10-15 m pixel size quality in images, which is why a large and destructive landslide such as the Danöfen earth flow (Figure 5a) is normally very well recognized in the overall geomorphological manifestation and geographical extension. On the other hand, for small and shallow earth flows (see Figure  5b) the degree of completeness and fragmentation of the classified pixels varies depending on several factors, such as: • the portion of denudation involved into the process; • the time span between the two images taken before and after the event; • the length of the transport trace; • the width of the landslide scar; • the degree of isolation of the single phenomena; • the suitable signature contrast between the land cover type and the landslide.

Data and Methods
The data used for this work were Terra-ASTER L2 products. The main advantage of the ASTER L2 Surface Radiance image in comparison to the L1 product is that the former is atmospherically and crosstalk corrected [28] for both the Visible Near-Infrared (VNIR) and Shortwave Infrared (SWIR) sensors ( Figure 6). The only drawback is that the last five TIR bands were not available and a better cloud masking algorithm, which needed the thermal bands [30], could not be performed. For this study the change analysis classification was applied to the images acquired on 17 September 2004 and 4 September 2005 only. The methodology used in this work is the standard workflow required for a pixel-based change detection analysis applied to optic satellite images, with the only exception of the final stage of the processing and for the whole post-processing. The entire procedure, illustrated in Figure 7 and computed in ArcGIS, can be subdivided in the following sub-steps: 1. Image pre-processing; 2. Image processing; 3. Image post-processing.

Image Pre-Processing
The initial stages of the ASTER image preparation require a series of standard routines which have the purpose to optimise the inter-comparison between different bands in space (same spatial resolution) and in time (different acquisitions) and to create masks of unwanted information such as clouds, snow and shadows which are caused by topography, vegetation and clouds. The first step consists of the radiometric correction of the bands, which is obtained by transforming the radiance L (expressing the at-surfaceradiance) into surface spectral reflectance R by calibrating the nine bands using Equation (1) [31]: where L is radiance, R is reflectance, d is the sun distance in astronomic units, Esun is the Exoatmospheric Solar Irradiances and Z is the sun elevation angle. The next step requires the resampling of the six bands from 30 m resolution, corresponding to the SWIR spectrum, into the finer spatial resolution of the VNIR bands (15 m). The method used was the nearest neighbourhood resampling technique. Furthermore, the georeferencing was preferred in comparison to the orthorectification in order to retain the original digital number pixel information useful for further change detection analysis. In fact, by means of the ArcGIS-10.5 georeferencing tools option "update georeferencing", a user can avoid to interpolate an image, whereas in performing orthorectification the original digital number is always resampled. The georeferencing is quite important because incorrected misalignment of two images can lead to false differential classification or incorrect validation in case of the aerial photos and sample inventory validation. The application of several Ground Control Points (GCPs) for each image in order to obtain a reasonable root-mean-square error (RMSE) of 15 m was required. This is not the case for the georeferencing of the SPOT-5 and Sentinel-2 (1C product) images because all single tiles have the same (i) centre position and (ii) constant viewing angle, and (iii) with a unique set of GCPs, it is possible to process a series of several images. Furthermore, a standard technique was used for snow identification. The Normalised Difference Snow Index (NDSI) based on a ratioing between Band 1 and Band 4 was adopted for this task, and by imposing a threshold of 0.4 [32] the masking of snow cover and hydrological features (such as lake and rivers) from the image was obtained. A further step required the creation of a False Colour Composite (FCC) 3-2-1 which gives results similar to traditional colour infrared aerial photography. This FCC is very useful for cloud and cloud shadow supervised classification. A further filter (Equation (2)) for the detection of cloud shadows [30] was applied by imposing the following threshold to NIR band (B3) and the blue band (B1) in the ArcGIS Raster Calculator:

Image Processing
A very good index for land cover change analysis between two images is the Normalised Differential Vegetation Index (NDVI). This index takes into account the extreme difference in reflection for the red band (B2) and NIR band (B3) for green vegetation areas and bare ground and is expressed by Equation (3): High index (near +1) means very vigorous vegetation, near -1 indicates water (ice) and near 0 bare soil, asphalt and built-up areas. The NDVI is very sensitive to rapid changes in land cover, and for this reason when a superficial landslide occurs a loss in vegetation vigor can be accurately detected. Another useful algorithm used in this work is the Tasseled Cap Transformation (TCT). The TCT is a tool able of compress spectral data into a few bands associated with physical scene characteristics. Both multi-and hyperspectral data have highly correlated bands [33]. As slope deformation by superficial landslides is a type of phenomenon influencing the denudation process in perennial (grass and woodland) and non-perennial (arable lands) vegetative areas, TCT was chosen as a tool for analyzing the time series ASTER images of the satellite-based study area. TCT indices give a measure of the brightness (soil), greenness (vegetation), or wetness (water) of a pixel and utilize a linear combination of nine ASTER frequency bands. A way to perform such a calculation is to use the tasseled cap algorithm expressed by Equation (4 where tascapi is the calculated tasselled cap index for brightness, greenness, or wetness depending on the coefficients used, the bands are the Top of the Atmosphere (TOA) reflectances and the coefficients are given by Yarbrough et al. [31] in Table 1. Successively, a method called image differencing [1], subtracting date one of the transformed remote sensing data to date two, was adopted. For this purpose, the three TCT components and the NDVI were combined to create two different image composites (Table 2): • The first with the three TCT differentiated images as RGB components; • The second with NDVI instead of change of surface greenness. A supervised classification (SC) method was adopted to classify the landslides evident on the second image composite (Table 2). More in detail, a number of 24 samples were chosen as landslide signatures and a group of 575 samples were labelled as nonlandslides. The supervised classification was run in ArcGIS Maximum Likelihood Classifier with "reject fraction" = 0.01 and "A priori probability weighting" = SAMPLE. In order to refine and improve the first classification, a further unsupervised classification (UC) algorithm of the first image composite in Table 2, imposing 15 classes and a 100 minimum pixel grouping, was adopted. A conditional statement on the selection of the most representative landslide classes in UC, which are matching the landslide recognised in SC, was imposed and a double classification (DC) was obtained by combining SC and UC.

Image Post-Processing
On the post-processing stage, the first step concerns the transformation of the rasterclassified image in vector as point excluding the pixels classified as non-landslides for the SC image and for the DC image. Two available DTMs (Digital Terrain Models) based on Airborne Laser Scanning (ALS) were used for the flat terrain correction: • a nationwide 10 m resolution ALS DTM; • a local ALS DTM Vorarlberg (for the core study area only) with 1-m resolution.
The local DTM was resampled to 5 m and then by using the publicly available building shape file for the Open Street Map Project [35] a correction of the artefacts present into the derived DTM slope model was performed. In this sense, the vector building shape was transformed in raster with a fixed value of three degrees for each pixel, and then this raster layer was fused to the DTM slope model. For the DTM 10 m, the same approach was applied. A further extraction to points of the slope values with a threshold greater than 6°, in order to remove unprovable candidate landslide pixels, was carried out. Furthermore, those two datasets of points were visually investigated using the pre-and post-disaster images as references and edited in order to eliminate misclassified groups of points, which were not correctly masked as clouds and shadows on the pre-processing stage. In this way, failed masked part of clouds and shadows that were incorrectly classified as landslides could be easily eliminated. Since false-positive landslides resulted from the widening of riverbeds [4], a buffer of 20 m around the river network was applied and those landslide candidates were classified as "Fluvial-torrent process". A further processing of the twopoint datasets required the compilation on the field value of a new class of code (Table 3).

Classification Method
The landslide classification accuracy was evaluated by means of the confusion matrix. A confusion matrix (Figure 8) is compound of a square matrix of numbers identifying the pixels in exams, which belong to a particular category. In our case, the two categories taken into account for are landslides and no change. A series of statistical measures of thematic accuracy are assessed by the matrix, such as the producer's accuracy, user´s accuracy and the Kappa coefficient (KHAT). The KHAT is a multivariate index which uses all the information in the confusion matrix and takes into account not only the overall accuracy but also the random accuracy. The random accuracy and the overall accuracy are calculated using the following equations (Equations (5) and (6) The KHAT, expressed by Equation (7) [36], for Landis and Koch [37] can be interpreted using the following main categorisations: • 0.81-1 indicates excellent agreement; • 0.61-0.8 indicates substantial agreement; • 0.41-0.6 indicates moderate agreement.

Results
The usual way to verify the effectiveness of the classification methods used is to create a confusion matrix in order to evaluate the classification accuracy reached. From our perspective it is also fundamental to evaluate how the detected landslides in different areas are representative in pixels' extension and position in comparison to our validation points identified as landslides. In order to judge the limitations and the strengths of the methods a nested approach, with subdivision of the whole image in four overlapping areas of interest (Figure 9), was proposed:  In Figure 10 the overall landslide detection map for the region of Vorarlberg, by taking into account the results obtained by applying the double classification method, is shown.

Landslide Classification Accuracy
All the pixels of the eight datasets in the four overlapping areas were assessed by means of the confusion matrix ( Figure 11). At this stage, it is important to evaluate the effectiveness of the DC process. In fact, apart from the satellite-based study area where the KHAT shows a moderate agreement with 48.2%, in the other study areas this statistic scored over 55%. The best scores were obtained in the central study area (KHAT = 79.3%) where the single classification with a KHAT of 70.9% also bore reliable accuracy. The SC for the slightly larger study area shown as well an almost substantial agreement (KHAT = 54.8%) with a producer accuracy for landslides of 81.7%. Another way to read the producer accuracy is to calculate the mutual omission errors. In larger areas the probability to miss small landslides is higher than in delimited areas, and, apart for case F, the single classification scores always better. For all of the four DC cases the producer accuracy is acceptable. Controversial instead it is the interpretation of the commission errors. The commission errors come primarily from the misclassification of fluvial-torrent process as common landslides and due to the similarity between landslide and fluvialtorrent process signatures.  Figure 9.
A further source of error is due to the time span between the first image and the second. The image was taken about at the same time one year before but the pixel appearance under different sun illumination effects, especially in high relief areas, the capacity of the classifier to assign the correct not-landslide class to the image. The consequence of the latter is that in high relief areas the labelled signature landslide can be misinterpreted with denudation processes or (in the majority of the cases) completely misclassified. Other sources of errors are new human activity detected which provoked a change in vegetation (forest clearance, new built up areas, forest-road, new infrastructures) and minor phenology variation affecting cultivars in specific crop seasons.
The good quality in central study area F is a consequence of the fact that there, a minor number of fluvial-torrent processes were documented in comparison to the rest of the study areas. A further remark concerns the river network extension which in this area it is not so preeminent. Those circumstances leaded to a lowered commission error. By removing the class fluvial-river process from the classifier ( Figure 12) several indexes were improved. For example, concerning the KHAT, it is possible to score for the three DC classified datasets almost an excellent agreement.
Furthermore, for a regional-scale classification (data G) this method demonstrates good potential for mapping landslides, obtaining a good 28% of omission errors and a substantial agreement performance with a KHAT of 59%, once the river processes are removed from the classifier.

Landslide Validation
The landslide validation has a duplex purpose: on one hand, to verify how many landslides the classifier was able to detect, and on the other hand, to quantify how well those pixels represent spatially specific landslide types. This stage of the study was conducted only on the three well-documented study areas.
From the original process inventory, a selection based on spontaneous landslides activated between 2005-2006 in the study area, in the aerial-photo of 2005, and the central study area was created. Furthermore, two new fields, called SAT_VAL and SAT_POS, were added and coded as from Tables 4 and 5: Table 4. ASTER validation assessment.

SAT_VAL-Satellite Validation Code
Class 0 No change detected 100 Landslide 200 Torrent process 300 Vegetation change 400 Anthropogenic change 500 Erosion (unspecific process) 600 Shadows 700 Clouds Table 5. ASTER position assessment. For this validation procedure, the code = 100 of SAT-VAL were assigned to the reclassified pixels which corresponded to the original inventory point data. Furthermore, some new points as new and reactivated landslides and torrent process, most of all on the areas not covered by the aerial photo of 2005, were added. In Figure 10, in fact, the biggest number of new and reactivated landslides was identified in the D areas where for almost the half of its surface only the 2006 aerial-photo was available. Here, between 50 and 110 new or reactivated landslides were detected by the used satellite-based approach. Despite those new findings the validated landslides are still under 50% with the DC datasets and under 70% for SC datasets (Figure 13). The second field SAT_POS was used to estimate the mapping accuracy and positioning of the pixels in the classified image. The majority of good quality data are found in the study areas D and B (Figure 14) but by looking at Figure 9 it can be clearly seen that area B, in terms of validation number and representative landslides, performs best considering as well the respectable KHAT. The results from this evaluation, which identify in each study site and for both datasets available accurate and precise landslides, are very important since they could be applied in the future for the creation of a landslide-susceptibility map of the same area.

Discussion
This study case focused on the 2005 catastrophic event in Vorarlberg shows some improvement on the producer and user accuracy, if only focused on the results delivered for the sub-areas (A, C and E), when compared to a previous attempt to adopt an objectbased image analysis method on landslide detection in different areas of Austria and Italy [38]. Moreover, another similar approach proposed by the same author for a regional case study [39], thanks to the higher revisit time of 5 days and to a multi-temporal ndvi profiling method, produced a higher kappa index score and most of the open/barren land could be discarded as landslides. In this sense, the use of new satellites with higher revisit times can solve the issue related to the high false positive rate documented in this study, especially on the regional scale attempt for the areas G and H. Considering that a similar high intensity precipitation episode occurred recently only on the 11th August of 2002 [40], it is ruled out that the false positives found could have been associated to another event that occurred between September 2004 and July 2005.
The work already mentioned [26] was addressed to the visual recognition and characterisation of different types of landslides in the region of Voralberg, being the focus of our group the localisation of the main scar for each landslide. This approach is generally driven by the aim of applying point source data to generate multi-variable landslide susceptibility maps at different scales [41].
In order to highlight visually the performance of our landslide detector, a series of figures will be used as showcases. Some explanation is required to enhance the reading of the following figures: • the image shown beside each aerial photos (in Figure 15, 17, 18 and 19) represents the RGB composites made of three index of spectral change (introduced in Table 2); • the bigger points with a black contour are the validation data whereas the smaller red points without a contour are the classified pixel as a landslide with our method; • the minimum mapping unit defined in our approach corresponds to two pixels (450 m 2 ) belonging to the same phenomenon, even if they are not contiguous.
The main classes established during the compilation of the landslide inventory were: 1. Debris slide: such as the two illustrated (with white points) in Figure 15 show higher contrast with the surrounding areas and shorter to nonexistant trail from the scar to the accumulation area. 2. Earth flow: such as the two illustrated (with yellow points) in Figure 15 show lower constrast with the surrounding areas but longer trail which indicates a flow-like movement over the topography. A distinctive characteristic of earth flows is the very unusual transport forms that they normally leave behind. 3. Mixed class: in the special case of a landslide difficult to distinguish due to the location of the process (forest land cover), to the quality of the aerial photo (haze, cloud or shadow) or to the complexity of the process (e.g., landslide initiated as a debris slide but evolved into an earth flow) a mixed class was adopted instead.
As already mentioned in Section 3, the major problem affecting landslide classification methods based on spectral signature analysis is the amount of water present in the mobilised debris at the time of the event. In our case by using, as verification for the model, only the three types of landslide defined above, we reduced the possibility of classifying debris flows and debris floods as gravitational phenomena affecting instable slopes.
Concerning the classification accuracy for each class of landslide, a preliminary conclusion can be drawn that, for example, for area C, 75% of the earth flows were identified with a good score (at least 300, see Table 5) in position accuracy and 87% of the total mapped landslides belong to this class. In area G, instead the percentuage of landslides classified in the inventory as earth flow is 54% of the total; moreover, here 80% of the correctly identified landslides correspond as well this class.
By summarising those outcomes in the two major sites investigated by this method and bearing in mind that only landlides bigger than two minimum mapping units (900 m 2 ) can be recognized, the majority of slope instabilities mapped were classified in the inventories as earth flows. Even though that consideration could be interpreted as a disadvantage, earth flows are rarely detectable one year after an event with aerial-photo visual identification, whereas the other two classes still maintain a distinct morphology at high resolution [26].
The most interesting example of the positive impact of ASTER images for landslide detection can be found in the area around the village Au that lies in the central study area (data E and Data F in purple) (see Figure 15). Here, a very conspicuous number of landslides were detected and positively validated (as shown in Figure 15c in the landslide inventory process legend). In this example the second biggest earth flow was detected, spatially, almost completely (77% of the total) and validated as extremely good (class 100 in Table 5) in position accuracy. Moreover 28 landslide correctly recognized here corresponds to the inventory class earth flows, whereas six corresponded to debris slide and only three to mixed class. If the double classification results are considered here, 62% of the total landslides are detected correctly and 82% of them score high in position accuracy.  Therefore, missed information retrieval from aerial-photo analysis can be overcome by optic satellite pixel-based landslide classification. Another interesting example ( Figure  17) is the area around Warth-Teschenberg on the river Lech where two landslides (in yellow) activated in 2005 were correctly identified. For those two earth flows, although the minimum mapping units are respected, the very low location accuracy prevents any further applications. On the other hand, the two old debris slides on the river bank were wrongly classified by the detector as active landslides; instead, both of them should have been categorised as river-debris floods. Figure 17. Landslide classification results in the area of Warth-Teschenberg on the river Lech (as seen in Figure 2). Figure 18 showcases the different degrees of detectability for a group of six landslides. In fact, in Hintereggalpe, out of six cases, only three types of landslides could not be completely detected. In two cases, highlighted with the red ellipse, their surfaces are too tiny (areas smaller than 900 m 2 which correspond to two minimum mapping units), whereas in the yellow circle the close proximity of two landslides allowed the classifier to identify only the one with a longer trail. A good outcome for the validation for this study case is that the three classified earth flows score very high in satellite position and accuracy scores, since for all of them the scar is very recognizable.  Figure 4b).
In Figure 19, the massive Danöfen earth flow represents the biggest landslide detected in the area by the method. In this case, it is possible to detect the whole extension and to assess extremely positively the position of the classified pixels. In fact, 71% of the total area of the landslide can be detected correctly and 80% of the pixels are inside the digitized polygon. In Figure 20a the results of the devastating earth flow and in Figure 20b its impressive perspective view from the valley floor on a field photo taken the day after the event can be seen.
At last, a consideration can be made concerning the impact of different land cover types and different lithology on the classification performance.
In fact, a distinction between study area D and G (the whole Vorarlberg) can be made. Althoght for both sites, the percentuage of grassland on the total remains relatively similar, forests decrease for area G, giving more room to bare rock and shrubs ( Figure 2) which are not suitable land cover types for landslide recognition based on NDVI changedetection techniques. Another interesting aspect is rock contrast of the metamorphic against sedimentary, with it being easier on lighter sediments and rocks to visually identify the landslide. The same statement cannot be made for the remote sensing method that relies on change in vegetation index and not in change of rock-vegetation contrast.

Conclusions
ASTER images were capable of detecting up to 82% of landslides mapped by visual cognitive interpretation of aerial-photos on a 30-times-finer resolution; more in detail, a combination of pixels classified with supervised classification SC and double classification DC (combination of SC and unsupervised classification) are particularly suitable for landslide detection at a sub-regional scale. It can be stated that the datasets for the 2005 aerial-photo extent B and central study area F can be entirely exploited, whereas in case of the whole study area (C-D) is better to rely on dataset C (with the KHAT and user accuracy in D not completely reliable). More in detail, if we combine dataset B with C and F we obtain a more reliable landslide classification for the whole study area (C-D). The conspicuous number of 110 landslides which are correctly identified and which are enough accurate in position and extension represents a very good dataset as input data for modelling landslide susceptibility maps or other hazard oriented maps. After having assessed DC pixels, another major achievement is that ASTER can be used for landslide recognition on a regional scale. In fact, by removing the class fluvial-river process from the statistics, on a regional scale an accuracy of 59% landslide recognition can be achieved. A large number of false negative pixels are detected on steep slopes also where slightly different sun elevations and viewing angles generate not-comparable shadowing effects. Those problems arise from the lack of a topographic correction performed on both scenes in exam. Sometimes the shadow and cloud masking were not efficient due to the lack of the thermal bands on 90 m resolution usually available on TERRA ASTER 1B products. The advanced ACCA cloud detection therefore could not be performed in this study. Some of those limitations could be overcome with the use of SPOT and SENTINEL-2 by applying topographic correction, better alignment of the scenes which helps to avoid geometric errors, and thanks to a higher temporal resolution (5-10 days) which allows the user in the post-processing stage to reduce the number of false positives. The feasibility of adoping a similar approach is already being investigated, based on high-resolution data, to detect landslides occurred over the course of other catastrophic event in Austria, such as Taxen

Funding:
This research was partially funded by the Österreichische Forschungsförderungsgesellschaft FFG in the Austrian Space Applications Program (ASAP 11) through the project "Land@Slide" (contract no: 847970)".

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author. The data are not publicly available because not yet being included completely in our landslide geo-database.

Conflicts of Interest:
The authors declare no conflict of interest.