On the Use of Single-, Dual-, and Quad-Polarimetric SAR Observation for Landslide Detection

: Remote sensing technologies, particularly with Synthetic Aperture Radar (SAR) system, can provide timely and critical information to assess landslide distributions over large areas. Most space-borne SAR systems have been operating in di ﬀ erent polarimetric modes to meet various operational requirements. This study aims to discuss how much detectability can be expected in the landslide map produced from the single-, dual-, and quad-polarization modes of observation. The experimental analysis of the characteristic changes of PALSAR-2 signals showed that quad-polarization parameters indicating signal depolarization properties revealed noticeable landslide-induced temporal changes for all local incidence angle ranges. To produce a landslide map, a simple change detection method based on characteristic scattering properties of landslide areas was proposed. The accuracy assessment results showed that the depolarization parameters, such as the co-pol coherence and polarizing contribution, can identify areas a ﬀ ected by landslides with a detection rate of 60%, and a false-alarm rate of 5%. On the other hand, the single- or dual-pol parameters can only be expected to provide half the accuracy with signiﬁcant false-alarms in areas with temporal variations independent of landslides. as detection and false alarm rates for two landslide sites in USA and Russia. They used pre-event optical and post-event dual-pol SAR data for detecting landslide areas. The detection rates varied from 0.48 to 0.87 with false alarms ranging from 0.001 to 0.003. Since they used small image chips for the evaluation and ancillary slope information for reﬁning detection results, it was not appropriate to directly compare the accuracies obtained previous study and this study. However, it was possible to examine the performance of the dual-pol parameter used in [10] in comparison with the result obtained in this study. single-, dual-, and quad-pol modes observations for deducing information about landslide areas.


Introduction
Landslides are one of the most widespread natural disasters. Producing a landslide map or event inventory map is an essential task to understand the extent and magnitude of landslides. Remote sensing techniques can accelerate the production of landslide maps for large spatial scales [1]. Since the landslide triggering event often hampers the acquisition of optical remote sensing data, the microwave remote sensing techniques including Synthetic Aperture Radar (SAR) can be very useful tools for near real-time landslides detection. The application of SAR images related to landslides has been mainly studied for monitoring slow-moving landslides through interferometric SAR technique based on phase information of SAR images [2]. However, the use of SAR amplitude information for the landslide application is not widespread.
Most studies on the detection of event landslides using SAR amplitude information have been carried out using polarimetric SAR data to overcome difficulties in discriminating landslide areas from land cover classes. Czuchlewski et al. [3] employed L-band airborne polarimetric SAR data (AIRSAR) and examined the usability of polarimetric parameters for detecting a large slope failure triggered by an earthquake in Taiwan. The scattering characteristics of the landslide region were able to be identified by the eigenvalue-eigenvector decomposition [4]. After the launch of space-borne polarimetric SAR systems, several studies have been conducted to analyze the availability of satellite SAR data for

Study Area and Data
A series of earthquakes occurred in Kumamoto, Kyushu Island, Japan in 2016 [11]. The first major shock (foreshock) occurred at 21:26 Japanese Standard Time (JMT) on 14 April 2016 with the moment magnitude of Mw 6.2. Subsequently, on 16 April at 01:26 JMT, a second strong shock (mainshock) occurred nearby with the magnitude of Mw 7.0. The foreshock occurred at a depth of 11 km. The focal mechanism exhibited right-lateral strike-slip faulting with a NNW-SSE tension axis. The epicenter of the mainshock is approximately 4.5 km northwest of the epicenter of the foreshock. Figure 1 shows the epicenter of the foreshock and mainshock. The focal mechanism of the mainshock is a right-lateral strike-slip fault with a tension axis in the NW-SE direction. It triggered numerous landslides in and around Minamiaso, which is located in the western part of Aso caldera [12].
terrain, geometric distortions of SAR image should be corrected before further analysis of SAR signals. To represent the SAR images as geometrically similar to the map coordinates, all SAR images were converted to the multilook covariance matrix format [14]. Then, the range-Doppler orthorectification was performed by using SRTM (Shuttle Radar Topography Mission) global 1 arc second digital elevation model (DEM). In addition, polarimetric speckle filtering using the IDAN filter [15] was applied to all PALSAR-2 data to reduce speckle.

Radar Scattering Characteristics of Landslide Area
In this Section, we recall the single-, dual-, and quad-pol SAR observables and evaluate their scattering characteristics of manually selected landslide affected areas. Among several observation modes of the ALOS-2 PALSAR-2 sensor, a fully polarimetric quad-pol mode data observed on 21 April 2016, about 5 days after the mainshock, was obtained in this study. In addition, a pre-event quad-polarization PALSAR-2 data observed on 3 December 2015 was obtained with the same observation mode. To evaluate landslide-induced changes of SAR signals, a study area (white rectangle in Figure 1) was selected in which several landslides and slope failures occurred by the Kumamoto earthquakes.
The PALSAR-2 images of the study area acquired before and after the earthquake are shown in Figure 2. Six landslides sites (marked in black in Figure 2b) for the examination of radar scattering responses were manually selected by the literature [12] and the aerial photos [13]. Three of the selected landslides (S1, S2, and S3) occurred in forested areas, while the other three (S4, S5, and S6) occurred in non-forested grasslands and shrublands. Since a landslide generally occurs in sloping terrain, geometric distortions of SAR image should be corrected before further analysis of SAR signals. To represent the SAR images as geometrically similar to the map coordinates, all SAR images were converted to the multilook covariance matrix format [14]. Then, the range-Doppler orthorectification was performed by using SRTM (Shuttle Radar Topography Mission) global 1 arc second digital elevation model (DEM). In addition, polarimetric speckle filtering using the IDAN filter [15] was applied to all PALSAR-2 data to reduce speckle. Figure 1. Location of the foreshock and mainshock of the 2016 Kumamoto earthquake and ALOS-2 PALSAR-2 data coverage. Boxed area in the PALSAR coverage is the Minamiaso study site in which many landslides occurred during the earthquake.

Radar Scattering Characteristics of Landslide Area
In this Section, we recall the single-, dual-, and quad-pol SAR observables and evaluate their scattering characteristics of manually selected landslide affected areas.

Radar Scattering Characteristics of Landslide Area
In this Section, we recall the single-, dual-, and quad-pol SAR observables and evaluate their scattering characteristics of manually selected landslide affected areas.

Single-Pol Scattering Characteristics
Most of conventional SAR systems have been operated in the single polarization channel, which transmit and receive either horizontally (H) or vertically (V) polarized electromagnetic wave. In this study, we will only consider the horizontally polarized signal. The horizontally polarized scattered field E s H is related the incident field E i H by the complex scattering amplitude S HH as where k is wavenumber and r is the distance between the target and the receiving antenna. The scattering amplitude depends on the dielectric and geometric properties of the scatterer in the direction of incident and scattered fields. The scattering properties of natural objects can be generally described by the backscattering coefficient σ 0 HH = S HH S * HH = |S HH | 2 , which is the averaged scattering intensity of a distributed scatterer. Figure 3a,b show temporal σ 0 HH histograms acquired over the selected landslide sites that occurred in forested areas and non-forested grasslands and shrublands, respectively. It is seen that σ 0 HH generally increases after the landslides occurred in non-forested areas, while there is no clear difference between the two histograms in forested areas. It is probably attributed to changes in the microwave scattering mechanisms, particularly the surface and double-bounce scattering mechanisms, in mountainous forest areas [16].
In the case of the sloping terrain, topography effects should be considered in the analysis of scattering response. Particularly, it was shown in the previous study [16] that the variations of local incidence angle due to local topography can lead to changes in the scattering processes. The local incidence angle θ l is defined as the angle between the line of sight direction and the surface normal as cos θ l = tan α rg sin θ + cos θ where α rg and α az are local slope angles in range and azimuth directions respectively, and θ is the incidence angle defined in the flat geometry. better explained by the local incidence angle. The temporal changes show different angular trends, particularly in forested areas. The σ 0 HH value increases after the landslide at the low local incidence angle region while it shows no apparent changes at the high local incidence angle regions.
scattering amplitude depends on the dielectric and geometric properties of the scatterer in the direction of incident and scattered fields. The scattering properties of natural objects can be generally described by the backscattering coefficient = 〈 * 〉 = 〈| | 〉 , which is the averaged scattering intensity of a distributed scatterer. Figure 3a,b show temporal histograms acquired over the selected landslide sites that occurred in forested areas and non-forested grasslands and shrublands, respectively. It is seen that generally increases after the landslides occurred in non-forested areas, while there is no clear difference between the two histograms in forested areas. It is probably attributed to changes in the microwave scattering mechanisms, particularly the surface and double-bounce scattering mechanisms, in mountainous forest areas [16]. In the case of the sloping terrain, topography effects should be considered in the analysis of scattering response. Particularly, it was shown in the previous study [16] that the variations of local incidence angle due to local topography can lead to changes in the scattering processes. The local incidence angle is defined as the angle between the line of sight direction and the surface normal as

Dual-Pol Scattering Characteristics
According to the observation strategy [17], PALSAR-2 primarily acquires data in the dual-pol mode. It employs a single polarization transmission (H) and dual (H and V) coherent reception. The dual-pol mode provides two complex scattering amplitude, S HH and S VH , which relate the incident and scattered signals as Then, the dual-pol covariance matrix [C 2 ] can be obtained by applying multilook processing to the dual-pol complex scattering vector The diagonal terms of the dual-pol covariance matrix provides additional scattering intensity in VH-polarization basis, such as σ 0 VH = S VH S * VH = |S VH | 2 . The off-diagonal term is the correlation between HH (co-pol) and VH (cross-pol) scattering amplitudes. Since there will be usually no coand cross-pol correlation term in natural reflection symmetric media, two scattering intensities or backscattering coefficients are the two important observables in the dual-pol mode. Figure 4a,b show temporal σ 0 VH values of landslide areas that occurred in forested and non-forested areas. As compared with the σ 0 HH case, it shows slightly better separation between pre-and post-landslide conditions at the high incidence angle region, while it shows no apparent differences at the low local incidence angle regions. The two polarization channels exhibit somewhat complementary angular trends, the ratio between VH-and HH-polarization intensities, such as the cross-pol ratio σ 0 VH /σ 0 HH , can be used to further characterize landslide-related backscatter changes. Figure 4c,d show the cross-pol ratio of landslide affected areas in relation to the local incidence angle. It shows a much weak angular dependency on the local incidence angle variations and emphasizes the landslide-induced changes of radar scattering properties in forested areas at either low or high local incidence angle region.
The scattering properties of natural media can be described by the quad-pol covariance matrix. In the monostatic backscattering case, when = by the scattering reciprocity, the quad-pol covariance matrix can be formed by the complex scattering vector = √2 , such as The Hermitian matrix is fully characterized by nine real parameters: three polarimetric intensities in the diagonal elements and three complex polarimetric correlations in the off-diagonal elements. As aforementioned, the co-and cross-pol correlation terms, 〈 * 〉 and 〈 * 〉, will be close to zero for the natural scatterers having reflection symmetry. On the other hand, the correlation coefficient between HH-and VV-polarization, such as 〈 * 〉/ 〈 * 〉〈 * 〉 , contains important information on the scatterers. The magnitude of the correlation, named co-pol coherence , indicates signal depolarization and varies from 0 for completely random signal to 1 for pure single scattering [18]. Since a landslide generally results in changes of land cover types from vegetation to bare surface, identification of signal depolarization properties will be important

Quad-Pol Scattering Characteristics
In general case, both the incidence wave and the wave scattered by an object can be adequately described by the two-dimensional orthogonal bases. The incidence wave of the quad-pol mode observation has both horizontal and vertical polarization component. Then, a scattering object can be characterized by the 2 × 2 complex scattering matrix [S], such as The scattering properties of natural media can be described by the quad-pol covariance matrix. In the monostatic backscattering case, when S HV = S VH by the scattering reciprocity, the quad-pol The Hermitian [C 3 ] matrix is fully characterized by nine real parameters: three polarimetric intensities in the diagonal elements and three complex polarimetric correlations in the off-diagonal elements. As aforementioned, the co-and cross-pol correlation terms, S HH S * HV and S HV S * VV , will be close to zero for the natural scatterers having reflection symmetry. On the other hand, the correlation coefficient between HH-and VV-polarization, such as S HH S * VV / S HH S * HH S VV S * VV , contains important information on the scatterers. The magnitude of the correlation, named co-pol coherence ρ HHVV , indicates signal depolarization and varies from 0 for completely random signal to 1 for pure single scattering [18]. Since a landslide generally results in changes of land cover types from vegetation to bare surface, identification of signal depolarization properties will be important in mapping landslide affected areas.
In addition to the polarimetric intensities and correlation terms in the [C 3 ] matrix, the general depolarization properties of a scatterer can be characterized by using the eigenvalue analysis [18]. By definition, [C 3 ] is a Hermitian matrix and can be decomposed into three real non-negative eigenvalues and three orthogonal unitary eigenvectors as: where λ i and → e i are eigenvalues and eigenvectors of [C 3 ], respectively. The eigenvalue decomposition expands [C 3 ] into the sum of three uncorrelated target matrices [C 3i ] where the eigenvalues λ i (λ 1 ≥ λ 2 ≥ λ 3 ≥ 0) are the statistical weights for the three components. The amount of polarizing contribution corresponds to the largest eigenvalue while the other eigenvalues are the depolarization components in the [C 3 ] matrix. Consequently, the relative contribution of the eigenvalues can be a good indicator of system depolarization properties. In order to interpret the eigenvalue contributions, we can define pseudo probabilities P i as The spread of probabilities can be represented by the single scaler measure, such as polarimetric entropy H defined as [19]: It indicates the scattering randomness which ranges from 0 (single nonzero eigenvalue) for the single scatter to 1 (three equal eigenvalues) for the random scatterer.
For a highly depolarizing media, the entropy is not a unique function of the eigenvalue ratios, and another eigenvalue parameter defined as the polarimetric anisotropy A defined as Equation (10) can be useful to describe two minor eigenvalues [20].
It indicates relative importance of the second and the third eigenvalues which also ranges from 0 to 1. High entropy and low anisotropy values imply strongly depolarizing systems with the presence of several scattering mechanisms.
In addition to the entropy and anisotropy parameters, each eigenvalue contribution in Equation (8) can be directly used as a measure of depolarization in a random media. Among them, the smallest eigenvalue P 3 has been used to infer the amount of vegetation. It is called Radar Vegetation Index (RVI) defined as [21] RVI = 4P 3 .
The factor 4 was introduced so that the RVI for a cloud of randomly oriented thin cylinders would be equal to 1. On the other hand, the amount of polarizing contribution among backscattered signals, which is represented by the largest eigenvalue P 1 [18], can be a useful parameter for identifying changes of signal depolarization levels caused by the landslide. P 1 has a maximum value of 1 when λ 2 = λ 3 = 0 and a minimum value of 1/3 when λ 1 = λ 2 = λ 3 . Therefore, as in the case of RVI, the normalized value of P 1 , named P pol , can be defined as: It indicates the polarizing contribution of a partially polarized system which ranges from 0 for the random noise process to 1 for the pure polarizing target. Figures 5 and 6 show the temporal changes of various polarimetric parameters of landslide areas that occurred in forested and non-forested areas, respectively. Among the diagonal components of the [C 3 ] matrix, only VV-polarization scattering intensity σ 0 VV = S VV S * VV is illustrated in Figures 5a  and 6a, because the HH-and HV-polarization scattering intensities were discussed in Figures 3 and 4. Similar to the single or dual-pol cases, it is difficult to distinguish signals from pre-and post-landside conditions in temporal σ 0 VV values, particularly in the forested areas.    It is seen that the polarimetric parameters defined by scattering amplitude and eigenvalue ratios in Figure 5b-f show better discriminability between pre-and post-landslide conditions of the selected sites. Among them, the polarimetric parameters indicating signal depolarization properties, particularly , entropy, and , exhibit significant landslide-induced changes for all local It is seen that the polarimetric parameters defined by scattering amplitude and eigenvalue ratios in Figure 5b-f show better discriminability between pre-and post-landslide conditions of the selected sites. Among them, the polarimetric parameters indicating signal depolarization properties, particularly ρ HHVV , entropy, and P pol , exhibit significant landslide-induced changes for all local incidence angle regions. The anisotropy and RVI parameters, which utilize the minimum eigenvalue, cannot clearly distinguish signals from pre-and post-landslide conditions, especially in non-forested areas as shown in Figure 6.

Landslide Detection
One of the goals of remote sensing data analysis for a landslide event is to generate the landslide map or landslide inventory map. A simple way to map landslide affected areas is to compare remote sensing data acquired at pre-and post-landslide conditions by image differencing. Consider two M × N polarimetric parameters, X 1 and X 2 , acquired at before and after the landslide, respectively. Figure 7 illustrates the difference image, ∆X = X 2 − X 1 , for the different polarimetric parameters. Figure 7a-d corresponds to the single-and dual-pol parameters, and Figure 7e-f corresponds to the quad-pol parameters. It also shows a manually generated reference image which illustrates regions of interests for landslide damaged areas and undamaged forest and agricultural areas. It is seen that changes in polarimetric parameters highlight the damaged areas except for the cross-pol ratio and the anisotropy. Three intensities σ 0 HH , σ 0 HV , and σ 0 VV show landslide-induced changes, but the direction of changes vary depending on the slope direction. The eigenvalue parameters and the co-pol coherence exhibit one-directional changes for landslide areas from high depolarization states at the pre-landslide condition to low depolarization states at the post-landslide condition. However, the amount of changes in the landslide areas and the level of general changes unrelated to landslides vary across the polarimetric parameters.
where the conditional probability density function of i th class is described by the mean and the standard deviation . Once the final parameter estimates of the three probability distributions are obtained, the optimal threshold values and can be derived according to the maximum a posterior decision rule, that is, solving following equations with respect to the variable ∆ : Since the EM-based threshold is applied exclusively to the pixel values of the difference image, the change detection usually provides noisy classification results. To further improve the detection After generating the difference image, the landslide map can be generated by a change detection approach. It aims to assign each pixel of the difference image to either changed class ω c or unchanged background class ω n . One of the simple and unsupervised solutions to this problem, which is important in emergency observation, is to select the global threshold in the difference image. In order to better reflect the characteristic changes of polarimetric parameters occurred in the landslide areas and other undamaged areas, change detection was carried out considering the negatively changed (ω c1 ) and the positively changed (ω c2 ) classes in this study. To find multiple thresholds from the difference image, we adopted the Expectation-Maximization (EM) thresholding method [22][23][24][25].
We assume that the probability distribution of the difference image p(∆X) is a mixture of three density components associated with the negatively and positively changed classes and unchanged background class, that is, p(∆X) = P(ω c1 )p(∆X|ω c1 ) + P(ω n )p(∆X|ω n ) + P(ω c2 )p(∆X|ω c2 ), where P(ω c1 ), P(ω c2 ), and P(ω n ) are the prior probabilities of negatively changed, positively changed, and unchanged classes, respectively, and p(∆X|ω c1 ), p(∆X|ω c2 ), and p(∆X|ω n ) are their conditional probability density functions. Both the unknown prior probabilities and the conditional probability density functions can be calculated iteratively by the EM procedure until convergence. If we further assume the Gaussian density functions, the prior and conditional probabilities of i th class where ω i ∈ {ω c1 , ω n , ω c2 } at (t + 1) iteration is given by: where the conditional probability density function of i th class is described by the mean µ i and the standard deviation σ i . Once the final parameter estimates of the three probability distributions are obtained, the optimal threshold values T 1 and T 2 can be derived according to the maximum a posterior decision rule, that is, solving following equations with respect to the variable ∆X: P(ω n )p(∆X|ω n ) = P(ω c1 )p(∆X|ω c1 ) and P(ω n )p(∆X|ω n ) = P(ω c2 )p(∆X|ω c2 ). (17) Since the EM-based threshold is applied exclusively to the pixel values of the difference image, the change detection usually provides noisy classification results. To further improve the detection results, spatial contextual information can be considered after obtaining three probability distributions and preliminary decisions from two threshold values. The Markov Random Field (MRF) model can provide a useful tool for characterizing contextual information [23][24][25][26]. The final change detection map after obtaining class labels for each pixel C l (m, n) ∈ {ω c1 , ω n , ω c2 } can be generated by minimizing the energy function U(∆X, C l ) such as: [U data (∆X(m, n), ω(m, n)) + U context (ω(m, n))] (18) where U data is the class conditional energy function. Under the Gaussian assumption, it can be written as: The contextual energy function U context (ω) describes spatial contextual information in a local spatial neighborhood system N given as: It counts the number of pixels in the neighborhood system assigned to the same class as pixel (m, n) with delta function and the β parameter tunes the importance of spatial context. In this study, we used the second-order neighborhood system containing eight neighboring pixels with β = 1.6.
After obtaining the class decision among the negatively and positively changed classes and unchanged background class for each pixel, the final landslide map can be generated accordingly with the characteristic changes of polarimetric parameters. Each pixel x(m, n) of the detection map is assigned to either 1 or 0, where 1 indicates landslide damaged areas and 0 indicates undamaged areas according to: Appropriate types of change ω ck for different polarimetric parameters (Table 1) were selected based on the analysis on the polarimetric scattering characteristics of landslide areas discussed in the previous Section. Table 1.
Selected changed class ω ck associated with the landslide area for different polarimetric parameters. Figure 8 illustrates the final binary decision results for landslide-related changes derived from the selected polarimetric parameters. Comparing the detection results with the reference landslide map in Figure 8i, it is seen that the landslide affected areas can be successfully identified by the automatic change detection from quad-polarization SAR observations. In order to evaluate the detection result quantitatively, several accuracy metrics were considered including the detection rate (Pd), false-alarm rate (Pfa), overall accuracy (OA), and Kappa coefficient (Kappa) [27]. Table 2 shows the accuracy analysis results for each binary classification result. The entries highlighted in bold correspond to the best and second-best polarimetric parameters of a specific accuracy metric. Detection rates of around 60% can be obtained using the co-pol coherence ∆ρ HHVV , the entropy ∆H, and the polarizing contribution ∆P pol . Changes of single-pol scattering intensities can also provide high detection rates but suffer from significant false alarms at the same time. The OA and Kappa provide an overall idea of detectability. Both accuracy metrics show that ∆ρ HHVV and ∆P pol are two best parameters for landslide detection. The OA shows much higher values than other accuracy metrics for all polarimetric parameters because of the low proportion of landslide areas in all pixels selected for accuracy analysis. In this case, the overall detection performance of different polarimetric parameters can be better evaluated by the Kappa. According to Kappa, we can expect about two times better performance by using quad-pol parameters with the significantly less false alarm errors than single-or dual-pol parameters.   Table 2 shows the accuracy analysis results for each binary classification result. The entries highlighted in bold correspond to the best and second-best polarimetric parameters of a specific accuracy metric. Detection rates of around 60% can be obtained using the co-pol coherence ∆ , the entropy ∆ , and the polarizing contribution ∆ . Changes of singlepol scattering intensities can also provide high detection rates but suffer from significant false alarms at the same time. The OA and Kappa provide an overall idea of detectability. Both accuracy metrics show that ∆ and ∆ are two best parameters for landslide detection. The OA shows much higher values than other accuracy metrics for all polarimetric parameters because of the low proportion of landslide areas in all pixels selected for accuracy analysis. In this case, the overall detection performance of different polarimetric parameters can be better evaluated by the Kappa. According to Kappa, we can expect about two times better performance by using quad-pol parameters with the significantly less false alarm errors than single-or dual-pol parameters.

Comparison with Previous Study
Although several papers have discussed landslide detection problems from single-, dual-, and quad-pol SAR data, there have been few studies that have provided accuracy analysis results of landslide detection. Plank et al. [10] analyzed the overall accuracy as well as detection and false alarm rates for two landslide sites in USA and Russia. They used pre-event optical and post-event dual-pol SAR data for detecting landslide areas. The detection rates varied from 0.48 to 0.87 with false alarms ranging from 0.001 to 0.003. Since they used small image chips for the evaluation and ancillary slope information for refining detection results, it was not appropriate to directly compare the accuracies obtained previous study and this study. However, it was possible to examine the performance of the dual-pol parameter used in [10] in comparison with the result obtained in this study.
The key parameter for the landslide identification in [10] was the dual-pol entropy H 2 . It can be derived from the two eigenvalues λ 1 and λ 2 of the dual-pol covariance matrix [C 2 ].
To evaluate the detectability of the H 2 parameter, the EM-MRF based landslide detection method was applied to the temporal change ∆H 2 shown in Figure 9a. In comparing ∆H 2 and the quad-pol entropy ∆H in Figure 8f, ∆H 2 provides more noisy changes in undamaged areas resulting in noisy detection result as shown in Figure 9b. of landslide detection. Plank et al. [10] analyzed the overall accuracy as well as detection and false alarm rates for two landslide sites in USA and Russia. They used pre-event optical and postevent dual-pol SAR data for detecting landslide areas. The detection rates varied from 0.48 to 0.87 with false alarms ranging from 0.001 to 0.003. Since they used small image chips for the evaluation and ancillary slope information for refining detection results, it was not appropriate to directly compare the accuracies obtained previous study and this study. However, it was possible to examine the performance of the dual-pol parameter used in [10] in comparison with the result obtained in this study.
The key parameter for the landslide identification in [10] was the dual-pol entropy . It can be derived from the two eigenvalues and of the dual-pol covariance matrix .
To evaluate the detectability of the parameter, the EM-MRF based landslide detection method was applied to the temporal change ∆ shown in Figure 9a. In comparing ∆ and the quad-pol entropy ∆ in Figure 8f, ∆ provides more noisy changes in undamaged areas resulting in noisy detection result as shown in Figure 9b.  Table 3 summarizes accuracy parameters of ∆ as compared with those of ∆ . The dual-pol entropy offers lower detectability with a significant amount of false alarms. According to the Kappa coefficient, the overall detection performance of the dual-pol entropy is only about half of the quad-  Table 3 summarizes accuracy parameters of ∆H 2 as compared with those of ∆H. The dual-pol entropy offers lower detectability with a significant amount of false alarms. According to the Kappa coefficient, the overall detection performance of the dual-pol entropy is only about half of the quad-pol entropy. Therefore, it is important to use quad-pol SAR data in the generation of the landslide map.

Influence of Slope and Land Cover on the Detectability
As discussed in Sections 3 and 4, the quad-pol parameters indicating depolarization properties exhibited clear landslide-induced temporal changes for all incidence angle regions, and provided much higher detectability than single-or dual-pol parameters. Nonetheless, there were some differences in the detection performance between signal depolarization parameters such as ρ HHVV , H, and P pol . In order to better understand the detection performance of depolarization parameters, the influence of environmental conditions, such as land cover types and slopes, on the detectability was further examined. Similar to Section 3, the regions of interest were divided into two categories, such as forested areas and non-forested grasslands and shrublands at the pre-landslide condition. The pixels belonging to each category were further subdivided based on the local incidence angle according to the local slope. Then, we recalculated the detection rates for the subdivided regions of interest. Figure 10 illustrates the effects of environmental conditions on the detectability of two selected depolarization parameters, ∆ρ HHVV and ∆H, that differed in detection accuracy. In Figure 10, the detection rates, which were recalculated for two different land cover types, were plotted with respect to the local incidence angle. It is seen that the detectability of polarimetric parameters can be affected by the local slope. The detection rates generally decrease with an increase of the local incidence angle. In the case of ∆H, which offers a slightly lower detection rate for the entire study area among depolarization parameters, the angular dependency of the detection rate is more evident regardless of the land cover condition. The detection rate of ∆ρ HHVV , however, is much less sensitive to the local incidence angle variation particularly in the non-forested area resulting in the higher overall detectability.
interest. Figure 10 illustrates the effects of environmental conditions on the detectability of two selected depolarization parameters, ∆ and ∆ , that differed in detection accuracy. In Figure 10, the detection rates, which were recalculated for two different land cover types, were plotted with respect to the local incidence angle. It is seen that the detectability of polarimetric parameters can be affected by the local slope. The detection rates generally decrease with an increase of the local incidence angle. In the case of ∆ , which offers a slightly lower detection rate for the entire study area among depolarization parameters, the angular dependency of the detection rate is more evident regardless of the land cover condition. The detection rate of ∆ , however, is much less sensitive to the local incidence angle variation particularly in the non-forested area resulting in the higher overall detectability. Figure 10. Variations of the detection rates according to the local incidence angle and the land cover type for the detection results from (left) ∆ and (right) ∆ .

Conclusions
Due to various triggering mechanisms, landslides are widespread in many parts of the world. With a significant triggering event, slope failures can be sparse and widely distributed across a large area. Remote sensing technologies, particularly with air-borne or space-borne SAR sensors, can

Conclusions
Due to various triggering mechanisms, landslides are widespread in many parts of the world. With a significant triggering event, slope failures can be sparse and widely distributed across a large area. Remote sensing technologies, particularly with air-borne or space-borne SAR sensors, can provide timely and critical information to assess landslide distributions over large areas. Nonetheless, studies on the systematic use of SAR data to produce event landslide maps are rare. In this context, this study discussed an automatic detection and mapping strategy of landslide areas from polarimetric SAR data. This study focused mainly on the usability of different polarimetric parameters of single-, dual-, and quad-pol modes observations for deducing information about landslide areas.
The experimental analysis on the characteristics of landslide-induced changes for different polarimetric parameters showed that, in general, SAR observations were largely influenced by the local slopes in which the landslide occurred. This is because the microwave scattering mechanisms can vary not only by landslide-induced land cover changes but also local incidence angle variations in the sloping terrain as reported in the previous study [25]. Consequently, single-or dual-pol SAR observations were not able to represent various landslide areas that occurred in different topographic conditions. However, quad-pol parameters indicating signal depolarization properties revealed noticeable landslide-induced temporal changes for all local incidence angle ranges.
To produce a landslide map, a simple change detection method based on EM thresholding and MRF contextual classification was also presented in this study. For an automatic identification of landslide areas, binary decision strategies based on the characteristic scattering properties of landslide areas were proposed. The accuracy assessment results showed that the depolarization parameters, such as the co-pol coherence and polarizing contribution, can identify areas affected by landslides in about 18.9 km 2 study area with a detection rate of 60%, a false-alarm rate of 5%, and a Kappa coefficient of 0.45. On the other hand, the single-or dual-pol parameters can only be expected to provide half the accuracy with significant false-alarms in areas with temporal variations independent of landslides. It is worth noting that the results were obtained exclusively from the polarimetric SAR data. The detection algorithm is unsupervised and fully automatic. Therefore, the proposed detection accuracy can be considered to be the minimum values that can be expected from the quad-pol scattering mechanism indicators if a rapid survey is required. There is enough room to improve accuracy through advanced change detection or supervised classification methods. In particular, considering the effects of environmental conditions on the SAR observations, it is considered that the detection accuracy can be improved by using ancillary information such as DEM and optical remote sensing data.
Author Contributions: Conceptualization, Sang-Eun Park and Sun-Gu Lee; methodology, software, validation, formal analysis and writing-Original draft preparation, Sang-Eun Park; writing-review & editing, Sang-Eun Park and Sun-Gu Lee; project administration, Sun-Gu Lee.