Exploring available landslide inventories for susceptibility analysis in Gipuzkoa province (Spain)

Similar to many mountainous regions of the world, landslides are a recurrent geological hazard in the Gipuzkoa province (northern Spain) that commonly cause damage to communication infrastructure, such as roads and railways. This geomorphological process also threatens buildings and human beings, albeit to a lesser degree. Over time, different institutions and academic research groups have individually collected crucial information on historic and ancient landslides in this region, resulting in various landslide inventories. However, these inventories have not been collectively assessed, and their suitability for landslide susceptibility modelling projects has often been assumed without comprehensive evaluation. In this study, we propose a simplified method to explore, describe, and compare the various landslide inventories in a specific study area to assess their suitability for landslide susceptibility modelling. Additionally, we present the results of an illustrative experiment that demonstrates the direct effect of using different inventories in landslide susceptibility modelling through a data-driven approach. We found that out of the five digitally available inventories in the study area, only three provide sufficient guarantees to be used as input data for susceptibility modelling. Furthermore, we observed that each individual inventory exhibited inherent biases, which directly influenced the resulting susceptibility map. We believe that our proposed methods can be easily replicated in other study areas where multiple landslide inventory sources exist, and that our work will induce other researchers to conduct preliminary assessments of their inventories as a critical step prior to any landslide susceptibility modelling project.


Introduction
Landslide assessment has been a recurrent research topic since the early 1980s (Caine 1980;Evans 1982;Cruden and Varnes 1996;Cendrero and Dramis 1996;Schuster et al. 2002;Hungr 2007;Petley 2012;Froude and Petley 2018), with the majority of studies using landslide inventories as baseline data. Spatial susceptibility modelling research is currently one of the most rapidly expanding fields within the realm of landslide studies (Fell et al. 2008;Guzzetti et al. 2012;Yong et al. 2022), for which landslide inventories are essential too. In recent decades, significant advancements have been made in the field of spatial data generation and computational processing (Maguire and Longley 2005;Moraga and Baker 2022). These advancements have also extended to the domain of automatic and semi-automatic landslide feature recognition (Mondini et al. 2021). As a result, there is now a wider range of approaches available for generating landslide susceptibility maps. Recently published reviews have cited up to 3000 articles on the topic of landslide susceptibility modelling. Furthermore, there has been a huge growth in the rate of publication, with 517 articles solely dedicated to landslide susceptibility modelling released in 2020 (Reichenbach et al. 2018;Yong et al. 2022).
Due to the extensive literature available, numerous methodologies have been developed to study landslide susceptibility over time. These methodologies encompass a variety of approaches, including knowledge-based methods, statistically-or machine-learning-based techniques, and physically-based approaches. Knowledge-based methods include expert evaluation of the factors that led to slope failures, and they frequently include numerical analysis matrices like in the analytic hierarchical process (Pourghasemi et al. 2012;Kayastha et al. 2013;Jun et al. 2017). These methods are highly dependent on the researcher's expertise, and they are often poorly reproducible (Corominas et al. 2014). Statistical and machine-learning approaches-also called data-driven methods-are based on the landslides that already occurred in a study area, and whose location are already known. Using Geographical Information System (GIS) tools, the spatial features of landslide occurrence locations are coded and mathematically processed in order to find a model that best fits what the landslide inventory represents (Bovis and Jakob 1999;Tien Bui et al. 2016;Kang et al. 2017;Chmiel et al. 2021). These techniques have the advantage of being easily reproducible and applicable to large areas, although they are completely dependent on the quality and completeness of the input data . Physically-based approaches, in contrast, aim to simulate the comprehensive behavior of slopes through numerical modelling. These approaches rely on specific slope conditions, including factors such as ground water conditions and the presence or absence of structural reinforcements Maleki and Mir Mohammad Hosseini 2022). However, a major scientific challenge associated with physically-based models is the need for high-quality and abundant input data (Hutter et al. 1994;Gariano and Guzzetti 2016;Park et al. 2019;Palacio Cordoba et al. 2020;Marin et al. 2021).
In the case of data-driven methods, the information derived from landslide inventories, which is utilized to train and validate the susceptibility models, is frequently accepted at face value. This is probably due the limited availability of data regarding past landslide occurrences. However, landslide databases often show systematic bias related to the spatiotemporal scale, the considered data sources, and the representation criteria, affecting to their completeness and representativeness (Valenzuela et al. 2017). In fact, some authors highlight the relevant role that landslide inventories play as baselines to calibrate datadriven models Jacobs et al. 2018;Bornaetxea et al. 2022), demanding to pay more attention to the quality and accuracy of landslide databases prior to their use as input data.
Literature review revealed some good examples about landslide inventory assessments by comparing various data sources for the same study area. For instance, Brardinoni et al. (2003) compared the information provided by a photo-interpretation-based inventory and a field-work-based inventory obtained for a densely forested river basin. Their findings showcased a significant loss of information if only the photo-interpretation-based inventory had been relied upon. Bernat Gazibara et al. (2019) also analyzed the differences of three inventories based on landslide statistics, frequency-area distribution, geographical discrepancy of landslides and landslide density maps, which allowed the authors to demonstrate the unreliability of the existing historical inventories. Also, other authors conducted similar studies in larger regions of central and south Italy (Galli et al. 2008;Pellicani and Spilotro 2015). Notwithstanding, all these works only examined the mutual differences among the existing inventories, and the specific result of using each of them to model landslide susceptibility remains unknown. To date, only Pokharel et al. (2021) studied the direct effect of using different inventories to model the landslide susceptibility in the same area. For that study, the authors used five landslide inventories undertaken by different researchers after an earthquake event in central Nepal Himalayas.

Purpose of study
Landslide data collection in the Gipuzkoa province (northern Spain, Fig. 1) has been ongoing since the 1980s. Over the years, various attempts have been made to develop regional and local landslide susceptibility models (see Sect. 3). However, it is important to note that each susceptibility model created thus far has been based on different landslide inventories, and their accuracy and representativeness have not been thoroughly examined. Hence, this study aims to conduct an exploratory analysis of the available data on landslide locations in the Gipuzkoa province. The goal is to use these data in controlled statistical simulations to assess the implications of their use in landslide susceptibility mapping. The specific objectives of this study are twofold: (1) to explore, describe, and compare the various landslide inventories utilized in the study area for generating susceptibility and landslide hazard maps, and (2) to visually illustrate the direct implications of employing different landslide inventories in relation to the final susceptibility map.
Unlike the approach proposed by Pokharel et al. (2021), this study presents comparison and assessment methods specifically tailored to a rainfall-induced landslide environment. These methods are designed to be applicable to various historical or geomorphological landslide data collections. To the best of our knowledge, this is the first direct comparison of landside inventories conducted in the Gipuzkoa province. The insights gained from this research will be valuable for future landslide data collection and susceptibility modelling projects in the area. Moreover, the experiments carried out in this study can be easily replicated in similar regions within Europe or worldwide. This will enable other researchers to assess the suitability of their landslide databases before using them for susceptibility modelling.
The paper is structured as follows: Sect. 2 provides a detailed description of the study area. Section 3 introduces the existing landslide inventories available in the study area. In Sect. 4, we outline the specific data employed to implement the methods presented in Sect. 5. The results obtained from the analysis are presented in Sect. 6, followed by a comprehensive discussion in Sect. 7. Finally, Sect. 8 summarizes the key conclusions drawn from this study.

Study area
The Gipuzkoa province (with an extension of approximately 2000 km 2 ) is located in the northern Iberian Peninsula, along the western end of the Pyrenees (Fig. 1), and represents the typical hilly and mountainous Atlantic landscape. This area presents generalized persistent precipitation values (average: 1597 mm/year) and two main rainy periods: November-January and April (Remondo et al. 2017). The study area comprises two main geological units: the Basque Arc, with Mesozoic and Cenozoic marine sedimentary materials, and the Paleozoic Massif of Cinco Villas, in the northeastern sector (EVE 2010). The heads of the basins are generally dominated by marls, limestone, and calcarenite, with some alternating detrital rock outcrops. In contrast, their middle and ending sectors present strata of marls, limestone, ophiolites, clays, gypsum, and thick detrital grain materials, usually covered by Quaternary deposits (Campos and García-Dueñas 1972;Campos et al. 1983). Our study area is prone to landslides due to this combination of morphological, climatic, and geological conditions. The most common mass movements are small to medium shallow slides (cf. Hungr et al. 2014), debris flows, and rockfalls triggered by either intense rainfall or changes in the slope condition caused by humanmade infrastructures. Larger landslides have also been identified (Remondo et al. 2003;Bonachea 2006;Mavrouli et al. 2019).

Landslide inventories for the Gipuzkoa province
Data on landslide occurrences in Gipuzkoa (Fig. 2) have been collected since the 1980s. Tamés et al. (1986) published a work entitled A Study of Natural Hazards of the Historical Territory of Gipuzkoa (Estudio de Riesgos Naturales del Territorio Histórico de Gipuzkoa) in 1986 ("Natural Hazard (NH) Risk Study" in Fig. 2). The project includes written reports and paper maps of different natural hazard phenomena, including landslide processes, at a scale of 1:5000. Not all municipalities were studied for this project; as a result, the provincial map compilation shows several gaps. To the best of our knowledge, this information has not yet been digitized, and therefore can only be accessed in hard copy.
The Spanish Geologic Survey (IGME for the Spanish acronym) published a synthesis map at national scale in 1987: Mapa de síntesis de movimientos de terreno en España ("IGME map" in Fig. 2;Ferrer Gijon and Ayala Carcedo 1987). Only 41 landslides mapped within our study area are available in this work. Salazar and Ortega (1990) updated the IGME inventory in 1990 including the coastal segment between the Igeldo and Mendizorrotz mountains (Fig. 1); the data are available only in hard copy.
The Gipuzkoa Provincial Council (DFG for the Spanish acronym) conducted a mapping project, simultaneously to the NH Risk Study, that was published in 1991: Geomorphological and Soil Study of Gipuzkoa (Estudio Geomorfológico y Edafológico de Gipuzkoa; "Geomorphologic Study" in Fig. 2; Diputación Foral de Gipuzkoa DFG 1991). The DFG produced printed maps and their associated reports for several municipalities in the province at a scale of 1:25 000. In this work, the presence of slope movements is represented as morphological units within the geomorphologic maps.
The public Spatial Data Infrastructure (IDE for the Spanish acronym) digitized and published those printed maps in 1999, resulting in the geoportal geoEuskadi ("geoEuskadi" in Fig. 2; IDE geoEuskadi 2022). Duque et al. (1990) studied the effects of intense rainfall events occurred on slopes in 1989. They exploited the collected data to assess the validity of the Natural Hazard (NH) Risk Study conducted in the early 1980s.
In 1995, the Basque Regional Council published a regional landslide inventory conducted by the private company INGEMISA (Gobierno Vasco 1995a). This inventory includes the three Basque Regional Council provinces (i.e., Gipuzkoa, Bizkaia, and Araba) and comprises three printed volumes (one per province) of individual fieldwork sheets with detailed information about each landslide and its coordinates. Ormaetxea and Sáenz de Olazagoitia (2017) contributed significantly by digitizing the landslide location points from the INGEMISA 1995 inventory.
In 2013, the Roads Management Department of the Provincial Council of Gipuzkoa conducted a geotechnical assessment of the entire road network aided by private company IKERLUR consultants (DFG 2013). The resulting PDF document contains detailed information about landslide points detected close to roads. Similarly, Mavrouli et al. (2019) designed and applied a rockfall risk assessment approach in a short segment of the road N-634 connecting Zarautz and Getaria municipalities (Fig. 1). Led by the Roads Fig. 2 Chronologic diagram of the different landslide data-collection initiatives for the Gipuzkoa province . Digitized databases are highlighted Management Department, they generated a rockfall inventory based on reports from the emergency services spanning nearly 10 years.
Some of the most relevant contributions to the recording of landslide events in Gipuzkoa came from the Universidad de Cantabria research group. Remondo (2001) generated a landslide inventory from aerial photos taken in 1991 and field visits. This inventory includes only four municipalities located in the Lower Deba valley (western Gipuzkoa; Fig. 1). Remondo et al. (2003) also conducted a second inventory for the same area including landslides that occurred between 1991 and 1997. In addition, since 2000, annual landslide inventories have been compiled for the municipality of Deba by Bonachea (2006), Bonachea et al. (2016), and Rivas et al. (2022), all of them using photointerpretation and fieldwork.
Lastly, Bornaetxea (2018) conducted two different landslide inventories covering the entire study area. The first is a fieldwork-based inventory obtained using field surveys during 2015 and 2016 ("Bornaetxea Field" in Fig. 2; Bornaetxea et al. 2018a). The second is based on a review of newspaper articles that mention landslide events spanning from 2006 to 2015 ("Bornaetxea News" in Fig. 2; Bornaetxea et al. 2018b).
All the abovementioned independent initiatives to collect and record landslide events suggest (i) a generalized interest from several institutions and individuals in recording landslide information due to the significant problems caused by landslides in the province and (ii) that the information is largely dispersed and sometimes inaccessible. These two factors are probably the motivations behind the recently published technical report Exposure of the Gipuzkoa Historical Territory to Landslides and Slope Instabilities under Climate Change Scenarios (Exposición del territorio histórico de Gipuzkoa frente a deslizamientos en masa y flujos bajo escenarios de cambio climático; DFG 2022). This work is part of the Strategy to Prevent Climate Change 2050 promoted by the provincial government of the Gipuzkoa province. According to the authors, the reason for not using data-driven approaches was "the existing unfavorable conditions in terms of quality, representativeness and quantity of the landslide inventory".
Many of the abovementioned inventories have been used in the past to produce landslide susceptibility maps, following both statistical and heuristic (experience driven) methodologies. The INGEMISA inventory was first used to generate a heuristic landslide susceptibility map of the entire Basque Region (Gobierno Vasco 1995b); then Ormaetxea and Sáenz de Olazagoitia (2017) used the same information to generate a new landslide susceptibility map following a semiquantitative approach. Remondo (2001) took advantage of his own inventory to produce a landslide susceptibility map of the Lower Deba area using, among other statistical methods, the favorability functions approach. Subsequently, Felicísimo et al. (2013) used that same inventory to obtain enhanced statistical models (multiple logistic regression, multivariate adaptive regression splines, classification and regression trees, and maximum entropy). The geoEuskadi inventory was adopted by a private company (CGS company) to run a landslide susceptibility model (LSM) for the Gipuzkoa province, applying the linear discriminant analysis method (Diputación Foral de Gipuzkoa DFG 2007). Finally, Bornaetxea et al. (2018a) exploited the information collected from fieldwork to produce several landslide susceptibility maps for the Gipuzkoa province using logistic regression as their modelling approach.
Besides the differences in the modelling approaches, these models were all implemented at different scales or resolutions and used different explanatory variables, thus preventing their direct comparison. However, since they all refer to the susceptibility of landslide occurrence (in its general definition) in the same geographic area, we believe the degree of agreement between them should be investigated.

Data
Among all the landslide inventories mentioned in Sect. 3, the present work focuses on the five digitized landslide inventories that cover, totally or partially, the province of Gipuzkoa, hereinafter referred to as INGEMISA, geoEuskadi, Lower Deba, Bornaetxea Field, and Bornaetxea News (Fig. 3). Specifically, Lower Deba comprises the inventories obtained by Remondo (2001), Remondo et al. (2003), and Bonachea (2006). Landslide occurrence updates mentioned in Bonachea et al. (2016) and Rivas et al. (2022) refer to one single municipality, so we decided not to include these data in the study to maintain homogeneity within the Lower Deba area. Consequently, all landslide occurrences included in Lower Deba predate the year 2000.
Descriptive data and references regarding these inventories are presented in Table 1. The period of time in which these were generated, the method used for data collection, and even their purposes are heterogeneous, challenging their direct comparison.
One of the mayor differences between inventories is the publication date. In an optimal scenario, the most recent inventories should contain all the information collected in previous databases. Another relevant variable is the length of time during which these databases were compiled: from one year, in the case of Bornaetxea Field, to nine years, in the case of Lower Deba.
The scale, directly related to the area originally covered in each project, is a factor that can affect the level of detail of the inventories; for instance, INGEMISA and geoEuskadi were originally conducted in a total area larger than the Gipuzkoa province, while Lower Deba focused in four municipalities within this province, a much smaller area (140 km 2 ).
The geometrical mapping approach chosen in each project can also affect the level of detail of the inventories. Three of the databases studied in the present work are represented spatially as single points, and the other two are mapped as polygons, providing information about the shape and size of each landslide. This is directly related to the data collection method followed in each case. For the INGEMISA inventory, for example, landslide occurrence information was obtained by reviewing the previous  100% -bibliographic references and then verifying this through field visits. In contrast, the Lower Deba inventory relies on the stereoscopic visual interpretation of aerial photographs followed by fieldwork observations. The specific approach followed for geoEuskadi remains unknown, while direct field observations were the main source of information for Bornaetxea Field. For Bornaetxea News, location and occurrence date of landslides were obtained from articles published in local newspapers. The type of landslides mapped in each inventory also varies. INGEMISA and geoEuskadi include every type of landslides-slides, rockfalls, flows, or complex movements-while Lower Deba and Bornaetxea Field focus only in shallow slides. Moreover, in Lower Deba, the observed landslides directly related to road-side works are systematically excluded. Bornaetxea News does not specify the typology of the movements due to lack of detailed information in the newspaper reports.
Regarding the purpose of the inventories, INGEMISA, Lower Deba, and Bornaetxea Field were produced specifically to investigate the landslide susceptibility of the territory; geoEuskadi was also used for this purpose even though its original aim was to describe the general geomorphological characteristics of the province ( Table 1). The aim of Bornaetxea News inventory was more specific: to study the triggering effect of rainfall in landslide occurrences.
Number of total observations varies significantly among datasets. The largest inventory is Lower Deba, with 2290 observations, despite covering a relatively small area compared with those covered in the other inventories. The second largest inventory, Bornaetxea Field, shows 542 observations, followed by INGEMISA, geoEuskadi, and Bornaetxea News with 481, 385, and 324 observations, respectively. Three of the inventories studied herein show a significant abundance of slide-type movements, whereas the other two show either a relatively small amount of slide typology (geoEuskadi; 22%) or no landslide type classification whatsoever (Bornaetxea News). The two inventories mapped as polygons indicate a large difference in size. In geoEuskadi, the mapped features range between 1144 m 2 and 6 772 478 m 2 , whereas in Bornaetxea Field, the polygons range between 7.27 m 2 and 11 620 m 2 . Therefore, landslides in the latter are more than three orders of magnitude smaller than those in the former. In fact, most mapped objects in geoEuskadi seem to represent sliding areas in general terms instead of the specific location of each single landslide (Fig. 3), therefore covering large extensions.

Methodology
Several examples of statistical and geospatial comparative analyses of different landslide inventories for a particular area were found in literature (Carrara et al. 1993;Galli et al. 2008;Pokharel et al. 2021;Bornaetxea et al. 2022). However, similar experiments have never been undertaken in the Gipuzkoa province. Furthermore, the diversity of the available landslide inventories in our study area forced us to adapt some of the experiments proposed in literature to the specifications of the data available. Hence, the methods described in this section were designed in order to highlight the discrepancies between multi-source inventories and to show their predictive capability when used to produce LSMs.
We performed pairwise comparisons that could be summarized in four successive analyses: 1 3

Computation of the pairwise spatial matching index
The study area was divided in slope units (SUs) following the approach suggested by Alvioli et al. (2016). As inputs, we used a digital elevation model with 5 × 5 m cell resolution, 150 000 m 2 as the minimum area for SU subdivision, a value of 0.2 as circular variance threshold, and a reduction factor of 5. We used a threshold of 25 000 m 2 to clean the SU map. Then, we counted the number of SUs containing at least one landslide. We decided to express the amount of SUs with landslides in relative form-i.e., in relation to the total number of SUs in the area covered by each inventory-given that the extension of the area covered by the Lower Deba inventory was smaller than that of those covered in the other databases. The use of SUs allowed the comparison between inventories with different geometry, such as points and polygons, and reduced the positional errors due to scale differences or spatial inaccuracies during digitation.
Next, we adapted the conventional error index proposed by Carrara et al. (1993) and used by other authors (Alvioli et al. 2018;Pokharel et al. 2021). We calculated the spatial matching index (SMI) between two inventories with the following expression: where A and B represent the number of SUs with landslides in each inventory and A ∩ B represents the number of SUs with landslides in inventory A that are also present in inventory B. The unity value (SMI = 100) meant that 100% of the SUs with landslides in inventory A coincided with those in inventory B, and vice versa. The goal of this test was to assess the degree of correspondence in displaying the same landslide locations between different inventories and to detect whether some inventories report landslide occurrences that others do not.

Assessment of the spatial homogeneity using the visibility class test
We followed the approach proposed in Bornaetxea et al. (2022) to assess the spatial homogeneity of the studied inventories. We simulated the visibility of a landslide of ~ 20 m 2 (the fifth percentile of the total landslide sizes among the two polygon-based inventories) from the main road network in Gipuzkoa, counted the number of landslides in each of the five visibility classes, and plotted them into normalized landslide count plots. The rationale behind the visibility class test was that a spatially homogeneous inventory should show a nearly regular distribution of landslides among the visibility classes ).

Exploration of the distribution of the explanatory variables
We selected slope, lithology, and land use as independent variables to perform distribution analyses. With the exception of Bornaetxea News, all inventories have been used in the past for landslide susceptibility modelling. Table 2 presents the spatial variables selected as the best predictors by the different authors who conducted LSMs in the past in Gipuzkoa. For this test, we selected only the explanatory variables present in all the precedent works, i.e., slope, lithology, and land use; these are the three most commonly used variables in literature (Reichenbach et al. 2018). We downloaded the thematic maps from the public spatial data infrastructure of the Basque Country (IDE geoEuskadi 2022) and rasterized them into a grid file with a resolution of 5 × 5 m. Since the original lithology and land use maps contained a large number of classes, these were reclassified following the expert criteria as specified in Bornaetxea (2018). More detailed description and the maps of the explanatory variables can be found in Appendix I. For land use and lithology, we divided the landslide frequency in each class by the total frequency of the same class in the entire study area to represent the distribution of these two variables in a relative manner. For slope, we sampled the slope values and then represented them in a frequency histogram with bins of five degrees each. To facilitate comparison between datasets, we decided to convert the two polygon-based inventories, i.e., geoEuskadi and Bornaetxea Field, into point location databases.

Landslide susceptibility modelling and model comparison
We performed LSMs using the different landslide inventories and following the logistic regression (LR) modelling approach. LR is one of the most commonly used statistical approaches for landslide susceptibility in literature (Reichenbach et al. 2018;Shano et al. 2020) and was useful and reliable in several studies (Goyes-Peñafiel and Hernandez-Rojas 2021; Xing et al. 2021; Rossi et al. 2022), and also for the study area of Gipuzkoa (Felicísimo et al. 2013;Bornaetxea 2018). The mathematical relationship between the dependent dichotomous variable (presence or absence of a landslide in the mapping unit) and the independent variables (i.e., slope, lithology, land use; X 1 , … , X n ) can be expressed as where p represents the probability of landslide occurrence and z is a linear combination of independent variables. The linear combination in Eq.
(1) reads as follows: where b 0 is the intercept of the linear model, b i (i = 0, 1, 2, … , n) represents the coefficient of the regression model, and x i (i = 0, 1, 2, … , n) represents the independent variable. A different landslide inventory was used as dependent variable for each LSM; slope, land use, and lithology were used as independent variables.
To ensure robust validation of our models, we employed the LAND-SUITE software ) for all landslide susceptibility assessments. Following the approach outlined in Rossi et al. (2022), we implemented a rigorous cross-validation strategy for each of the four experiments. Notably, the training and validation subsets were independently divided for every experiment.
In the training phase, we randomly selected 70% of the inventory to train the models, ensuring a representative dataset for learning. Subsequently, the remaining 30% of the dataset was used as an independent validation set to evaluate the models' performance against real observations.
To compare the models' performance, we calculated the area under the receiver operating characteristic (ROC) curve (AUC). This metric allowed us to objectively compare and assess the predictive capabilities of each model. Additionally, we examined the spatial similarity among the resulting four susceptibility maps using pairwise scatter plots and the Pearson correlation index.
(   Fig. 4. We observe a moderate spatial coincidence between INGEMISA and Lower Deba as these presented the same information in 32% of SUs. A slightly smaller matching level was found between Bornaetxea Field and Lower Deba (20% of spatial coincidence). The third larger SMI corresponds to the combination of INGEMISA and Bornaetxea Field (10% of spatial coincidence). Both INGEMISA and geoEuskadi showed low reciprocal SMI (6.5%), even though these two inventories were conducted almost in the same period of time covering the same area. Bornaetxea News shows very low values of spatial coincidence compared with those of the rest of the inventories, and the mismatch is the highest with respect to the geoEuskadi inventory (2.4% of coincidence).
Therefore, the group composed by INGEMISA, Lower Deba, and Bornaetxea Field showed moderate spatial matching between them, whereas geoEuskadi and Bornaetxea News showed poor coincidence with any of the other inventories.

Analysis 2
The results from the spatial homogeneity assessment are shown in Fig. 5. INGEMISA (40% in class 1), Bornaetxea Field (60% in class 1), and Bornaetxea News (67% in class 1) show a high accumulation of landslides in the areas most visible from the roads (class 1) and a pronounced monotonic decreasing trend over the visibility classes.
The case of Lower Deba is more complex. In Fig. 5c we can observe a slight reduction of the collected information related to the visibility classes. However, this behavior is much less pronounced compared with the three abovementioned inventories. Indeed, the difference between the first class (32% of the landslides) and the rest of the classes is not as large as in Fig. 5a, d, and e. Moreover, the amount of data within the less visible classes (23% in class 4 and 5) is much more numerous.
In opposition, the geoEuskadi inventory shows a completely different distribution (Fig. 5b). Landslides are evenly distributed among visibility classes 2-5, although in the first class, a relevant abundance of observations can be found, similarly to the other inventories. Figure 6 presents the relative distributions of land use and lithology, along with the frequency distribution of slope, for each of the inventories. We decided to convert the two polygon-based inventories, i.e., geoEuskadi and Bornaetxea Field, into point location databases to facilitate the comparisons. Based in the chosen methodology and the size of the mapped landslides, we assumed Bornaetxea Field to fairly represent the boundaries of the mass movements; therefore, the points were sampled at the higher position of the polygon to simulate the location of the primary displacement scar. For geoEuskadi, polygon centroids were sampled, since the large areas covered by many of the features suggest an inexact individualization of the mapped landslides.

Analysis 3
The relative distribution of land use (Fig. 6a) shows some similarities between INGEMISA, Lower Deba, and Bornaetxea Field. For these three inventories, the most important land use classes are crops and grassland. However, anthropic land use and rock outcrops are the third and fourth most relevant classes, respectively, for INGEMISA, while these are almost irrelevant for the Lower Deba inventory. For the latter, scrubs or hedges are more significant, whereas forest land use shows a moderate presence in the landslide locations. In Bornaetxea Field, scrubs and hedges are also the third class, closely followed by anthropic land use. Forest is the least present land use, and rock outcrop is completely absent.
In geoEuskadi and Bornaetxea News, the distribution of land use classes shows substantial differences with respect to the other three inventories. The former presents rock outcrops as the most relevant class, while the latter locates the most part of the landslides Relative distributions among the lithology classes are more heterogeneous (Fig. 6b). For the INGEMISA inventory, the most relevant lithology is the magmatic rock, followed by marls and limestone. Clay and detrital rock, together with surface deposits, show a moderate relative frequency, while slate is clearly the less important lithology class. For geoEuskadi, marls are the most frequent lithology, followed by limestone, surface deposits, and clay and detrital rock. In this case, slate is also the least relevant class, while magmatic rock, the most important for INGEMISA, is the second least important lithology. For Lower Deba, the abundance of landslides in marls is remarkable over the rest of the classes; clay and detrital rock, limestone, and magmatic rock maintain a moderate relative frequency, while surface deposits are less relevant. The complete absence of slate responds to the real absence of this lithology in the area analyzed for Lower Deba. In contrast, differences between lithology classes are much less pronounced in Bornaetxea Field. Magmatic rocks, marls, and limestone show similar relative frequencies, and slate is more relevant than for any of the other inventories. Surface deposits remains the least important lithology. For Bornaetxea News, however, the majority of landslides fall in the surface deposits lithological class. Limestone is the second most frequent class, followed by magmatic rock, marls, and clay and detrital rock. Again, the least abundant lithology is slate.
Slope frequency distributions are shown in Fig. 6c, where INGEMISA, Lower Deba, and Bornaetxea Field present again some similitude. In these three cases, frequency follows a normal distribution, with a peak slope angle in the range of 30-35º, although a higher presence of landslides can be found with lower slope values in INGEMISA than in Lower Deba and Bornaetxea Field. The geoEuskadi inventory also follows a Gaussian distribution, but the peak can be found in much lower slope values, i.e., in the range of 15-25°. Once again, the inventory that differs the most is Bornaetxea News, for which a flattish distribution can be observed, suggesting that landslides can be observed at almost any slope inclination from 0 to 40º. Figure 7 shows the LSMs calibrated with INGEMISA, geoEuskadi, Lower Deba, and Bornaetxea Field. In view of the evident differences resulting from the previous tests, we decided not to include Bornaetxea News inventory in this experiment. Considering the spatial patterns shown by the LSMs, we highlight some remarkable differences between inventories. We observe distinctive zones of very low susceptibility in the east sector of both INGEMISA and geoEuskadi maps ( Fig. 7a and b); a susceptibility class that does not appear in the Bornaetxea Field map (Fig. 7d). This area corresponds to the only slate outcrop in the province (see Fig. 10b), so the maps agree well with the relative distribution of the lithology observed in Fig. 6b, where slate shows the lowest relative frequency for INGEMISA and geoEuskadi, but not for Bornaetxea Field. In Fig. 7a, the very low susceptibility class is scarce, as is the very high class. Thus, the majority of the province territory is classified between low and high susceptibility. Conversely, Fig. 7d shows very polarized susceptibility values, with the very low class covering a large sector of the territory, and some very high-class zones distributed throughout the study area. The smooth spatial distribution of the susceptibility obtained with INGEMISA and Bornaetxea Field contrasts with the evident concentration visible in geoEuskadi map and, in a less pronounce manner, in the Lower Deba map (Fig. 7b and c). In both these maps, the valley bottom areas are classified as very low susceptibility and the highest values are defined by marls lithology type, suggesting that susceptibility is highly controlled by lithology in these two models. However, Fig. 7b shows complete absence of very high susceptibility values.

Validation of the models
Each susceptibility model presented above are calibrated with a random sample of 70% among the respective inventories. Computed susceptibility values corresponding to the remaining 30% of the observations are extracted to obtain ROC curves for each. The resulting areas under the curves (AUC) are shown highlighted in red in Fig. 8. The figures obtained suggest that the model calibrated with the geoEuskadi inventory performs the worst (AUC = 0.65); the model calibrated with the INGEMISA model performs slightly better (AUC = 0.68); and the models calibrated with Lower Deba and Bornaetxea Field perform the best (AUC = 0.73 and AUC = 0.85, respectively). Therefore, our results indicate that Lower Deba and Bornaetxea Field fit the data used for calibration satisfactorily, whereas INGEMISA and geoEuskadi show a lower predictive capacity.
Clearly different susceptibility maps are obtained by the visual interpretation of map evidence in each of the four inventories. However, the pairwise scatter plots and their numerical correlations observed in Fig. 8 offer a more objective way to compare these results. The first map, INGEMISA, shows a very low coincidence with geoEuskadi, with which it shares a correlation coefficient value of 0.38 of; the correlation with the Lower Deba map increases to 0.45; while the largest similarity, a correlation coefficient of 0.68, is shown with the Bornaetxea Field map. The geoEuskadi inventory shows the largest correlation coefficient with Lower Deba (0.79); however, its results agree with neither INGEMISA nor Bornaetxea Field (with which the correlation is the lowest: 0.33). Lower Deba and Bornaetxea Field show a significant correlation of 0.67, which indicates that these two maps are moderately in agreement, at least within the Lower Deba area boundaries.

Discussion
The purpose of this study was to investigate, describe, and compare different landslide inventories used for creating susceptibility and hazard maps in the Gipuzkoa province (northern Spain). Additionally, we aimed to visually demonstrate how the use of different landslide inventories directly affects the final susceptibility map. By comparing and analyzing multiple landslide inventories, we sought to advance the understanding of landslide susceptibility modelling and improve the accuracy of the resulting maps. We stress that the objective of this paper is not the finding of the best landslide susceptibility model possible for this region. To do so, apart from the need of a previously assessed and reliable inventory, we would also need to experiment with multiple modelling approaches and an extended list of explanatory variables (Zeng et al. 2023). Instead, in this paper we want to propose some specific methods that can be used to compare the existing inventories in an area, and to highlight the direct effect of using those different inventories on the final susceptibility maps. In this section, we discuss the results of our research and their implications, considering their alignment with existing literature, while addressing their limitations.
According to the results showed in Sect. 6.1, Lower Deba present a very large number of SUs with landslides. This responds to the number of records in the database due to the detailed scale in which this inventory was conducted. However, the second larger inventory, in terms of data quantity, is Bornaetxea Field, which shows the second lowest percentage in Table 3. The use of direct field observations to develop this inventory implies that a relevant part of the study area remained unvisited. As a result, observations are only concentrated in some areas, so more than one landslide record may be located in the same SU. In contrast, Lower Deba was developed using photointerpretation at a higher resolution combined with fieldwork, which allows a more homogeneous treatment of the study area. Surprisingly, the geoEuskadi inventory covers a relevant part of the SUs (1009 SUs with landslide), despite the small number of records in this dataset (385 landslides; Table 1), indicating that many landslide polygons cover more than one SU. This is likely the effect of large and very large polygons that represent sliding zones-rather than specific locations of single landslides-that cross SU borders covering various units at the same time. This presumably overestimates the real extension of individual landslides.
The relevant spatial mismatch observed in Fig. 4 cannot be attributed to the spatial inaccuracy of the mapped landslides because we used SUs for the comparison of landslide locations. Thus, the absence of landslides in some inventories should be attributed to other factors. A possible reason is the large spatial incompleteness of information in some cases, which restricted data collection to specific areas, leaving other parts of the study area unsurveyed. As mentioned above, this issue is often linked to the methodology and the type of data source considered . The interpretation of aerial photographs allows a large percentage of the study area to be covered, while fieldwork is often limited by accessibility, time constrains, or economic reasons (Brardinoni et al. 2003). When press archives are used, the majority of the records are referred to landslide causing some damage, which implies a bias towards the infrastructure network and the most populated areas (Stanley and Kirschbaum 2017). The chosen methodology also determines the type of landslide inventoried, since a recent database may ignore certain types of landslides previously registered in older inventories because different objectives were pursued.
So, these tests confirm that the five inventories analyzed provide a considerably different data, which was something to be expected. Firstly, because differences in the working spatial scale inevitably affect the amount of collected data (Malamud et al. 2004;Galli et al. 2008;van Westen et al. 2013); but also because of the different duration of the data collection time periods. While the Lower Deba database has been almost continuously updated for a period of 9 years, data in Bornaetxea Field was compiled between 2015 and 2016. This may be relevant since the geomorphic features of landslides in steep terrain swiftly disappear due to the action of weathering agents, growth of vegetation, or changes in land use. Thus, a database covering a period of several decades that is frequently updated is more likely to register a larger number of landslides. In contrast, a short period of data collection or the irregular updates of an inventory can result in the loss of information between the successive landslide mapping campaigns (Yang et al. 2018;Tanyas et al. 2021).
Regarding Sect. 6.2, Bornaetxea et al. (2022) highlighted that, inventories created using field surveys and/or historic legacy data may suffer from compromised uniformity in landslide mapping. This is often reflected in the high concentration of landslides in the areas most visible from the roads (class 1 in Fig. 5), and in a pronounced monotonic decreasing trend, indicating significant heterogeneity in completeness. In fact, INGEMISA, Bornaetxea Field, and Bornaetxea News (Fig. 5a, d and e) exhibited the expected behaviour, suggesting that these inventories may be considered complete only within specific regions, such as those proximal to roads and urban areas, but not across the entire study area.
In the case of data collection for Lower Deba, a different methodology was employed, although a slight, but still monotonic, descending pattern can also be observed in Fig. 5c. This can be attributed to the inclusion of landslides identified during fieldwork near the dense road network prevalent in the study area. Furthermore, the presence of roads could be a contributing factor to the higher abundance of inventoried landslides in the more visible classes (Bruschi et al. 2013), potentially influenced by unstable cut slopes, altered water runoff patterns due to drainage works, and other factors (Taylor et al. 2020). Notably, a relevant percentage of the road network within the study area runs through the hillsides of the valleys, where the presence of steeper slopes is a clear conditioning factor for the occurrence of landslides. Thus, we cannot reject the hypothesis that a greater number of landslides actually occur in close proximity to roads in Gipuzkoa. Indeed, class 1 is the most prevalent column in Fig. 5 for all cases. Notwithstanding, the decreasing presence of landslides in classes 2-5 indicates that all inventories, except geoEuskadi, exhibit some bias toward areas visible from roads. Among them, Lower Deba displays a relatively lower impact of this bias, as evidenced by the considerable abundance of observations in less visible zones (classes 4 and 5). Conversely, the nonmonotonic distribution observed in the geoEuskadi inventory (Fig. 5b) is likely attributable to inaccurately located landslides rather than the inventory's completeness.
Concerning the distributions of explanatory variables among the inventories (see Sect. 6.3), we observe that INGEMISA, Lower Deba, and Bornaetxea Field show similar behavior regarding land use and slope. In the three inventories, crops and grassland types are dominant, and forests are among the less represented classes. This is likely due to the generalized underestimation of landslide occurrences in forested areas, which is most common in inventories obtained by means of traditional approaches (Brardinoni et al. 2003;Schmaltz et al. 2017). The comparable slope distributions also indicate that landslides can be expected in similar slope conditions, according to these databases. Furthermore, considering that Lower Deba and Bornaetxea Field collected exclusively shallow slide types of movements, they are consistent with the almost complete absence of rock outcrop, since these mass movements require, by definition, the presence of well-formed soil or altered regolith over the base rock (Cruden and Varnes 1996;Hungr et al. 2014).
Conversely, geoEuskadi and Bornaetxea News show relevant differences with respect to other inventories. The former presents rock outcrops as the most relevant class, while the latter locates most landslides in relation to anthropic land use. This behavior can be directly related to the data collection methods applied in each inventory. The large polygons mapped in geoEuskadi likely include large unstable bedrock walls outcropping within them. On the other hand, the information obtained through the newspaper chronicles for Bornaetxea News inevitably concentrates data around urbanized areas (Pereira et al. 2014;Taylor et al. 2015). Moreover, in densely urbanized areas and along roadsides, the spatial resolution of the digital elevation model utilized in this research may not be accurate enough to correctly consider the abrupt slope changes. This is probably the reason behind differences also observed in slope distributions that must necessarily result in significantly different susceptibility map. In addition, the generalized mismatch between inventories with respect to lithology indicates these databases fail to display, in fully representative manner, landslide occurrence in the Gipuzkoa province.
Taking into consideration the factors mentioned above, we made the decision to exclude Bornaetxea News from the susceptibility modelling experiment (Sect. 6.4). More specifically, we found three key reasons that influenced this choice: (i) based on the findings presented in Sects. 6.1-6.3, the accuracy of landslide locations in this inventory is not guaranteed, and there is a notable bias towards urbanized areas; (ii) the absence of information about landslide types in this inventory introduces uncertainty in the interpretation of the results; and (iii) the original purpose of this inventory was not specifically geared towards susceptibility modelling. In contrast, while there are valid arguments to also exclude geoEuskadi from the experiment, we decided to include this inventory due to its previous use in a landslide susceptibility modelling project (Diputación Foral de Gipuzkoa DFG 2007). It is nevertheless interesting to observe the susceptibility map that can result using this dataset, and to compare it with those derived from the other inventories, in a controlled experiment.
By observing Fig. 7, it becomes evident that in statistical landslide susceptibility modelling projects, if we only change the input landside inventory, while keeping all other parameters constant, the resulting maps exhibit significantly different spatial patterns; even if these inventories correspond to the same study area and were obtained for the same purpose. This matter highlights the extraordinary capacity of the statistical and machine learning approaches (such as the Logistic Regression) to find very tailored patterns that fit the input data (Goetz et al. 2015;Kavzoglu et al. 2019;Akinci and Zeybek 2021;Mehrabi 2022). We believe that this is a great opportunity for landslide susceptibility modelling in cases where there is a fully representative and reliable landslide inventory available. However, when the origin of the available data and its completeness are not verified, this may become a real drawback, as we have seen in this illustrative experiment.
As an example, in INGEMISA, Lower Deba, and Bornaetxea Field, data collection methods are well defined, and field verifications were performed in all cases. Therefore, we assume that location accuracy has been assessed, and thus the statistics obtained through these datasets present low spatial uncertainty. However, one of the reasons for the low performance of the model calibrated with INGEMISA (AUC = 0.68) could be the heterogeneity of the landslide typologies contained in the inventory, where 21% of the points represent landslides different to slide type (Table 1). Instead, in Lower Deba (AUC = 0.73) and Bornaetxea Field (AUC = 0.85), 100% of the dataset refers to the same landslide type; this specificity in the event typology probably facilitates the finding of common spatial features. In fact, correlation coefficients in Fig. 8 show a moderate agreement between these three maps, and their model fitting performance resulted in satisfactory AUC metrics, according to the conventions in literature (Zhao and Chen 2020). On the other hand, the susceptibility maps ultimately obtained reveal that the models calibrated with these inventories are heavily influenced by their inherent biases, such as a bias towards roads and urban areas (see Sect. 6.2). This results in relevant differences in the spatial distribution of the areas with extreme levels of landslide susceptibility, both very low and very high. Indeed, the highly polarized susceptibility map obtained with Bornaetxea Field inventory (Fig. 7d) probably indicates an overfitting situation, challenging the accuracy of this model if validated with future landslide occurrences.
The Lower Deba inventory represents the opposite scenario. The combination of photointerpretation and periodic field visits, together with the long observational time-span resulted in an abundant and accurate dataset, which ensures the detection of landslides terrain with limited ground visibility or in hardly accessible zones, so reducing the bias toward roads and urban area observed in the other inventories. The performance of the susceptibility model is also good (Fig. 8). Unfortunately, the small study area in which this inventory is applied is not fully representative of the entire Gipuzkoa province, and hence Lower Deba, by itself, is not applicable to our entire study area.
Although a larger variety of landslide types (only 22% of slide type; Table 1) is contained in geoEuskadi, we ignore the data collection methodology chosen for this inventory. Therefore, we cannot assume that the accuracy of the data location has been assessed by field verifications, for example, unlike the other three inventories. In addition, the location accuracy error obtained by considering the large polygons in this inventory as single centroid points represents another relevant issue that makes this database unsuitable for generating LSMs. We thus believe that too much uncertainty is implicit in the results obtained by the geoEuskadi inventory, even though its susceptibility map shows moderate to low prediction capacity (AUC = 0.65).
The data verification tests proposed in this research have been demonstrated to adapt correctly to any kind of landside inventory type, no matter the geometry or the collection methodology, and they allowed us to compare not only event-based inventories, but also historical and geomorphological inventories. These tests lead us to conclude that currently the only valid inventories for statistical susceptibility modelling in Gipuzkoa are INGEMISA, Lower Deba, and Bornaetxea Field. Furthermore, we believe that the Fig. 9 a LSMs obtained using LR modelling and calibrated with the combination of INGEMISA, Lower Deba, and Bornaetxea Field inventories. b ROC curve computed with the 30% random sample reserved for validation. c Empirical cumulative density function (ECDF) curve of the modelled susceptibility values for each inventory combined use of these three inventories could compensate for the inherent bias of each, resulting in a more homogeneous and representative susceptibility model. To verify this hypothesis, we conducted one last susceptibility model.
We extracted only the points corresponding to slide mass movements from INGEMISA and combined them with the Lowe Deba and Bornaetxea Field inventories. Subsequently, we conducted another logistic regression (LR) model following the exact settings applied in all previous experiments. The resulting combined landslide susceptibility map presented in Fig. 9, exhibited a prediction performance of 0.77 AUC, which can be considered good, albeit not optimal. Comparatively, the model calibrated with only the Bornaetxea Field inventory demonstrated better performance (Fig. 8). However, as noted in the literature (Reichenbach et al. 2018;Steger et al. 2021), high AUC values do not necessarily imply geomorphic causation in data-driven LSMs, and models with lower AUC may be more reliable and valuable from a geomorphological perspective. Moreover, the three inventories (individually) show a good agreement with the susceptibility map. The empirical cumulative density function (ECDF) curves displayed in Fig. 9c indicate that approximately 55% of the landslides present in each of the three inventories show susceptibility values greater than 0.5. Furthermore, this value increases up to approximately 70% when we reduce the threshold to a 0.4 of susceptibility. These results support the idea that the combined use of different inventories (once their suitability has been confirmed) can contribute to the generation of more representative susceptibility models.

Conclusions
The present work highlights the exceptional efforts made by different institutions and private individuals to compile sixteen landslide inventories (eleven in hard copy and five in electronic form) with different purposes within the Gipuzkoa province. The statistical analyses of five digitized landslide inventories-four of them used in the past to generate susceptibility maps-revealed remarkable differences between them. Such differences are mainly related to the characteristics of the single/multiple data sources used to obtain the information, in addition to the objectives of each database which condition the size of the study area, the temporal period covered, and the type of landslides studied, among other aspects. Notwithstanding, despite their methodological disparities, we observe a moderate spatial coincidence between some inventories and a great variability between others.
The results from our comparison tests (Sects. 6.1-6.3) indicate that the Bornaetxea News inventory does not offer a spatially representative picture of landslide occurrence in the study area. The spatial coincidence with respect to the other inventories analyzed is minimal, as shown in Figs. 5e and 6, and data are highly concentrated in urbanized areas.
Our results indicate the geoEuskadi inventory has minimal spatial accuracy, with absence of documentation detailing the precise methods followed to generate this dataset hindering our conclusions. A visual analysis of the inventory map (Fig. 3) and the results obtained from comparisons with other inventories suggest that geoEuskadi offers an extremely generalized landslide location pattern, and therefore it is hardly suitable for specific landslide analysis projects, such as susceptibility, hazard, or risk assessment.
Hence, out of the five landslide inventories digitally available in our study area, we conclude that only INGEMISA, Lower Deba, and Bornaetxea Field present enough guarantees to perform susceptibility analysis. However, we must acknowledge that specific limitations of each of these inventories are directly propagated into the final products, as shown in Fig. 7. Therefore, the combined use of all three databases offers more plausible and generalizable results according to the AUC and ECDF tests, and to our own judgement (Fig. 9).
The methods presented in this work allowed the comparison of heterogenous landslide databases and can be easily reproduced in other study areas where multiple landslide inventory sources exist. The preliminary assessment of such a relevant input should become a critical step in any landslide susceptibility modelling project.
In addition, we believe that the continued effort of landslide data collection is key to advance on the research of susceptibility, hazard, and risk models; and a cost-benefit approach to this laborious activity should always be a priority, especially for large regions such as the Gipuzkoa province. In Appendix II, we propose some general recommendations for the compilation of landslide inventories that will hopefully be useful to motivate future landslide mapping projects in either the Gipuzkoa province, or any other national or foreign territory.

Appendix I
Description of the explanatory variables used in the experiments mentioned in Sects. 5.3 and 5.4. See Fig. 10.

Land use
It is the corresponding part to the study area of the Spanish National Forest Inventory of the 2010 (for its acronym in Spanish, IFN version (4), which is an update of IFN version (3), carried out in 2005. The update was conducted by aerial photo-interpretation of images with 25 cm of spatial resolution obtained in 2009. The original 32 classes were re-classified by expert criteria as it is shown in Bornaetxea (2018). The source layer can be found under the code INV_FORESTAL_2010_10000_ETRS89 in the geoportal geoEuskadi (IDE geoEuskadi 2022).

Lithology
Lithology, or rock typology, constitutes the geological substrate of the territory classified according to its composition and its physical-chemical behavior. The original layer, available in the geoportal geoEuskadi (IDE geoEuskadi 2022) and named CT_LITOLOG-ICO_25000_ETRS89, was undertaken in 1999 based on the contributions of the Geological Map of Euskadi developed by the EVE (Ente Vasco de Energía). According to this map, 21 different rock typologies can be found in study area. Knowing that the resistance against the shear tension, as well as the water infiltration capacity could be considered similar in some of those lithologies, the original 21 classes were re-classified into 6 classes following the expert criteria specified in Bornaetxea (2018).

Slope
The slope was digitally calculated as a derived product of the DEM. By means of Q-GIS software, the maximum rate of change in value from each cell to its neighbours was calculated, assuming that the maximum change in elevation over the distance between the cell and its eight neighbours identified the steepest downhill descent from the cell. For this reason, the calculations were carried out over the original DEM, and then the resulting layer was clipped in order to fit the study area. This way, border errors were avoided. This variable represents the existing angle between the terrain surface and the horizontal plane in degrees, and it shows a direct relation with the tangent and normal cutting stress of the surface formation. In addition, the slope influences very importantly the water flow velocity and distribution.
sponding to each observation must be properly registered to facilitate the decision of its use by future potential users. 6. The concentration of records in the vicinity of the road network is attributed to both natural and anthropogenic causes, but this can also be magnified depending on the type of data source used or the subsequent methodology. Both overrepresentation and underrepresentation of landslide records in these areas can affect to the calibration of landslide susceptibility models. The initial exclusion of presumed human-induced landslides from data collection is discouraged; their use can be discarded in subsequent phases of the modelling process.