Intense rainfall events often produce a great number of shallow landslides events, which in many cases can hits large areas or an entire regional territory. These slope instabilities cause damage to many roads, buildings, and infrastructures and often human loss. In these conditions, it is useful to refine shallow landslides susceptibility maps at regional scale progressively more reliable and efficacy. To take the highlighted goal it is opportune to promote the use of a circular approach that can considers knowledge (data, methods, models, solutions, etc.) constantly upgraded. To achieve this aims we propose a method that introduces structurally in a possible circular approach (progressive better results with constantly upgraded knowledge) the use of a comprehensive geo-database of shallow landslide events and related implemented through a collection and analysis of numerous sources, including published inventory maps, scientific literature, technical reports and newspapers, integrated by a multi-temporal interpretation of remote sensing images and several field surveys. The method is applied referring to the Calabria region, which is largely affected by this landslide category. The refined geo-database realized includes 22,028 shallow landslides, occurred between 1951 and 2017. The relationship between spatial pattern of the shallow landslides and the analyzed predisposing factors (lithological units, fault density, land use, drainage density, slope gradient, TWI, SPI and LS) showed that the high values of slope gradient, LS factor and drainage density, coupled to low values of TWI, displayed a strong control on the shallow landslide occurrence. The efficacy of the geo-database realization proves their usefulness in order to estimate and validate shallow landslide susceptibility map, which was optimally obtained applied a simple bivariate statistical method. The susceptibility map was classified into five classes and about 26% of the study area falls in high and very high susceptible classes and most of the shallow landslides mapped (76%) occur in the same classes. The AUC value of the prediction rate curve was 0.81, indicating a good prediction capability of the susceptibility map. The interaction between shallow landslide susceptibility map and road network map highlighted that the 20% of the roadways of the region area falls in high and very high ARTICLE HISTORY