Detailed and accurate geomorphological historical landslide inventory maps
are an invaluable source of information for many research topics and
applications. Their systematic preparation worldwide has been advised by
many researchers as it may foster our knowledge on landslides, their spatial
and temporal distribution, their potential interaction with the built
environment, their contribution to landscape dynamics, and their response to climate change in the past. Due to the extreme variability of the
morphological and radiometric elements that can reveal historical
landslides, geomorphological historical landslide inventory maps are
produced by expert interpretation, which makes it a time-consuming and expensive process, which often discourages wide-area mapping activities. In this paper we present a new geomorphological historical landslide inventory
map for a 1460 km
Landslides are widespread natural hazards that occur worldwide, posing threats to populations, structures and infrastructures and causing relevant economic damage to society (Donnini et al., 2017; Froude and Petley, 2018; Petley, 2012). Availability of information on landslide locations, types, and sizes is key basic information for many activities, from research to emergency support and land planning. Nevertheless, expert landslide mapping is not an investment priority in many countries and institutions, and hence homogeneous wide-scale landslide maps are rare (Guzzetti et al., 2012; Bucci et al., 2021), although the availability of a landslide inventory map is fundamental to the production of landslide hazard and risk maps (Thiery et al., 2020). As a result, several applications spanning from land management to susceptibility modelling that would benefit from more homogeneous and complete datasets rely on often poor or incomplete/inconsistent data.
Any landslide mapping activity assumes (i) that landslides leave discernible features on the territory, i.e. a “landslide signature” made by radiometric and morphologic elements (Fiorucci et al., 2018; Guzzetti et al., 2012; Santangelo et al., 2022), (ii) that, in the case of use of remote sensing, the data used (e.g. images) portray them, (iii) that the technique adopted is adequate to the data used, and (iv) that the operator (or the automatic or semi-automatic classification system) is able to exploit those features to detect and map landslides. Therefore, any landslide mapping activity is affected by a degree of incompleteness and inconsistency that stems from these basic assumptions. Landslides can remain unnoticed since (i) features may be cancelled by erosion, other landslides, and/or human activity (Malamud et al., 2004), (ii) the type of image or its spectral and spatial resolution may be inadequate for some types or sizes of landslides or for some locations (e.g. slope orientation), and (iii) interpreters or (semi-)automatic systems may be unable to map some landslides using certain images or techniques (e.g. lack of experience, poor or incorrect system training).
Landslide inventory maps (LIMs) are the simplest tools to represent landslide distribution in a territory. They usually display landslide locations preferably as polygons rather than points or lines, which should be used only for scale representation issues. LIMs usually store landslide attributes such as type of movement, estimated time of occurrence, estimated depth, activity at the date of the observation if available, and other more detailed data depending on availability and on the scale of the work (Bucci et al., 2021; Guzzetti et al., 2012; Santangelo et al., 2015).
Different types of inventories can be prepared. Event inventories (E-LIMs, Antonini et al., 2002; Ardizzone et al., 2012; Donnini et al., 2017; Fiorucci et al., 2018; Harp and Jibson, 2017; Mondini et al., 2012; Santangelo et al., 2022) report landslides that were triggered by specific events (i.e. rainfall, rapid snowmelt, volcanic eruption, earthquake). Geomorphological inventories (G-LIMs) report landslides that can be recognised by geomorphologists, usually from the expert interpretation of stereoscopic aerial photographs (Bucci et al., 2021; Cardinali et al., 2001; Guzzetti et al., 2012; Santangelo et al., 2015), but lidar-derived images are also widely used (Niculiţă et al., 2016; Razak et al., 2013; Schulz, 2004; Van Den Eeckhaut et al., 2007). They can be seen as the result of many landslide events over thousands of years (Malamud et al., 2004). Multi-temporal landslide inventories (M-LIMs) include the information of G-LIMs and also report landslides that occurred in the last tens of years as visible on historical images available at (possibly) regular time steps of several years (Galli et al., 2008; Guzzetti et al., 2012; Zumpano et al., 2020).
G-LIMs can be prepared for wide areas and provide a fundamental source of information about landslides that occurred in the last tens of thousands of years. This particular type of inventory is an invaluable source of data but is time-consuming and requires a high level of training and experience to get accurate and reliable results (Guzzetti et al., 2012), which usually hampers the systematic preparation of G-LIMs at regional or even national/continent scales.
This paper presents a new G-LIM prepared for the Daunia mountains in the Apulia (Puglia) region, southern Italy. For this area, historically affected by diffuse slow-moving large landslides, the regional government required the preparation of
a new landslide inventory to overcome the inhomogeneity among the existing
inventories in terms of the spatial distribution of landslides and of techniques adopted and working and publication scales. The landslide inventory map was
produced through the interpretation of two sets of stereoscopic aerial
photographs, taken in 1954/55 and 2003, supplemented by targeted field
checks. The new landslide inventory is an entirely original dataset openly
accessible at
The study area extends for 1460 km
The southern Apennines are the result of compressive tectonic dynamics characterised by an oblique collision of the Calabrian forearc with the Apulian margin (Filice and Seeber, 2019) that evolved from an antecedent subductive phase of thrust emplacement (Vitale et al., 2011; Vitale and Ciarcia, 2013). Currently, the external (eastern) margin preserves the compressive tectonics, while the internal sector is characterised by an extensive dynamic with tension faults that dissect the pre-existent fold and thrust structure (Brozzetti et al., 2009; Schiattarella et al., 2003).
The western sector of the Daunia Apennines is characterised by a wide variety of formations, mainly clay-rich flysch lithologies with different mechanical properties highly affected by folds and faults (Cotecchia et al., 2020; Losacco et al., 2021; Pellicani et al., 2014b), resulting in medium to high relief, with the elevations ranging between 55 and 1152 m a.s.l., where the highest peak is represented by Monte Cornacchia. Morphologically, it is strongly influenced by both the prolonged tectonic stresses and by the differences in composition and erodibility, where selective erosion has alternatively produced gentler and rounded slopes on the more erodible clayey formations and steep slopes cut on harder rocks. In this sector, the orography associated with the widespread presence of clay-rich materials and the intense deformation is the main cause of landslides (Ciarcia et al., 2003; Losacco et al., 2021; Spalluto et al., 2021; Wasowski et al., 2010; Zumpano et al., 2020). The eastern sector is characterised by the sub-Apennine clay formation and alluvial deposits and by gentle slopes slightly dipping towards the north-east. In this sector, several orders of terraced fluvial deposits overlay the sub-Apennine clay formation. Here, mass movements are mainly concentrated along the scarps of the terraces.
The Daunia Apennines are influenced by a Mediterranean sub-humid climate with mild and often wet winters and usually hot and dry summers. In general, inter-annual temperature variations are significant and the total yearly precipitation seldom exceeds 1000 mm (Wasowski et al., 2012, 2010). In the area, torrential streams flow NE–SW, draining the sediments from the Daunia Apennines through the Tavoliere plain towards the Adriatic Sea.
In response to the lithostructural and morphoclimatic setting of this area, landslides are widespread. They are mostly slow-moving slide-type and flow-type movements, involving soil and earth, and secondarily rock material. Rapid-moving (i.e. debris flows) and fast-moving (i.e. rockfalls) landslides are less abundant (Wasowski et al., 2010, 2012; di Lernia et al., 2022; Zumpano et al., 2020). The widespread presence of landslides in the area is a major cause of damage to urban settlements, to the road network, and to the stability of buildings. These features determine a direct effect on the development and the economy of the area (Pellicani et al., 2014b, a; Zumpano et al., 2020).
Location map of the study area. Base map: digital terrain model (DTM) at 8 m resolution made available by the Regione Puglia administration.
The procedure adopted to prepare the landslide inventory map is illustrated in Fig. 2. The overall procedure consists of three main phases: (i) data collection, (ii) photo interpretation, and (iii) editing. The entire work was carried out by a “mapping team” made up of four interpreters and one supervisor and an “editing team” made up of three GIS editors.
The first stage of the work consisted in the collection of data available
for the study area (i.e. ancillary data, Fig. 2; complete list in Appendix A) and the pre-processing of the aerial photographs to produce oriented stereo-models in absolute coordinates to be used for the visual
interpretation in a photogrammetric GIS environment. Such photogrammetric
pre-processing included the interior and exterior orientations of each pair of aerial photographs. For the interior orientation, a non-metric camera
model was adopted for which the focal length of the camera, the flight
altitude, and the pixel size were required, while for the representation of the ground-to-image geometry (exterior orientation), ground control points (GCPs) were used. The planimetric coordinates (
The second phase of the work is the photo interpretation (Fig. 2) of the aerial photographs. At the very initial stage of the work, the entire mapping team defined a legend through an expeditious photo interpretation in representative geomorphological settings and according to well-established landslide classification schemes (Cruden and Varnes, 1996; Hungr et al., 2014; WP/WLI, 1993). Then, divided into pairs, for each stereo-pair the mapping team prepared a preliminary interpretation (first-level G-LIM in Fig. 2) that was then reviewed by the supervisor and discussed with the mapping team. After this step, an updated version of the map was produced (second-level G-LIM in Fig. 2). Each of the two mapping team pairs performed their preliminary photo interpretation in alternating strips along the flight plan. The side lap between subsequent strips was therefore common between the different pairs of interpreters, and it was used to compare, correct, and homogenise the interpretation made by different sub-teams independently. Such a continuous interaction among geomorphologists made it possible to best define the characteristics of the identified landslides, limiting interpretation inconsistencies among operators. Finally, at significant advancement steps, extensive field checks were performed by the entire mapping team to check, validate, and correct the inventory based on field evidence. After the field check stage, the interpretation phase was concluded and the landslide mapping considered quasi-definitive; i.e. occasional changes might be needed if inconsistencies were observed at the third stage (e.g. consistency with contour maps and hydrography).
Flowchart of the procedure to prepare the geomorphological landslide inventory map. Pink rectangles indicate contributions of different teams/operators.
The third and last step of the procedure (editing) consists in the preparation of the geographical database. The interpreted features drawn as
polylines by the mapping teams were verified on the digital topography,
i.e. contour lines at scale
The G-LIM of the Daunia Sub-Apennines was prepared through systematic visual interpretation of a set of 270 b/w stereoscopic aerial photographs taken in 1954/55 at 1 : 33 000 scale and of 384 b/w stereoscopic aerial photographs taken in 2003 at 1 : 30 000 scale. Aerial photographs were provided as 800 dpi scanned copies. All the ancillary data used for interpretation are listed and described in detail in Appendix A.
Descriptive statistics of landslides by relative age.
For the interpretation phase we used the digital stereoscopic vision 3D
PluraView System, composed of two monitors for the stereoscopic vision of
digital images, and a 45
The inventory covers an area of 1460 km
In this section the different groups of elements comprising the inventory will be analytically described. Section 4.1 to 4.3 describe landslides classified according to their relative age, type classification, and estimated depth. Section 4.4 describes the widespread landslides (i.e. landslides represented as groups of failures due to their small size and high spatial frequency). Section 4.5 gives the descriptive statistics of the landslide inventory, and Sect. 4.6 presents the geomorphological elements (i.e. elements that are considered in relation to slope evolution and that can provide useful information for landslide identification and mapping).
Based on their appearance in aerial photographs, landslides were classified according to their relative age, namely (i) very old relict landslides, (ii) very old landslides, (iii) old landslides, and (iv) recent landslides (Fig. 3). The four relative age levels are the result of a major subdivision, referring to the general morphologic appearance of landslides (Keaton and DeGraff, 1996), and are based on the assumption that evidence of landslides becomes less obvious with increasing age, due, for example, to erosion processes, vegetation growth, and occurrence of other landslides. Hence, older landslides are more difficult to detect than more recent ones.
Landslide modifications within the two epochs (1954 and 2003) were negligible for very old relict and very old landslides, whereas they were more frequent for old landslides (i.e. pre-2003). Therefore, the final landslide delineation considers by definition the latest evidence within the 49-year time interval.
In our G-LIM, the most represented landslides in the study area are old
landslides (Table 1, Fig. 3c), corresponding to 14 793 landslides covering
a total landslide area of 339 km
Spatial distribution of
Figure 4b–e show the inventory disaggregated according to the four relative age classes recognised. The sequence of the four panels represents a timeline where landslides of each time step are shown with their own colour, the same as Fig. 4a, whereas landslides belonging to antecedent time steps are represented in white.
In addition to the four mentioned relative age classes, a generation principle based on crosscut relationships defines minor sub-divisions within each age class: younger landslides cover the older ones. For the overall inventory, this criterion allowed us to detect up to (i) two landslide generations within the relict landslides, (ii) two generations within the very old landslides, (iii) four generations within the old landslides, and (iv) two generations within the recent landslides. Such minor sub-divisions are applicable only to landslides that overlap; i.e. they are not applicable to landslides that do not overlap. For instance, in Fig. 4d the age class “old” records up to three overlapping classes, indicated, respectively, in light orange (first failure “old”), orange (second failure “old”), and red (third failure “old”).
Examples of the relative age classes in the area of the Pietramontecorvino municipality (see Fig. 1 for reference).
Landslides classified as very old relict (Fig. 4b) generally affect entire slopes, where the energy of the relief is greatest. They are usually controlled by the geological and lithological structure and the presence of main faults (Bucci et al., 2016a). Very old relict landslides are deeply dismantled by the erosive action of water courses and often reshaped by recurring gravitational phenomena. These morphological modifications are often related to the regional morphological and tectonic evolution, which has determined considerable base-level variations in the study area. As a consequence, these landslides are often suspended with respect to the present-day base level and are totally or partially isolated from the recent evolution of the drainage network. Therefore, very old relict landslides are considered to have occurred under geomorphological and climatic conditions different from the present day (WP/WLI, 1993).
Landslides classified as very old (Fig. 4c) show evidence of active and passive interaction with river dynamics, being often partly eroded from them and locally modifying their path. Also, morphological modifications of very old landslides can be induced by the occurrence of other mass movements over time. Very old landslides are generally large and in agreement with the present-day river network: they are mainly distributed where the local relief is high, and they often occur in or near very old relict landslides as local reactivations.
Landslides classified as old (Fig. 4d) present recognisable evidence on both the 1954/55 and 2003 aerial photograph epochs. These landslides generally present morphological characteristics typical of landslides which are not modified by erosion (e.g. concave–convex shape of the slope, presence of steps, and back slopes in the deposit area). Old landslides tend to cluster spatially, as is evident from the several generations of mass movements that were recognised. Spatial clustering of old landslides is even more evident in the vicinity of very old or very old relict landslides. Old landslides are hypothesised to have occurred in more recent times (as yet undefined) compared to very old and very old relict landslides.
Landslides classified as recent (Fig. 4e) are represented by landslides that show evidence of having occurred close to the date of the 2003 aerial image. The diagnostic elements for the identification of recent landslides refer not only to the morphological evidence, but above all to the spectral elements (photographic tone and contrast), which can be perceived even with a bi-dimensional (2D) view of the images. Recent landslides are mainly shallow failures essentially involving the soil horizons or the alteration of the debris cover for a few metres in depth. Their spatial distribution is ubiquitous and only limited to the areas where the triggering event induced the landslides. Therefore, they were found both in previously unaffected slopes or within pre-existing landslides.
As opposed to the recent landslides, whose evidence is stronger, some uncertainty characterises the delineation of older landslide borders, whose boundaries can be affected or even dismantled by different degrees of erosion or covered by more recent slope failures. General considerations on the relative age criteria of classification (Keaton and DeGraff, 1996) support our own experimental observations on landslide morphological appearance that the uncertainty degrees in landslide mapping increase with the increasing landslide age (Bucci et al., 2021). However, we do not have information on the absolute ages of the landslides in our study area and only hypothesise that the majority of the mapped landslides occurred in the last 10–20 Kyr, as suggested by findings of recent studies on landslide age in the southern Apennines (Gioia et al., 2011) and elsewhere (Niculiţă et al., 2016; Pánek et al., 2014).
The types of landslides are assigned based on the morphological and radiometric signatures of each landslide on the aerial photographs. Landslide types included in the legend of the G-LIM are (i) slide-type landslides (which include deep-seated rock slides or earth slides, shallow soil slides, or earth slides), (ii) earth flow, (iii) slide-earth flow, (iv) debris flow, (v) rockfall, and (vi) sackung (Table 2). In Table 2, widespread landslides refer to polygons representing groups of landslides whose size is smaller than the smallest feature that can be represented at the publication scale of the final map.
Descriptive statistics of landslides by type.
According to the type of movement, the most represented landslides in the
study area are slides (9867 landslides, Fig. 5a) covering a total
landslide area of 287 km
The main features of the three more represented landslide types are illustrated in Fig. 6. Inspection of the figure highlights the main cartographic differences (Fig. 6a, c, e) between the different landslide types (Fig. 6b, d, f).
Figure 6a, b illustrate the slide-type landslides. Such failures present a well-defined scarp which can be semi-circular (rotational slides) or angular (translational slides). The slide deposit is convex, with a morphologically depressed head of the deposit characterised by a centripetal drainage and local counter-slopes and a toe characterised by a typical upward bulge (Fig. 6b). It is reasonable to expect a volume balancing between escarpment and deposit areas since slides are usually characterised by low mobility of landslide material. This is clear in Fig. 6a, b, which shows a planimetric and altimetric shift of about 25 m of the displaced material, with no evidence of chaotic rearrangement and/or volumetric changes.
Slide-earth flows (Fig. 6c, d) start as slides and then evolve into flows. Therefore, they show the characteristics of slides (most commonly rotational) in the escarpment area and in the head of the deposit, where local endorheic conditions can promote seasonal swamps with diffuse organic soil development (Fig. 6c, d). On the other hand, the transport zone and accumulation zone are more similar to earth flows, characterised by the typical lobe shape of the accumulation zone (Fig. 6c, d).
Spatial distribution of
Earth flows (Fig. 6e, f) are characterised by an overall elongated planar shape, with the median part narrower than the detachment and accumulation zones. The mobility of earth flows is generally higher than slides, which is also evident in their more elongated shape. A typical earth flow is usually bounded by lateral streams draining the deposits into the main river valley, which typically experiences narrowing of the valley section and erosion on the opposite side of the landslide toe (Fig. 6f).
Debris flows and rockfalls are not statistically represented in our inventory because the Daunia Apennines lack the environmental conditions favouring their development, such as steep rocky slopes and sub-vertical cliffs. However, local geomorphological conditions (e.g. steep edges of suspended fluvial terraces) can promote the formation of coalescing phenomena of small debris flows and small rockfalls, which cannot be mapped individually but which were included within areas with widespread landslides (Fig. 5d, Sect. 4.4). Finally, sackungs are poorly represented in our study area.
Examples on the map
The estimated depth of a landslide is assigned based on the main morphological characteristics of the landslide, such as slope height, and the extent of curvatures along the slope (convexity and concavity). Based on these morphological characteristics, landslides were classified as deep-seated or shallow.
Deep-seated landslides are predominantly represented by failures that cover
considerable extents, even exceeding 1 km
Shallow landslides are generally small (
The areas of widespread landslides (Table 2, Fig. 5d) are represented by landslides whose sizes are smaller than the smallest feature that can be represented at the publication scale of the final map. In the map we distinguished between (i) shallow earth/soil slides and earth flows, (ii) debris flows, and (iii) rockfalls. Their spatial distribution is mainly related to locally steep slopes, for instance in the vicinity of badlands, in the crown areas of pre-existing large landslides, or along the edges of ancient suspended fluvial terraces (Fig. 7). Widespread shallow slides and flows typically involve weathered clayey soils, while widespread debris flows and rockfalls are commonly fed by loose debris or poorly cemented conglomerates and breccias. Finally, it was not possible to assign a well-defined age to widespread landslides, although their morphological evidence suggests seasonal reactivations.
Descriptive statistics of landslide number and size are shown in Fig. 7. Box plots in Fig. 7a, b show the distribution of landslide areas for the landslide classified based on the expected relative age (Fig. 7a) and the type of movement (Fig. 7b).
Inspection of the plot in Fig. 7a confirms that slides, slide-earth flows, and flows are the most represented landslides of the inventory. Among these three classes, the slide-type landslides are characterised by the largest size variability, while the slide-earth flows are on average slightly larger than the others. In the plot of Fig. 7a, landslides (on the left) are separated by areas of widespread landslides (on the right) because the latter are not related to a single feature and cannot be directly compared to single landslides. This is the reason why the areas of widespread landslides are on average larger than the single landslides. In particular, the areas of widespread soil slides and earth flows are larger than the others because they typically involve edges of fluvial terraces characterised by strong lateral continuity.
The plot in Fig. 7b shows the size distribution of the four relative age classes of landslides. The clear separation of the age classes according to their median values is a statistically robust indication of the reliability of this landslide classification strategy, which is under-used in the landslide-oriented international literature.
Comparing Fig. 7a and b, it is worth noting that the median area of all
the landslide types is around 10
Figure 7c shows landslide count (represented by a colour gradient and labels) and cumulated area (proportional to circle sizes), grouped by relative age class and landslide type. Visual inspection of the plot reveals that relict and very old landslides are not represented by small landslides, and in particular by debris flows and rockfalls, which tend to be small and easily obliterated by erosion and subsequent failures. Further evidence is that a large portion of the total landslide area is occupied by few relict and very old kilometre-scale slope failures, whereas a decreasing trend in both cumulated area and the total number of landslides is evident through time for all landslide types within the old landslide age class. Such evidence suggests a size threshold effect of previous landslides on subsequent slope movements, probably controlled by slope-scale morphological and hydrological perturbations induced by the occurrence of the first failure.
Plots summarising landslide statistics.
In the landslide inventory map, some geomorphological elements are also reported in addition to the mass movements (Fig. 8). Such elements can be considered in relation to slope evolution and can provide useful information for landslide identification and mapping. In particular, our map includes alluvial deposits, fluvial terraces, and alluvial fans. Alluvial deposits are flat and always located in the lowest portions of the Tavoliere delle Puglie plain and along the main water courses draining the hilly and mountainous areas of the Daunia Apennines. Older alluvial/fluvial deposits are nowadays suspended over the present-day base level and are organised in fluvial terraces at different elevations, recording the ongoing deepening of the drainage network during the Quaternary. Alluvial fans were recognised and mapped into two relative age classes according to their appearance. They were classified as relict if dismantled and dissected due to the incision of the present-day river network. They were instead classified as recent if well preserved.
The dataset is available at
The landslide inventory map presented in this paper is a new catalogue of landslide phenomena in the Daunia Apennines. It provides new data uniformly distributed over the study area at an unprecedented cartographic detail.
The new landslide inventory map was produced upon request of the authorities of the Apulia region to solve the heterogeneity and limitations in the currently available landslide inventories (Pellicani et al., 2014b; Pellicani and Spilotro, 2015), mainly due to (i) partial coverage of the territory, (ii) heterogeneity in the degree of cartographic detail, (iii) lack of metadata and other methodological information, and (iv) geospatial inconsistencies where the different products overlap.
The new inventory represents a significant improvement upon earlier mapping
because (i) it presents a new homogeneous landslide coverage for the Daunia Apennines, (ii) it is acquired with a uniform method and based on rigorous and reproducible criteria, (iii) it keeps the same cartographic detail
throughout the entire area, and (iv) it presents a systematic classification
of landslides by type, relative age, and depth, hence providing basic
information at
Analysis of the data revealed that the great majority of landslide volume was mobilised by relict and very old landslides, whereas the landslide maximum size decreased over time and tended to cluster around pre-existing failures. The results have relevance for determining the statistics of landslide size (Malamud et al., 2004) and for triggering detailed studies on the types, patterns, and distributions of landslides in relation to geology (Bucci et al., 2016b), structure (Bucci et al., 2013), and tectonic (Bucci et al., 2014) and climatic forcing (Schiattarella et al., 2017). The presented dataset also documents the relationships between landslides and other geomorphological elements, in particular with ancient and recent alluvial deposits and alluvial fans. These fluvially related landforms can be used as input data to study Quaternary depositional (Mirabella et al., 2018) and erosional events (Mancini et al., 2020) and, in combination with landslide information, to study the interplay of gravitational and fluvial processes (Santangelo et al., 2013).
In the G-LIM presented, the vast majority of landslides are slow-moving, whereas rapid and fast landslides are rare (debris flows and rockfalls). Such evidence should be further confirmed by in-depth studies (e.g. multi-temporal inventories or susceptibility models) as it may indicate a relatively low exposure of human life compared to other mountain areas of the southern Apennines. On the other hand, it must be taken into account that the increasing frequency of extreme rainfall events may also trigger rapid- and fast-moving landslides in this area, differently from what was observed in the past. This aspect further shows that systematic landslide mapping to prepare new inventories and update the existing ones is crucial, particularly in the ongoing climate change scenario.
We acknowledge that scarcity of data represents a limitation and a challenge for the landslide community. Most regional-scale studies are less effective (or analysed) than expected due to the lack of easy sharing of landslide information in digital format. Our work is also intended as a contribution to the broader geomorphological community, promoting the sharing of a landslide dataset at regional scale from various environments.
Finally, it is worth remembering that a landslide map, even if accurate, does not answer all questions regarding hazards or risk at any scale. A landslide map is informative about the area covered by landslides, but nothing can be
said about the remaining territory. The map dataset should be used as an
indication to optimise resources and to plan investigations aimed at
determining landslide hazard and risk at a larger scale. In our case, the
map dataset is published at the
In the first phase of the work, information and cartographic/thematic products were collected to support the preparation of the geomorphological landslide
inventory map of the Daunia Apennines. Some of the listed products have been
provided by Regione Puglia and others can be consulted on the web through cartographic portals or web servers.
Regional technical maps in digital format at scale Geological maps at scale Geological maps at scale Hydrogeomorphological map at scale Digital elevation model (DEM, cell 8 DEM made by lidar (1 Urbanised and land use maps provided by Regione Puglia ( Orthophoto maps provided by Regione Puglia ( Landslide inventory maps of Apulia from the following:
Official archives: AdBP photo interpretation: landslides from the photo interpretation study
conducted by AdBP (Basin Authority of the Apulia region) Province of Foggia: this deals with the landslides surveyed by the Province of Foggia within the framework of the agreements signed with the AdBP concerning the “Activity of collection and classification of data on
landslide phenomena in the Province of Foggia” and “Conduct of studies for
the deepening of the aspects related to the classified landslide phenomena”.
We declare that all the authors' contributions in preparing both the dataset and the manuscript are equivalent to that of a first author. This is why all the authors appear in alphabetical order.
The contact author has declared that none of the authors has any competing interests.
Publisher's note: Copernicus Publications remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
This work was supported by the Civil Protection of the Apulia region in the framework of the project “Integrated assessment of geohydrological instability phenomena in the Apulia region, interpretative models and definition of rainfall thresholds for landslide triggering” funded by the P.O.R. Puglia 2014-2020, Asse V – Azione 5.1 (project identification no. B82F16003840006).
This paper was edited by Achim A. Beylich and reviewed by Louise M. Vick and one anonymous referee.