Spatial implementation of frequency ratio, statistical index and index of entropy models for landslide susceptibility mapping in Al-Balouta river basin, Tartous Governorate, Syria

Landslide vulnerability prediction maps are among the most important tools for managing natural hazards associated with slope stability in river basins that affect ecosystems, properties, infrastructure and society. Landslide events are among the most hazardous patterns of slope instability in the coastal mountains of Syria. Thus, the main goals of this research are to evaluate the performance of three different statistical outputs: Frequency Ratio (FR), Statistical Index (SI) and Index of Entropy (IoE) and therefore map landslide susceptibility in the coastal region of Syria. To this end, we identified a total of 446 locations of landslide events, based on the preliminary inventory map derived from fieldwork and high-resolution imagery surveys. In this regard, 13 geo-environmental factors that have a high influence on landslides were selected for landslide susceptibility mapping. The results indicated that the FR method outperformed the SI and IoE models with a high AUC of 0.824 and better adaptability, followed by the SI with 0.791. According to the SCAI values, although the FR model achieved the best reliability, the other two models also showed good capability in determining landslide susceptibility. The result of FR-based modelling showed that 18.51 and 19.98% of the study area fall under the high and very high landslide susceptible categories, respectively. In the map generated by the SI method, about 36% of the study area is classified as having high or very high landslide sensitivity. In the IoE method, whereas 14.18 and 25.62% of the study area were classified as “very high susceptible” and “high susceptible,” respectively. The relative importance analysis demonstrated that the slope aspects, lithology and proximity to roads effectively motivated the acceleration of slope material instability and were the most influential in both the FR and SI models. On the other hand, the IoE model indicated that the proximity to faults and roads, along with the lithology factor, were important influences in the formation of landslide events. As a result, the statistical bivariate models-based landslide mapping provided a reliable and systematic approach to guide the long-term strategic planning procedures in the study area.


Introduction
Landslides are classified as one of the most important physical hazards affecting human life, infrastructure and sustainable development (Alsabhan et al. 2022;Nwazelibe et al. 2022). Landslide event occurs when the shear strength of the material, that forms the slope, is greater than gravity and other types of shear stress within a path (Reichenbach et al. 2018a, b;Shi et al. 2020;Karaman et al. 2022). Unlike other natural hazards, the landslide risk assessment process is described as complex due to the difficulty of landslide inventorization and the spatial interaction between the motivating factors (Naceur et al. 2022;Jana et al. 2019;Valiante et al. 2021;Hateffard et al. 2021). Globally, landslides cause a real threat to life by total or partial destruction of infrastructure projects (Emberson et al. 2021;Jamir et al. 2022).
Landslides are among the greatest significant issues that governments strive to reduce their destructive impacts on lives, properties and infrastructure, especially in mountainous areas (Dikshit et al. 2020;Nakileza and Nedala 2020;Batar and Watanabe 2021;Rahman et al. 2022). Many studies confirmed that increasing the landslide events are motivated by the spatial integration between physical and anthropological factors (Razavizadeh et al. 2017;Akinci and Yavuz Ozalp 2021;Yamusa et al. 2022;Jaafari et al. 2022). The first fundamental step related to efficient spatial management of landslide risk is the preparation of landslide susceptibility mapping (Das et al. 2022;Guo et al. 2021;Skrzypczak et al. 2021). A landslide susceptibility map, however, delineates zones vulnerable to future landslide events within a given area. the high-quality generation of landslide susceptibility map depends on the effectiveness of the fieldwork, access to the most suitable data, the overall determination of the spatial conditioning factors and the appropriate selection of modelling and simulation methods (Mersha and Meten 2020;Wubalem 2021).
As such, landslide susceptibility models can be classified into three categories: machine learning-based models, empirical approaches-based models and statistical-driven models. However, each of the aforementioned landslide susceptibility models has its own set of benefits and disadvantages (Reichenbach et al. 2018a, b). The machine learning-based landslide susceptibility models offer superior flexibility and adaptability, but they are limited by the model parameters used and the quantity of the training dataset (Zhou et al. 2018). In spite of the fact that the empirical models make use of past information and experience, the analysis results may differ significantly from the natural conditions (Ghosh et al. 2011). Whilst statistical-based landslide susceptibility models can reflect the relation between input conditional factors and output assessment outcomes, the linear model is very simple and subject to aberrations (Reichenbach et al. 2018a, b). As a result, prior studies demonstrate that the efficacy of established susceptibility models differs depending on various conditioning factors, and no one method is preferable in all settings (Zhou et al. 2018;Argyriou et al. 2022).
According to the current literature, the Mediterranean terrain is one of the areas most affected by the landslides risk (Argyriou et al. 2022;. The landslide is considered one of the most direct impacts on slope instability in Mediterranean environments (Ullah et al. 2022). Hence, the prediction of landslides is one of the most important pivots of geological and geomorphological studies in these environments, such as Isparta−Antalya highway (D-685), Turkey (Hepdeniz 2020), Ionian Islands, Greece (Mavroulis et al. 2022), Mila town, Algeria (Bounemeur et al. 2022) and the prefecture of Chania, Crete (Psomiadis et al. 2020).
Due to fragile physical characteristics and the acceleration of human activities, the coastal mountain area (CMA) is considered one of the most vulnerable areas to geomorphological hazards in Syria (Alsafadi et al. 2022;Abdo 2018;Mohammed et al. 2020aMohammed et al. , b, 2021. A literature review of geomorphological hazards in the eastern Mediterranean revealed an almost complete absence of studies related to landslide susceptibility in CMA which is prone to annual landslide events. For example, landslides in the winter of 2019 caused severe consequences regarding casualties and the partial destruction of the infrastructure in the Tartous Governorate. In this regard, the importance of this study can be justified by the urgent need to conduct more landslide assessment studies in a highly vulnerable area such as CMA. In addition, the abundance of related spatial data is considered a significant challenge in light of the consequences of the current war in Syria and spatial data availability (Chaaban et al. 2022). Thus, bivariate statistical methods (BSM) provide a reliable assessment with constructive results, especially in areas with scarce data. Moreover, the relevant literature has demonstrated the flexibility and performance quality of BSM, especially in the Mediterranean mountain environment (Karim et al. 2019;Karaman et al. 2022;Akter and Javed 2022).
The main goals of the present assessment are: (1) to digitize the current landslide events; (2) to map landslide conditioning factors; (3) to map landslide susceptibility using frequency (FR), Statistical Index (SI) and index of entropy (IOE); and (4) to assess the accuracy performance and outputs. The major importance of this analysis is to conduct an accurate landslide susceptibility assessment for Al-Balouta river basin through the combination and comparison of those models. This study provides spatial insights of high importance to national planners in terms of landslide hazard managing landslide risk in the study area, especially during the current war period in Syria, which caused a great gap in relevant studies.

Study area
The Al-Balouta river basin is sited in northwestern Syria, between Latitude 34° 57ʹ-35° 04ʹ North and longitudes 36° 01ʹ-36° 17ʹ East with an area of 116 km 2 (Fig. 1). This basin boarded by Al-Khawabi river basin to the west, Al-Ghab river basin to the west, Qays River basin to the south and Marqya River basin to the north. Geologically, the lithological structure of the study area varies from the Jurassic to the Quaternary (Ponikarov et al. 1967). Jurassic and Cretaceous formations consist of limestone, dolomite, marls, ophiolites, limy marl and sandy limestone. Neogene structure consists of basalt and sedimentary stones. Quaternary consists of Pleistocene and Holocene formations with fluviatile gravels, boulders, and deposits in riverbeds. The study area can be categorized into two major geomorphological regions based on the phase of terrestrial development (Abdo 2020). The hills area (110-400 m) consists of mainly upland parts distinguished by relatively sloping valleys. Dissected mountains area (400-1133 m) is featured faulted walls, steep slopes, narrow valleys and a diversity of geomorphological processes, especially karstification. The study basin is mainly subjected to the Mediterranean climate pattern: mild and rainy winter and long, dry and hot summer. The average summer temperature is 23.6° while in the winter is 10.3°. The annual precipitation is between 1000 and 1300 mm, with maximum precipitation recorded of about 320 mm in January. The study area is located in a wet climate Csa according to the Köppen−Geiger climate classification (Beck et al. 2018;Mohammed et al. 2020a).

Data
The assessment of landslide susceptibility was carried out in the study area based on a set of data from various sources (Fig. 2). Digital elevation model (DEM) obtained from shuttle radar topographic mission (SRTM) was used for mapping the topographical factors i.e., slope angle, slope aspect, curvature, plan curvature, profile curvature, elevation, streams and TWI. Fault and lithological structure data were acquired from the Ministry of Oil and Mineral Resources, Geology Directorate−Lattakia. Vegetation data was extracted from Landsat-8 OLI sensor satellite data collected from the United States Geological Survey (USGC). Rainfall data was collected from the Directorate of Meteorology, Tartous governorate. Road network data was obtained from the Ministry of Transport and Communications, Directorate of Transport and Public Roads. However, the details of the data sources have illustrated in Table 1. Figure 2 illustrated the methodology applied in this study.

Landslide inventory map
The landslide inventory map generally represents current and historical landslide events of an assigned area using various data sources, including GPS-based fieldwork, previous reports, interviews with locals and satellite image interpretation (Eitvandi et al. 2022;Abu El-Magd et al. 2021;Ali et al. 2021Ali et al. , 2020Abu El-Magd et al. 2021). The literature review reveals that a landslide inventory map can be generated either by collecting past evidence of different landslide patterns or using high-resolution satellite data combined with field surveys or by recognizing specific sites of landslide events in Google Earth by high-resolution satellite imageries. Conventionally, the inventory map was created using direct field investigation to confirm the actual locations. However, at present with the development of spatial technology, this task turns out to be easier to obtain accurate and rapid results (Chen et al. 2017a;Ali et al. 2020;Sachdeva et al. 2020;Pham et al. 2021b). It is frequently suggested that equivalent numbers of non-landslide locations should also be chosen for preparation and validation (Chen et al. 2017b;Abu El-Magd et al. 2021;Ali et al. 2021). In the current study, the inventory map was prepared according to the following strategy: (1) create a preliminary inventory map derived from fieldwork, information from local authorities and interviews with locals, (2) use high-resolution images (interpretation of Google Earth images (Google Earth Pro tools) to complete the digitizing of landslide events, (3) verify the final inventory map quality by conducting extensive fieldwork in various parts of the basin. This strategy was relied upon due to the lack of national inventory maps that could assist in conducting landslide assessment studies. Thus, this strategy can be used in other areas of CMA. However, the Geostatistical Analyst tools in ArcGIS (Geostatistical Analyst-Subset features) were used to divide the landslide events into a training dataset and a test dataset (Zhu et al. 2021;He et al. 2021;Pham et al. 2021a, b). A total of 446 locations of landslide events, however, are determined in study area. In this context, 70% (312 events) of landslides were selected randomly for the training dataset, whereas the rest of 30% (134 events) data was utilized for model validation (Fig. 1). The fieldwork revealed the diversity of patterns and types of landslides as a result of the complex integration between the driving geographical characteristics in the study area (

Landslide susceptibility index Frequency ratio (FR)
The frequency ratio is one of the well-known and widely used BSM which have been frequently applied for hazard susceptibility mappings, such as flood and landslide (Ali et al. 2020;Kincal and Kayhan 2022). FR is a widely utilized method for understanding the potential linking between current landslide events and causative geofactors (Rahman et al. 2022). However, FR can be defined through a spatial relationship between dependent and independent variables, where landslide inventory is the dependent variable and landslide conditioning factors are the independent variable. The weight of frequency ratio was estimated by dividing the pixels containing landslide points in each conditioning factors class and the total pixels of the considered area. The final weight in the frequency ratio for each conditioning factor is estimated using the following equation (Eq. 1): where Xpixel Li is No. of pixels containing landslide points in class X , Ypixel Ti is the No. of total pixels covering in class X over the study area, Xpixel Li is the sum of pixels containing landslide points in class X and Ypixel Ti is the sum of pixels covering in class X over the study area.
The value of frequency ratio >1 indicates there is a positive relationship between training points and each class of landslide conditioning factor and high landslide susceptibility, whereas a value <1 directs negative relation and low landslide susceptibility (Mind'je et al. 2020).

Statistical index (SI)
A statistical index (SI) is also BSM commonly used for landslide susceptibility assessment (Thapa and Esaki 1970). Using the SI, the particular class of a conditioning factor can be weighted based on the pixel concentration of landslide points of specific criteria and the total pixel concentration of landslide points across the whole study area (Bourenane et al. 2021). The weight of landslide factors computed using SI can be expressed as Eq. (2) where ld ij is the landslide density for i class of j factor, td is the total landslide density of the whole study region, L ij is the number of landslides in i class of j factor, L tn is the total number of landslides in the whole study region, pix ij is the number of pixels in i class of j factor, and pix t is the total pixels of the whole study region. The positive and negative value calculated using SI shows the presence and absence of a link between each class of landslide causative factors and landslide current events of landslide, respectively (Razavizadeh et al. 2017).

Index of entropy (IoE)
Index of entropy (IOE) is another BSM used for landslide susceptibility mapping, where determined on the basis of certain variables which calculate the weight of each variable. In this method, entropy depicts the level of uncertainty, imbalance, disorder and instability (Youssef et al. 2015). The role of different conditioning factors on the occurrence of landslides is represented by the entropy which provides an index system (Sahana et al. 2020) So, the value of entropy is useful for calculating the factor's weight . The equations which are used for calculating the information coefficient and W f expressing the weight of factors as a whole are as follows (Eqs. 3-8) (4) (P rs ) = P rs L s s=1 P rs where p and q are the area and percentage of landslide, respectively, P rs is the density of possibility, E s and E smax are the entropy values, I s and W f are information coefficient and factor's weight as a whole, respectively.

Landslide-controlling criteria
A landslide is a physical phenomenon that occurs due to the synergic action of many environmental factors (Dikshit et al. 2020). Therefore, to accurately estimate the spatial risk of landslides, it is very critical to select the triggering factors which have a high influence on the landslides. However, it should be noted that these were selected depending on the initial data availability, type and pattern of landslides, geographical features of the study area and Mediterranean landslide literature (Yu et al. 2022;Senouci et al. 2021;Ullah et al. 2022). At the present assessment, 13 landslide conditioning factors were selected (i.e. slope degrees, slope aspect, curvature, plan curvature, profile curvature, altitude, distance from faults, distance from streams, distance from roads, Normalized Difference Vegetation Index (NDVI), precipitation, geology and Topographic Wetness Index (TWI).

Slope angle
Slope angle is a main terrestrial motivating factor influencing slope debris stabilization (

Curvature
Curvature represents another morphometric landslide predictor, derived from DEM, which has a high influence on landslide susceptibility by shifting the values of slope angle or aspect. Curvature effects throw given by the control that the values of the curvature have over the phenomenon of soil erosion and the rapid flow of water on the slopes (Costache and Tien Bui 2020). In this regard, negative curvature shows an area of land that is concave, while a positive curvature highlights a convex surface of the ground. In this assessment, the values of the curvature were divided into the following classes: <− 0.05 (negative curvatures), − 0.05-0.05 (flat) and >0.05 (positive curvatures) as shown in Fig. 3c.

Plan curvature
The plan curvature, which is derived also from DEM in GIS environment, is an important landslide conditioning factor because its values can indicate the areas with a convergent or a divergent runoff (Costache and Tien Bui 2020) Plan curvature affects the dramatic change of water channels distance by the slope flowing (Ullah et al. 2022). Figure 3d shows the plan curvature map that was extracted from the DEM in GIS environment and categorised into three groups: <− 0.05 (negative plan curvatures), − 0.05-0.05 (flat) and >0.05 (positive plan curvatures) as presented in Fig. 3d.

Profile curvature
Profile curvature is an important morphometric landslide influencing factor because it shows the areas with accelerated soil erosion and runoff (Costache 2019). Thus, profile curvature essentially triggered the slope debris movement (Di et al. 2019) Profile curvature values in this study were calculated by using DEM in GIS environment and classified in the following classes <− 0.05 (negative plan curvatures), − 0.05-0.05 (flat) and >0.05 (positive plan curvatures) as showed in Fig. 3e. Abdo et al. Geoscience Letters (2022)

Elevation
The elevation conditioning factor is a highly used parameter in landslide susceptibility studies due to the fact that this factor controls many other climatic as well as geomorphological parameters (Costache and Tien Bui 2020). It is often considered that the susceptibility to landslides is higher in areas with high altitudes and vice versa (Chen et al. 2018;Costache and Tien Bui 2020). In this study, the elevation factor map was extracted using DEM and divided into five classes: <200 m, 200-400 m, 400-600 m, 600-800 m and >800 m as shown in Fig. 3f.

Proximity to faults
The proximity to faults is an essential geological factor triggering landslide sensitivity since the likelihood of landslides will increase as the distance to faults decreases (Fang et al. 2021). Primarily, the landslide phenomenon will occur along the faults that enhance slope instability. The proximity to faults map, in the current evaluation, was calculated using the Euclidean Distance tool in GIS software as shown in Fig. 3g.

Proximity to rivers
Similar to the case of the previous factor, the proximity to rivers is a significant triggering criterion for landslide susceptibility (Ma et al. 2020). In fact, a watercourse can destabilize the stabilization of the slope by flow energy and erosion capacity, especially down the slope (Sahana et al. 2020). In the case of the current investigation, the proximity to rivers was grouped into five classes as shown in Fig. 3h.

Proximity to roads
The proximity to roads is an anthropogenic factor that influences directly triggers slope materials instability (Karlsson et al. 2017). The road construction process and the weight and traffic of vehicles could frequently lead to landslide occurrence. The road network was derived from the Directorate of Transportation of Tartous Governorate (DTTG) and buffered with a Euclidean distance of 100 m in the GIS environment to paper the proximity to the road map (Fig. 3i).

NDVI
The Normalized Difference Vegetation Index (NDVI) is a wide landslide conditioning factor used in the majority of the previous landslide assessment literature (Hussain et al. 2022;Yousefi et al. 2022;Huang and Zhao 2018;Wu et al. 2020). NDVI map was created by utilizing Landsat 8 imagery and Eq. 9: where NIR is the near-infrared band (band 4, 0.76-0.90 μm) and RED is the red band (band 3, 0.63-0.69 μm). Figure 3j showed the classes of NDVI values in the study area.

Rainfall
Rainfall is a principal climatic criterion that motivates slope materials instability by the heavy rainfall intensifications, the high kinetic energy of the raindrops and generated runoff (Youssef et al. 2015;Bui et al. 2019;Sahana et al. 2020). In this regard, landslide events accelerate in the Mediterranean mountainous regions due to the orographic precipitation pattern which produces heavy rainstorms with great peaks of runoff (Mohammed et al. 2020a). The rainfall map in the study area was delineated depending on the five rainfall stations data  obtained from the Directorate of Meteorology in Tartous governorate. The inverse distance weighted (IDW) method, however, was utilized for mapping (Mohammed et al. 2020a;Valiante et al. 2021). Figure 3k depicted the five spatial domains of rainfall: >1000 mm, 1000-1100 mm, 1100-1200 mm, 1200-1300 mm and >1300 mm.

Lithology
Lithology is a landslide conditioning factor that can pro- where α is the cumulative up slope area draining through a point (per unit contour length) and tan β is the slope (10) TWI = ln α tan β angle at the point. In the present assessment, the TWI values were divided into four spatial domains as: <3) low wetness), 3-6 (moderate wetness), 6-9 (high wetness) and >9 (very high wetness), as shown in Fig. 3m.

Landslides correlation with conditioning criteria Application of FR
As per the equation of FR (Eq. 1), if the value is >1, it signifies that there is a positive correlation between training points and a particular class of landslide conditioning factors and high landslide susceptibility. The result of the FR model showed that the slope with 20°-25° and >25° have a ratio value of 1.31 and 1.26, respectively. The slope aspect facing east, southeast and south, all have a ratio value >1. Convex curvature with ratio value 1.12, 600-800 m and >800 m elevation with ratio value 1.22 and 1.23, respectively, NDVI ranges between 0.3 and 0.6 with a ratio value 1.19, proximity to a fault between 100 and 200 m with a ratio value of 1.07, proximity to rivers between 100 and 200 m and <100 m with ratio value 1.33 and 1.11, respectively, rainfall >1300 with 1.33 ratio value, among lithology J2 and J1-2 with ratio value 1.81 and 1.82, respectively and greater value of TWI (>9) has a ratio value 1.07, indicating that these all are a more relevant class of the selected criteria having a significant role in landside occurring and positive correlation with landslide susceptibility ( Table 3).

Application of SI
In the case of a SI, a positive and negative value of each class of landslide conditioning factors indicates the presence and absence of a relationship with landslide susceptibility. The result of SI is more similar to FR. Table 3 shows that those classes of landslide conditioning factors have an FR value of >1, having a positive SI value. The result of the calculated SI value shows that the slope is between 20° and 25° and >25° have SI values of 0.27 and 0.23, respectively. Slope aspect facing east (SI = 0.08), southeast (SI = 0.52) and south (SI = 0.32), all have a positive value. Convex curvature (SI = 0.11), elevation with 600-800 m and >800 m (SI = 0.20 and 021, respectively), NDVI ranges between 0.3 and 0.6 (SI = 0.18), proximity to a fault between 100 and 200 m (SI = 0.07), proximity to river between 100 and 200 m and <100 m (SI = 0.28 and 010, respectively), rainfall >1300 (SI = 0.32), among lithology J2 and J1-2 (SI = 0.59 and 0.60, respectively) and the value of TWI with >9 (SI = 0.07), representing important classes causing landslides in the stud area.

Application of IoE
Based on Eqs. (3-8), the individual factor's weight has been calculated for preparing the landslide susceptibility index. A higher value of the index of entropy (IOE) indicates more causative for landslide occurring. The result of IOE weight reveals that proximity to faults, lithology and proximity to roads having an IOE value of 1.26, 1.22 and 1.19, respectively, followed by proximity to river (IOE = 0.99), rainfall (IOE = 0.91) and TWI (IOE = 0.74). Hence, these are more important factors for evaluating landslide susceptibility in the study area out of the thirteen selected. On the other hand, curvature (IOE = 0.004), plan curvature (IOE = 0.01), slope (IOE = 0.02), profile curvature (IOE = 0.04), elevation (IOE = 0.04) and slope aspect (IOE = 0.07) are less important factors for landslide susceptibility assessment in the study area with lower IOE value (Table 4).

Landslide delineation and assessment
In this analysis, three different statistical approaches were tested in modelling the landslide susceptibility in Al-Balouta river basin. In this regard, data derived from augmented fieldwork and remote sensing in a GIS environment were used in the modelling process. Figures 4, 5 and 6 show the results of the multi-criteria modelling process after being categorized using the Natural Brecks method in a GIS environment into five degrees of severity: very low, low, moderate, high and very high.
Moreover, Fig. 7 indicates the classification of landslide pixels for each susceptibility degree in the study area. In the FR-based modelling, 18.51% and 19.98% of the study area fall under high and very high landslide susceptibility, respectively. While moderate, low and very low covered 29.09%, 24.63% and 7.78% of the area under investigation, respectively. In the map generated by the SI method, about 36% of the study area is classified under high and very high landslide sensitivity. Whilst, the remaining 25%, 23.57% and 15.30% of landslides are classified as moderate, high and very high landslide susceptibility, respectively. In the IOE method, 14.18% and 25.62% of the study area were classified as very high and high landslide susceptible, respectively. While moderate, low and very low landslide susceptibility covered 19.87%, 30.04% and 10.29% of the study area, respectively.

Validation of landslides susceptibility map
Landslide vulnerability prediction maps are among the most important tools for managing natural hazards associated with slope stability in river basins. Statistical models that spatially linking between landslide incidents and causative factors provide a reliable and constructive approach in the landslide sensitivity mapping process. However, evaluating the accuracy of the conducted modelling process outputs is a necessary final procedure in order to verify the applied models performance. At the present study, landslide capability maps acquired by the FR, SI and IOE were evaluated using the verifying data sets excluded from the modelling process. The area under curve (AUC) of the receiver-operating characteristics (ROC) was utilized for the standard accuracy evaluation of both landslide sensitivity outputs and the three models performance. The forecasting rate curves (Fig. 8a) showed that the applied models performance in landslide susceptibility modelling is greater for FR with 0.841 of AUC followed by SI and IOE with 0.821 and 0.788 of AUC, respectively. Moreover, Fig. 8b illustrated that the AUC of performed models success rate was larger for FR with 0.824, followed by SI and IOE with 0.791 and 0.776 of AUC. In this sense, prior validation outcomes illustrate that the used curative factors have generated constructive maps with a great rate of accuracy. Although the FR method has achieved the highest performance accuracy, the accuracy of the SI and IOE models is considered constructive in assessing the landslides susceptibility in the study area.

Importance of implemented models and key parameter
In the present investigation, the Seed Cell Area Index (SCAI) method was used to assess the significance of each applied model. SCAI enables proportional calibration between the test dataset with derived landslide susceptibility zones by each applied model (Süzen and Doyuran 2004;Li et al. 2021a, b;Rehman et al. 2022).
In the context of SCAI values, higher values in very low susceptibility and lower values in very high susceptibility indicate higher reliability of the model (Rahman et al. 2022). Although the FR model achieved the best reliability according to the SCAI results (Table 5), the other two models also showed a good capability in determining landslide susceptibility in the study area. Additionally, the relative importance of each key parameter in generating landslide events was calculated depending on the weights calculated in Tables 3 and 4. Figure 9a shows that the slope, lithology, proximity to roads and elevation parameters were the most influential in the FR model. Similarly, lithology, slope and proximity to roads effectively motivated the acceleration of slope materials instability in SI model (Fig. 9b). in IoE model, Fig. 9c depicts the important influence of proximity to faults, proximity to roads and proximity to rivers factors in the creation of landslide events.

Discussion
Spatial assessment of landslide susceptibility is a critical basis for creating safe spatial development, especially in areas with insufficient spatial data. Thus, many studies aim to produce landslide susceptibility maps using different modelling methods in a GIS environment (Chowdhuri et al. 2020;Tesfa 2022;Zhang et al. 2022). Population and infrastructure are exposed to frequent landslide events in CMA as a result of the complex spatial interaction between a set of physical and human geographical factors. Hitherto, there is a large gap in the national literature concerned with conducting landslide prediction studies. Hence, there is an urgent need to provide  in-depth spatial assessments of landslide risk that assist in the management and mitigation process. In this evaluation, the performance of three bivariate statistical methods was tested and compared in landslide susceptibility mapping in Al-Balouta river basin. The inventorization process of landslide events is the primary step in landslide susceptibility modelling. This process involves multiple procedures in areas without documented records of current landslide events. Despite the high frequency of landslide events in most of the watersheds in CMA, there is a loss of historical records for the landslide locations. In the current study, a technical-field strategy was based on producing an inventory map of 446 landslide locations. However, this strategy can be relied on in most areas of CMA in generating similar inventory maps for the different types of slope material movement (Es-Smairi et al. 2022). This inventory map enabled a perfect modelling process with thirteen driving factors. The selection of causative factors is a critical task in the context of the most reasonable extraction of landslides using bivariate statistical models (Wubalem et al. 2022). The fieldwork carried out showed the diversity of landslide patterns, including rainfall-induced rapid mass movements, rock falling, bouncing and rolling, shallow landslides and debris and mud flow. These patterns represent a direct reflection of the complex interaction  Fig. 7 The proportion (%) of landslide pixels for each susceptibility degrees between natural and human characteristics in the study area, Bivariate statistical methods are flexible in implementation with constructive spatial outcomes that assist in managing the landslide risk Kincal and Kayhan 2022). In this investigation, the FR model provided the best performance in comparison with the SI and IoE models. In addition to the performance accuracy, the FR model was also identified as the best model in terms of landslide classification capacity. The FR model, however, is based on a direct spatial correlation through the application of an independently relativistic mathematical structure between current landslide events and classifications of causative factors. Thus, direct monitoring of the spatial sensitivity of these classifications through quantitative discrimination of landslide events generation. However, this result allows reducing limitations in this assessment and consistents with landslide assessment studies conducted worldwide (Babitha et al. 2022;Alsabhan et al. 2022;Akter and Javed 2022). Despite the better performance of the FR model, the SI and IoE models also provided reliable performance in landslide susceptibility mapping. In the context of each key factor influencing, slope, lithology and proximity to roads factors were the most influential in landslide events creation. This result is consistent with observations of fieldwork that emphasized the influence of lithological structure and human activity as actual driving factors. Similar results were reported by Tesfa (2022), Senouci et al. (2021) and Yamusa et al. (2022).
In this regard, the outputs of the modelling process indicated that the high and very high landslide susceptibility zones were mainly concentrated in the eastern and northeastern parts with some middle parts along the riverbeds. The fieldwork revealed that the terrain of these areas is characterized by structural fragility and steep slopes with more than 50° in some locations. Also, these areas are characterized by high-intensity rainstorms that cause more landslide events as a result of crossing rainfall thresholds for landslide ). These  characteristics are integrated with the great acceleration of the karstification process, which causes many slope instability patterns, including karst rockfalls and slides of rocks and soil mixture. Karst land systems are characterized by their significant susceptibility to a combination of physical processes and human activities (Chen et al. 2022). Many studies have indicated a great correlation between the karstification process and landslides, especially in the Mediterranean environment (Devoto et al. 2021;Pisano et al. 2022). Moreover, the acceleration of human activity, especially infrastructure projects, has contributed to the instability of slopes. Unfortunately, the spatial distribution of landslide events is not considered when planning these projects. For instance, the road network is built on steep slopes with large loads on slope aspects. These spatial interpretations are consistent with Jamir et al. (2022), Das et al. (2022) and Yamusa et al. (2022). In this evaluation, FR, SI and IoE models provided a satisfactory performance for a comprehensive assessment of landslides susceptibility in Al-Balouta river basin. The outputs of this study are critical for planners and decision-makers in light of the paucity of relevant geographical data. Moreover, this study represents a constructive contribution in the context of enhancing the national landslides studies, especially during the current war conditions in Syria.

Conclusion
In this research, landslide hazard maps were produced based on three models (FR, SI and IOE). Although there are advantages and limitations to applying the landslide susceptibility models, they provided a reliable and constructive approach to landslide sensitivity mapping and proved that there is no statistical model that suits all geo-environmental variables. Three alternative statistical approaches were investigated in order to increase the accuracy and flexibility of landslide susceptibility mapping within a GIS framework. To demonstrate the efficiency of the suggested models, a case study of landslide susceptibility mapping in the Al-Balouta river basin in northern Syria is done. The output of this research can be summarized as follows: 1. For FR and SI models, the relevant classes for the selected criteria that have a direct impact on landslide susceptibility were: (1) slopes with 20°-25° and >25°; (2) the slope aspect facing east, southeast and south; (3) convex curvature; (4) elevations between 600 and 800 m and >800 m and >800 m; (5) NDVI ranges between 0.3 and 0.6; (6) proximity to faults between 100 and 200 m; (7) proximity to rivers between 100 and 200 m and <100 m; (8) rainfall >1300 mm, (9) lithology J2 and J1-2; and (10) a greater value of TWI (>9). 2. For IOE model, the most important criteria were: (1) proximity to faults; (2) lithology; (3) proximity to roads; (4) proximity to river; (5) rainfall; and (6) TWI. 3. Based on the three applied models, the landslide hazard maps indicate that 14-18% of the study area was categorized as very high susceptibility to landslide, while 20-25% of the study area was classified into high susceptibility to landslide. . FR method has achieved the highest performance accuracy, however, the accuracy of the SI and IOE models is considered constructive in assessing the landslides susceptibility in the study area. 5. The relative importance analysis demonstrated that the slope aspects, lithology and proximity to roads effectively motivated the acceleration of slope material instability and were the most influential in both the FR and SI models. On the other hand, the IoE model indicated that the proximity to faults and roads, along with the lithology factor, were important influences in the formation of landslide events.