The study of land cover and land use dynamics are fundamental to understanding the radical changes that human activity is causing locally and globally and to analyse the continuous metamorphosis of landscape. In Europe, the Copernicus Program offers numerous territorial monitoring tools to users and decision makers, such as Sentinel data. This research aims at developing and implementing a land cover mapping and change detection methodology through the classification of Copernicus Sentinel-1 and Sentinel-2 satellite data. The goal is to create a versatile and economically sustainable algorithm capable of rapidly processing large amounts of data, allowing the creation of national-scale products with high spatial resolution and update frequency for operational purposes. Great attention was paid to compatibility with the main activities planned in the near future at the national and European level. In this sense, a land cover classification system consistent with the European specifications of the EAGLE group has been adopted. The methodology involves the definition of distinct sets of decision rules for each of the land cover macro-classes and for the land cover change classes. The classification refers to pixels’ spectral and backscatter characteristics, exploiting the main multi-temporal indices while proposing two new ones: the NDCI to distinguish between broad-leaved and needle-leaved trees, and the Burned Index (BI) to identify burned areas. This activity allowed for the production of a land cover map for 2018 and the change detection related to forest disturbances and land consumption for 2017–2018, reaching an overall accuracy of 83%.

Soil sealing is the destruction or covering of natural soils by totally or partially impermeable artificial material. ISPRA (Italian Institute for Environmental Protection Research) uses different remote sensing techniques to monitor this process and updates yearly a national-scale soil sealing map of Italy. In this work, for the first time, we tried to combine soil sealing indicators as additional parameters within a landslide susceptibility assessment. Four new parameters were derived from the raw soil sealing map: Soil sealing aggregation (percentage of sealed soil within each mapping unit), soil sealing (categorical variable expressing if a mapping unit is mainly natural or sealed), urbanization (categorical variable subdividing each unit into natural, semi-urbanized, or urbanized), and roads (expressing the road network disturbance). These parameters were integrated with a set of well-established explanatory variables in a random forest landslide susceptibility model and different configurations were tested: Without the proposed soil-sealing-derived variables, with all of them contemporarily, and with each of them separately. Results were compared in terms of AUC ((area under receiver operating characteristics curve, expressing the overall effectiveness of each configuration) and out-of-bag-error (estimating the relative importance of each variable). We found that the parameter “soil sealing aggregation” significantly enhanced the model performances. The results highlight the potential relevance of using soil sealing maps on landslide hazard assessment procedures.

This study focuses on the July-August 2019 eruption-induced wildfires at the Stromboli island (Italy). The analysis of land cover (LC) and land use (LU) changes has been crucial to describe the environmental impacts concerning endemic vegetation loss, damages to agricultural heritage, and transformations to landscape patterns. Moreover, a survey was useful to collect eyewitness accounts aimed to define the LU and to obtain detailed information about eruption-induced damages. Detection of burnt areas was based on PLÉIADES-1 and Sentinel-2 satellite imagery, and field surveys. Normalized Burn Ratio (NBR) and Relativized Burn Ratio (RBR) allowed mapping areas impacted by fires. LC and LU classification involved the detection of new classes, following the environmental units of landscape, being the result of the intersection between CORINE Land Cover project (CLC) and local landscape patterns. The results of multi-temporal comparison show that fire-damaged areas amount to 39% of the total area of the island, mainly affecting agricultural and semi-natural vegetated areas, being composed by endemic Aeolian species and abandoned olive trees that were cultivated by exploiting terraces up to high altitudes. LC and LU analysis has shown the strong correlation between land use management, wildfire severity, and eruption-induced damages on the island.

Land consumption is the increase in artificial land cover, which is a major issue for environmental sustainability. In Italy, the Italian Institute for Environmental Protection and Research (ISPRA) and National System for Environmental Protection (SNPA) have the institutional duty to monitor land consumption yearly, through the photointerpretation of high-resolution images. This study intends to develop a methodology in order to produce maps of land consumption, by the use of the semi-automatic classification of multitemporal images, to reduce the effort of photointerpretation in detecting real changes. The developed methodology uses vegetation indices calculated over time series of images and decision rules. Three variants of the methodology were applied to detect the changes that occurred in Italy between the years 2018 and 2019, and the results were validated using ISPRA official data. The results show that the produced maps include large commission errors, but thanks to the developed methodology, the area to be photointerpreted was reduced to 7300 km² (2.4% of Italian surface). The third variant of the methodology provided the highest detection of changes: 70.4% of the changes larger than 100 m² (the pixel size) and over 84.0% of changes above 500 m². Omissions are mainly related to single pixel changes, while larger changes are detected by at least one pixel in most of the cases. In conclusion, the developed methodology can improve the detection of land consumption, focusing photointerpretation work over selected areas detected automatically.

Geological maps convey different and multifaceted information including lithology, age, tectonism and so on. This complex information is not fully exploited in landslide susceptibility (LS) studies, as a single parameter is usually derived from the geological map of the study area (e.g. the area is divided into lithological or lithostratigraphic or geological units). The aim of this work is testing different approaches to extract significant information from geological maps, creating different parameterizations, and analyzing the sensitivity of a LS model to these variations.Our test site is a 3100 km2 wide area in Tuscany (Italy) characterized by a very complex geological setting. A 1:10000 scale geological map subdivides the area into 194 different lithostratigraphic units. This map was reclassified according to different criteria, creating 6 different parameters derived from the same geological map: lithology (6 lithological classes), age of deposition (the area was subdivided into 6 chronological units), paleogeography (6 units were differentiated on the basis of their environment of formation), genesis of the bedrock (5 classes accounted for the mechanism of formation of the outcropping rock/terrain), broad tectonic domain (the mapped elements were grouped into 5 broad structural units accounting for their tectonic history), detailed tectonic domain (same as before but with a more detailed subdivision into 10 classes).Some of these parameters have already been used in LS studies, others have been used here for the first time; however, all of them have some connections with landslide predisposition. These parameters were used (one by one and altogether) to run seven times a landslide susceptibility model based on the widely used random forest machine learning algorithm. The model configurations and resulting maps were evaluated in terms of AUC(Area Under Curve) and OOBE(out of bag error): while the former expresses the forecasting effectiveness of each configuration, the latter expresses, among a single configuration, the importance of each input parameter.We discovered that the results are very sensitive to the approach used to consider geology in the susceptibility assessment, with AUC values ranging from 63.5% (using chronological units) to 70.0% (using genetic units) and 75.2% (using all the geology-derived parameters simultaneously). These results are in line with OOBE statistics, which showed a similar relative importance of the geologically-driven parameters.These outcomes can to assist future landslide susceptibility studies for different reasons:(i)at least in our study area, lithology, which is commonly used in LS, did not provide the best results;(ii)as geological maps provide multifaceted information, a single classification approach cannot fully grasp this complexity; therefore, the best results can be obtained using different geology-based parameters simultaneously, because each of them can account for specific features connected to landslide predisposition (to our knowledge, a similar approach has never been attempted before in LS literature).(iii)When using thematic maps to feed LS models, it is important to fully understand the nature and the meaning of the information provided by the geology-related maps: results are very sensitive to this kind of information and the interpretation of the results should take it into account.