2022 | 75/1 | 3–16 | 9 Figs. | 6 Tabs. | www.geologia-croatica.hr Journal of the Croatian Geological Survey and the Croatian Geological Society 1. INTRODUCTION Landslides have been in the focus of many researchers in the field of engineering geology, geotechnics, and geomorphology for many years, especially because they present a geohazard event that can cause extensive damage and casualties (VARNES & IAEG, 1984; BELL, 2003; HAQUE et al., 2016; MIHALIĆ AR- BANAS et al., 2017), and significantly contribute to landform evolution (CROZIER, 2010; BELL et al., 2012). For landslide analysis related to both of these subjects, preparation of a land- slide inventory map for a particular area is an initial step (GUZ- ZETTI et al., 2012), as it is believed that the locations of past and present landslides are the key to predicting future ones (VARNES & IAEG, 1984). A landslide inventory primarily records the land- slide location and, when known, the type of movement, activity, date of occurrence, and other landslide characteristics (MALA- MUD et al., 2004; FELL et al., 2008). Techniques of landslide mapping are constantly developing. In particular, detailed information about surface topography pro- vided by airborne laser scanning (ALS), also known as airborne light detection and ranging (LiDAR), improved the investigation of the Earth’s surface processes (JABOYEDOFF et al., 2012; TAROLLI, 2014). In landslide studies, the use of LiDAR has grown greatly during the last decade and has become a promising method for landslide mapping, monitoring, and modelling (DER- RON & JABOYEDOFF, 2010). Considering the LiDAR-based landslide inventories, many authors have highlighted the advan- tages of surface interpretation using bare earth high-resolution (HR) LiDAR digital terrain models (DTM) over conventional methods, i.e., geomorphological field mapping and visual inter- Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod, Croatia Iris Bostjančić*, Radovan Avanić, Tihomir Frangen and Mirja Pavić Croatian Geological Survey, Sachsova 2, 10000 Zagreb, Croatia; (*corresponding author: iris.bostjancic@hgi-cgs.hr) doi: 10.4154/gc.2022.03 Abstract A preliminary analysis of landslide spatial distribution and their geometric characteristics is pre- sented for the area of Slavonski Brod, located in the northeastern part of Croatia and belonging to the Pannonian Basin System. A landslide inventory for the study area of 55.1 km2 is accom- plished for the first time, based on the visual interpretation of a high resolution LiDAR digital ter- rain model. In total, 854 landslide polygons are delineated, corresponding to an average den- sity of 15.5 landslides per square kilometre. The average landslide area is 839 m2, and most of the landslides can be classified as small landslides (76 %). The spatial relationship between landslides and geological units is analysed and expressed as a landslide index. The Late Pan- nonian sands with silts and gravel interlayers and Pliocene clay, sands, gravels, and coal are determined as the units that are most susceptible to landslide processes. The majority of land- slides (85 %) are concentrated within these two units, for which a detailed analysis is performed, determining the morphometric parameters (slope and relief) and drainage network. The param- eters’ classes that create favourable preconditions to slope instabilities are defined, based on the landslide density within individual classes. Besides, the geometric characteristics of land- slides (size and shape) within these two units are compared. The results serve as the basis for further investigations. They help to foresee the area of future landslides through landslide sus- ceptibility maps, and offer a better understanding of the influence of fluvial-denudation and slope processes on recent landscape evolution and form. pretation of stereoscopic aerial photographs (e.g., VAN DEN EE- CKHAUT et al., 2007; JABOYEDOFF et al., 2012; GUZZETTI et al, 2012; GÖRÜM, 2019). The application of visual interpreta- tion techniques are tested for different landslide types and under a wide range of environmental conditions (e.g., ARDIZZONE et al., 2007; VAN DEN EECKHAUT et al., 2007; PETSCHKO et al., 2016; BERNAT GAZIBARA et al., 2019; PAWLUSZEK, 2019; GÖRÜM, 2019; JAGODNIK et al., 2020). Along with the increas- ing availability of HR digital elevation models, advances in Geo- graphical Information System (GIS) have improved quantitative analysis of the earth’s surface, including landslide zonation studi es related to inventory, and susceptibility, hazard or risk assessment (VAN WESTEN, 2000; CHACON et al., 2006). Until now, in Croatia, LiDAR-based landslide inventories exist for smaller areas of the City of Zagreb and the Vinodol Val- ley. In the populated area of the City of Zagreb, 24 km2 (MIHALIĆ ARBANAS et al., 2016), and 21 km2 (BERNAT GAZIBARA et al., 2019) were examined. In the Vinodol Valley, 18.55 km2 were studied in the central part of the Dubračina river basin (TOŠEVSKI, 2013), and 0.48 km2 in the Slani Potok gully (JAG- ODNIK et al., 2020). The study area is one of the six pilot areas selected in the framework of the safEarth project (“Transnational advanced management of land use risk through landslide susceptibility maps design”), which was carried out from 2017 to 2019 under the Interreg IPA Cross-border Cooperation Programme Croatia- Bosnia and Herzegovina-Montenegro 2014-2020. The main pro- ject goals were related to landslide susceptibility maps at different scales. For the project implementation, an ALS was performed, Article history: Manuscript received December 08, 2020 Revised manuscript accepted July 02, 2021 Available online February 28, 2022 Keywords: Landslide inventory, LiDAR, geological units, landslide spatial distribution, North Croatian Basin G eo lo gi a C ro at ic a 4 Geologia Croatica 75/1 which resulted in HR DTMs for an area of the app. 300 km2 in total. It gave the opportunity to accomplish landslide inventories over a wide area and increased the knowledge of landslide chara- cteristics and their spatial distribution. The pilot area presented in this paper (55.1 km2) is located in the northeastern part of Croatia, in the hinterland of Slavonski Brod. It represents the southwestern part of the Pannonian Basin System (HORVÁTH & ROYDEN, 1981; ROYDEN, 1988) and belongs to the North Croatian Basin (PAVELIĆ & KOVAČIĆ, 2018). In this area, mostly Miocene, Pliocene, and Quaternary sediments are exposed on the surface (ŠPARICA et al., 1979a; ŠPARICA, 1986). Most of the study area is covered by forest, and the most populated area, the centre of Slavonski Brod, is outside of the study area. However, the expansion of the settlements towards the north resulted in the highest city elevations of approx. 250 m. Landslide occurrence and damage to property have been recently documented after a period of heavy precipitation in March 2018. While severe damage occurred to private houses and local road infrastructure, a natural disaster was declared for the hilly areas of Slavonski Brod. This has led to increased interest in landslide processes and related problems by the local authorities. Here, a LiDAR-based landslide inventory for the area of Slavonski Brod has been completed for the first time. The pre- liminary results presented aim to define the landslide spatial dis- tribution and geometric characteristics within different geological units (GUs). This was accomplished through several steps: (i) identifying and mapping landslide polygons using HR LIDAR DTM derivatives, (ii) correlating the landslide distribution with GUs identified within the study area, (iii) defining morphometric characteristics within the two most representative GUs, and (iv) examining in more detail the relationship between morphometric parameters and landslide occurrence, as well as differences in landslide geometric characteristics within these units. 2. STUDY AREA 2.1. Geographical and geomorphological settings The study area is located in the northeastern continental part of Croatia (Fig. 1a) and covers an area of 55.1 km2. The area is situ- ated northwest of the City of Slavonski Brod, the main centre of Brod-Posavina County (Fig. 1b), and covers the area of three mu- nicipalities, Slavonski Brod, Podcrkavlje, and Sibinj (Fig. 1c). Slavonski Brod is situated at the junction of two spatial units. The larger, southern part of the city area belongs to the lowland area along the Sava river, and the smaller, northern part, to the wider hilly belt, in the hinterland of which is Dilj Mt. with the highest peak of 461 m. Several suburban settlements are situated within the study area, where 11 % (app. 7,100 inhabitants) of the popu- lation of Slavonski Brod (BPŽ, 2019) reside in Podvinje and Brod- sko Vinogorje. In these settlements, landslides were reported to Figure 1. Location map of the study area: Brod-Posavina County in Croatia (a); the study area northwest of the City of Slavonski Brod (b); elevation map derived from LiDAR DTM and drainage network obtained from the 1:25.000 topographic map (c). G eologia C roatica 5Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod local authorities in 2010 and 2018, when a large amount of pre- cipitation was recorded. According to the Climate Atlas of Croatia (ZANINOVIĆ et al., 2008), the climate of the Slavonski Brod area is characterized as a temperate continental climate, with harsh winters, short springs, and warm and humid summers. Based on a 30 year mea- surement period (1971–2000), and data recorded at the meteoro- logical station „Slavonski Brod“, in the study area the average annual air temperature is 10.7 °C. Average monthly temperatures rise until July when they reach a maximum of 21.0 °C. The avera ge annual amount of precipitation is 748.1 mm. Concerning the annual course of monthly precipitation, the main maximum occurs in July, and the main minimum in February. Still, certain years show a significant deviation in monthly amounts from the average precipitation conditions. Based on geomorphological regionalization (BOGNAR, 2001), the study area belongs to the mega-macro-geomorphologi- cal region of the Pannonian Basin, within which the hills of Dilj Mt. as a micro-geomorphological region stands out. In the study area, the elevation ranges from 94 m to 353 m, and most of the area (almost 73 %) is within the elevation class of 100-200 m (Fig. 1c). The elevation class of 300-350 m represents less than 1 % of the study area. The slopes vary from gentle to steep, up to 83°. The southern slopes of Dilj Mt. are very branched, expanding into a dense network of numerous gullies and valleys. Pediments and hills are formed by slope and fluvial-denudation processes (BOG- NAR, 1996). According to the area zonation by the degree of sta- bility, on the Engineering Geological Map of SFRY in a scale of 1:500.000 (ČUBRILOVIĆ et al., 1967) a large part of the study area is defined as a “prevailing unstable area under natural con- ditions and under the man’s activity mostly unstable”. 2.2. Geological setting In the investigated area eight GUs are recognized and distin- guished (Fig. 2): GU1 – Limestones and marls (Badenian); GU2 – Clayey limestones and marls (Early-Late Pannonian); GU3 – Marls with sandy interlayers (Late Pannonian); GU4 – Sands with silts and gravel interlayers (Late Pannonian); GU5 – Clay, sands, gravel, and coal (Pliocene); GU6 – Loess (Pleistocene); GU7 – Deluvial sediments (Holocene); and GU8 – Alluvial gravel, sands, and clay (Holocene). GU1 – Limestones and marls (Badenian) are located in the northwestern part of the study area. Marls are structureless and contain Lower Badenian (Orbulina universa, O. suturalis, Glo- bigerinoides bisphaericus, Globoquadrina altispira, and Uvige- rina macrocarinata), and Upper Badenian (Uvgerina venusta ve- nusta, U. venusta liesigensis, Bolivina dilatata, and Bulimina pyrula) planktonic and benthic foraminifera (ŠPARICA et al., 1979b). Limestones contain algal fragments (Lithothamnium), and mollusc shells (Pecten fuchsi stiryacus, Chlamys auensis zol- likoferi and Osrea sp. ŠPARICA et al., 1979b). Marls (Vejalnica Fm.) represent suspension deposits of the offshore (AVANIĆ et al., 2018a, b). Limestones with lithothamnium algae and mollusc shells indicate the nearshore environment of the Vrapče Fm. (AVANIĆ et al., 2018a, b). GU2 – Clayey limestones and marls (Early-Late Pannonian) are located in the northwestern part of the investigated area. They unconformably overlie Badenian limestones and marls showing an erosional contact. Clayey limestones contain typical macro- fossils including Radix croatica, Limnea extensa, and Gyraulus praeponticus (ŠPARICA et al., 1979b), while the marls are chara- cterized by the ostracods Hungarocypris auriculata, H. Hiero- gliphica, and Cyprideis heterostygma obesa (ŠPARICA et al., 1979b). Clayey limestones were deposited in the littoral area (Croatica Fm.), and marls are deposited in the deeper part of the lake (Medvedski Breg Fm.) (KOVAČIĆ et al., 2011; FILJAK et al., 2016; AVANIĆ et al., 2018a, b). GU3 – Marls with sandy interlayers (Late Pannonian) are lo- cated in the northern part of the study area and continuously over- lie the Lower-Upper Pannonian marls. A rich mollusc and ostra- cod fauna was recognized within these sediments (Paradacna abichi, P. lenzi, Didacna otiophora, Congeria digitifera, C. Za- grabiensis, Amplocypris recitulata, Candona (Caspiocypris) alta, C. (Caspiocypris) labiata, C. (Camtocypria) balcanica, Hemicy- theria pejinovicensis, H. Marginata, and P. croatica) (ŠPARICA et al. 1979b, 1987). Marls with sandy interlayers of the Andraševec Fm. (FILJAK et al., 2016), were deposited in the prodelta environ- ment (KOVAČIĆ et al., 2011; AVANIĆ et al., 2018a, b). GU4 – Sands with silts and gravel interlayers (Late Panno- nian) are situated in the central part of the studied area and con- tinuously overlie the marls with sandy interlayers. These sedi- ments contain a rich macrofauna (Congeria rhomboidea, C. Markovići, Dreissensia auricularis, Limnocardium petersi, L. (Euximicardium) ocetophorum, Didacna otiophora, Paradacna okrugići) (ŠPARICA et al., 1979b, 1987). Sands with silts and gravel interlayers (Late Pannonian) can be identified within the Nova Gradiška Fm. (FILJAK et al., 2016), and were deposited in the delta front environment (KOVAČIĆ et al., 2011; AVANIĆ et al., 2018a, b). GU5 – Clay, sands, gravel, and coal (Pliocene) are located in the southern part of the study area, and unconformably overlie sands with silts and gravel interlayers. The deposits are represented by a typical mollusc assemblage (Viviparus bifracinatus, V. dez- manianus, V. hoernesi, V. sturi, V. oncoforus, V. noihus, V. ornatus, V. neumayeri, V. avelanus, V. stricturatus, V. sadleri, Melanopsis hastata, M. hybostoma, M. pterochilia, and Unio sturzae) (ŠPARICA et al., 1979b). Gravel and sands are characteristic for alluvial deposits, and clay and coal accumulated in swamps, which can be identified as the Vrbova Fm. (AVANIĆ et al., 2018a, b). GU6 – Loess (Pleistocene) extends along the southern boundary of the investigated area and overlies Pliocene sedi- ments. Its granulometric composition is dominated by silt (50- 70%), with up to 30% sand and clay (ŠPARICA et al., 1979b). It contains the typical loess malaco-fauna (Trichia hispida, Suc- cinea oblonga, Pupilla muscorum, P. sterri, Vallonia tenuilabris, V. enniensis, and Arianta arbustorum) (ŠPARICA et al., 1979b, 1987). Pleistocene loess is aeolian sediment. GU7 – Deluvial sediments (Holocene) are located in the northeastern part of the study area and are mostly redeposited marls with sandy interlayers, formed after a short period of trans- port (ŠPARICA et al., 1987). GU8 – Alluvial sediments (Holocene) fill the stream valleys in the central and northeastern parts of the area and are composed of gravel, sand, and silt. These deposits originate from older se- diments of the area, from the Badenian to the Pleistocene. Based on the described GUs a simplified evolution of the study area can be proposed. During the Badenian in the deeper marine part (offshore) marls (Vejalnica Fm.), and in the shallow nearshore part lithothamnia limestones (Vrapče Fm.) were de- posited. In the area of Dilj Mt. the Sarmatian is characterized by marine sedimentation in an environment with reduced salinity (Dolje Fm.). In contrast to the northern part of the Dilj Mt., in the study area, no sedimentation occurred during the Sarmatian. This G eo lo gi a C ro at ic a 6 Geologia Croatica 75/1 is probably due to local tectonics during that time. Later, due to the Early Pannonian transgression, clayey limestones of the lit- toral environment begin to accumulate, followed by the Early and Late Pannonian marls of the Medvedski Breg Fm. The marls were deposited in a relatively deeper part of the newly formed brack- ish Pannonia lake (AVANIĆ et al., 2018a, b). During the Late Pan- nonian, marls with sandy interlayers accumulated in the prodelta environment (Andraševec Fm.), and sands with layers of silt and gravel of the delta front (Nova Gradiška Fm.), and during the Pliocene, sands, gravels, clays, and coals related to rivers and swamps. Sedimentation from the Lower to the Upper Pannonian has characteristics of a transgressive sequence, and from the Late Pannonian to the Pliocene of a regressive sequence. In the trans- gressive sequence, there was a deepening from the littoral to the lake basin, and in the regressive sequence, there was a shallow- ing or propagation of the river clastic system, when after the lake sediments, the deposits of prodelta, delta-front, alluvial, and swamp plains accumulated. The final terrestrial sedimentation occurred during the Quaternary and is represented mostly by loess, as well as by deluvial and alluvial sediments. 3. MATERIALS AND METHODS 3.1. Input data The HR DTM used in this study is a result of an ALS, which was performed in the early spring of 2018 by Slovenian company Fly- com Technologies d.o.o., fulfilling the requirements for at least 20 points per square meter and assuring the accuracy for each point of at most 10 cm in all directions. As a final product of raw LiDAR data post-processing, DTM with a 0.5 m resolution is ob- tained. At the same time, topographic scanning was performed resulting in HR orthophotos with a 10 cm resolution. Based on geological reconnaissance and existing basic geo- logical maps (ŠPARICA et al., 1979a; ŠPARICA, 1986), a geo- logical map of the investigated area was compiled and individual GUs (FILJAK et al., 2016) were interpreted. 3.2. Landslide mapping Noticeable morphological signs that remain after landslide oc- currence can be recognized and mapped through the interpreta- tion of digital topographic surface model (GUZZETTI et al., 2012). The methodological approach for landslide recognition is based on visual interpretation and geomorphic feature extraction (SCAIONI et al., 2014) using HR LiDAR DTM derivatives. For that purpose, several DTM derivatives most commonly used for landslide inventory preparation were extracted, using standard ArcGIS tools. These are: (i) hillshade maps (according to SCHULZ (2004) different sun azimuth angles of 315°, 135°, and 45°, and a constant altitude angle of 45° were used); (ii) slope map; (iii) profile curvature map; (iv) contour line map with elevation equidistance of 1 m (Fig. 3). DTM derivatives were used indi- vidually or in combination to delineate landslide polygons. Indi- Figure 2. Geological map of the study area modified according to 1:100,000 scale basic geological maps (ŠPARICA et al., 1979a; ŠPARICA, 1986), with landslide in- ventory. G eologia C roatica 7Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod vidual landslides are identified and digitized as polygons at a large scale close to 1:500, comprising the completely affected area from the landslide crown to the toe. When more different land- slide generations were recognized, each landslide polygon was delineated separately whenever possible in order to obtain the real and representative landslide geometric characteristics. As most of the study area is covered with forest, the use of orthopho- tos was limited, although in some locations the landslide main characteristics were clearly recognized (e.g., main scarp, crown cracks, flank scarps) (Fig. 4). All these morphological signs are generally sharper, better developed and more clearly visible for recently active landslides and are modified through time by different processes of weather- ing, erosion, and deposition (MCCALPIN, 1984). As the land- slide ages, the more its morphology is reshaped and the less rec- ognisable the landslide features become in the field (MCCALPIN, 1984; KEATON & DEGRAFF, 1996; BELL et al., 2012; PETSCHKO et al., 2014). Thus, the recognition of typical land- slide characteristics is evaluated by several authors titled as land- slide freshness (VAN DEN EECKHAUT et al., 2007), confidence (BURNS & MADIN, 2009), certainty of mapping (PETSCHKO et al., 2016), or certainty of landslide identification (BERNAT GAZIBARA et al., 2019). In this study, the landslide outline confidence level (LOCL) is presented as the degree of landslide contour certainty based on the visibility and expressiveness of the landslide main features and is expressed with scores from 1 to 10. Each of the following landslide features is rated with a confidence score from zero to two: (i) main scarp, (ii) lateral flanks, (iii) toe, (iv) internal land- slide morphology (e.g. hummocky or undulating topography, in- ternal cracks, minor scarps), and (v) the ability of discernment of depletion and accumulation areas. Based on the confidence scores, four LOCL classes are defined (Tab. 1). A landslide inventory was produced as the basis for the analy- sis of landslide spatial distribution and geometric characteristics. All of the analyses were performed for landslides with a confi- dence score higher than two, i.e., for LOCL classes II, III, and IV. In addition, the minimum landslide area for mapping was set to Figure 3. Examples of landslide polygon delineation for high landslide outline confidence level (A) and moderate landslide outline confidence level (B), using DTM derivatives: contour line map over the hillshade map (a); hillshade map with the transparency of 60 % over the slope map (b); hillshade map with the transparency of 60 % over the profile curvature map (c); landslide polygon on the hillshade map with movement direction (d). Figure 4. Visibility of landslide characteristics on HR orthophotos: landslide main scarp with several minor scarps (a); landslide main scarp and crown cracks (b). Table 1. Description of landslide outline confidence level (LOCL). LOCL class Confidence score Description I 1-2 Possible landslides with questionable delineated boundaries II 3-4 Landslides with unreliable delineated boundaries in some parts III 5-7 Landslides with imprecise delineated boundaries in some parts IV 8-10 Landslides with high precise delineated boundaries G eo lo gi a C ro at ic a 8 Geologia Croatica 75/1 50 m2, given the scale of the input DTM. In the analysis per- formed, no distinction was made between the different landslide types, or according to their depth and relative ages. The LiDAR-based landslide inventory was verified during a field engineering geological inspection, conducted in March 2019. 3.3. Landslide inventory analysis Preliminary analysis of the landslide inventory was performed, which aimed to define landslide size, shape, and spatial distribu- tion in the study area. Geometric parameters used to analyse landslide size and shape are presented for every landslide polygon with several fea- tures that were determined automatically (area), and manually (length, width) using standard ArcGIS tools, or calculated by equation (aspect ratio). Based on the landslide area (expressed in square metres), landslides were grouped into five classes (ŠESTANOVIĆ, 2001): very small (<100), small (100-1,000), moderate (1,000-10,000), large (10,000-50,000), and very large (>50,000). The landslide aspect ratio (LA) represents the ratio of land- slide length to width. According to the classification based on the landslide aspect ratios proposed by TIAN et al. (2017), landslides were grouped into four classes: transverse (LA≤0.8), isometric (0.83). The relationship between GUs and landslide occurrence within the study area is evaluated and represented by the land- slide index. The landslide index, expressed in percentages, is cal- culated as a ratio between landslide area in a given GU and the area of that unit. The landslide index of each GU is compared to the landslide index of the entire study area, also expressed in per- centages and calculated as a ratio between the total landslide area and total study area (BARTELLETTI et al., 2017). In terms of statistics, the correlation between landslides and morphometric parameters (slope, relief, distance to streams) of two major GUs (Fig. 2) is examined and expressed through land- slide density for each parameter class. Individual landslide den- sity for each parameter class is compared to the landslide density of each GU, which is calculated as a ratio between landslide num- ber within an individual GU and an area of that unit. For that pur- pose, centroids of landslide polygons are created. The slope angle was automatically derived from DTM resa- mpled to 5 m resolution. In order to observe landslide distribution from the slope angle perspective, the terrain was classified ac- cording to DEMEK (1972), whose classification is based on the dominant processes that activate on a slope of a particular angle. Six slope classes (expressed in degrees) are used: 0-2 (plain ter- rain), 3-5 (slightly sloping terrain), 6-12 (sloping terrain), 13-32 (very sloping terrain), 33-55 (very steep terrain), and > 55 (cliff). Relief is expressed as a difference between the maximum and minimum elevation in a certain area and is calculated using the ArcGIS tool Focal statistics. The radius of the circle taken around each cell to calculate the elevation range is set to 250 cells, i.e., 125 m. The six classes (expressed in metres) are defined au- tomatically using natural breaks and are: 0-18 (flattened area), 19-34 (very low relief), 35-47 (low relief), 48-60 (moderate re- lief), 61-75 (high relief), and 76-109 (very high relief). To calculate the distance to streams, a detailed stream net- work was automatically extracted from the DTM resampled to 5 m resolution, using the Arc Hydro tools (v2.0). The standard pro- cedure included filling sinks, computing the flow direction, and flow accumulation. Stream definition is based on a flow accumu- lation and a specified threshold value, which represents the num- ber of cells indicating the start of a stream. The different thresh- old areas were tested and a representative network was chosen based on visually cross-checking the main streams using a topo- graphic map at a scale of 1:25,000. However, strong topographic diversity influenced the selection of different threshold areas that are associated with the various GUs. The stream order was auto- matically calculated using the method proposed by STRAHLER (1956). A distance to stream map is created as buffer zones around streams and is classified into six classes (expressed in metres): <20, 20-40, 40-60, 60-80, 80-100, and >100. The landslide inventory and spatial analysis were performed using ESRI ArcGIS 10.2. The official Croatian coordinate projec- tion, HTRS96/TM, is used for all maps presented herein. The visualization of graphs was aided by Daniel’s XL Toolbox add-in for Excel, version 7.3.2 (KRAUS, 2014). 4. RESULTS 4.1. Landslide inventory In the study area, a total of 854 landslide polygons were deline- ated, corresponding to an average landslide density of 15.5 land- slides per square kilometre. The total mapped landslide area is 0.72 km2 (Tab. 2), which represents 1.31 % of the total area. The boundaries for three landslides exceed the study area, but the polygons are completely delineated within the DTM buffer zone, and their actual sizes are used for landslide area statistics. Thus, the total study area affected by landslides is, therefore, smaller by 0.004 km2. The landslide areas range from 50.43 m2 to 17,855.32 m2, with an average value of 838.93 m2 (Tab. 2). The distribution of landslide area is presented in more detail via the frequency den- sity and a cumulative number of landslides (Fig. 5). It is shown that 3.5 % of landslides are smaller than 100 m2, and 20.5 % are larger than 1,000 m2. Most landslides (76 %) range between 100 m2 and 1,000 m2, of which 54.4 % are in the range of 100 m2 to 500 m2. Within this range, the highest frequency density of 1.48 m2 is presented for landslide areas ranging from 200 m2 to 300 m2. Only five landslides are larger than 10,000 m2. Considering different classes of LOCL, the number of land- slides decreases as LOCL increases (Tab. 2). Thus, almost 60 % of landslides mapped are described with the LOCL class of II, and a little over 40 % of landslides with LOCL class of either III or IV. As LOCL gives an overview of the expression of landslide features, these results show that in the study area, all typical land- Table 2. Statistics of landslide size in the study area, considering different land- slide outline confidence level (LOCL). LOCL (Confidence score) II (3, 4) III (5, 6, 7) IV (8, 9, 10) All landslides (3-10) No. of landslides (#) 501 327 26 854 No. of landslides (%) 58.67 38.29 3.04 100.00 Landslide density (#/km2) 9.10 5.94 0.47 15.50 Landslide area (m2) 286,486.50 358,124.49 71,832.25 716,443.24 Smallest landslide area (m2) 50.43 69.86 374.03 50.43 Largest landslide area (m2) 12,878.26 17,855.32 15,618.94 17,855.32 Average landslide area (m2) 571.83 1,095.18 2,762.78 838.93 Median landslide area (m2) 322.36 552.91 1,451.41 408.48 G eologia C roatica 9Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod slide features are not completely visible for most of the landslides (60 %) and the boundaries are particularly hard to determine un- ambiguously. It can be interpreted as either landslide activation occurred a long time previously explaining why their sharp boundaries and morphological signs have been masked over time, or the landslides are younger but their main features are lost due to fast and strong erosion processes. In addition, the average landslide area increases as LOCL increases (Tab. 2), which is es- pecially a result of a great variability of the smallest landslides within each LOCL class (the smallest landslide for LOCL IV is almost 7.5 times larger than the smallest landslide for LOCL II). The landslide inventory is verified in the field by checking delineated polygons and identifying the morphological features typical of landslides. The verification is conducted for landslides with different confidence scores (from 3 to 8), and a range of ar- eas between approximately 70 m2 and 7,600 m2, covering the area within three different GUs (GU3, GU4, and GU5). From 97 land- slides visited, 94 were verified, which represents 11 % of the Li- DAR-mapped landslides. The remaining three landslides are all located in the forest and are marked with LOCL of II, i.e. confi- dence score of 3, but in the field, typical landslide morphology was not confirmed. 4.2. Landslides and geological units It is recognized in many landslide studies that the geological set- ting greatly influences landslide occurrence. As one of the main predisposing environmental factors, geological data are used as one of the main input layers in landslide susceptibility assessment (VAN WESTEN et al., 2008). Thus, for the landslide spatial dis- tribution analysis, first, the relationship between landslide occur- rence and GUs is investigated. The relationship is expressed with the landslide index (Fig. 6), and the landslide statistics for each GU are presented (Tab. 3). The landslide index of each GU was compared to the landslide index of the entire study area (BAR- TELLETI et al., 2017), enabling determination of which units are more correlated to landslide source areas. Considering the landslide area distribution for different GUs, it can be seen that landslides are mainly concentrated within GU4 and GU5, where 85 % of landslides are present (Fig. 6). Within GU4, 423 landslide polygons are delineated, and within GU5, 303 landslide polygons, which corresponds to landslide densities of 31.66 and 18.73 landslides per km2, respectively. These GUs refer to Late Pannonian sands with silts and gravel interlayers (GU4) and Pliocene clay, sands, gravels, and coal (GU5). The landslide index in these units is higher than the landslide index of the study Figure 5. Frequency density and cumulative number of landslides considering landslide area. Figure 6. Landslide area distribution in different geological units (GU1 – Limestones and marls (Badenian); GU2 – Clayey limestones and marls (Early - Late Panno- nian); GU3 – Marls with sandy interlayers (Late Pannonian); GU4 – Sands with silts and gravel interlayers (Late Pannonian); GU5 – Clay, sands, gravels, and coal (Pli- ocene); GU6 – Loes (Pleistocene); GU7 – Deluvial sediments (Holocene); GU8 – Alluvial gravels, sands, and clay Holocene). The landslide index for each geological unit and for the study area is presented. G eo lo gi a C ro at ic a 10 Geologia Croatica 75/1 area (1.30 %) at 2.03 % and 2.11 %, respectively. At the same time, as already previously stated, these units are most represented in the study area, covering almost 30 km2 (54 % of the study area). However, the highest landslide index (4,13 %) is calculated for GU1 (Badenian limestones and marls), representing only 21 land- slides. Consequently, the area covered by GU1 is only 0.6 % of the study area. For all other GUs, the landslide index is less than 1, i.e., below the landslide index of the study area (Tab. 3). According to the average landslide area, the largest land- slides occur in GU5, and the smallest in GU4 (Tab. 3). The dif- ferences in landslide average width are not so noticeable, with values ranging from 20.05 m to 25.61 m. The values of landslide average length, compared to others, for GU5 and GU6 stands out, with values of 51.55 m and 51.86 m, respectively. 4.3. Landslides within the predominant geological units 4.3.1. Morphometric analysis The diversity of GU4 and GU5 environments is represented by their morphometric features, including slope angle, relief, and drainage network (Fig. 7). Although the range and average values of the slope angles within GU4 and GU5 are similar, the slope angle class distribu- tion shows differences (Fig. 7a). Almost 90 % of the area of GU4 is within three classes in the range of 6° to 55°, with the most rep- resentative class of 13-32° (54 %). Most of the area of GU5, almost 95 %, is within three classes in the range of 3° to 32°, with the most representative class of 6-12° (49 %). Slope angle class with values above 55° is the least represented, with an area below 1 % for both GUs. The values of relief for GU4 and GU5 range from 5 m to 109 m, and from 9 m to 100 m, respectively (Tab. 4). For GU4, 81 % of the area is within three relief classes in the range of 48 m to 109 m, with the most representative class of 61-75 m (34 %). For GU5, 87 % of the area is within three relief classes in the range of 19 m to 60 m, with the most representative class of 35-47 m (44 %). For both GUs the relief class with values below 19 m is least represented covering approx. 1 % in both GUs (Fig. 7b). Concerning the classes defined to analyse the distance to streams, 65 % of the area of GU4 is within the first three classes, with distances to streams less than 60 m. Within this unit, most of the area (25 %) is in the stream buffer zone of 0-20 m. For GU5, all classes are almost equally represented, with significant repre- sentation by class (43 %) where the distance to stream is more than 100 m (Fig. 7c). These results indirectly indicate that the stream network within GU4 is more branched and denser than within GU5. This is also shown by the results of the analysis of stream length for different stream orders (Tab. 4). Since the LiDAR poly- gon does not include the entire basin, only streams that are com- pletely developed (streams of a certain order and all lower orders) in a single GU, are considered. There are half as many 1st order streams within GU4 than GU5, although the total length of 1st order streams does not differ much between these two units. Con- sequently, there is a big difference in their average stream lengths, 87.70 m within GU4 and 164.06 m within GU5. Thus, GU4 can be described as a unit with a short but dense series of 1st order streams, and GU5 as a unit where 1st order streams are rarer but much longer. The same trend is presented for 2nd order streams. Furthermore, the difference in the stream number and stream lengths decreases for 3rd order streams, and is even less for 4th order streams. 4.3.2. Landslide spatial distribution In two zones characterized by different geological and geomor- phological conditions, landslide spatial distribution is analysed. Classes of parameters that create favourable preconditions to slope instability, presented in the previous sub-chapter, are de- fined based on landslide density (Fig. 8). Considering the slope angle, landslides within GU4 are mainly concentrated within classes of 13-32° (52 % landslides) and 33-55° (45 % landslides). The highest landslide density is also related to these classes, with the highest peak of 115 landslides per km2 within the class of 33-55° (Fig. 8a). Within GU5, a large number of landslides is also related to classes of 13-32° (64 % landslides) and 33-55° (12 % landslides). Still, landslides are also concentrated on somewhat milder slopes, within a class of 6-12° (23 % landslides). The highest landslide density is related to the class of 33-55°, with 79 landslides per km2 (Fig. 8a), same as for GU4. The differences in the relationship between slope angle and Table 3. Landslide statistics for each geological unit. Geological unit 1 2 3 4 5 6 7 8 No of landslides 21 39 38 423 303 7 0 23 No of landslides (%) 2.46 4.57 4.45 49.53 35.48 0.82 - 2.69 No of landslides LOCL II 11 26 21 271 157 2 - 13 No of landslides LOCL III 10 13 17 142 130 5 - 10 No of landslides LOCL IV 0 0 0 10 16 0 - 0 Landslide index 4.13 0.95 0.34 2.03 2.11 0.18 - 0.25 Landslide density (# / km2) 61.93 11.80 4.30 31.66 18.73 2.01 - 2.69 Average landslide area (m2) 667.53 802.74 793.41 642.24 1,127.97 883.18 - 927.88 Average landslide length (m) 31.90 35. 33 31.50 34.10 51.55 51.86 - 39.39 Average landslide width (m) 24.95 24.26 25.61 20.05 22.77 22.14 - 22.13 Table 4. Number and length of streams considering different stream orders within geological units 4 and 5. Geological unit Stream order No of streams analysed Total stream length (m) Average stream length (m) 4 1 380 33327.74 87.70 2 69 6258.184 90.70 3 11 739.37 67.22 4 2 115.51 57.75 5 1 183 30022.92 164.06 2 39 6773.36 173.68 3 8 574.84 71.85 4 1 134.70 134.7 G eologia C roatica 11Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod Figure 7. Spatial distribution of parameters’ classes within geological units 4 and 5, concerning: slope (a); relief (b); and distance to streams (c). G eo lo gi a C ro at ic a 12 Geologia Croatica 75/1 landslide occurrence can be expressed by the average values of slope angle within landslide polygons and the rest of the area. Thus, for GU4 the average slope angle within areas affected by landslides is 27°, and for GU5 is 18°. For the area not affected by landslides, the average slope angle for GU4 is 17° and GU5 is 12°. In a view of relief, 96 % of landslides within GU4 are con- centrated within classes with relief values higher than 48 m, with the highest concentration within the class of 61-75 m (45 % land- slides). The highest landslide density is also related to these classes (Fig. 8b). Within GU5, the highest landslide concentration is related to classes with relief values from 35 m to 75 m (64 % landslides). The highest landslide density is, however related to the class of >76 m (Fig. 8b). The differences in the relationship between relief and landslide occurrence expressed by the average values of relief within landslide polygons and the rest of the area are more significant compared to the slope angle. For GU4 the average relief within areas affected by landslides is 65 m, and for GU5 it is 51 m. For the area not affected by landslides, the ave- rage relief for GU4 is 61 m and for GU5 it is 45 m. Considering the distance to streams, 92 % of landslides within GU4 are concentrated within the classes 0-20 m (62 % landslides) and 20-40 m (30 % landslides). The highest landslide density is also related to these classes, with the highest peak of 78 landslides per km2 within the class of 0-20 m (Fig. 8c). Within GU5, 81 % of landslides are concentrated within classes with a distance below 60 m. Still, the highest landslide density is related to the class of 0-20 m, with a value of 52.6 landslides per km2 (Fig. 8c). 4.3.3. Landslide geometric characteristics Several landslide characteristics are considered, including area, length, width, and aspect ratio (Table 5), in order to compare landslide size and shape within GU4 and GU5. Comparing landslide sizes, the average landslide area in GU4 (642.24 m2) is almost half the size of the average landslide area in GU5 (1,127.97 m2). The range of landslide areas within GU4 and GU5 is quite different (see box plots, Fig. 9a) though the outliers are not represented, thus the maximum landslide areas are not shown. The difference in landslide width is not as significant as land- slide length, as was the case with the analysis of landslide inven- tory for each GU (Tab. 3). The average landslide lengths for GU4 Figure 8. Landslides and landslide density (number of landslides per area of each parameter class) for geological units 4 and 5, concerning: slope angle (a); relative relief (b); and distance to streams (c). Landslide density for each geological unit (number of landslides per area of the geological unit) is also presented. G eologia C roatica 13Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod and GU5 are 34.10 m and 51.55 m, respectively (Tab. 5). This dif- ference influences not only the landslide area but also the land- slide shape. The average aspect ratio is thus slightly larger for landslides within GU5. The average aspect ratio for GU4 and GU5 is 1.88 and 2.27, respectively (Tab. 5). The aspect ratio classes indicate that for both GUs most landslides occur within two classes and landslides are mainly longitudinal. Further, it can be stated that an isometric shape is more associated with land- slides within GU4 and an elongated shape with landslides within GU5 (Fig. 9b). 5. DISCUSSION Landslide inventories are the basis for assessing landslide sus- ceptibility, hazard, and risk assessment (SOETERS & VAN WESTEN, 1996; VAN WESTEN et al., 2008). Without know- ledge about landslide locations, the methods for susceptibility models that predict landslides based on past conditions remain inapplicable. The landslide inventory accomplished for the study area includes landslides that occur over time (periods of tens, hundreds, or even more years), and are not associated with a trig- ger event. Thus, this inventory can be interpreted as a historical geomorphological inventory (MALAMUD et al., 2004; GUZ- ZETTI et al., 2012). The landslide inventories are widely used to analyse land- slide statistics, including their size, shape, and location (MAL- AMUD et al., 2004). Considering the landslide size and shape, geometricl parameters and their ratios can be further analysed to characterize landslides and assess their volumes, motion mecha- nisms, runout distances, and hazard (TIAN et al., 2017). Consi- dering the landslide distribution, it has been proven that landslides mostly occur at preferred locations. Among others, geological- lithological and morphometric parameters have been used as the main environmental predisposing factors to analyse the effects on landslide spatial distribution (VAN WESTEN et al., 2008; REICHENBACH et al., 2018). The distribution of landslides con- cerning the main factors controlling shallow landslides (litho- type, regolith thickness, slope angle, and land use) is discussed for the Picentino river basin in southern Italy (LUPIANO et al., 2019). CONFORTI & IETTO (2020) highlighted that the litho- logy, fault density, slope gradient, and relief have an important role in determining landslide occurrence. GUZZETTI et al. (2008) examined the relationship between landslide distribution and local morphological and lithological settings, in the Upper Tiber River basin in central Italy. Their analysis confirmed that geological conditions influence the abundance of landslides in the catchment. JACOBS et al. (2016) also concluded that landslide processes largely depend on the prevailing lithology. In their study, slope angle is distinguished as the main controlling topo- graphic factor for landslides. BARTELLETTI et al. (2017) inves- tigated the influence of geological-morphological and land use settings on the occurrence of shallow landslides. They used the landslide index to identify significant parameters correlated with landslide source areas. In lithologically homogeneous areas, (as each of the presented Gus may be considered), COCO & BUC- COLINI (2015) indicated that the slope morphometry, especially angle, height, and form were the most influential factors on land- slide occurrence. They also confirm that landslides are coupled with a drainage network. In this study, the geological setting is proven to be the sig- nificant influencing factor that controls landslide occurrence and abundance. Besides analysing the number of landslides within individual GUs, the correlation between GUs and landslide oc- currence is estimated and expressed with the landslide index. The highest landslide index is defined for the unit represented by the Badenian limestones and marls (GU1). Since the area covered by GU1 represents only 0.6 % of the study area where only 21 land- slides are mapped, further investigations are needed in areas where this unit is much more representative. The next two units for which the landslide index is above the landslide index of the study area refer to Late Pannonian sands with silts and gravel in- terlayers (GU4) and Pliocene clay, sands, gravels, and coal (GU5). As 85 % of the mapped landslides are concentrated within these units, they can be considered to be the most susceptible to land- slide processes in the study area. When taken into consideration the uniform geographical, climatic, and tectonic settings for GU4 and GU5, it can be rea- sonably assumed that the geological setting is the prevailing fac- tor in the evolution of this area. Significant differences in mor- phometric features are strongly correlated to the geological Table 5. Descriptive statistics for landslide characteristics for geological units 4 and 5. Parameter Geological unit Min Max Average Median Area (m2) 4 50.77 17,855.32 642.24 368.40 5 63.14 15,618.94 1,127.97 475.36 Length (m) 4 9.00 233.00 34.10 30.00 5 8.00 375.00 51.55 37.00 Width (m) 4 6.00 115.00 20.05 16.00 5 6.00 101.00 22.77 19.00 Aspect ratio 4 0.48 4.71 1.88 1.78 5 0.40 8.90 2.27 2.09 Figure 9. Area (a) and aspect ratio (b) for landslides within geological units 4 and 5. G eo lo gi a C ro at ic a 14 Geologia Croatica 75/1 setting, i.e. lithological features and material characteristics, which are clearly visible on Fig. 7. Still, it is necessary to empha- size the scale of the input of geological data (Basic geological map in scale 1:100,000), which may affect the presented results in the boundary areas of GUs. This can also explain the fact that 2.69 % of landslides is placed within GU8 (alluvial gravel, sands, and clay), as presented in Table 3. The unified results of detailed analysis performed for GU4 and GU5, elaborating the relationship between morphometric fea- tures and landslide source areas, and geometrical characteristics of landslides, are presented in Table 6. Observing the classes of morphometric parameters analysed, it is evident that the predominant classes for slope and relief within these two units differ for one to two class categories (Tab. 6). Considering the distance to streams, the dominant classes dif- fer even more. Significant morphometric characteristics within GU4 refer to high relief, steep slopes, the dense stream network, and branched gullies. Dominant distance from stream classes are less than 60 m. Morphometric characteristics within GU5 refer to low relief, less steep slopes, and a less dense stream network, when compared to GU4, with dominant distance from stream classes greater than 100 m. Regardless of these differences, parameters’ classes defined as the most influenced on landslide spatial distribution are almost equal for both of these units. The classes for which the landslide density of parameter class is higher or very close to the landslide density of the entire GU are singled out. Based on these results, landslides are predominantly occurring on very steep and very sloping terrain, with very high, high, and moderate relief, and are strongly coupled to streams (Tab. 6). The fluvial-denudation pro- cesses, together with slope processes have the dominant influence on the development of local landscape. The geomechanical prop- erties of sands are a dominant factor for current landscape for- mation in GU4. Their drainage capabilities and low cohesion causes the development of a dense gully network due to strong erosion processes. Consequently, the majority of the landslides are associated with the steep gully banks. Typical landslides are small and shallow, often without clearly visible landslide ele- ments in the foot area where the landslide material is removed by permanent or intermittent streams. In contrast to GU4, low drain- age properties of GU5 govern landscape shapes in their respec- tive area, where the soil complex comprises impermeable clay interlayers or lenses. Landslides are predominantly triggered by high groundwater pore pressures. Landslides are slightly larger and more elongated, when compared to landslides within GU4, which implies the flow of material when the landslide is activated. The differences in landslide geometric characteristics within these units implies that geological characteristics and consequent morphometric features influence not only the landslide density and spatial distribution but also the landslide size and shape. 6. CONCLUSION In this paper, preliminary analysis of the initial LIDAR-based landslide inventory for the area of Slavonski Brod is presented for the first time. Using HR LiDAR DTM derivatives (hillshade, slope, curvature, and contour line maps), the geomorphological historical landslide inventory for an area of 55.1 km2 was com- pleted. In total, 854 landslide polygons were delineated, corre- sponding to an average landslide density of 15.5 landslides per square kilometre. The landslide area ranges from 50.43 m2 to 17,855.32 m2, with an average value of 838.93 m2, according to which they can be considered as small landslides. Given the strong relationship between landslide occurrence and geological setting in many studies, first, the spatial landslide distribution was evaluated, according to eight GUs defined in the study area. The correlation between GUs and landslide occur- rence is estimated and expressed with the landslide index. Ac- cording to differences in the landslide index, the geological set- ting can be considered as a significant feature in landslide susceptibility assessment. Late Pannonian sands with silts and gravel interlayers (GU4) and Pliocene clay, sands, gravels, and coal (GU5) are those units where the majority of landslides are concentrated. 85 % of landslide polygons mapped are present within these units. Besides, these units are predominantly repre- sented in the study area, covering an area of almost 30 km2. The landslide index for both of these units is above the landslide index of the study area. According to all the above, GU4 and GU5 were chosen for detailed investigation, in order to analyse the landslide spatial distribution and compare landslide geometrical character- istics. Defining the relationship between landslide distribution and environmental parameters enables a better understanding of landslide characteristics and occurrence, and can help to predict the future landslides. Parameters analysed in this study (slope, relief, and distance to stream) are proven to influence landslide spatial distribution, and classes with a stronger correlation with landslide density are defined. Table 6. Main morphometric and landslide characteristics within geological units 4 and 5. Geological unit Late Pannonian sands with silts and gravel interlayers (GU4) Pliocene clay, sands, gravels, and coal (GU5) Main morphometric characteristics – Sloping to very steep terrain (dominant class of 13-32°) – Moderate to very high relief (dominant class of 61-75 m) – Dense stream network with large number of 1st order streams (average length of 87.70 m) – Branched gullies with steep side walls – Slightly to very sloping terrain (dominant class of 6-12°) – Very low to moderate relief (dominant class of 35-47 m) – Rarer stream network with less number of 1st order streams (average length of 164.06 m) – Gullies with less steep side walls Landslide spatial distribution – Landslide density of 31.66 per km2 within entire unit – Highest landslide density for very steep (33-55°) and very sloping (13-32°) terrain – Highest landslide density for high (61-75 m), very high (76-109 m), and moderate (48-60 m) relief – Highest landslide density at a distance of 0-20 m and 20-40 m from streams – Landslide density of 18.73 per km2 within entire unit – Highest landslide density for very steep (33-55°) and very sloping (13-32°) terrain – Highest landslide density for very high (76-109 m), high (61-75 m), and moderate (48-60 m) relief – Highest landslide density at a distance of 0-20 m, 20-40 m, and 40-60 m from streams Main landslide characteristics – Very small to moderate landslide size, according to average area (642.24 m2) small landslides – Longitudinal to isometric shape – Dominant mapping landslide outline confidence level of II – Very small to large landslide size, according to average area (1,127.97 m2) moderate landslides – Longitudinal to elongated shape – Dominant mapping landslide outline confidence level of II and III G eologia C roatica 15Bostjančić et al.: Spatial distribution and geometric characteristics of landslides with special reference to geological units in the area of Slavonski Brod The results of landslide spatial distribution compared to ge- ological features, slope morphometry, and stream network, pro- vide useful information for further studies related to landslide susceptibility and hazard. In that sense, the basic information of significant predisposing factors controlling landslide occurrence in the study area is given. Together with the landslide inventory created, a large proportion of relevant input data is prepared. Those data are of great importance and can be used to produce statistically-based landslide susceptibility maps and to predict areas of higher landslide susceptibility. It would allow the local authorities to direct detailed engineering geological investiga- tions and financial resources toward landslide-prone areas, highly relevant for the hilly area of Slavonski Brod. ACKNOWLEDGEMENT The basis for landslide inventory and all the analysis presented in this study was a high resolution LiDAR-based digital terrain model. It was obtained within the implementation activities in the framework of the safEarth project, which was co-financed by the Interreg IPA Cross-border Cooperation Programme Croatia – Bosnia and Herzegovina – Montenegro 2014-2020. REFERENCES ARDIZZONE, F., CARDINALI, M., GALLI, M., GUZZETTI, F. & REICHENBACH, P. (2007): Identification and mapping of recent rainfall-induced landslides using elevation data collected by airbone Lidar.– Natural Hazards and Earth System Sci- ences, 7, 637–650. doi: 10.5194/nhess-7-637-2007 AVANIĆ, R., KOVAČIĆ, M., PAVELIĆ, D. & PEH, Z. (2018a): The Neogene of Hrvat- sko zagorje.– In: TIBLJAŠ, D., HORVAT, M., TOMAŠIĆ, N., MILEUSNIĆ, M. & GRIZELJ, A. (eds.): 9th Mid-European Clay Conference, Conference book, Field Trip Guide book, Zagreb, 128–129. AVANIĆ, R., ŠIMUNIĆ, AN. & PEH, Z. (2018b): Geology of the Croatian Zagorje Re- gion.– In: RMAN, N., MARKOVIĆ, T. & BRENČIĆ, M. (eds.): 5. Slovenski geološki kongres, Post congress field trip book, Ljubljana, 35. BARTELLETTI, C., GIANNECCHINI, R., D’AMATO AVANZI, G., GALANTI, Y. & MAZZALI, A. (2017): The influence of geological–morphological and land use settings on shallow landslides in the Pogliaschina T. basin (northern Apennines, Italy).– Journal of Maps, 13/2, 142–152. doi: 10.1080/17445647.2017.1279082 BELL, F.G. (2003.): Geological Hazards: Their Assessment, Avoidance and Mitigation.– Spoon Press, London, 648 p. BELL, R., PETSCHKO, H., RÖHRS, M. & DIX, A. (2012): Assessment of landslide age, landslide persistence and human impact using airborne laser scanning digital terrain models.– Geografiska Annaler: Series A, Physical Geography, 94, 135–156. doi: 10.1111/j.1468‐0459.2012.00454.x BERNAT GAZIBARA, S., KRKAČ, M. & MIHALIĆ ARBANAS, S. (2019): Landslide inventory mapping using LiDAR data in the City of Zagreb (Croatia).– Journal of Maps, 15/2, 773–779. doi: 10.1080/17445647.2019.1671906 BOGNAR, A. (1996): Tipovi klizišta u Republici Hrvatskoj i Republici Bosni i Herce- govini – geomorfološki i geoekološki aspekti. [The main types of landslides in the Republic of Croatia and Republic of Bosnia and Hercegovina – geomorphological and landscape ecological aspects – in Croatian].– Acta Geographica Croatica, 31, 27–39. BOGNAR, A. (2001): Geomorfološka regionalizacija Hrvatske [Geomorphological re- gionalisation of Croatia – in Croatian].– Acta Geographica Croatica, 34, 7–29. BPŽ, BRODSKO-POSAVSKA ŽUPANIJA (2019): Procjena rizika od velikih nesreća: Identifikacija, analiza, vrednovanje i obrada rizika od velikih nesreća za područje Brodsko-posavske županije [Risk assessment of major accidents: Identification, analysis, evaluation, and treatment of the major accident risk in the area of Brod- Posavina County – in Croatian]. –Brodsko-posavska županija, 301 p. BURNS, W.J. & MADIN, I.P. (2009): Protocol for inventory mapping of landslide de- posits from Light Detection and Ranging (LIDAR) imagery.– Special paper 42, Oregon Department of Geology and Mineral Industries, 30 p. CHACON, J., IRIGARAY, C., FERNANDEZ, T. & EL HAMDOUNI, R. (2006): Engi- neering geology maps: landslides and geographical information systems.– Bulletin of Engineering Geology and the Environment, 65, 341–411. doi: 10.1007/s10064- 006-0064-z COCO, L. & BUCCOLINI, M. (2015): The Effect of Morphometry, Land-Use and Li- thology on Landslides Susceptibility.– In: SCHWECKENDIEK, T., VAN TOL, A.F., PEREBOOM, D., VAN STAVEREN, M.TH. & COOLS, P.M.C.B.M. (eds): Geotechnical Safety and Risk V, IOS Press BV, 779–784. doi: 10.3233/978-1- 61499-580-7-779 CONFORTI, M. & IETTO, F. (2020): Influence of Tectonics and Morphometric Features on the Landslide Distribution: A Case Study from the Mesima Basin (Calabria, South Italy).– Journal of Earth Science, 31, 393–409. doi: 10.1007/s12583-019- 1231-z CROZIER, M.J. (2010): Landslide geomorphology: An argument for recognition, with examples from New Zeland.– Geomorphology, 120/1–2, 3–15. doi: 10.1016/j.geo- morph.2009.09.010 ČUBRILOVIĆ, P., PALAVESTRIĆ, LJ. & NIKOLIĆ, T (1967): Inženjerskogeološka karta SFR Jugoslavije u mjerilu 1:500000 [Engineering Geological Map of SFRY in scale 1:500.000 – in Croatian].– Savezni geološki zavod, Beograd. DEMEK, J. & IGU (1972): Manual of Detailed Geomorphological Mapping.– IGU Co- mission on Geomorphological survey and mapping, Chechoslovak Academy of Science, Prague, 344 p. DERRON, M.-H. & JABOYEDOFF, M. (2010): Preface “LIDAR and DEM techniques for landslides monitoring and characterization”.– Natural Hazards and Earth Sys- tem Sciences, 10, 1877–1879. doi: 10.5194/nhess-10-1877-2010 FELL, R., COROMINAS, J., BONNARD, C., CASCINI, L., LEROI, E. & SAVAGE, W.Z. (2008): Guidelines for landslide susceptibility, hazard and risk zoning for land use planning.– Engineering geology, 102/3–4, 85–98. doi: 10.1016/j.eng- geo.2008.03.022 FILJAK, R., PIKIJA, M., AVANIĆ, R., BAKRAČ, K. & MIKNIĆ, M. (2016): Osnov- na geološka karta Republike Hrvatske mjerila 1 : 50000 – List Slavonska Požega 4 [Basic Geological Map of Republic of Croatia 1:50000, Slavonska Požega 4 sheet – in Croatian].– Hrvatski geološki institut, Zavod za geologiju, Zagreb. GÖRÜM, T. (2019): Landslide recognition and mapping in a mixed forest environment from airborne LiDAR data.– Engineering Geology, 258, 105155. doi: 10.1016/j. enggeo.2019.105155 GUZZETTI, F., ARDIZZONE, F., CARDINALI, M., GALLI, M., REICHENBACH, P. & ROSSI, M. (2008): Distribution of landslides in the Upper Tiber River basin, central Italy.– Geomorphology, 96, 105–122. doi: 10.1016/j.geomorph.2007.07.015 GUZZETTI, F., MONDINI, A.C., CARDINALI, M., FIORUCCI, F., SANTANGELO, M. & CHANG, K.-T. (2012): Landslide inventory maps: New tools for an old problem.– Earth-Science Reviews, 112, 42–66. doi: 10.1016/j.earsci- rev.2012.02.001 HAQUE, U., BLUM, P., DA SILVA, P.F., ANDERSEN, P., PILZ, J., CHALOV, S.R., MALET, J.-P., AUFLIČ, M.J., ANDRES, N., POYIADJI, E., LAMAS, P.C., ZHANG, W., PESHEVSKI, I., PÉTURSSON, T.K., DOBREV, N., GARCÍA- DAVALILLO, J.C., HALKIA, M., FERRI, S., GAPRINDASHVILI, G., ENG- STRÖM, J. & KEELLINGS, D. (2016): Fatal landslides in Europe.– Landslides, 13/6, 1545–1554. doi: 10.1007/s10346-016-0689-3. HORVÁTH, F. & ROYDEN, L.H. (1981): Mechanism for the Formation of the Intra- Carpathian Basins: A Review.– Earth-Science Reviews, 3–4, 307–316. JABOYEDOFF, M., OPPIKOFER, T., ABELLÁN, A., DERRON, M.-H., LOYE, A., METZGER, R. & PEDRAZZINI, A. (2012): Use of LIDAR in landslide investi- gations: a review.– Natural Hazards, 61, 5–28. doi: 10.1007/s11069-010-9634-2 JACOBS, L., DEWITTE, O., POESEN, J., MAES, J., MERTENS, K., SEKAJUGO, J. & KERVYN, M. (2016): Landslide characteristics and spatial distribution in the Rwenzori Mountains, Uganda.– Journal of African Earth Sciences, 134, 917–930. doi: 10.1016/j.jafrearsci.2016.05.013 JAGODNIK, P., JAGODNIK, V., ARBANAS, Ž. & MIHALIĆ ARBANAS, Ž. (2020): Landslide types in the Slani Potok gully, Croatia.– Geologia Croatica, 73/1, 13–28. doi: 10.4154/gc.2020.04 KEATON, J.R. & DEGRAFF, J.V. (1996): Surface observation and geologic mapping.- In: TURNER, A.K. & SCHUSTER, R.L. (eds.): Landslides – Investigation and Mitigation, National Research Council, Washington, DC, 178–230. KOVAČIĆ, M., HORVAT, M., PIKIJA, M. & SLOVENEC, Da. (2011): Composition and provenance of Neogene sedimentary rocks of Dilj gora Mt. (south Pannonian Basin, Croatia).– Geologia Croatica, 64/2, 121–132. doi: https://doi.org/10.4154/ GC.2011.10 KRAUS, D. (2014): An open-source add-in for the Microsoft Excel® spreadsheet soft- ware. Downloaded from https://www.xltoolbox.net/ (21.09.2020) LUPIANO, V., RAGO, V., TERRANOVA, O.G. & IOVINE, G. (2019): Landslide in- ventory and main geomorphological features affecting slope stability in the Picen- tino river basin (Campania, southern Italy).– Journal of Maps, 15/2, 131–141. doi: 10.1080/17445647.2018.1563836 MALAMUD, B.D., TURCOTTE, D.L., GUZZETTI, F. & REICHENBACH, P. (2004): Landslide inventories and their statistical properties.– Earth Surface Processes and Landforms, 29, 687–711. doi: 10.1002/esp.1064 MCCALPIN, J. (1984): Preliminary age classification of landslides for inventory map- ping.– Proceedings of the Annual Symposium on Engineering Geology and Soil Engineering, 21, Moscow, Idaho, 99–120. MIHALIĆ ARBANAS, S., KRKAČ, M. & BERNAT, S. (2016): Application of advanced technologies in landslide research in the area of the City of Zagreb (Croatia, Eu- rope).– Geologia Croatica, 69/2, 231–243. doi: 10.4154/gc.2016.18 MIHALIĆ ARBANAS, S., SEČANJ, M., BERNAT GAZIBARA, S., KRKAČ, M., BEGIĆ, H., DŽINO, A., ZEKAN, S. & ARBANAS, Ž. (2017): Landslides in the G eo lo gi a C ro at ic a 16 Geologia Croatica 75/1Geologia Croatica 74/2 Dinarides and Pannonian Basin – from the largest historical and recent landslides in Croatia to catastrophic landslides caused by Cyclone Tamara (2014) in Bosnia and Herzegovina.– Landslides, 14, 1861–1876. doi: 10.1007/s10346-017-0880-1 PAVELIĆ, D. & KOVAČIĆ, M. (2018): Sedimentology and stratigraphy of the Neogene rift-type North Croatian Basin (Pannonian Basin System, Croatia): A review.– Ma- rine and Petroleum Geology, 91, 455–469. doi: 10.1016/j.marpetgeo.2018.01.026 PAWLUSZEK, K. (2019): Landslide features identification and morphology investiga- tion using high-resolution DEM derivatives.– Natural Hazards, 96, 311–330. doi: 10.1007/s11069-018-3543-1 PETSCHKO, H., BELL, R. & GLADE, T. (2014): Relative Age Estimation at Landslide Mapping on LiDAR Derivatives: Revealing the Applicability of Land Cover Data in Statistical Susceptibility Modelling.– In: SASSA, K., CANUTI, P. & YIN, Y. (eds.): Landslide Science for a Safer Geoenvironment, Springer, 337–343. doi: 10.1007/978-3-319-05050-8_53 PETSCHKO, H., BELL, R. & GLADE, T. (2016): Effectiveness of visually analyzing LiDAR DTM derivatives for earth and debris slide inventory mapping for statisti- cal susceptibility modelling.– Landslides, 13, 857–872. doi: 10.1007/s10346-015- 0622-1 REICHENBACH, P., ROSSI, M., MALAMUD, B.D., MIHIR, M. & GUZZETI, F. (2018): A review of statistically-based landslide susceptibility models.– Earth- Science Reviews, 180, 60–91. doi: 10.1016/j.earscirev.2018.03.001 ROYDEN, L.H. (1988): Late Cenozoic Tectonics of the Pannonian Basin System.– In: ROYDEN, L.H. & HORVATH, I.F. (eds.): The Pannonian Basin: A study in Basin Evolution. AAPG Memoir, 45, 27–48. doi: 10.1306/M45474C3 SCAIONI, M., LONGONI, L., MELILLO, V. & PAPINI, M. (2014): Remote Sensing for Landslide Investigations: An Overview of Recent Achievements and Perspec- tives.– Remote Sensing, 6/10, 9600–9652. doi: 10.3390/rs6109600 SCHULZ, W.H. (2004): Landslides mapped using LIDAR imagery, Seattle, Washington. U.S. Geological Survey Open–File Report 2004–1396, 11 p. SOETERS, R. & VAN WESTEN, C.J. (1996): Slope instability recognition, analysis, and zonation.– In: TURNER, A.K. & SCHUSTER, R.L. (eds.): Landslides: Inves- tigation and Mitigation. Transportation Research Board, 247, National Academy Press: Washington DC, 129–177. STRAHLER, A.N. (1956): Quantitative Slope Analysis.– Geological Society of Ameri- ca Bulletin, 67, 571–596. doi: 10.1130/0016-7606(1956)67%5B571:QSA%5D2. 0.CO;2 ŠESTANOVIĆ, S. (2001): Osnove geologije i petrografije [Fundamentals of geology and petrography – in Croatian].– Sveučilište u Splitu, Građevinski fakultet, 236 p. ŠPARICA, M., JURIŠA, J., CRNKO, J., ŠIMUNIĆ, AN., JOVANOVIĆ, Č. & ŽIVANOVIĆ, D. (1979a): Osnovna geološka karta SFRJ 1:100.000, list Nova Ka- pela [Basic Geological Map of SFRY 1:100.000, Nova Kapela sheet – in Croa- tian].– Institut za geološka istraživanja, Zagreb, Geološki zavod, Sarajevo, Savez- ni geološki zavod, Beograd. ŠPARICA, M., JURIŠA, J., CRNKO, J., ŠIMUNIĆ, AN., JOVANOVIĆ, Č. & ŽIVANOVIĆ, D. (1979b): Osnovna geološka karta SFRJ 1:100.000. Tumač za list Nova Kapela [Basic Geological Map of SFRY 1:100000, Geology of the Nova Kapela sheet – in Croatian].– Institut za geološka istraživanja, Zagreb, Geološki zavod, Sarajevo, Savezni geološki zavod, Beograd. ŠPARICA, M. (1986): Osnovna geološka karta SFRJ 1:100.000, list Slavonski Brod [Basic Geological Map of SFRY 1:100.000, Slavonski Brod sheet – in Croatian].– Geološki zavod, Zagreb, Savezni geološki zavod, Beograd. ŠPARICA, M., BUZALJKO, R. & MOJIĆEVIĆ, M. (1987): Osnovna geološka karta SFRJ 1:100.000. Tumač za list Slavonski Brod [Basic Geological Map of SFRY 1:100.000, Geology of the Nova Kapela sheet – in Croatian].- Institut za geološka istraživanja, Zagreb, Geološki zavod, Sarajevo, Savezni geološki zavod, Beograd. TAROLLI, P. (2014): High-resolution topography for understanding Earth surface pro- cesses: opportunities and challenges.– Geomorphology, 216, 295–312. doi: 10.1016/j.geomorph.2014.03.008 TIAN, Y., XU, C., CHEN, J., ZHOU, Q. & SHEN, L. (2017): Geometrical characteris- tics of earthquake-induced landslides and correlations with control factors: a case study of the 2013 Minxian, Gansu, China, Mw 5.9 event.– Landslides, 14, 1915– 1927. doi: 10.1007/s10346-017-0835-6 TOŠEVSKI, A. (2013): Inventar klizišta u slivu rijeke Dubračine [Slide inventory in Du- bracina river basin – in Croatian].– Rudarsko-geološko-naftni zbornik, 27, 1–16. VAN DEN EECKHAUT, M., POESEN, J., VERSTRAETEN, G., VANACKER, V., NYSSEN, J., MOEYERSONS, J., VAN BEEK, L.P.H. & VANDEKERCKHOVE, L. (2007): Use of LIDAR-derived images for mapping old landslides under for- est.– Earth Surface Processes and Landforms, 32, 754–769. doi: 10.1002/esp.1417 VAN WESTEN, C.J. (2000): The modelling of landslide hazards using GIS.– Surveys in Geophysics, 21, 241–255. VAN WESTEN, C.J., CASTELLANOS, E. & KURIAKOSE, S.L. (2008): Spatial data for landslide susceptibility, hazard, and vulnerability assessment: An overview.– Engineering Geology, 102/3–4, 112–131. doi: 10.1016/j.enggeo.2008.03.010 VARNES, D.J. & IAEG (1984): Landslide hazard zonation: a review of principles and practice.– IAEG Commission on Landslides and Other Mass Movements, Natural Hazards 3, UNESCO, Paris, 63 p. ZANINOVIĆ, K., GAJIĆ-ČAPKA, M., PERČEC TADIĆ, M., VUČETIĆ, M., MILKOVIĆ, J., BAJIĆ, A., CINDRIĆ, K., CVITAN, L., KATUŠIN, Z., KAUČIĆ, D., LIKSO, T., LONČAR, E., LONČAR, Ž., MIHAJLOVIĆ, D., PANDŽIĆ, K., PATARČIĆ, M., SRNEC, L. & VUČETIĆ, V. (2008): Klimatski atlas Hrvatske: 1961-1990, 1971-2000 [Climate atlas of Croatia: 1961-1990, 1971-2000 – in Cro- atian].– Državni hidrometeorološki zavod, Zagreb, 200 p.