Articles | Volume 23, issue 4
Research article
14 Apr 2023
Research article |  | 14 Apr 2023

PERL: a dataset of geotechnical, geophysical, and hydrogeological parameters for earthquake-induced hazards assessment in Terre del Reno (Emilia-Romagna, Italy)

Chiara Varone, Gianluca Carbone, Anna Baris, Maria Chiara Caciolli, Stefania Fabozzi, Carolina Fortunato, Iolanda Gaudiosi, Silvia Giallini, Marco Mancini, Luca Paolella, Maurizio Simionato, Pietro Sirianni, Rose Line Spacagna, Francesco Stigliano, Daniel Tentori, Luca Martelli, Giuseppe Modoni, and Massimiliano Moscatelli

In 2012, the Emilia-Romagna region (Italy) was struck by a seismic crisis characterized by two main shocks (ML 5.9 and 5.8) which triggered relevant liquefaction events. Terre del Reno is one of the municipalities that experienced the most extensive liquefaction effects due to its complex geostratigraphic and geomorphological setting. This area is indeed located in a floodplain characterized by lenticular fluvial channel bodies associated with crevasse and levee clay–sand alternations, related to the paleo-Reno River. Therefore, it was chosen as a case study for the PERL project, which aims to define a new integrated methodology to assess the liquefaction susceptibility in complex stratigraphic conditions through a multi-level approach. To this aim, about 1800 geotechnical, geophysical, and hydrogeological investigations from previous studies and new realization surveys were collected and stored in the PERL dataset. This dataset is here publicly disclosed, and some possible applications are reported to highlight its potential.

1 Introduction

In the last few years, an increasing number of source data is publicly disclosed, allowing for wider access to research activities. Key examples are the huge amount of free satellite imagery (i.e., Sentinel, Landsat) provided by the main space agencies and the cutting-edge tools and procedures integrated in widely known and open-source EO platforms such as Google Engine. A multitude of algorithms and codes are available for all the fields of knowledge concerning natural hazards, while their application is made easier by the increasing number of open-access inventories of natural phenomena (i.e., Martino et al., 2014; Guarino et al., 2018; Tanyaş et al., 2022). However, only a few examples of datasets of in situ investigations and related parameters are publicly disclosed, and this is a gap to be filled. With regard to macro types of investigations (i.e., geological, geophysical, geotechnical, hydrogeological, etc.), some databases are currently available worldwide (i.e., Orgiazzi et al., 2017; Kmoch et al., 2021; Geyin et al., 2021; Minarelli et al., 2022), as well as in the Italian national territory. An example is provided by Vannocci et al. (2022), which includes geotechnical and hydrological soil parameters for shallow landslide modeling.

However, there are only a few examples of freely available products which integrate different macro typologies of in situ investigations in a unique database, especially with reference to the Italian territory (Gaudiosi et al., 2021). In light of the above, the aim of the authors is to make a dataset of about 1800 geological, geophysical, geotechnical, and hydrogeological in situ investigations and related parameters collected in the Terre del Reno municipality freely available (Emilia-Romagna region, Italy). The study area is affected by severe seismic hazards and prone to seismically induced effects, as extensively documented by the 2012 seismic sequence which was characterized by more than 2000 earthquakes (Facciorusso et al., 2016). Two main shocks were recorded during the crisis: the first one on 20 May with ML 5.9 and the epicenter in Finale Emilia and the second one on 29 May with ML 5.8 and the epicenter in Medolla, both in Modena Province.

As widely reported in the bibliography, the propagation of seismic waves through the upper portion of the soil can be modified by local site conditions (i.e., Bozzano et al., 2017; Fabozzi et al., 2021; Falcone et al., 2020, 2021; Gautam, 2017; Luo et al., 2020; Meza-Fajardo et al., 2019) and can determine the triggering of earthquake-induced effects at ground surface (i.e., Forte et al., 2021; Martino et al., 2017, 2019; Giannini et al., 2022; Paolella et al., 2022). In Terre del Reno, these earthquakes triggered several earthquake-induced effects (Chini et al., 2015; Papathanassiou et al., 2015), among which linear and punctual liquefaction effects were the most prominent. These effects may occur when saturated granular deposits are shaken by a seismic action, and their magnitude depends on the combination of earthquake intensity and soil condition. Literature reports plenty of liquefaction events happened in complex geological conditions and were triggered by earthquakes with various magnitudes such as for instance the Gorkha, Nepal (Gautam et al., 2017), Christchurch, New Zealand (Maurer et al., 2019), Urayasu, Japan (Baris et al., 2021), 2008 Wenchuan, China (Zhou et al., 2020), and 2019 Dürres, Albania (Mavroulis et al., 2021) earthquakes. Liquefaction effects in Terre del Reno were mainly related to the complex sedimentological and stratigraphic setting of the areas (i.e., Stefani et al., 2018; Tentori et al., 2022), characterized by multiple and alternate sandy and silty-sandy packing, hosting local (shallow) and regional (deep) aquifers (Regione Emilia-Romagna, 1998). Several authors (i.e., Ecemis, 2021; Jain et al., 2022) have highlighted that the presence of a tiny alternation of silt and sand seems to influence the liquefaction occurrence, while other studies focused on the role of silty sands and soil packing condition on liquefaction triggering (i.e., Naeini and Baziar, 2004; Stamatopoulos, 2010; Gobbi et al., 2022a). To overcome the difficulties related to heterogeneously complex soil conditions, integrated approaches are applied to predict the occurrence of liquefaction by combination of numerical and experimental methods (Gobbi et al., 2022b; Rios et al., 2022; Paolella et al., 2022). Further steps toward this direction were made for the Terre del Reno case study by pursuing two main objectives of the PERL project in order to (i) define a new integrated methodology to assess the liquefaction susceptibility in complex stratigraphic settings through a multi-level approach and (ii) perform the seismic microzonation of the municipality for land and civil protection planning purposes. This project allowed for the collection and analysis of the abovementioned in situ investigations and the elaboration of thousands of related parameters that were stored in a harmonized and standardized dataset (named PERL) conceived to guarantee interoperability with existing ICT (information and communication technologies) solutions and data models. The availability of such a dataset of surveys, catalogued and processed according to shared standards, makes Terre del Reno one of the best-characterized municipalities in Italy in terms of seismic hazard and earthquake-induced effects. This flexible dataset can be manipulated and combined to tackle different problems and represents a powerful resource for the scientific community, for those who cannot set up and manage a living laboratory or directly perform on-site investigations.

For these reasons, the authors provide complete access to the dataset through the supplementary materials and present two different applications herein used as references to highlight the potential of the PERL dataset.

2 Geological setting

2.1 Structural and stratigraphic setting

The study area is located within the southern portion of the Po alluvial plain, which represents the sedimentary cover of the Po Basin infill (Fig. 1). The geological substrate of the study area, which lies along the northern sectors of the Apennine chain, shows complex fold and thrust structures with arcuate geometry associated with strongly asymmetrical foredeep basins. Although the Po Basin represents both the Alpine retro-foreland basin and the Apennine foredeep, its Cenozoic structural evolution was mainly driven by the northeast migration of the external front of the northern Apennines, which consists of four arcuate fold-and-thrust systems: the Monferrato Arc, the Emilia Arc, the Ferrara Arc, and the Adriatic Arc (Pieri and Groppi, 1981; Royden et al., 1987; Scrocca et al., 2007). These systems that are buried beneath the present Po plain have been active since the late Miocene (Fig. 1) and are still considered seismogenic (Boccaletti et al., 2011; Ghielmi et al., 2013). In particular, the movement of a segment of the Ferrara Arc thrust system (i.e., the Mirandola thrust system) was responsible for the 2012 Emilia seismic events (ISIDe Working Group, 2010), which triggered numerous co-seismic effects associated with liquefaction phenomena in the provinces of Ferrara, Modena, and Bologna. In the study area, the shallowest Quaternary sedimentary fill consists of marine deposits (marine Quaternary in Fig. 1b) and 100 kyr spaced transgressive–regressive cycles constituted by nearshore sands and alluvial deposits, formed during interglacial and glacial periods, respectively (continental Quaternary in Fig. 1b). The stratigraphic framework of the topmost late Pleistocene to Holocene Po Basin succession (at 0–40 m depth from the ground surface) documents a succession of tabular-shaped fluvial sands (i.e., glacial) overlain up-section by the Holocene's poorly drained and mud-rich floodplain and swamp and marsh succession with subordinate lenticular fluvial sandy channel bodies associated with crevasse and levee clay–sand alternations fed by the paleo-Reno River (i.e., interglacial) (Bruno et al., 2021; Stefani et al., 2018; Tentori et al., 2022). The Reno River's modern drainage basin extends for about 2500 km2 in the northern Apennines. Owing to the low topographic gradients in the area, the paleo-Reno River experienced fast aggradation and frequent avulsion episodes during recent and historical times (see Tentori et al., 2022 and references therein).

Figure 1(a) Tectonic sketch map of the Po Plain (northern Italy) showing the main buried faults of the northern Apennines and southern Alps and epicenters of the two mainshock from the 2012 Emilia sequence (red stars) (modified from Tentori et al., 2022 and Bruno et al., 2021). The black-lined rectangle encloses the study area. (b) Simplified stratigraphic cross section and major tectonic structures along the trace A–A in panel (a) (modified after Boccaletti et al., 2004). The tectonic structures are enucleated into Mesozoic to Paleogene carbonate successions and largely controlled the sedimentary evolution of the terrigenous basins during the Neogene (Ghielmi et al., 2013; Rossi et al., 2015; Ricci Lucchi, 1986). The Pliocene–Pleistocene boundary records the transition from turbiditic sedimentation to marine clay deposition, whereas the Quaternary sedimentary fill consists of marine deposits, nearshore sands, and alluvial deposits (see text for details).

2.2 Hydrostratigraphic setting

The hydrostratigraphic architecture reflects the depositional and tectonic evolution of the southern Po sedimentary basin from the Pleistocene to Holocene (Molinari et al., 2007; Emilia-Romagna Region and ENI-AGIP, 1998). The aquifers from the most superficial hydrostratigraphic group (e.g., Group A), consist of six lower-order hydrostratigraphic units belonging to the Quaternary fluvio-deltaic and alluvial depositional systems. In the study area, Group A aquifers consist of the sandy fluvial bodies deposited during glacial periods, separated by the muddy-dominated intervals of transgressive alluvial facies (aquitards) deposited during interglacial periods. The more surficial composite aquifer system named A0 by Molinari et al. (2007) consists of two sandy-dominated aquifer units hosted within the late Pleistocene–Holocene channelized bodies and encased by alluvial floodplain muds. Based on the piezometric level dating back to the summer of 2012, Calabrese et al. (2012) placed the groundwater level of the shallower semi-confined aquifer at about 3–4 m depth below the levee and about 1–2 m in the floodplain.

3 Materials and methods

The PERL dataset was obtained by merging three databases provided by different institutions. Additional 17 geotechnical investigations were specifically performed in the framework of the PERL project.

The three existing databases are the following:

  • Modello Unico Digitale per l'Edilizia – Unique Digital Model for Building (MUDE database). The MUDE database consists of 384 records including punctual and linear in situ investigations. Data were extracted from a series of technical reports produced to plan the reconstruction works of buildings collapsed during the 2012 seismic crisis. Since the digital formats of these investigations were originally not available, geolocalization, key information, and measured parameters were obtained from the digital scans of technical and geological reports.

  • Regione Emilia-Romagna - Emilia-Romagna Region database (RER database). The RER database is composed of 906 geolocalized, punctual records, associated with a set of key information (typology, date, coordinates, and maximum depth) and a scan of the investigation sheet. Parameters were extracted from investigation sheets, as they are not available in a digital format. This database is available at (last access: 12 September 2022).

  • Seismic Microzonation Studies (SM database). The SM database is composed of 1284 records, including punctual, as well as linear in situ investigations. These investigations are geolocalized and organized in a standardized structure according to Commissione tecnica per la microzonazione sismica (2015). The key information (typology, date, coordinates, etc.) of each investigation is stored in a dedicated table, while all the measured parameters are reported in chained tables. This database is available at (last access: 12 September 2022).

The first problem faced when merging these databases was the presence of duplicate information. To avoid duplicates, a methodology to discern and verify the uniqueness of an investigation was elaborated.

This methodology is based on the implementation of a series of multiple, progressive true/false (TF) controls applied to various control parameters (CP) relative to all the investigations included in the pertinence area. The latter was defined as a circle with a radius equal to 200 m centered in the correspondence of the considered investigation. The progressively considered CP (Fig. 1) are (CP1) the absence of another investigation within the area of pertinence, (CP2) unmatching of the investigation typology, (CP3) unmatching of the date of the survey, and (CP4) matching of the maximum depth reached by the investigation. Each CPm (m= 1, 2, 3, 4) is checked in a dedicated TF test (TFn with n= 1, 2, 3, 4). Starting from TF1, an investigation that verifies CP1 is moved to TF2 for CP2 verification up to TF4. Each time a CPm in a TFn is not verified, the investigation is defined as “unique”. If an investigation verifies all the control parameters, it is defined as “redundant” and removed from the database. The application of this methodology allowed us to identify and remove 32 % of the investigations, obtaining a final dataset composed of 1805 unique investigations (Fig. 2).

Figure 2Synthetic workflow of the method used to merge the MUDE, RER, and SM databases, as well as the new realization investigations, into the PERL dataset.


4 Data description

The PERL dataset consists of two shapefiles implemented into a GIS system and an associated geodatabase. The two shapefiles are named ind_pc and ind_ln and correspond to punctual and linear investigations, respectively (EPSG:32633). The associated attribute tables contain the main information of each investigation:

  • ID. Unique identification number for each investigation.

  • Investigate. Investigation typology.

The complete set of investigations, and the related measured parameters are reported in an Excel file following this structure:

  • ID. Unique identification number of each investigation.

  • Type_par. Parameter typology (see list of parameters for legend).

  • Value. Parameter value.

  • Depth_top. Depth (m) of the layer top to which the parameter refers.

  • Depth_bottom. Depth (m) of the layer bottom to which the parameter refers.

The PERL dataset includes 1805 records corresponding to as many investigations and guarantees an average density of 35 punctual (Fig. 3a) and/or linear measurements (Fig. 3b, c) per square kilometer of the Terre del Reno municipality (about 51 km2). Focusing on investigation typologies (see the list of typologies and codes and Fig. 4a), the dataset consists of 71 % of penetrometer tests (CPT, CPTU, CPTE, SCPT, SPT, DN), 16 % of boreholes and trenches (S, T, SP, SC, SD), 12 % of punctual and linear geophysical investigations (CH, DH, HVSR, MASW, ESAC_SPAC), and 1 % of laboratory and hydrogeological tests (CR, CI, CD, SM, ED, TD, LF; see Appendix A for abbreviation definitions) (Fig. 3).

Figure 3Spatial distribution of the in situ punctual (a) and linear (b, c) investigations composing the PERL dataset. The digital elevation model (DEM) was retrieved from Regione Emilia-Romagna (2015).

Figure 4PERL dataset characteristics. (a) Classes of in situ investigations and (b) depth reached by penetrometer tests, geophysical investigations (CH and DH), boreholes, and trenches.


Penetrometer tests, geognostic boreholes, trenches, and borehole geophysical tests are characterized by a depth of investigation ranging from a few meters to more than 100 m (maximum depth: 265 m) (Fig. 4b). About 90 % of them reach a maximum depth of investigation of 35 m. Thus, the most represented depth classes are 30–35 and 10–15 m with 330 (21 %) and 310 (20 %) investigations, respectively. Penetrometer tests are characterized by depths ranging between 5 and 50 m, with the 30–35 m class being the most represented. On the contrary, boreholes and trenches cover the entire spectrum of the dataset depth. However, it is worth noticing that about 60 boreholes and trenches reach a depth higher than 55 m, which is the most represented classes together with the 10–15 m class. Penetrometer tests, boreholes, trenches, and geophysical tests are characterized by investigation depths ranging from a few to some hundred meters with a maximum of 265 m. About 90 % of them reach a maximum depth of investigation equal to 35 m. Most investigations are carried out up to a depth of 30–35 m and 10–15 m (330 (21 %) and 310 (20 %) investigations for each class, respectively). On the contrary, boreholes and trenches cover the whole spectrum of depth classes. However, it is worth noticing that about 60 boreholes and trenches reach a depth higher than 55 m, which is the most represented class together with the 10–15 m class.

5 Examples of applications

To address some of the conceptual points discussed in the Introduction and to better highlight the uniqueness and potential of the PERL dataset, we present two different applications. In the first case study, we take advantage of the PERL database to represent the complex geology beneath the San Carlo alluvial plain. The second case history focuses on a statistical inference of the PERL geophysical data to obtain soil dynamics when experimental information is missing.

5.1 Stratigraphic reconstruction of liquefiable layer thickness in the San Carlo subsoil

The PERL database includes several sedimentological, geotechnical, geophysical, and hydrogeological data which can be used to reconstruct the stratigraphic architecture of the Terre del Reno subsurface and provide a reliable geological framework for future studies devoted to earthquake-induced hazard mitigation. The position and the thickness of the liquefiable portion within the subsoil are key information for liquefaction risk assessment and mitigation. The possibility of an automatically built three-dimensional subsoil model with advanced procedures represents a current topic of the applied technological research. Here, a combination of these two approaches is presented to spotlight the potential of the PERL dataset.

As an example, Fig. 5 shows the geostatistical interpolation, performed with the ordinary kriging, of the cumulated thickness of the liquefiable layers (CTL) in the district of San Carlo. In particular, the CTL has been manually extracted from 33 boreholes and automatically obtained on 148 CPTs by applying the procedure proposed by Spacagna et al. (2022). The obtained map has been overlayed on the map of liquefaction evidence that occurred after the Emilia-Romagna 2012 seismic sequence, showing a good match between the liquefaction-induced surficial manifestations and the CTL distribution.

Figure 5Geostatistical maps of the cumulated thickness of the liquefiable layer.

5.2 Statistical analysis of shear waves variability with depth

As widely represented by literature data, the amount of available investigations progressively decreases with soil depth. Thus the uncertainty in subsoil characterization increases from the ground surface to the deepest layers of the soil. To overcome this issue, an example of statistical inference of soil parameters is presented. Based on the available data, a correlation analysis between the values of Vs (m s−1) and depth (m) was carried out to infer information when the depth of investigation does not guarantee a correct soil parameterization. Considering the high geological and stratigraphic complexity of the case study, a lithological-based statistical inference was performed. Specifically, for each lithology L, the scatter plot of the value of Vs (m s−1) as a function of depth, D (m), and the corresponding linear regression were calculated.

The linear regression model is defined by the following relationship:


where a and Vs0 correspond to the slope and intercept of the model line, respectively.

The PERL database can rely on 164 Vs profiles mainly identified from penetrometer tests (SCPT), down hole (DH), MASW, ESAC_SPAC, SDMT, and cross hole (CH) tests. Each of these Vs profiles was discretized with a step size of 1 m in depth and, through an automated procedure, each meter of depth was associated with lithological (L) information extracted from proximal boreholes.

For each lithology L, a statistical analysis was performed, and the linear regression models were calculated. The results obtained for MH and SP lithologies are here presented as case examples (Fig. 6). The corresponding values of coefficients a and b and the coefficient of determination R2 are reported in Table 1.

Figure 6Scatter plot and regression line model for lithology MH (a) and SP (b). Equations and R2 are also reported.


Table 1Coefficients (a= slope, Vs0= intercept) of the regression model for two lithologies (L) and related coefficient of determination (R2).

Download Print Version | Download XLSX

As expected, the obtained a and b values highlight a positive slope as the depth mean value increases together with the mean of the dependent variable Vs. At the same time, these values allow us to quantify how much the Vs value changes per meter of depth, showing an SP variation rate greater than that characterizing MH. Moreover, the coefficient of determination R2 obtained for MH and SP lithologies are characterized by a 0.59 and 0.79 value, respectively, proving the reliability of the fitting. When experimental data are lacking, for depths that fall within the variability range of the available data, the regression models allow us to obtain Vs by interpolation, while an extrapolation can be applied for greater depths.

Results may be used in the future for comparison with other Italian estimates (Romagnoli et al., 2022) or combined with ambient vibration measurements to define the thickness of the resonant sedimentary layers (D'Amico et al., 2008; Giannini et al., 2021).

6 Conclusions

As part of PERL project, a considerable number of investigations were collected in the Terre del Reno municipality, Emilia-Romagna region (Italy). This area experienced the most extensive liquefaction effects during the 2012 Emilia-Romagna seismic crisis and remains exposed to severe seismic hazards and seismically induced effects due to its complex geological setting.

Thanks to this study, complete and free access to the PERL dataset, which includes 1805 punctual and linear in situ investigations consisting of geological, geotechnical, geophysical, and hydrogeological data, is provided. The database is composed of 71 % of penetrometer tests, 16 % of boreholes and trenches, 12 % of geophysical investigations, and 1 % of laboratory and hydrogeological tests.

Two applications of the PERL dataset are presented to highlight its potentiality and to show that a high-quality large dataset could be critical to infer information in areas and/or portions of the soil characterized by poor or sparse data. The first examples pointed out the database potentials in overcoming problems due to the uneven distribution of surveys across the territory, while the second one spotlights its capability to provide information at subsoil depth not reached by investigations.

Other major outcomes of the PERL project in the Terre del Reno municipality included (i) a detailed reconstruction of the subsoil geology to examine the stratigraphic control on earthquake-induced liquefaction (Tentori et al., 2022); (ii) data-driven and automatic subsoil characterization through the analysis of CPT-based soil behavior type (SBT) and soil behavior type indexes (Ic), combining geostatistical and artificial intelligence genetic approaches (Baris et al., 2022); and (iii) a third-level seismic microzonation study of the Terre del Reno municipal area (Varone et al., 2022) to mitigate seismic and seismically induced hazards through sensible urban planning.

Ongoing studies currently focus on the definition of a comprehensive methodology to quantify liquefaction susceptibility in areas dominated by complex geostratigraphic conditions by applying a multi-level approach laying on simplified models and to promote the identification of potentially liquefiable granular bodies, thus mitigating earthquake-induced hazards to allow for a sustainable development of this urban area.

Appendix A: List of investigation typologies and codes
CH Cross-hole test
CPT Cone penetration test
CPTE Electrical cone penetration test
CPTU Piezocone penetration test
CR Resonant column test
DH Down-hole test
DMT Dilatometric test
DN Dynamic cone penetration test
ESAC_SPAC Seismic array elaborated by ESAC/SPAC methods
HVSR Ambient noise measurements elaborated by HVSR technique
LF Lefranc test
MASW Multichannel analysis of surface waves
PA Borehole (water well)
S Non-destructive borehole
SC Borehole with collection of samples
SCPT Seismic cone penetration test
SD Borehole
SDMT Seismic dilatometer Marchetti test
SM Laboratory test on soil sample
SP Piezometer
SPT Standard penetration test
T Trench
TD Shear strength test
Appendix B: List of parameters and codes
AR Clay classification obtained from laboratory test (%)
C Effective cohesion (MPa)
CAM Number of sample (–)
CU Undrained cohesion (MPa)
E1 Index of voids (–)
F1 Effective soil friction angle ()
FR Resonance frequency (Hz)
FS Skin friction (MPa)
G Shear modulus (MPa)
GH Gravel classification obtained from laboratory test (%)
IP Plasticity index (–)
K Classification obtained from laboratory test (m s−1)
L Layer lithology (–)
LID Lithology of the hydrolayer (–)
LM Silt classification obtained from laboratory test (%)
PT Number of SPT blows (–)
PTM Number of DN blows (–)
PV Weight of the unit of volume (kN m−3)
QC Tip resistance (MPa)
SA Sand classification obtained from laboratory test (%)
SG Water table level (m)
U Hydrostatic pressure (MPa)
VP Compressional waves velocity (m s−1)
VS Shear-waves velocity (m s−1)
W Water content (%)
Appendix C: List of lithologies (L) and codes
CH Inorganic clays with high plasticity
CL Inorganic clays with low-medium plasticity, gravelly or sandy, silty clays
GP Clean gravel poorly graded, mixture of gravel and sand
MH Inorganic silts, fine sands, micaceous or diatomitic silts
ML Inorganic silts, silty or clayey fine sands, silts clayey sands with low plasticity
OH Organic clays of medium-high plasticity, organic silts
OL Organic silt, silty clays with low plasticity
PT Peats and peaty soils
RI Anthropogenic filling
SC Clayey sands, mixture of sand and clay
SM Silty sands, mixture of sand and silt
SP Sands poorly graded
SW Sands well graded, gravelly sands
Data availability

The authors confirm that the data supporting the findings of this study are available in the Supplement.


The supplement related to this article is available online at:

Author contributions

CV, GC, and FS conceptualized the methodology. CV, GC, AB, MCC, SF, CF, IG, SG, LP, MS, PS, RS, and DT curated the data. GC implemented the algorithm. CV, AB, DT, and RS curated the examples of applications. LM, GM, and MaMo conceptualized and administrated the project. CV wrote the original draft of the manuscript. CV, GC, AB, MCC, SF, CF, IG, SG, MaMa, LP, MS, PS, RS, FS, DT, LM, GM, and MaMo reviewed and edited the manuscript.

Competing interests

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.

Special issue statement

This article is part of the special issue “Earthquake-induced hazards: ground motion amplification and ground failures”. It is a result of the EGU General Assembly 2021, 19–30 April 2021.


The authors wish to thank two anonymous referees for their valuable suggestions that improved the quality of the manuscript. This work was carried out in the framework of the project Accordo di collaborazione tra Regione Emilia-Romagna, Amministrazione Comunale di Terre del Reno, Consiglio Nazionale delle Ricerche – Istituto di Geologia Ambientale e Geoingegneria, Università degli Studi di Cassino e del Lazio Meridionale – Dipartimento di Ingegneria Civile e Meccanica, finalizzato a definire una strategia multilivello per valutare il rischio da liquefazione in presenza di argini e situazioni geologiche e morfologiche complesse. The authors wish to thank Olga Mantovani (municipality of Terre del Reno) for its support to the project.

Review statement

This paper was edited by Giovanni Forte and reviewed by two anonymous referees.


Baris, A., Spacagna, R.L., Paolella, L., Koseki, J., and Modoni, G.: Liquefaction fragility of sewer pipes derived from the case study of Urayasu (Japan), B. Earthq. Eng., 19, 3963–3986,, 2021.  

Baris, A., Caciolli, M. C., Fabozzi, S., Gaudiosi, I., Mancini, M., Martelli, L., Modoni, G., Moscatelli, M., Paolella, L., Razzano, R., Simionato, M., Spacagna, R. L., Stigliano, F., Tentori, D., and Varone, C.: Automatized CPT-based soil profile characterization for liquefaction susceptibility assessment, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-9779,, 2022. 

Boccaletti, M., Bonini, M., Corti, G., Gasperini, P., Martelli, L., Piccardi, L., Severi, P., and Vannucci, G.: Carta sismotettonica della Regione Emilia-Romagna, scala 1:250.000, Con note illustrative, Regione Emilia-Romagna–SGSS, CNR-IGG, SELCA, Firenze, 2004. 

Boccaletti, M., Corti, G., and Martelli, L.: Recent and active tectonics of the external zone of the Northern Apennines (Italy), Int. J. Earth Sci., 100, 1331–1348,, 2011. 

Bozzano, F., Buccellato, A., Coletti, F., Martino, S., Marra, F., Rivellino, S., and Varone, C.: Analysis of the seismic site effects along the Ancient Via Laurentina (Rome), Ann. Geophys., 60, S0435,, 2017. 

Bruno, L., Amorosi, A., Lugli, S., Sammartino, I., and Fontana, D.: Trunk river and tributary interactions recorded in the Pleistocene–Holocene stratigraphy of the Po Plain (northern Italy), Sedimentology, 68, 2918–2943,, 2021. 

Calabrese, L., Martelli, L., and Severi, P.: Stratigrafia dell'area interessata dai fenomeni di liquefazione durante il terremoto dell'Emilia (maggio 2012), Atti 31 Convegno del Gruppo Nazionale di Geofisica della Terra Solida, Potenza, 20–22 novembre 2012, sess. 2.2, 119–126, 2012. 

Commissione tecnica per la microzonazione sismica, Microzonazione sismica, Standard di rappresentazione e archiviazione informatica, Versione 4.0b, 2015. 

Chini, M., Albano, M., Saroli, M., Pulvirenti, L., Moro, M., Bignami, C., Falcucci, E., Gori, S., Modoni, G., Pierdicca, N., and Stramondo, S.: Coseismic liquefaction phenomenon analysis by COSMO-SkyMed: 2012 Emilia (Italy) earthquake, Int. J. Appl. Earth Obs., 39, 65–78,, 2015. 

D'Amico, V., Picozzi, M., Baliva, F., and Albarello, D.: Ambient Noise Measurements for Preliminary Site-Effects Characterization in the Urban Area of Florence, Italy, B. Seismol. Soc. Am., 98, 1373–1388,, 2008. 

Ecemis, N.: Experimental and numerical modeling on the liquefaction potential and ground settlement of silt-interlayered stratified sands, Soil Dyn. Earthq. Eng., 144, 106691,, 2021. 

Emilia-Romagna Region and ENI-AGIP: Book: Riserve idriche sotterranee della Regione Emilia-Romagna, edited by: Regione Emilia-Romagna, Servizio Geologico, Sismico dei Suoli, ENI – AGIP, 1998. 

Fabozzi, S., Catalano, S., Falcone, G., Pagliaroli, A., Peronace, E., Porchia, A., Romagnoli, G., and Moscatelli M.: Stochastic approach to study the site response in presence of shear wave velocity inversion: application to seismic microzonation studies in Italy, Eng. Geol., 280, 105914,, 2021. 

Facciorusso, J., Madiai, C., and Vannucchi, G.: The 2012 Emilia earthquake (Italy): Geotechnical characterization and ground response analyses of the paleo-Reno river levees, Soil Dyn. Earthq. Eng., 86, 71–88,, 2016. 

Falcone, G., Mendicelli, A., Mori, F., Fabozzi, S., Moscatelli, M., Occhipinti, G., and Peronace, E.: A simplified analysis of the total seismic hazard in Italy, Eng. Geol., 267, 105511,, 2020. 

Falcone, G., Acunzo, G., Mendicelli, A., Mori, F., Naso, G., Peronace, E., Porchia, A., Romagnoli, G., Tarquini, E., and Moscatelli, M.: Seismic amplification maps of Italy based on site-specific microzonation dataset and one-dimensional numerical approach, Eng. Geol., 289, 106170,, 2021. 

Forte, G., Verrucci, L., Di Giulio, A., De Falco, M., Tommasi, P., Lanzo, G., Franke, K. W., and Santo, A.: Analysis of major rock slides that occurred during the 2016–2017 Central Italy seismic sequence, Eng. Geol., 290, 106194,, 2021. 

Gaudiosi, I., Romagnoli, G., Albarello, D., Fortunato, C., Imprescia, P., Moscatelli, M., and Stigliano, F.: Verso un approccio totalmente probabilistico alle stime di pericolosità sismica: studio della variabilità delle curve del modulo secante normalizzato, G/G0 (γ) e del rapporto di smorzamento, D (γ) con la deformazione di taglio, associate alle unità geologico-tecniche degli Studi di Microzonazione Sismica, 39 Convegno Nazionale GNGTS – Gruppo Nazionale di Geofisica della Terra Solida (GNGTS), 22–24 June 2021,, 2021 (in Italian). 

Geyin, M., Maurer, B. W., Bradley, B. A., Green, R. A., and van Ballegooy, S.: CPT-based liquefaction case histories compiled from three earthquakes in Canterbury, New Zealand, Earthq. Spectra, 37, 2920–2945,, 2021. 

Ghielmi, M., Minervini, M., Nini, C., Rogledi, S., and Rossi, M.: Late Miocene–Middle Pleistocene sequences in the Po Plain–Northern Adriatic Sea (Italy): the stratigraphic record of modification phases affecting a complex foreland basin, Mar. Petrol. Geol., 42, 50–81, 2013. 

Giannini, L. M., Varone, C., Esposito, C., Scarascia Mugnozza, G., and Schilirò, L.: The potential of spatial statistics for the reconstruction of a subsoil model: A case study for the Firenze-Prato-Pistoia Basin, Central Italy, J. Appl. Geophys., 194, 104466,, 2021. 

Giannini, L. M., Varone, C., Esposito, C., Marmoni, G. M., Scarascia Mugnozza, G., and Schilirò, L.: Earthquake-induced reactivation of landslides under variable hydrostatic conditions: evaluation at regional scale and implications for risk assessment, Landslides, 19, 2005–2019,, 2022. 

Gobbi, S., Reiffsteck, P., Lenti, L., Santisi d'Avila, M. P., and Semblat, J. F.: Liquefaction triggering in silty sands: effects of non-plastic fines and mixture-packing conditions, Acta Geotech., 17, 391–410,, 2022a. 

Gobbi, S., Santisi d'Avila, M. P., Lenti, L., Semblat, J. F., and Reiffsteck, P.: Liquefaction assessment of silty sands: Experimental characterization and numerical calibration, Soil Dyn. Earthq. Eng., 159, 107349,, 2022b. 

Guarino, P. M., Santo, A., Forte, G., De Falco, M., and Niceforo, D. M. A.: Analysis of a database for anthropogenic sinkhole triggering and zonation in the Naples hinterland (Southern Italy), Nat. Hazards, 91, 173–192,, 2018. 

Gautam, D.: Unearthed lessons of 25 April 2015 Gorkha earthquake (MW7.8): geotechnical earthquake engineering perspectives, Geomat. Nat. Haz. Risk, 8, 1358–1382,, 2017. 

Gautam, D., Santucci de Magistris, F., and Fabbrocino, G.: Soil liquefaction in Kathmandu valley due to 25 April 2015 Gorkha, Nepal earthquake, Soil Dyn. Earthq. Eng., 97, 37–47,, 2017. 

Iside Working Group: Italian seismological instrumental and parametric database, (last access: 23 January 2023), 2010. 

Jain, A., Mittal, S., and Kumar Shukla, S.: Liquefaction proneness of stratified sand-silt layers based on cyclic triaxial tests, Journal of Rock Mechanics and Geotechnical Engineering, in press,, 2022. 

Kmoch, A., Kanal, A., Astover, A., Kull, A., Virro, H., Helm, A., Pärtel, M., Ostonen, I., and Uuemaa, E.: EstSoil-EH: a high-resolution eco-hydrological modelling parameters dataset for Estonia, Earth Syst. Sci. Data, 13, 83–97,, 2021. 

Luo, Y., Fan, X., Huang, R., Wang, Y., Yunus, A. P., and Havenith, H. B.: Topographic and near-surface stratigraphic amplification of the seismic response of a mountain slope revealed by field monitoring and numerical simulations, Eng. Geol., 271, 105607,, 2020 

Martino, S., Prestininzi, A., and Romeo, R. W.: Earthquake-induced ground failures in Italy from a reviewed database, Nat. Hazards Earth Syst. Sci., 14, 799–814,, 2014. 

Martino, S., Bozzano, F., Caporossi, P., D'Angiò, D., Della Seta, M., Esposito, C., Fantini, A., Fiorucci, M., Giannini, L. M., Iannucci, R., Marmoni, G. M., Mazzanti, P., Missori, C., Moretto, S., Rivellino, S., Romeo, R. W., Sarandrea, P., Schilirò, L., Troiani, F., and Varone, C.: Ground effects triggered by the August 24th 2016, Mw 6.0 Amatrice (Italy) Earthquake: Survey and Inventorying to update the CEDIT catalogue, Geogr. Fis. Din. Quat., 40, 1–20,, 2017. 

Martino, S., Bozzano, F., Caporossi, P., D'Angiò, D., Della Seta, M., Esposito, C., Fantini, A., Fiorucci, M., Giannini, L. M., Iannucci, R., Marmoni, G. M., Mazzanti, P., Missori, C., Moretto, S., Piacentini, D., Rivellino, S., Romeo, R. W., Sarandrea, P., Schilirò, L., Troiani, F., and Varone, C.: Impact of landslides on transportation routes during the 2016-2017 Central Italy seismic sequence, Landslides, 16, 1221–1241,, 2019. 

Maurer, B. W., Green, R. A., van Ballegooy, S., and Wotherspoon, L.: Development of region-specific soil behavior type index correlations for evaluating liquefaction hazard in Christchurch, New Zealand, Soil Dyn. Earthq. Eng., 117, 96–105,, 2019. 

Mavroulis, S., Lekkas, E., and Carydis, P.: Liquefaction Phenomena Induced by the 26 November 2019, Mw = 6.4 Durrës (Albania) Earthquake and Liquefaction Susceptibility Assessment in the Affected Area, Geosciences, 11, 215,, 2021. 

Meza-Fajardo, K. C., Varone, C., Lenti, L., Martino, S., and Semblat, J. F.: Surface wave quantification in a highly heterogeneous alluvial basin: Case study of the Fosso di Vallerano valley, Rome, Italy, Soil Dyn. Earthq. Eng., 120, 292–300,, 2019. 

Minarelli, L., Amoroso, S., Civico, R., De Martini, P. M., Lugli, S., Martelli, L., Molisso, F., Rollins, K. M., Salocchi, A., Stefani, M., Cultrera, G., Milana, G., and Fontana, D.: Liquefied sites of the 2012 Emilia earthquake: a comprehensive database of the geological and geotechnical features (Quaternary alluvial Po plain, Italy), B. Earthq. Eng., 20, 3659–3697,, 2022. 

Molinari, F. C., Boldrini, G., Severi, P., Duroni, G., Rapti-Caputo, D., and Martinelli, G.: Risorse idriche sotterranee della Provincia di Ferrara, Regione Emilia Romagna, edited by: MAP, D. B., Firenze, 61 pp., 2007. 

Naeini, S. A. and Baziar, M. H.: Effect of fines content on steady-state strength of mixed and layered samples of a sand, Soil Dyn. Earthq. Eng., 24, 181–187,, 2004 

Orgiazzi, A., Ballabio, C., Panagos, P., Jones, A., and Fernández-Ugalde, O.: LUCAS Soil, the largest expandable soil dataset for Europe: a review, Eur. J. Soil Sci., 69, 140–153,, 2017. 

Paolella, L., Modoni, G., Spacagna, R. L., and Baris, A.: A generalized severity number to predict damage with lateral spreading, Géotechnique, 22 pp.,, 2022. 

Papathanassiou, G., Mantovani, A., Tarabusi, G., Rapti, D., and Caputo, R.: Assessment of liquefaction potential for two liquefaction prone areas considering the May 20, 2012 Emilia (Italy) earthquake, Eng. Geol., 189, 1–16,, 2015. 

Pieri, M. and Groppi, G.: Subsurface Geological Structure of the Po Plain (Italy), CNR – Progetto Finalizzato Geodinamica, 414, 1–13, 1981. 

Regione Emilia-Romagna: DTM 5x5, Regione Emilia-Romagna [dataset], (last access: 22 September 2022), 2015. 

Ricci Lucchi, F.: Oligocene to recent foreland basins of northern Apennines, in: Foreland Basins IAS Special Publication, edited by: Allen, P. and Homewood, P., Blackwell Scientific, London, 8, 105–139, 1986. 

Rios, S., Millen, M., Quintero, J., and da Fonseca, A. V.: Analysis of simplified time of liquefaction triggering methods by laboratory tests, physical modelling and numerical analysis, Soil Dyn. Earthq. Eng., 157, 107261,, 2022. 

Romagnoli, G., Tarquini, E., Porchia, A., Catalano, S., Albarello, D., and Moscatelli, M.: Constraints for the Vs profiles from engineering-geological qualitative characterization of shallow subsoil in seismic microzonation studies, Soil Dyn. Earthq. Eng., 161, 107347,, 2022. 

Rossi, M., Minervini, M., Ghielmi, M., and Rogledi, S.: Messinian and Pliocene erosional surfaces in the Po Plain-Adriatic Basin: Insights from allostratigraphy and sequence stratigraphy in assessing play concepts related to accommodation and gateway turnarounds in tectonically active margins, Mar. Petrol. Geol., 66, 192–216, 2015. 

Royden, L., Patacca, E., and Scandone, P.: Segmentation and configuration of subducted lithosphere in Italy: an important control on thrust-belt and foredeep-basin evolution, Geology, 15, 714–717, 1987. 

Scrocca, D., Carminati, E., Doglioni, C., and Marcantoni, D.: Slab retreat and active shortening along the central–northern Apennines. Thrust Belts and Foreland Basins, Springer, 471–487, ISBN: 1863-4621, 2007. 

Spacagna, R. L., Baris, A., Paolella, L., and Modoni, G.: Data-driven soil profile characterization using statistical methods and artificial intelligence algorithms, in: Cone Penetration Testing 2022, CRC Press, 708–714,, 2022. 

Stamatopoulos, C. A.: An experimental study of the liquefaction strength of silty sands in terms of the state parameter, Soil Dyn. Earthq. Eng., 30, 662–678,, 2010. 

Stefani, M., Minarelli, L., Fontana, A., and Hajdas, I.: Regional deformation of late Quaternary fluvial sediments in the Apennines foreland basin (Emilia, Italy), Int. J. Earth Sci., 107, 2433–2447,, 2018. 

Tanyaş, H., Görüm, T., Fadel, I., Yıldırım, C., and Lombardo, L.: An open dataset for landslides triggered by the 2016 Mw 7.8 Kaikōura earthquake, New Zealand, Landslides, 19, 1405–1420,, 2022. 

Tentori, D., Mancini, M., Varone, C., Spacagna, R., Baris, A., Milli, S., Gaudiosi, I., Simionato, M., Stigliano, S., Modoni, G., Martelli, L., and Moscatelli, M.: The influence of alluvial stratigraphic architecture on liquefaction phenomena: a case study from the Terre del Reno subsoil (southern Po plain, Italy), Sediment. Geol., 440, 106258,, 2022. 

Vannocci, P., Segoni, S., Masi, E. B., Cardi, F., Nocentini, N., Rosi, A., Bicocchi, G., D'Ambrosio, M., Nocentini, M., Lombardi, L., Tofani, V., Casagli, N., and Catani, F.: Towards a National-Scale Dataset of Geotechnical and Hydrological Soil Parameters for Shallow Landslide Modeling, Data, 7, 37,, 2022. 

Varone, C., Baris, A., Caciolli, M. C., Fabozzi, S., Fortunato, C., Gaudiosi, I., Giallini, S., Mancini, M., Martelli, L., Modoni, G., Moscatelli, M., Paolella, L., Simionato, M., Sirianni, P., Spacagna, R. L., Stigliano, F., Tentori, D., and Razzano, R.: An integrated approach for engineering - geological modelling in view of seismic microzonation, EGU General Assembly 2022, Vienna, Austria, 23–27 May 2022, EGU22-8586,, 2022. 

Zhou, Y., Xia, P., Ling, D., and Chen, Y.: Liquefaction case studies of gravelly soils during the 2008 Wenchuan earthquake, Eng. Geol., 274, 105691,, 2020. 

Short summary
In 2012, Italy was struck by a seismic crisis characterized by two main shocks and relevant liquefaction events. Terre del Reno is one of the municipalities that experienced the most extensive liquefaction effects; thus it was chosen as case study for a project devoted to defining a new methodology to assess the liquefaction susceptibility. In this framework, about 1800 geotechnical, geophysical, and hydrogeological investigations were collected and stored in the publicly available PERL dataset.
Final-revised paper