GIS-Based Landslide Susceptibility Modelling in Urbanized Areas: A Case Study of the Tri-City Area of Poland

: This paper presents the results of landslide prediction modelling for young glacial areas performed using statistical methods. The area in question is urbanized and therefore mass wasting activity is a matter of interest to both the local community and the authorities. The analysis was based on the 2011 ‘Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdansk’ and the 2012 incomplete ‘Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdynia’. The research took into account geological, geomorphological, hydrological, hydrogeological, and anthropogenic conditions. The landslide susceptibility map was created using the statistical landslide index. The calculated indices were used to create a map of Gdansk’s landslide susceptibility. In Gdansk, 84.50% of the total diagnosed landslide area belongs to the high susceptibility class, 14.25% to the moderate susceptibility class, and only 1.25% to the low or very low susceptibility class. After extrapolation, the data was also used to create a susceptibility map for the remaining parts of the Tri-City area, Sopot and Gdynia. The difﬁculty of extrapolating landslide data for neighboring urban areas was indicated. In Gdansk, which had been covered by geological mapping, the best modelling results were obtained with a large number of causal factors. In Gdynia and Sopot, for which the statistical landslide index value was extrapolated from Gdansk, the best results were obtained when selected causal factors were considered. In Sopot and Gdynia, 81.6% of the landslide area belongs to the high susceptibility class, 15.1% to the moderate class, and 3.3% to the low susceptibility class. These results emphasize a different role of some causal factor classes in the occurrence of landslides in neighboring urban areas. The resultant maps show the areas in which mass wasting is the most probable in the future.


Introduction
Mass wasting processes, controlled by gravity, hydrologic response, and strong-ground motion conditions, are important geodynamic hazards that can significantly affect local economies. The damaging nature of landslides is particularly dangerous in densely populated urban areas. The Tri-City area of Gdansk, Gdynia, and Sopot is an example of such an area, located in a terrain with very diverse topography. The negative impacts of mass wasting may lead to enormous damage to property. In 2001 alone, after violent rains, damage caused by landslide activation in Gdansk was estimated at EUR 780,000 [1]. Currently, many landslides in Gdansk are still active, as confirmed by terrestrial laser scanning techniques monitoring of selected shallow landslides carried out in 2014-2022 [2]. However, awareness of these types of hazards in the Tri-City area is low, as comprehensive geological research on the area's landslides has been carried out only recently [1][2][3][4][5][6][7].
Risk management involves identifying the extent of risk, introducing appropriate preventative methods to mitigate the risk, and decreasing the possibility of the harmful incidence or minimizing the vulnerability or exposure to the risk factors [8]. One of the many methods of managing landslide risk, commonly used worldwide, is to determine GeoHazards 2022, 3 509 locations susceptible to mass wasting using geographical information systems (GIS) tools and statistical methods. During the last decades, the ever-growing availability of earth observation data and the well-established use of GIS has led to the development of automated workflows based on GIS toolboxes [9,10].
Landslide susceptibility maps are an instrument used to support land management. Landslide susceptibility assessment is a valuable tool for land use planning and reducing the costs of landslides [11,12]. Landslide susceptibility assessment can also be considered as an important tool in early warning system techniques, as well as the initial step towards the landslide hazard and risk assessment [13]. The literature provides many examples of GIS being harnessed in mass wasting research (e.g., References [14][15][16][17][18]). Many diverse techniques are used to study landslide susceptibility, including the empirical, statistical, and deterministic approach used in large-scale engineering geology projects [15,19].
In Poland, geoinformation analyses of landslides have been performed only since the early 2000s, mainly in the mountainous areas of southern Poland [20][21][22][23] and for the remaining part of Poland in the areas characterized by intensive human activity [6,7,24]. The importance of landslide susceptibility assessments for environmental policy and decisionmaking processes in Europe is set forth in the framework of the European Union's Thematic Strategy for Soil Protection [25]. The relevance of landslide zoning through spatial susceptibility assessments in Europe is additionally recognized in the European Commission's approach to natural and man-made disaster prevention [26]. Additionally, landslides are among the natural hazards that the European Union considers particularly important and had a high priority in the EU forest strategy for 2030 [27]. Therefore, many studies evaluate landslide susceptibility in the context of multi-hazard assessments [18].
This paper presents the results of shallow landslide prediction modelling in the urbanized Tri-City area of Poland using the statistical landslide index (SLI) method. The bivariate SLI method has been used by various researchers [6,7,17,18,[28][29][30][31]. So far, however, whether the susceptibility indices developed by this method would be applicable in neighboring areas has not been validated. This article presents the results of the extrapolation and, in addition, a comparison of the two susceptibility maps for the cities of Gdansk [6] and Gdynia [32] produced using the SLI method. The article also compares the role of individual passive factors for these two neighboring areas.
This type of extrapolation of data for neighboring areas has not been discussed extensively in the literature. According to Smith et al. [33], spatial analysis is a process that in many instances follows a number of (often iterative) stages. Therefore, some researchers carry out multiple iterations of susceptibility modelling of the same area to improve prediction performance and comparison [34,35].
The primary source of information is the 2011 'Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdansk'. In all, 96 landslides were identified in Gdansk at that time. They account for 0.12% of the city's area (Jurys et al., 2014) [1]. The results of the landslide inventory in Gdansk were then extrapolated to the neighboring cities Gdynia and Sopot. The extrapolation results were verified by landslides occurring in 2012 for selected areas of Gdynia. Extrapolation for Gdynia and Sopot was carried out in 2013, but the results have not been published so far. In 2015, a complete landslide inventory was prepared for Gdynia [5], which was then updated as part of the author's doctoral thesis [32]. On the basis of this research, susceptibility maps were produced for Gdynia [7,32].

Location and General Features of the Study Area
The Tri-City consists of three metropolitan centres: Gdansk, Gdynia, and Sopot. The area is greatly diversified in terms of its geomorphology and geology ( Figure 1). altitude ranges from −8.5 to 205.23 m above sea level. The very location of the Tri-City, characterized by the hypsometric diversity of the terrain surface, affects the high variability of atmospheric precipitation. According to the Polish Climate Atlas, normal precipitation in the time period between 1971 and 2000 ranged from about 550 mm in the area of the coastal zone to about 700 mm in the eastern part of the moraine plateau of the Kashubian Lakeland [36]. The mean annual maximum air temperature observed in the period between 1951 and 1998 was 11 °C [37]. As late as the 19th century, buildings in Gdansk and Sopot were located on flat coastal and river areas [38][39][40]. In the Tri-City, for the last few decades, the slopes of uplands and the slopes of valleys have also been included in the development [1]. In the Tri-City between 1997 and 2016 the greatest changes took place in arable land, which altogether decreased by 17%, and forest land (13% in total). All this was to the benefit of grassland (which increased by 23%) and built-up areas (increasing by 10% in total) (Kwoczyńska, 2021) [41]. The land cover change caused an increase in the landslide occurrence, causing more severe property damage (Małka, 2018) [32].
In terms of geomorphology, the study area is mainly located within the moraine plateau ( Figure 1). It was formed during the Weichselian glaciation as a result of the accumulation of moraine deposits, i.e., tills, sands, and gravels. During deglaciation, the moraine The local climate is influenced by elevation and the distance from the Baltic Sea. The altitude ranges from −8.5 to 205.23 m above sea level. The very location of the Tri-City, characterized by the hypsometric diversity of the terrain surface, affects the high variability of atmospheric precipitation. According to the Polish Climate Atlas, normal precipitation in the time period between 1971 and 2000 ranged from about 550 mm in the area of the coastal zone to about 700 mm in the eastern part of the moraine plateau of the Kashubian Lakeland [36]. The mean annual maximum air temperature observed in the period between 1951 and 1998 was 11 • C [37].
As late as the 19th century, buildings in Gdansk and Sopot were located on flat coastal and river areas [38][39][40]. In the Tri-City, for the last few decades, the slopes of uplands and the slopes of valleys have also been included in the development [1]. In the Tri-City between 1997 and 2016 the greatest changes took place in arable land, which altogether decreased by 17%, and forest land (13% in total). All this was to the benefit of grassland (which increased by 23%) and built-up areas (increasing by 10% in total) [41]. The land cover change caused an increase in the landslide occurrence, causing more severe property damage [32].
In terms of geomorphology, the study area is mainly located within the moraine plateau ( Figure 1). It was formed during the Weichselian glaciation as a result of the accumulation of moraine deposits, i.e., tills, sands, and gravels. During deglaciation, the moraine plateau was incised by a glacio-marginal outflow system [42]. The edges of the moraine plateau were cut by erosional valleys [43]. The Tri-City's western upland, located within the boundaries of the peripheral part of the Kashubian Lakeland which developed in the Pleistocene, is the most prone to landslides. Furthermore, also threatened by landslides are the slopes of the Kashubian Coast upland, covering the eastern part of Sopot and Gdynia, along with a relatively small north-western part of Gdansk [44]. The eastern boundary of the Kashubian Coast comes in direct contact with the Gulf of Gdansk to create a cliff shoreline, which is the most susceptible to landslides of all of the Tri-City areas (Figure 2a-d). plateau was incised by a glacio-marginal outflow system (Mojski 2005) [42]. The edges of the moraine plateau were cut by erosional valleys (Woźniak 2014) [43]. The Tri-City's western upland, located within the boundaries of the peripheral part of the Kashubian Lakeland which developed in the Pleistocene, is the most prone to landslides. Furthermore, also threatened by landslides are the slopes of the Kashubian Coast upland, covering the eastern part of Sopot and Gdynia, along with a relatively small north-western part of Gdańsk (Kondracki 1998) [44]. The eastern boundary of the Kashubian Coast comes in direct contact with the Gulf of Gdańsk to create a cliff shoreline, which is the most susceptible to landslides of all of the Tri-City areas (Figure 2a-d). The mass wasting which occurs there is initiated and aggravated by abrasion. The cliffs began to form approximately 6000 years ago and their general shape remains similar to this day (Uścinowicz 2003;Gałka et al., 2018) [46,47]. The earliest information on threats related to mass movements in the area of present-day Gdynia dates from the end of the 16th century and relates to cliff coasts (Małka 2019) [48]. The first geological survey for these cliffs was carried out in the early 20th century (Sonntag 1910; Figure 2a,b) [45].
At present, the development of Baltic cliff shorelines is the most influenced by the eustatic factor, with its activity regionally modified by the neotectonic factor with minimal influence from glacial isostatic movements (Zawadzka-Kahlau 1999) [49]. In the Tri-City area, negative (subsiding) vertical movements of the Earth crust surface were detected at a value of (−1.1) mm/year, to additionally increase the relative growth of sea level in the area (Wyrzykowski 1971) [50].
From a geological point of view, the Tri-City area features Quaternary formations of significant thickness, reaching as much as 145 m [51], with Neogene lacustrine formations locally exposed in situ and in the form of glacial float [45,[51][52][53]. For the most part, the The mass wasting which occurs there is initiated and aggravated by abrasion. The cliffs began to form approximately 6000 years ago and their general shape remains similar to this day [46,47]. The earliest information on threats related to mass movements in the area of present-day Gdynia dates from the end of the 16th century and relates to cliff coasts (Małka 2019) [48]. The first geological survey for these cliffs was carried out in the early 20th century (Figure 2a,b) [45].
At present, the development of Baltic cliff shorelines is the most influenced by the eustatic factor, with its activity regionally modified by the neotectonic factor with minimal influence from glacial isostatic movements [49]. In the Tri-City area, negative (subsiding) vertical movements of the Earth crust surface were detected at a value of (−1.1) mm/year, to additionally increase the relative growth of sea level in the area [50].
From a geological point of view, the Tri-City area features Quaternary formations of significant thickness, reaching as much as 145 m [51], with Neogene lacustrine formations locally exposed in situ and in the form of glacial float [45,[51][52][53]. For the most part, the Tri-City area is covered by sandy-gravelly fluvioglacial sediments and tills of the morainic upland, with silty-clayey marginal lake sediments occurring in places [54]. The lowered terrain is filled with diverse Holocene sedimentary facies [52,53]. In the Tri-City's urbanized sections, a clear relationship can be observed between the geological structure and its resultant land morphology and the occurrence of shallow landslides. Pronounced drops in terrain, varied slope angles, and the conducive geological structure produce good conditions for landslides. The presence of non-homogenous glacial, marginal lakes and fluvioglacial sediments, as well as multiple glaciotectonic deformations, is highly significant. Moreover, an important role is played by formations with aleurite grain size, which are prone to being washed away and liquefied. They are exposed at the fringe of the Kashubian Lakeland's morainic upland [1].
Land development methods and human activity important factors initiating the occurrence and development of shallow landslides in the Tri-City. Many present-day landslides originated from anthropogenic changes in land use. Undercutting, over-steeping of slopes, and construction burden, coupled with faulty water and sewage management, have an especially negative influence [1]. The problem of landslides occurring in the urbanized areas of the Tri-City has been known since the second half of the 20th century and is still ongoing, causing significant damage to infrastructure (Figure 3a-d).
In the Tri-City's urbanized sections, a clear relationship can be observed between the geological structure and its resultant land morphology and the occurrence of shallow landslides. Pronounced drops in terrain, varied slope angles, and the conducive geological structure produce good conditions for landslides. The presence of non-homogenous glacial, marginal lakes and fluvioglacial sediments, as well as multiple glaciotectonic deformations, is highly significant. Moreover, an important role is played by formations with aleurite grain size, which are prone to being washed away and liquefied. They are exposed at the fringe of the Kashubian Lakeland's morainic upland (

Input Data
The research material for landslide susceptibility maps of the Tri-City area was provided by Poland's National Geodetic and Cartographic Resource, with additional source data from the Polish Geological Institute, National Research Institute, and public domain materials ( Table 1). The following data have been used for the purpose of the GIS analysis: the digital elevation model (DEM) from airborne laser scanning (ALS), also described as light detection and ranging data (LIDAR data) [55] from the Information System of Country Protection Against Extraordinary Hazards (ISOK) project; the Polish Topographic object database (BDOT database); maps of sediment type at a depth of 1 and 4 m b.g.l. from the 'Geological and engineering atlas of the Tri-City of Gdansk-Sopot-Gdynia Urban Area' [54]; a map of the groundwater table level; the Urban Atlas database [56]; the 'Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdansk' [3]; and the 'Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdynia' [4]. The accuracy of the layers used in the analysis corresponds to the cartometric accuracy of a 1:10,000 map. The resolution of all the raster maps used in the analysis is 5 × 5 m.
Thematic layers were collected, sorted, and prepared to showcase the causal factors. The research considered geological, geomorphological, hydrological, hydrogeological, and anthropogenic conditions. Selected causal factors visualized in the form of maps were used to model the landslide susceptibility of Gdansk and Gdynia (Table 2) [6,7]. Slope aspect has an influence on solar radiation intensity, which controls the local temperature and evaporation [58]. The slope aspect also affects soil moisture, rainfall intensity, and vegetation density on the hillsides [29]. River erosion and sea abrasion as well as drainage density determine the likelihood of a landslide in young morainic areas [6,7]. Close river proximity induces instability at the foot of a slope due to the concentration of groundwater flow; further, river and sea undercutting also have a destabilizing effect [19].
9 ground type at a depth of 1 m b.g.l. (G 1)

Digital engineering geological data
Influence on the weathering of the bedrock and the resulting soils [62] 10 ground type at a depth of 4 m b.g.l.

Digital engineering geological data
Deep geological structure has an impact on landslides in young glacial areas as it illustrates changes in permeability of sediments and potential slip zones [7].

groundwater table level (GW) Digital engineering geological data
Well-sorted, fine-grained deposits in which the groundwater table is within 3 m may be subject to possible severe liquefaction [63].
12 elevation (E) * DEM Elevation relates to landslides trigger, i.e., heavy rainfall [7] 13 profile curvature (P) * DEM Divergent or convergent curvature has an influence on the water run-off concentration [64] * Factor used only in landslide susceptibility mapping of Gdynia.
The selection of geospatial data for the Tri-City area was determined by relevance, reliability, and compatibility of scales. Although the Tri-City has a substantial amount of geospatial data, many of these sources (e.g., detailed geological maps of Poland, geoengineering data) were rejected in the analysis as they did not meet the above-mentioned criteria.
The five selected factors are shown in Figure 4. GeoHazards 2022, 3, x FOR PEER REVIEW 8 of 21

Research Methodology
The landslide inventory was carried out with a scale of 1:10,000, in the framework of the LCS (SOPO) project, using standardized methods (Grabowski et al., 2008) [65]. The "Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdańsk" (Jurys et al., 2011) [3] was the basic source of information for the analysis and data validation of the susceptibility map of Gdańsk.
The analysis considered the western part of the city of Gdansk, as the eastern part is flat and not affected by landslides ( Figure 1). As Gdańsk only has a small number of active landslides, all the landslides (96) were included to achieve reliable geoinformation analysis using statistical methods. For susceptibility modelling, the entire landslide area (landslide body) was considered.
Due to the already-performed comprehensive geological charting of Gdańsk area landslides [1,3], a statistical analysis was performed for this area first. The set of inventoried landslides in Gdańsk was divided into two spatially diversified subsets of 66 and 30 (Małka 2015) [6]. The larger subset and parameter maps were used to develop the model, while the second subset was used to verify it (Małka 2015) [6]. In modelling the landslide susceptibility of the area under research, 11 causal factors visualized as maps were used (Małka 2015, Table 2) [6]. The topographic attributes, computed directly from the digital elevation model, included the slope angle, slope aspect, relief (computed in the preset vicinity of 5 m from each raster cell), topographic wetness index (TWI), stream power index (SPI), and the LS sediment transport capacity index. The primary and secondary

Research Methodology
The landslide inventory was carried out with a scale of 1:10,000, in the framework of the LCS (SOPO) project, using standardized methods (Grabowski et al., 2008) [65]. The 'Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdansk' [3] was the basic source of information for the analysis and data validation of the susceptibility map of Gdansk.
The analysis considered the western part of the city of Gdansk, as the eastern part is flat and not affected by landslides ( Figure 1). As Gdansk only has a small number of active landslides, all the landslides (96) were included to achieve reliable geoinformation analysis using statistical methods. For susceptibility modelling, the entire landslide area (landslide body) was considered.
Due to the already-performed comprehensive geological charting of Gdansk area landslides [1,3], a statistical analysis was performed for this area first. The set of inventoried landslides in Gdansk was divided into two spatially diversified subsets of 66 and 30 [6]. The larger subset and parameter maps were used to develop the model, while the second subset was used to verify it [6]. In modelling the landslide susceptibility of the area under research, 11 causal factors visualized as maps were used (Table 2) [6]. The topographic attributes, computed directly from the digital elevation model, included the slope angle, slope aspect, relief (computed in the preset vicinity of 5 m from each raster cell), topographic wetness index (TWI), stream power index (SPI), and the LS sediment transport capacity index. The primary and secondary topographic attributes were calculated directly from the high-resolution DEM obtained by airborne laser scanning (ALS) from 2014. The point density in the ALS-based point cloud was 12 points/m 2 . The Jenks natural break algorithm was used for grouping slope angles, elevation, relief energy, TWI, SPI, and LS [6,[66][67][68]. The remaining parameter layers include: land use and land cover, distance from watercourses and reservoirs, maps of ground type at a depth of 1 and 4 m b.g.l., and the groundwater table level, computed on the basis of borehole data from geological and engineering documentation by means of Euclidean allocation (Table 2) [6].
The data obtained during field work (MOTZ map) and the maps which visualize the environmental factors were used to calculate the landslide susceptibility indexes by means of the landslide susceptibility [14,17,31]. The statistical landslide index (SLI) method used in the analysis is based on the following formula: ln W i -the weight (landslide susceptibility index) given to a certain parameter class; Densclas-the landslide density within the parameter class; Densmap-the landslide density within the entire map; N pix(Si)-the number of pixels which contain landslides in a certain parameter class; N pix(Ni)-the total number of pixels in a certain parameter class.
The main premise of this method is the 'crossing' of a landslide map with a certain parameter map. The calculated landslide susceptibility indexes were used to create a landslide susceptibility map of Gdansk with a scale of 1:10,000 [6].
After extrapolation, the data was also used to create a susceptibility map with a scale of 1:10,000 for Sopot and Gdynia, the remaining areas of the Tri-City. To assign the SLI values to individual classes the intervals were selected by the Jenks (1967) [66] natural break algorithm. Determined in this way, the intervals made it possible to obtain the model's optimum predictive capacity, while at the same time narrowing down the susceptible areas, which is of great significance to the Tri-City.
The susceptibility map produced for Sopot and Gdynia was validated with the incomplete 'Register of landslides and areas prone to mass movements with a scale of 1:10,000 for the city of Gdynia' [4].
The article also compares the results of the susceptibility assessment for Gdansk and Gdynia, performed using the same methods for the entire landslide area. Comparison of the two susceptibility maps for neighboring urbanized areas was performed on the basis of the landslide susceptibility map of Gdansk, developed using the landslide susceptibility index (LSI) and eleven parameters ( Table 2) [6], and the landslide susceptibility map of Gdynia was created using LSI and thirteen parameters ( Table 2) [32].
The calculated weights (ln W i ) allowed the graphical representation of the role of the individual causal factors. The graphs were developed using Microsoft Excel 2010; for the sake of readability, the description of the individual factor classes was omitted in the graphs.
The spatial analysis of the landslide susceptibility maps of the Tri-City was performed using ArcGIS 10.2 software (https://www.esri.com (accessed on 23 August 2022)). Global Mapper 14.2 software (https://www.bluemarblegeo.com (accessed on 23 August 2022)) was used in an auxiliary capacity for data conversion and to change the resolution of the DEM data.

Results
Landslide susceptibility is determined by considering the interplay between landslide occurrences and impacts of factors that induce instability [7]. In this analysis, the use of the statistical landslide index and eleven thematic layers has revealed the importance of many environmental factors. Various combinations of instability factors were used in the modelling. Figure 5 presents the landslide susceptibility map of Gdansk with a scale of 1:10,000 developed using the SLI method, eleven variables, and all landslides. the importance of many environmental factors. Various combinations of instability factors were used in the modelling. Figure 5 presents the landslide susceptibility map of Gdańsk with a scale of 1:10,000 developed using the SLI method, eleven variables, and all landslides. In the next stage, the LSI index in Gdańsk was used to create a landslide susceptibility map for the neighboring counties of Sopot and Gdynia, which have a similar geological structure and equivalent geomorphological conditions. Not all the layers computed for the Gdańsk area were applicable to the rest of the Tri-City. The limitation was the different scopes of the relief map and the LS sediment transport capacity index map for Gdańsk and the rest of the Tri-City. For this reason, these data were not included in the secondary analysis. Greater slope angle values were found in the area of Sopot and Gdynia (from 0 to 63.4°) compared to Gdańsk (from 0 to 48.1°). However, there was a growing dependence between the values of the lnWi index and the slope angle in Gdańsk (Table 3), which is why the slope angle was applicable in the predictive model. In the next stage, the LSI index in Gdansk was used to create a landslide susceptibility map for the neighboring counties of Sopot and Gdynia, which have a similar geological structure and equivalent geomorphological conditions. Not all the layers computed for the Gdansk area were applicable to the rest of the Tri-City. The limitation was the different scopes of the relief map and the LS sediment transport capacity index map for Gdansk and the rest of the Tri-City. For this reason, these data were not included in the secondary analysis. Greater slope angle values were found in the area of Sopot and Gdynia (from 0 to 63.4 • ) compared to Gdansk (from 0 to 48.1 • ). However, there was a growing dependence between the values of the lnWi index and the slope angle in Gdansk (Table 3), which is why the slope angle was applicable in the predictive model. Table 3. Spatial relationship between the selected important conditioning factors and all landslides using the SLI for the entire landslide area in Gdynia and Gdansk. The LSI values (lnWi) were calculated on the basis of the updated landslide inventory for Gdynia from 2018 [32] and for Gdansk on the basis of the complete landslide inventory from 2011 [3]; values with different signs (plus/minus) are highlighted in yellow.  At the next stage, similar causal factor maps as for Gdansk were developed for the remaining part of the Tri-City (Gdynia and Sopot). Computed for the area of Gdansk by means of the SLI method, the landslide susceptibility indexes (lnWi) were then extrapolated onto the areas of Sopot and Gdynia. As a result of multiple iterations, a number of landslide susceptibility models were produced with the use of various parameter layers.
Five parameter layers were used in the final modelling of landslide susceptibility in Sopot and Gdynia: slope angle, slope aspect, land use, distance from watercourses and reservoirs, and the groundwater table level ( Figure 6). The other parameter layers were not included in analysis, as they did not enhance the model's predictive capacity.
As a result, a landslide susceptibility map with a scale of 1:10,000 was produced for Sopot and Gdynia. Figure 7 presents the landslide susceptibility map of Sopot and Gdynia with a scale of 1:10,000, developed using the SLI method and five variables.
Based on the statistical analysis, a spatial layout of landslide-susceptible areas in the Tri-City was produced. The landslide susceptibility maps were verified in the case of Gdansk by means of the second landslide subset [6] and in in the case of Gdynia by means of the landslides documented during the incomplete 2012 cartographic field work for selected areas of Gdynia. Consequently, it was found that 84.50% of the landslide area diagnosed in Gdansk belongs to the high susceptibility class, 14.25% to the moderate susceptibility class, and only 1.25% to the low or very low susceptibility class [6]. High and very high landslide susceptibility classes are found mainly in the north-east (the districts of Oliwa, VII Dwór, Brętowo, Wrzeszcz Górny) and in the east (Orunia-Św. Wojciech-Lipce) of Gdansk ( Figure 5). Based on the statistical analysis, a spatial layout of landslide-susceptible areas in the Tri-City was produced. The landslide susceptibility maps were verified in the case of Gdańsk by means of the second landslide subset (Małka 2015) [6] and in in the case of Gdynia by means of the landslides documented during the incomplete 2012 cartographic field work for selected areas of Gdynia. Consequently, it was found that 84.50% of the total area diagnosed in Gdańsk belongs to the high susceptibility class, 14.25% to the moderate susceptibility class, and only 1.25% to the low or very low susceptibility class (Małka 2015) [6]. High and very high landslide susceptibility classes are found mainly in the north-east (the districts of Oliwa, VII Dwór, Brętowo, Wrzeszcz Górny) and in the east (Orunia-Św. Wojciech-Lipce) of Gdańsk ( Figure 5). In Sopot and Gdynia, 81.6% of the landslide area belongs to the high susceptibility class, 15.1% to the moderate class, and 3.3% to the low susceptibility class. High landslide susceptibility is related to the cliff in the Oksywie and Redłowo Hillocks (Kępa Oksywska, Kępa Redłowska) and the area of the Tri-City Landscape Park.
Comparison of the two susceptibility maps for neighboring urbanized areas was performed on the basis of the author's previous research [6,7,32]. The application of the SLI method and many variables (11 for Gdansk and 13 for Gdynia) allowed for the determination of the importance of predisposing factors for Gdansk and Gdynia (Table 2) [6,7,32].
In the modelling of the landslide susceptibility of Gdansk using the SLI method and many causal factors, the most significant geo-environmental variables were the topographic attributes such as slope angle, relief energy, and LS; noteworthy also were TWI and SPI (Figure 8). Similar causal factors had the greatest impact on landslide formation in Gdynia ( Figure 9) [7,32]. Furthermore, analogous classes had an impact on landslide formation in both cities. However, there was some differences in the individual classes of these factors for the two cities (Figures 8 and 9; Table 2). In Sopot and Gdynia, 81.6% of the landslide area belongs to the high susceptibility class, 15.1% to the moderate class, and 3.3% to the low susceptibility class. High landslide susceptibility is related to the cliff in the Oksywie and Redłowo Hillocks (Kępa Oksywska, Kępa Redłowska) and the area of the Tri-City Landscape Park.
Comparison of the two susceptibility maps for neighboring urbanized areas was performed on the basis of the author's previous research [6,7,32]. The application of the SLI method and many variables (11 for Gdańsk and 13 for Gdynia) allowed for the The highest (>2) landslide susceptibility index values were determined for the LS index, slope angle, and relief energy. Most landslide areas occur at the angles of 10-29 • . The value of the landslide susceptibility index increases with the greater angles both in Gdansk and Gdynia. Additionally, the interdependence between the increase in LSI index value and the intensities of slope failures is evident. Geological structure, especially the lithological diversity of strata, is an important triggering factor for mass wasting processes in Gdansk and Gdynia. The LSI index values show that landslides develop in loose (permeable) sediments (in Gdansk, for soils lnWi = 0.46, for sand and gravels lnWi = 0.04; in Gdynia, for soils 1.38) that occur above cohesive (impermeable/poorly permeable) ones (in Gdansk, for glacial till lnW i = 0.06; in Gdynia-for glacial till lnWi = 0.17). Hydrological and hydrogeological conditions are the next factors controlling mass wasting in Gdansk and Gdynia. The slopes in the vicinity (up to 50 m) of watercourses and reservoirs (in Gdansk, lnWi = 0.84, in Gdynia, lnWi = 0.88) are especially at risk. In Gdansk, the SLI reach high values for the TWI interval of −0.6-2.6 (lnWi = 1.84). In Gdynia, the highest values of LSI (lnWi = 1.12) were observed for the TWI interval of 0-3, accounting for 43.33% of landslides [7]. This indicates the link between predominantly shallow mass movements and soil moisture parameters in Gdansk and Gdynia. A shallow groundwater table is frequently associated with the surface of rupture and demonstrates its role in landslide susceptibility (in Gdansk, for the groundwater table interval of 1-2 m lnWi = 0.24; in Gdynia, for the groundwater table interval of 0.5-1 m lnWi = 0.88).  Table 2) [6,7,32].
In the modelling of the landslide susceptibility of Gdańsk using the SLI method and many causal factors, the most significant geo-environmental variables were the topographic attributes such as slope angle, relief energy, and LS; noteworthy also were TWI and SPI (Figure 8). Similar causal factors had the greatest impact on landslide formation in Gdynia (Małka 2018(Małka , 2021; Figure 9) [7,32]. Furthermore, analogous classes had an impact on landslide formation in both cities. However, there was some differences in the individual classes of these factors for the two cities (Figures 8 and 9; Table 2).    Table 2) [6,7,32].
In the modelling of the landslide susceptibility of Gdańsk using the SLI method and many causal factors, the most significant geo-environmental variables were the topographic attributes such as slope angle, relief energy, and LS; noteworthy also were TWI and SPI (Figure 8). Similar causal factors had the greatest impact on landslide formation in Gdynia (Małka 2018(Małka , 2021; Figure 9) [7,32]. Furthermore, analogous classes had an impact on landslide formation in both cities. However, there was some differences in the individual classes of these factors for the two cities (Figures 8 and 9; Table 2).

Discussion
Susceptibility maps were made for both Gdansk [6] and Gdynia [7,32] on the basis of a detailed landslide inventory [3,4,32,36]. For the last 16 years, inventories of landslides in Poland have been compiled under a nationwide project, the Landslide Counteracting System [65,69]. However, these inventories are very expensive and are not yet produced for the entire area of Poland. It is therefore appropriate to consider to what extent landslide data from neighboring areas can be used for areas for which we do not have such information.
However, the landslide susceptibility map prepared for Gdynia and Sopot on the basis of data extrapolated from Gdansk can only be indicative. Repeated modelling of susceptibility for Gdynia and Sopot using different variants of eleven causal factors and extrapolated landslide susceptibility index values for Gdansk and validation showed that such a wide set of causal factors does not produce reliable results. After multiple iterations, limiting the analysis to the five selected factors (Figure 6) yielded the best results, as shown by the validation results with the landslides in Gdynia, documented in 2012. On the other hand, the complete inventory of landslides for Gdynia and the susceptibility assessment performed on this basis using the landslide susceptibility index clearly show that sometimes different classes of causal factors are of different importance in the formation of landslides in Gdansk and Gdynia (Table 3). Particular attention should be paid to classes where this significance is quite different. The northern slope aspect is susceptible to landslides in Gdansk (lnWi = 0.43), while in Gdynia it is an unfavorable factor (lnWi = −0.17). On the other hand, the eastern slope aspect in Gdynia is the most favorable factor (lnWi = 0.49), whereas in Gdansk it is an unfavorable factor (lnWi = −0.10). In the upland area of Gdansk, landslides occur mainly along the entire slope lengths or at their bottom parts, with north, north-west and north-east aspects. The north-facing slopes retain snow longer than the south-oriented slopes. They also have an increased infiltration of precipitation and snowmelt. In the Tri-City, valley slopes with a northern orientation show higher relative humidity values compared to the southern ones, sometimes up to 3% [70]. Increased infiltration of the northern slopes may affect soil moisture parameters and the liquefaction of fine-grained and silty sediments. As a result, shallow mass movements may be activated. In the case of the slope aspect of Gdynia, eastern slopes (lnWi = 0.49) are the most predisposed to landslides, which should be associated with the orientation of coastline (Figure 1; Table 3). In the young morainic area of Gdynia, coastal abrasion is crucial for the reactivation of mass movement.
The landslide susceptibility index in Gdansk reaches higher values in forests (lnWi = 0.68; Table 3). The high landslide susceptibility of forested areas should be treated as a consequence, rather than a cause. Forests are a slope stabilizing factor as vegetation, through interception, retains some of the water soaked up by the earth, while capillary phenomena prevent the soil's saturation with water, on top of which deep roots contribute to soil drainage and lower pore pressure [71]. The higher landslide susceptibility of forested areas determined by the analysis is related to anthropogenic land use. Most of the Tri-City forests are located in areas with steep terrain and high relief energy (Figures 4 and 6). Centuries of human impact on the environment have caused the deforestation of theŻuławy Lowland and the Vistula Spit, while the hypsometrically diversified fringes of the upland have remained largely undeveloped.
In the case of the Tri-City area, cutting down forests on varied lithology slopes and in the areas susceptible to being washed away is certainly a landslide activating factor. In Gdynia, forested areas are not predisposed to landslides, (ln W = −0.13; Table 3). Additionally, the cutting down of a forest in Sopot Górny (Upper Sopot) in connection with the planned construction of a ski jump caused the occurrence of a new landslide. The landslide developed in the upper part of the slope with tills and multiple Miocene sandy-clayey xenoliths with lignite, and deeper situated fluvioglacial sands and gravels. Likewise, in Oliwa (Gdansk), the devastation of the forest activated mass wasting in the form of earth flow and earth creep on the slopes of Pachołek Hill.
The situation is similar for the other classes of causal factors (Table 3). However, factor slope does not show such a large variation. In addition, the slope factor is easy to count in most geoinformation programs such as ArcGIS, QGIS, and SAGA. For this reason, the slope factor can be most suitable for the development of susceptibility maps in young glacial areas, where there is no information on landslides.
Most landslides in the Tri-City occur in areas with steep slopes > 14 • (Table 3). These results differ from those obtained in similar GIS analyses performed for the Carpathian Mountains (Poland) where, for example, in the Lower Beskid (Beskid Niski) range the greatest susceptibility was observed on slopes of 9-14 • , whereas greater gradients display a clear decrease in the lnWi index [22]. The neighboring area near Szymbark (Lower Beskid), in turn, displays the greatest susceptibility with angles of 5-7 • and 18-20 • [20]. Similar results were obtained by Wojciechowski near Nowy Sącz, where landslides occupy the largest areas on slopes of 7 • -11 • (2009) [72]. The reason may lie in a different geological structure, because young glacial areas contain cohesive tills. They form steep slopes, especially in cliffs, with some of them angled at >80 • , while earth fall predominates in mass wasting processes. The high gradient of the Tri-City slopes with landslide occurrences can be traced back to a different role of the relief-forming processes which shaped the slopes. The processes of the last substage of the Vistulian Glaciation played a decisive role in shaping the major relief features [43]. Most slopes were formed by erosion over a relatively short period of time, from 10,000 to 12,000 years ago. In the case of the Carpathians, gentle slopes were formed on low-strength shale sediments. These sediments have poor resistance to degrading factors and are the most susceptible to landslides [73]. Furthermore, highresolution LIDAR data was used in this analysis to calculate slope angles, which may also have influenced the results. More methodologically correct detailed angle comparisons for slopes with landslide occurrences in various regions of Poland require the same kinds of source data. According to Burdziej and Kunz [74], the DEM acquisition method has a significant influence on DEM's accuracy in high resolutions (from 1 to 25 m), with the most serious errors corresponding to the forms with the greatest slope angles.
Not all the causal factors which affect mass wasting can be visualized as a parameter map; in the case of young glacial areas, the greatest difficulties are caused by crustal deformation associated with glacial thrusting (glaciotectonic deformation; Figure 2b,c) and the presence of perched groundwater, which are some of the major factors in landslide occurrence. Quaternary deposits are also a problem as the available digital geological data depicts only the general regularities of Pleistocene and Holocene formations. These sediments are of a highly variable character in terms of type and method of deposition, even over a distance of several meters. Such a detailed digital visualization of geological data is currently not possible for large areas of cities.

Conclusions
The digital elevation model, originating from airborne laser scanning (LIDAR data) proved very useful in landslide susceptibility studies. Based on the DEM, primary and secondary topographic attributes can be calculated. The TWI and the SPI are well suited for modelling landslides in which earthflow is the predominating mass movement. These types of landslides often occur in unconsolidated Quaternary formations (sands, gravels, silts).
The performed landslide susceptibility analysis for the Tri-City demonstrated relationships between landslides and causal factors: the LS index, slope angle, slope aspect, TWI, SPI, lithology, distance from watercourses and reservoirs, the groundwater table level, and land use. Potential landslide areas have been indicated, particularly in the Tri-City Landscape Park, the fringes of the morainic upland, and a number of erosionaldenudational valleys of the Kashubian Lakeland morainic upland. The SLI analysis for Gdansk and Gdynia has shown increased landslide susceptibility on slopes with a 10-29 • gradient, developed in sediments with lithological diversity in their stratigraphic profile, characterized by the presence of interbedded silty sands, silts, and clays. The co-occurrence of loose sediments above cohesive ones is conducive to the infiltration of precipitation and to the initiation of slip surfaces. A close proximity of surface water is significant due to the major role played in mass wasting by fluvial erosion, as well as the drainage function.
The new aspect presented in this article is the comparison of two susceptibility maps created using the same method for neighboring urbanized areas. The difficulty of extrapolating landslide data for neighboring urban areas was indicated. In the case of Gdansk, which had been covered by comprehensive geological mapping, the best modelling results were obtained when a large number of causal factors were taken into account. In Gdynia and Sopot, in turn, for which the SLI value was extrapolated from Gdansk, the best results were obtained when only five selected causal factors were included. This can be related to a different role of some causal factor classes in the occurrence of landslides in these areas, as well as the major significance of marine abrasion in landslide occurrence on the cliffs of Oksywie and Redłowo Hillocks.
The comparison of the two landslide susceptibility maps of Gdansk [6] and Gdynia [7] made using the SLI method also confirms the different importance of causal factors for neighboring urban areas.