Landslide Susceptibility Mapping in Brazil: A Review

: Landslide susceptibility studies are a common type of landslide assessment. Landslides are one of the most frequent hazards in Brazil, resulting in signiﬁcant economic and social losses (e.g., deaths, injuries, and property destruction). This paper presents a literature review of susceptibility mapping studies in Brazil and analyzes the methods and input data commonly used. The publications used in this analysis were extracted from the Web of Science platform. We considered the following aspects: location of study areas, year and where the study was published, methods, thematic variables, source of the landslide inventory, and validation methods. The susceptibility studies are concentrated in Brazil’s south and southeast region, with the number of publications increasing since 2015. The methods commonly used are slope stability and statistical models. Validation was performed based on receiver operating characteristic (ROC) curves and area under the curve (AUC). Even though landslide inventories constitute the most critical input data for susceptibility mapping, the criteria used for the creation of landslide inventories are not evident in most cases. The included studies apply various validation techniques, but evaluations with potential users and information on the practical applicability of the results are largely missing.


Introduction
Natural hazards are frequent in Brazil, and landslides are considered one of the most common hazards in the country [1][2][3]. They are responsible for causing losses almost every year, resulting in many economic and social problems. The high frequency of landslides in Brazil is due to natural characteristics such as geomorphology, the tropical climate, and frequent rainfall [3].
Most of rainfall-triggered fatal landslide events in South America are concentrated in Brazil (37%) [4]. Several high-magnitude events have occurred in the country in the last few decades, most in the south and southeast region [2]. All high-magnitude events were triggered by rainfall and resulted in fatalities and infrastructure damage.
The south and southeast regions of Brazil are the most affected by landslides due to the presence of mountainous areas [3] (Figure 1). The Serra do Mar is the main mountain range near the coast in the south and southeast. Mass movements frequently occur on its steep slopes, resulting in economic and social losses [2]. Important highways cross the area, connecting the plateau to the coastal cities. Thus, landslide susceptibility studies in the Serra do Mar are critical and may help to avoid future losses. The maximum elevation of the Serra do Mar varies from 900-1200 m in São Paulo up to 1800 m in Paraná [5]. Despite having morphological conditions that are favorable to mass movements, the Serra do Mar also presents high values of rainfall. The annual average precipitation for the area is between 1000 and 2000 mm [3], eventually reaching 5000 mm when associated with tropical cyclones [1].
Landslide studies in Brazil address several topics related to the process and try to understand under which conditions landslides occur, offering potential mitigation actions to avoid future losses. One of the first steps in landslide analysis is the construction of landslide inventories. These inventories can be created by means of different methods (manual, semi-automatic, and automatic) for various purposes, such as the documentation of recent events on the regional or local scale, and as the first step in susceptibility, vulnerability, and risk analysis [6]. Landslide inventories in Brazil are usually created manually through visual image interpretation, without considering standards and common guidelines, and their methodological basis is not addressed frequently in publications [7]. Therefore, guidelines for landslide inventorying are crucial and can influence the accuracy of subsequent analyses based on these inventories. The susceptibility analysis of landslides is one of the main procedures designed to understand which conditions were responsible for past landslide occurrences, which could generate new landslides in the future [6,8]. Landslide susceptibility refers to the likelihood of landslides to occur in certain conditions (predisposing and triggering factors) on a local or regional scale [9,10], or the tendency of an area to generate landslides [11]. Its determination aims to provide spatial information that helps to prevent future loss of human lives and infrastructure. The analysis is based on the characteristics of the environment and strongly relies on the reliability, accuracy, and completeness of the landslide inventories used. According to Aleotti and Chowdhury [8] the analysis can be qualitative or quantitative, and depends on the aim of the study, the study area, scale, and the data that are available.
Qualitative approaches are based on expert experience and fieldwork [8], whereby the expert is responsible for the selection and mapping of the factors that influence the occurrence of landslides and the weighting of each factor based on personal experience [12,13].
In the quantitative approaches, homogeneous zones are created based on overlay maps and weighted factors with the application of statistical analysis, physically-based models, and machine learning techniques, rather than personal experience [14][15][16][17][18][19]. In general, quantitative approaches are used when detailed geotechnical information is not available for the whole study area, so the models recognize patterns through computation [8]. Statistical approaches can compare the spatial distribution of landslides to predisposing factors (e.g., morphological and geological factors) and thus identify the most decisive factors for indicating potentially unstable areas. Physically-based models evaluate a local analysis (i.e., of specific sites and slopes) of geotechnical characteristics and the calculation of safety factors (SFs). On the other hand, machine learning techniques rely on algorithms that react to training data, which must be capable of learning and create scenarios of susceptibility [20].
Landslide susceptibility studies address several aspects: (1) location-susceptibility analysis is undertaken for areas with previous landslide events, since past events indicate that a particular locality can be susceptible; (2) thematic variables-certain conditions can facilitate or influence landslide occurrence, such as morphological, geological, and geotechnical factors; (3) validation-the final susceptibility map needs to be validated to ascertain if the mapping can be used or not in prevention and mitigation activities.
This paper aims to analyze the methods and input data commonly used in landslide susceptibility mapping in Brazil, via a bibliographic analysis, considering six main criteria: location of the study areas; year of publication and where it was published; susceptibility assessment methods; thematic variables; origin of the landslide inventory; and validation methods.

Materials and Methods
We applied a systematic search to find publications about susceptibility assessments in Brazil (Figure 2), using one of the main and largest bibliographic repositories commonly used in several geoscience studies [21][22][23][24], the Web of Science™(WoS), as of April 2021. Journals and other publication media indexed in the WoS need to fulfill several quality criteria; thus, by relying on the WoS it can be expected that the identified articles correspond to a certain scientific standard. Publications not listed in the WoS were not considered in this study since this would have required an extensive manual search, wherein a complete and objective coverage of studies would have hardly been possible.
The keywords applied were "Brazil", "Susceptibility", and "Landslide" for the title, keywords, and abstract of the available publications. The analysis considered only peerreviewed articles and proceedings papers written in English and Portuguese that applied a susceptibility analysis to landslides in the Brazilian territory. Based on our search strategy and criteria, the publications were identified within the WoS platform, and a manual check was performed to verify whether all publications fitted the criteria (Figure 2).
We analyzed the following aspects from the compiled publications: (1) Location of study areas-we collected the name of the federal state of each study area and checked the number of studies per region. (2) Year of publication and where it was published-we identified the year, if the publication was published in a journal or proceedings, and the name of the journal and the event (conference). (3) Susceptibility assessment methodswe collected all the methods and models for the creation of the susceptibility mapping.
(4) Thematic variables-we collected all variables used for the assessment and organized them into six categories (climatic, geological, geotechnical, hydrological, land cover, and morphological). (5) Origin (source) of the landslide inventory, and if it was constructed specifically for the susceptibility assessment or if the authors used existing inventories. (6) Validation methods-we identified the method used to evaluate the model performance.

Results and Discussion
The search on the WoS platform resulted in 67 publications. Only 33 were directly related to landslide susceptibility mapping in Brazil. These 33 publications, published in international, in Brazilian journals that publish articles in the Portuguese language as well as in conference proceedings, were considered for further analysis. Several Brazilian journals are not indexed in the WoS, so articles from these journals were not considered in our study.
The primary reason for the concentration of Brazilian landslide susceptibility studies in the southern part of the country is the presence of the Serra do Mar mountain range, located between the states of Rio de Janeiro and Santa Catarina [5]. This region is frequently affected by mass movements (Table 1 and Figure 4), especially the states of Rio de Janeiro and São Paulo, which have together registered more than 3200 deaths due to landslides since 1928 [2].

Year of Publication and Where It Was Published
The number of publications related to landslide susceptibility mapping in Brazil has increased significantly since 2015 ( Figure 5). The first two publications were from 2006 [39,47]. Between 2015 and April 2021, 27 articles were published, comprising about 82% of the database. According to Reichenbach et al. [54] there has been an increasing global trend of new publications on the susceptibility analysis of landslides starting in the year 2000. In Latin America, the number of studies on landslide hazards has been increasing since 2011 [55]. Brazil is following this regional and global trend ( Figure 5).
In the WoS database, there are two types of publication: articles (26) and proceedings papers (7). The proceedings papers were published in English and were from four different conferences: the International Symposium on Landslides (4) [13,33,50,51]; the Lusophone Conference on Sciences and Technologies of Geographic Information (1) [52], the International IAEG Congress (1); and theInternational Conference on Evolution, Monitoring, Simulation, Management and Remediation of the Geological Environment and Landscape (1) [47]. The articles were published in 15 different geosciences journals ( Figure 5). The journals address several topics, such as geology, geomorphology, and remote sensing. We found that 73% of the articles were published in international journals in the English language. According to Carrera et al. [55], in Latin America 48.8% of all publications were in English, followed by Spanish (43.1%). These numbers show that Brazilian researchers tend to publish in English rather than in Portuguese (the native language). Natural Hazards is the journal with the highest number of publications about susceptibility assessments in Brazil (≈27%), and it is among the five most common journals with publications related to landslides in Latin America [55]. According to Reichenbach et al. [54], Natural Hazards is also one of the journals with the highest number of publications about landslide assessments worldwide. The remaining 27% of articles were published in Brazilian journals: GEO UERJ, Anais da Academia Brasileira de Ciências (Annals of the Brazilian Academy of Sciences), and Confins and Revista Brasileira de Geomorfologia (Brazilian Journal of Geomorphology). Only the publication from Anais da Academia Brasileira de Ciências was in English [40]; the other articles were in Portuguese. Revista Brasileira de Geomorfologia is the Brazilian journal with the highest number of publications.

Susceptibility Assessment Methods
The publications applied 21 different susceptibility model types ( Figure 6). About 37% were slope stability models: the Shallow Landsliding Stability Model (SHALSTAB) (20.9%), Transient Rainfall Infiltration and Grid-Based Regional Slope Stability (TRIGRS) (9.3%), Stability Index Mapping (SINMAP) (4.7%) and Safety Factor FIORI (SF FIORI) (2.3%) models. In general, slope stability models use hydrological, geomorphic and geotechnical input parameters, computation of safety factors (SFs), and in some cases they enable the construction of susceptibility scenarios based on rainfall thresholds [49]. SHALSTAB was the most frequently applied model [26][27][28]31,38,[48][49][50][51]. It is a physically-based digital terrain model for mapping slope stability potentially related to landslides [56]. Two publications compared some of these models and assessed their performance: Ávila et al. [15] used TRIGRS and SF FIORI and concluded that TRIGRS leads to better results. Vieira et al. [28] compared SHALSTAB and TRIGRS, and concluded that both models were effective in predicting landslides. The first publication applying a slope stability model was published in 2010 [32], and the most recent was published in 2021 [15]. Statistical models are also commonly used, such as bivariate analysis, logistic regression, and fuzzy analysis. According to Reichenbach et al. [54], logistic regression analysis counts for 18.5% of all statistical models applied to landslide susceptibility assessments worldwide. In Brazil, the most frequently applied statistical method is "informative value", a bivariate statistical approach used in landslide susceptibility studies [57], whereas logistic regression analysis has a share of 4.7%. Bivariate statistical analysis calculates the weight of the classes of predisposing factors (e.g., the weight of each rock type in the study area) by landslide density in each class [8]. The first paper that used statistical models for landslide susceptibility mapping in our database was from 2019 [16]. The authors compared several statistical approaches (Bayesian model, informative value, weights of evidence, logistic regression, and discriminant analysis), and the most recent papers that applied these methods were published in 2021 [25,40]. One publication compared the SHALSTAB model and the informative value approach in São Paulo state. The results suggested that the informative value approach shows better performance than SHALSTAB [26]. However, further analysis should be carried out in other regions of Brazil for comparison.
The application of slope stability models to landslide susceptibility assessments is more common in Brazil, whereas statistical models seem to be becoming more prominent, following the worldwide trend [54].

Thematic Variables
The 60 identified input thematic variables were grouped into six thematic clusters: climatic, geological, geotechnical, hydrological, land cover, and morphological ( Figure 7). Some named thematic variables that were considered synonyms were grouped; for example, slope angle, slope gradient, and slope were grouped into slope.
Reichenbach et al. [54] analyzed the thematic variables that are most frequently used worldwide and found that slope was used most often, followed by geo-lithology, aspect, curvature, and land cover. The authors only analyzed papers that applied statistical approaches. Thus, it is impossible to compare the mechanical factors used in Brazil and their worldwide use.

Origin of the Landslide Inventory
Compiling an inventory of landslides is the first step in a susceptibility assessment [6]. Landslide inventories are primary data commonly used to trace the patterns of conditions that may cause landslides in the future. They can be organized based on different perceptions and aims [58,59] and can be used for several purposes, such as landslide documentation on the local or regional scale, or susceptibility, vulnerability, and risk analysis [6]. According to Dias et al. [7] there are no standard guidelines for the construction of landslide inventories in Brazil. Most inventories are constructed manually through visual image interpretation (using aerial photographs, unspecified satellite images, and Google Earth Pro images). Dias et al. [7] identified only six publications about landslide inventory mapping in Brazil. Information on the inventories' origin, completeness, accuracy, and up-to-dateness are essential for further analysis, such as landslide susceptibility mapping. However, this type of information is often unavailable, even though the criteria applied to landslide recognition and mapping are essential. The accuracy and reliability of susceptibility mapping results strongly depend on the quality and completeness of the inventory; an inventory with errors can influence the susceptibility analysis, creating distortions [60]. We found that 32% of the landslide susceptibility studies in our database used an inventory of landslides that was based on previous work (Figure 8). The authors generated a new inventory as an input for susceptibility mapping in 56% of the articles, and 12% did not mention any inventory in their analysis.
Only a few of the studies that created a new landslide inventory (eight of 19) defined the criteria for the recognition and mapping in the article [13,16,28,30,32,34,40,49]. Table 2 shows that there is no standard method for the recognition of landslides in the Brazilian literature. Each publication applied different criteria, and in several cases, the criteria are not clear. The terms "source of the landslide", "where the landslide started", and "rupture area" are examples of criteria that are not clear to the reader, and are subjective. Two publications applied criteria based on the literature [16,40], but they did not specify them. Visual interpretation was applied by five out of eight publications, whereas the other three did not describe how they carried out the recognition. The images used for landslide mapping varied between aerial photographs, orthophotos, satellite images (often without providing specifications), and free high-resolution Google Earth images. König et al. [27] is one example of publication that described the satellite image used: an image from the IKONOS satellite, dated April 2000, two months after the landslide event. Verification of the features during fieldwork was also commonly carried out by the authors [13,15,16,45,51].

Authors Criteria Data
Vieira et al., 2010 [32] Lighter patches with more contrast in texture, and in areas without vegetation (polygons).
Michel et al., 2014 [49] Only the source of the landslides. The transport and deposition areas were not analyzed.
Nery and Vieira, 2014 [30] Scar geometry, absence of vegetation, position on the slope, contour lines, and texture analysis.
Visual analysis.  [61]. Polygons and points representing centroids in the depletion zones, and fieldwork.
Google Earth Pro images and Digital Terrain Models (DTM).
Canavesi et al., 2020 [34] The scars were mapped along their entire length with polygon geometry, and also as points positioned where the slide started.
Google Earth images (visual interpretation), scientific papers, event reports provided by Brazil's National Center for Monitoring and Early Warning of Natural Disasters (CEMADEN). Rosa et al., 2021 [40] Based on Rogers and Doyle (2003) [62].
Orbital and aerial images, and fieldwork.

Validation Methods
The validation of the landslide susceptibility mapping results is the final step in the analysis and enables the comparison of the mapping with other studies in Brazil and worldwide. Nine methods of validation were identified (Figure 9). Sometimes different names were used for the same validation procedure by the authors. We found that 54.1% of all publications applied area under the curve (AUC) (27%), receiver operating characteristic curve (ROC) (16.2%), cross tabulation (8.1%) and success curve (2.7%) analyses.
The first paper that applied statistical validation methods was published in 2015 [52]. Between 2015 and 2021, 27 articles were published and 15 out of them applied cross tabulation [36,52], success curve [26], AUC [14,16,25,35,40,42], ROC [43,46], or AUC/ROC [15,34,44,45] to validate the susceptibility maps. Two publications mentioned validation curves, but did not provide details about them [13,50]. These results demonstrate a growing trend towards the use of these validation methods. Saboya et al. [39] considered only expert opinions to validate the susceptibility mapping, such as those provided by geotechnical engineers and geologists. Publications that did not mention any landslide inventories in their analysis carried out landslide susceptibility validations based on expert assessments [39] or using the risk potential index [31]; four publications did not mention any method of validation [12,33,37,47].
Some of the publications highlight the importance of the developed work and the application of landslide susceptibility studies in Brazil. The authors discuss how the work can be used by research centers, universities, government, and other institutions interested in the prevention, monitoring, and forecasting of landslides [12,15,29,31,[33][34][35][36][37]41,43,44,46,50,52]. Vieira et al. [28] mention the importance of models to provide relevant information about areas susceptible to landslides in a cost-efficient manner. However, there was no mention of practical validation with users and local or regional authorities.

Conclusions
This literature review and analysis enabled us to understand and identify the methods and input data commonly used for landslide susceptibility mapping in Brazil. The results revealed a considerable heterogeneity of methods, thematic variables, the origins of the landslide inventories used, and the validation methods used to assess the final susceptibility maps.
The susceptibility studies were concentrated in Brazil's south and southeast regions due to the very high landslide activity of these areas. Although landslides are also common in the proximity of big cities of the northeast region, such as Fortaleza and Salvador, this region was not studied in the identified publications. A reason that could explain this is the absence of high-magnitude events in these regions. The number of publications about landslide susceptibility has increased since 2015, and most of these have been published in international journals in English. The most commonly used methods were slope stability models (e.g., SHALSTAB) and statistical models (e.g., informative value), and the model validation was accomplished using ROC curves and AUC analysis in most publications. The thematic variables vary, but morphological factors were the most common. Slope was used in most publications, followed by aspect and elevation. In more than half of the studies, the landslide inventory was created particularly for the susceptibility analysis. Even though the landslide inventory constitutes the most essential input data for susceptibility mapping, the criteria for the recognition and mapping of landslides were not clarified in most cases. Moreover, the publications that mentioned the landslide inventory mapping criteria presented no standardized guidelines, confirming the findings of Dias et al. [7]. The number of Brazilian studies about landslide susceptibility mapping has been increasing over the last decade, following the global trend. The reviewed studies applied various validation techniques, but there was a lack of validation together with stakeholders and decision-makers to assess the practical applicability of the generated landslide susceptibility maps.
This review lacks studies in the Portuguese language published in journals not indexed in the WoS. Unfortunately, a manual search to identify all existing studies is hardly feasible and is difficult to implement. However, by relying on the WoS database, we believe that we were able to provide a comprehensive overview of landslide susceptibility mapping in Brazil. It is recommended that more Brazilian journals seek international indexing in reputable citation databases, in order to promote the dissemination of Brazilian research.
For future studies, it is suggested that landslide susceptibility analysis should also be carried out for other regions of the country, and not only in the south and southeast regions. It is also recommended that studies should provide details about the landslide inventory used and its construction, since the resulting susceptibility maps strongly rely on these input data.

Abbreviations
The following abbreviations are used in this manuscript: