Detection of Cliff Top Erosion Drivers through Machine Learning Algorithms between Portonovo and Trave Cliffs (Ancona, Italy)

Rocky coastlines are characterised by steep cliffs, which frequently experience a variety of natural processes that often exhibit intricate interdependencies, such as rainfall, ice and water run-off, and marine actions. The advent of high temporal and spatial resolution data, that can be acquired through remote sensing and geomatics techniques, has facilitated the safe exploration of otherwise inaccessible areas. The datasets that can be gathered from these techniques, typically combined with data from fieldwork, can subsequently undergo analyses employing/applying machine learning algorithms and/or numerical modeling, in order to identify/discern the predominant influencing factors affecting cliff top erosion. This study focuses on a specific case situated at the Conero promontory of the Adriatic Sea in the Marche region. The research methodology entails several steps. Initially, the morphological, geological and geomechanical characteristics of the areas were determined through unmanned aerial vehicle (UAV) and conventional geological/geomechanical surveys. Subsequently, cliff top retreat was determined within a GIS environment by comparing orthophotos taken in 1978 and 2022 using the DSAS tool (Digital Shoreline Analysis System), highlighting cliff top retreat up to 50 m in some sectors. Further analysis was conducted via the use of two Machine Learning (ML) algorithms, namely Random Forest (RF) and eXtreme Gradient Boosting (XGB). The Mean Decrease in Impurity (MDI) methodology was employed to assess the significance of each factor. Both algorithms yielded congruent results, emphasising that cliff top erosion rates are primarily influenced by slope height. Finally, a validation of the ML algorithm results was conducted using 2D Limit Equilibrium Method (LEM) codes. Ten sections extracted from the sector experiencing the most substantial cliff top retreat, as identified by DSAS, were utilised for 2D LEM analysis. Factor of Safety (FS) values were identified and compared with the cliff height of each section. The results from the 2D LEM analyses corroborated the outputs of the ML algorithms, showing a strong correlation between the slope instability and slope height (R2 of 0.84), with FS decreasing with slope height.