Landslide, a significant natural hazard in the high mountain area, has drawn worldwide attention due to increasing awareness of its socioeconomic impacts, as well as, the increasing pressure of urbanization on the mountain environment  . Landslide study establishes a significant constraint to development and urban planning. Landslide susceptibility mapping can be the preliminary step in mitigating the damages  . Moreover, landslide susceptibility assessment is an important process for prediction and management of natural disasters. It is also a necessary step for integrated watershed management, hazard mitigation, and urban planning in government policies   . A landslide susceptibility or hazard map depicts areas likely to have landslides in the future by correlating some of the principal factors that contribute to landslides with the past distribution of slope failures. Recently, DMR has been assigned to study about landslide hazard zonation in Na Haew District, Loei Province, NE Thailand.
2. Study Area
The Na Haew District is located between 17˚16'43'' and 17˚35'43''N latitudes and 100˚49'50'' and 101˚107'36''E longitudes. It covers a total area of approximately 625 km2. The study area is represented by the rocks of Phu Khat and Khao Ya Puk Formations.
3. Data and Method
In current study, a landslide inventory map was prepared using field surveys, local information, and satellite interpretation obtained from Google Earth software. The landslide scars found in Na Haew District were shown in Figure 1(a).
In this study, nine independent factors were chosen as potentially contributing to landslide susceptibility, slope, elevation, aspect, lithology, distance from drainage, distance from lineaments, NDVI, SPI and land-use. The availability of thematic data varies widely, depending on the type, scale, and method of data acquisition.
For landslide susceptibility mapping within Na Haew we choose statistical approach and frequency ratio methods. The landslide inventory map was compared with the various factor maps using a selected statistical method. By using statistical analysis, the factor and/or combination of factors that have resulted in slope instability in the past are determined. Quantitative predictions can then be made for areas where no landslides are currently present, but which have similar
Figure 1. Landslide inventory map (a) and landslide susceptibility map (b) of Na Haew District.
conditions to those where landslides exist. In frequency ratio ananysis the Frequency Ratio value (Fr) for each factor’s range was calculated by dividing the landslide occurrence ratio by the area ratio. Then, the frequency ratios were summed to calculate the landslide susceptibility index (LSI) as shown in Equation (1).
Once a landslide susceptibility map is created, it is necessary to divide this map into different susceptibility classes. The simplest method is to assign the categories using expert opinion as done by  . The most common method for this purpose depends on the optimum bandwidth classification of the histograms of various parameters  . However, this is not straight forward as there are no statistical rules which can categorize continuous data automatically  . Recently, mathematical methods for data classification have become available in GIS software. In literature, there are many methods used    . These are based on manual or natural breaks, equal intervals, or statistical consideration, and can be described as follows. For the susceptibility classification for this study, the natural break method was used, and the susceptibility classes were divided into 5 classes (very low, low, moderate, high, very high) as shown in Figure 1(b).
This work is supported by DMR of Thailand and is a contribution to UNESCO/ IUGS/IGCP 679 project.