Search

Europe has experienced over the last years earthquakes of moderate magnitude (Mw 5-6), yet destructive, reminding us of the seismogenic potential of slowly deforming regions. Among them, the 2020 Mw 6.4 Petrinja earthquake ruptured the Petrinja-Pokupsko Fault (PPKF) in Central Croatia, about 50-km southeast of Zagreb, a region in which the caracterisation of seismogenic faults had been insufficiently studied before that event. Understanding the strain accommodation through time and space is critical for accurate assessment of the regional seismic hazard.Using field observations and high-resolution topographical data derived from airborne LiDAR (~10 cm resolution) and tri-stereo satellite images (Pléiades, resolution 50 cm), we accurately mapped the fault trace, underlined at several sites by geomorphic markers such as valleys, terrace risers, and alluvial fans that have recorded cumulative displacements ranging from 5 to > 50 m and potentially up to ~180 m. Along the studied section, our fault mapping is composed of a clear NW-SE-trending 10-km-long strand between Donja and Cepelis, and of 1-4-km-long right-stepping segments marked by a non-negligible vertical component. The southern strand is composed of 2-3 sub-parallel segments that accommodate the deformation within a < 500 m wide fault zone.We have identified several sites on the main southern strand where offsets have been accurately measured and where displaced markers have been sampled for cosmogenic nuclide exposure dating and radiocarbon datings. This will allow to estimate the slip-rate for this fault at different sites and over several time spans.The mapped fault appears very discontinuous with the deformation absorbed by a series of small fault sections rather than on a single fault strand. This likely reflects a recent transpressive deformation, with immature faults, in agreement with the source parameter of the 2020 Petrinja earthquake derived from seismology.Finaly, the 2020 coseismic surface ruptures affected the northern section of the PPKF, while the mapped cumulative displacements appears more prominent along the southern section. A better knowledge of the seismic history of this entire fault system is thus crucial for seismic hazard assessment of this area.

The Val d’Agri (VA) oilfield in the Lucanian Apennines (southern Italy), represents the largest onshore in Europe. Since the 1990's, hydrocarbons are produced from a fractured carbonate reservoir with an average extraction rate of 7*104 barrels/day of oil and 3*106 Smc/day of gas. Part of the wastewater has been re-injected since 2006 into a marginal portion of the reservoir by a high-rate well (Costa Molina 2, CM2). Charged by the Italian oil and gas safety authority, the National Institute of Geophysics and Volcanology (INGV) monitors the VA industrial hydrocarbon operations through the research activity of a dedicated working group (CMS, Centro di Monitoraggio del Sottosuolo) and according to the governmental monitoring guidelines. The CMS operates the real-time acquisition and offline analyses of seismic data recorded at 56 seismic stations associated with public and private local seismic networks. The principal aim of the CMS is to investigate the risk associated with industrial activities that can induce or trigger seismic events by producing stress changes within the upper crustal volume. Previous works have highlighted a spatio-temporal relationship between micro-seismicity (ML ≤ 2.2) and wastewater injection, delineating a NE-dipping back-thrust near the CM2. Part of the microseismicity recorded in the southwestern portion of the VA has also been associated with the water level changes of the Pertusillo lake. One of the main challenges is to define an accurate structural setting of the VA to understand the potential of earthquakes in the area and investigate the presence of active faults. The VA consists of a Quaternary extensional tectonic basin and it is one of the areas of highest seismic hazard in Italy (Basilicata, 1857, M7 earthquake). The basin is bounded by two parallel and oppositely dipping normal fault systems: the Monti della Maddalena Fault System (MMFS) on its western side and the Eastern Agri Fault System (EAFS) on the eastern one. The characterization of the ongoing tectonic activity of the MMFS and EAFS, and their hierarchical relationship is still generating debate among the scientific community. We adopt a multidisciplinary approach based on detailed geological-structural, geophysical and seismic analyses, and electrical resistivity tomography, aimed at reconstructing the subsurface geology of the area and recognizing and characterizing the active and capable faults, and the associated potential for local seismic hazard. We present and discuss the results of this work, focusing on the relative location of seismic events that occurred between March and June 2022. The outcomes allow inferring interesting geologic constraints, highlighting the relationships between the distribution of local seismicity and the structural setting of the area in the uppermost crust (depth < 6 km).

This paper provides a new contribution to the construction of the complex and fragmentary mosaic of the Late Holocene earthquakes history of the İznik segment of the central strand of the North Anatolian Fault (CNAF) in Turkey. The CNAF clearly displays lower dextral slip rates with respect to the northern strand however, surface rupturing and large damaging earthquakes (M > 7) occurred in the past, leaving clear signatures in the built and natural environments. The association of these historical events to specific earthquake sources (e.g., Gemlik, İznik, or Geyve fault segments) is still a matter of debate. We excavated two trenches across the İznik fault trace near Mustafali, a village about 10 km WSW of İznik where the morphological fault scarp was visible although modified by agricultural activities. Radiocarbon and TL dating on samples collected from the trenches show that the displaced deposits are very recent and span the past 2 millennia at most. Evidence for four surface faulting events was found in the Mustafali trenches. The integration of these results with historical data and previous paleoseismological data yields an updated Late Holocene history of surface-rupturing earthquakes along the İznik Fault in 1855, 740 (715), 362, and 121 CE. Evidence for the large M7 + historical earthquake dated 1419 CE generally attributed to this fault, was not found at any trench site along the İznik fault nor in the subaqueous record. This unfit between paleoseismological, stratigraphic, and historical data highlights one more time the urge for extensive paleoseismological trenching and offshore campaigns because of the high potential to solve the uncertainties on the seismogenic history (age, earthquake location, extent of the rupture and size) of this portion of NAFZ and especially on the attribution of historical earthquakes to the causative fault.

This is the peer reviewed version of the following article: Smedile, A., Molisso, F., Chague, C., Iorio, M., De Martini, P.M., Pinzi, S., Collins, P.E.F., Sagnotti, L. & Pantosti, D. 2019, "New coring study in Augusta Bay expands understanding of offshore tsunami deposits (Eastern Sicily, Italy)", Sedimentology, , which has been published in final form at doi:10.1111/sed.12581. This article may be used for non-commercial purposes in accordance with Wiley Terms and Conditions for Use of Self-Archived Versions. © 2019 The Authors. Sedimentology © 2019 International Association of Sedimentologists.

We present a 1:10,000 scale map of the coseismic surface ruptures following the 26 December 2018 Mw 4.9 earthquake that struck the eastern flank of Mt. Etna volcano (southern Italy). Detailed rupture mapping is based on extensive field surveys in the epicentral region. Despite the small size of the event, we were able to document surface faulting for about 8 km along the trace of the NNW-trending active Fiandaca Fault, belonging to the Timpe tectonic system in the eastern flank of the volcano. The mapped ruptures are characterized in most cases by perceivable opening and by a dominant right-oblique sense of slip, with an average slip of about 0.09 m and a peak value of 0.35 m. It is also noteworthy that the ruptures vary significantly in their kinematic expression, denoting locally high degree of complexity of the surface faulting.

International audience; After the 2004 Indian Ocean (IOT) and the 2011 Tohoku-oki tsunamis, new research in tsunami-related fields was strongly stimulated worldwide and also in the Mediterranean. This research growth yields substantial advancements in tsunami knowledge.Among these advancements is the “Paleotsunami” research that has marked particular progress on the reconstruction of the tsunami history of a region. As an integration of the historical documentation available in the Mediterranean and the Gulf of Cadiz areas, geological and geoarchaeological records provide the insights to define the occurrence, characteristics, and impact of tsunamis of the past. Here, we present the recent advancements done for both the onshore and offshore realms.As for the onshore, we discuss case studies dealing with recent high-resolution works based on: a) direct push in situ sensing techniques, applied to identification and characterization of typical paleotsunami deposits features; b) combined XRF-X-CT approach, implemented for the identification of fine-scale sedimentary structures useful for the definition of the causative flow dynamics; c) the geoarchaeological “new field” contribution, with the development of specific diagnostic criteria in search for tsunami impact traces in archaeological strata; d) comparison of multiple dating methods and of different modeling codes for the definition of the potential source for the displacement of boulders of exceptional dimension, identified by 3D size calculation.As for the offshore advancements, we present case studies focusing on the recognition of tsunami deposits and their sedimentary traces in the geological record from the nearshore, thanks to diver-operated coring equipment, down to the continental slope, by means of vibracorer and long gravity core sampling in deeper areas. The examples provided show a multiproxy approach with a high potential of retrieving a complete record of paleotsunami traces at least during the Holocene. This is based on the combination of multidisciplinary approaches including X-ray imaging, high-resolution measurement of physical properties, X-ray fluorescence data, grain-size analysis, micropaleontology, palynological content, isotopic and optically stimulated luminescence dating methods.

The 29 December 2020, Mw 6.4 Petrinja earthquake nucleated at a depth of ~10 km in the Sisak-Moslavina County in northern Croatia, ~6 km WSW of the Petrinja town. Focal mechanisms, aftershocks distribution, and preliminary Sentinel-1 InSAR interferogram suggest that the NW-SE right-lateral strike-slip Pokupsko-Petrinja fault was the source of this event.The Croatian Geological Survey, joined by a European team of earthquake geologists from France, Slovenia and Italy, performed a prompt systematic survey of the area to map the surface effects of the earthquake. The field survey was guided by geological maps, preliminary morphotectonic mapping based on 1:5,000 topographical maps and InSAR interferogram. Locally, field mapping was aided by drone survey.We mapped unambiguous evidence of surface faulting at several sites between Župić to the NW and Hrastovica to the SE, in the central part of the Pokupsko-Petrinja fault, for a total length of ~6.5 km. This is probably a minimum length since several portions of the fault have not been explored yet, and in part crossing forbidden uncleared minefields. Surface faulting was observed on anthropic features (roads, walls) and on Quaternary sediments (soft colluvium and alluvium) and Miocene bedrock (calcarenites). The observed ruptures strike mostly NW-SE, with evidences of strike-slip right-lateral displacement and zones of extension (opening) or contraction (small pressure ridges, moletracks) atlocal bends of the rupture trace. Those ruptures are interpreted as evidences of coseismic surface faulting (primary effects) as they affect the morphology independently from the slope direction. Ground failures due to gravitational sliding and liquefaction occurrences were also observed, mapped and interpreted as secondary effects (see Amoroso et al., and Vukovski et al., this session). SE of Križ, the rupture broke a water pipeline with a right-lateral offset of several centimetres. Measured right-lateral net displacement varies from a few centimetres up to ~35 cm. A portion of the maximum measured displacement could be due to afterlisp, as it was mapped several days after the main shock. Hybrid surface ruptures (shear plus opening and liquefaction), striking SW-NE, with cm-size left-lateral strike-slip offsets were mapped on the northern side of the Petrinja town, ~3 km NE of the main fault.Overall, the rupture zone appears discontinuous. Several factors might be inferred to explain this pattern such as incomplete mapping of the rupture, inherited structural discontinuities within the Pokupsko-Petrinja fault system, or specific mechanical properties of the Neogene-Quaternary strata

Earthquakes and related coseismic effects at the surface, both primary and secondary, such as liquefaction and lateral spreading, can impact humans due to induced economic or social disruptions (e.g. slope, bridge and building foundation failures, flotation of buried structures). In this respect, it results of primary interest to map liquefaction induced evidences soon after an earthquake. On the 29th December 2020, a major earthquake (Mw 6.4) occurred in Croatia, close to Petrinja, 45 km south of Zagreb, generating widespread liquefaction and lateral spreading phenomena in a radius of approximately 20 km from the epicentre. A European team of researchers (geologists and engineers), in strict collaboration with the Croatian Geological Survey, performed field reconnaissance campaigns with the aim to provide a detailed identification and characterization of the primary and secondary geological and geotechnical coseismic effects induced by the Croatian earthquakes. Specifically with reference to the liquefaction phenomena, the Working Group integrated the data collected directly in the field with those from remote survey by drone aerial photos acquired in the post-event immediate. The adopted process allowed the collection of the liquefaction record with the highest possible completeness both in terms of pattern and distribution of the phenomena. The database includes several detailed case studies typified by the following characteristics: (1) liquefaction occurring on alluvial plain sites (Kupa river, Sava river and Glina river); (2) blows made by sand and/or gravel with local presence of shells and armored mud balls; (3) lateral spreading phenomena along road and river embankments; (4) sand ejecta of different grain size and matrix, even at the same site; (5) sand and/or gravel ejecta along fault traces. The characteristics of these features are discussed with reference to the alluvial setting and tectonic context. All together, the detailed survey of these recent liquefaction features will assist to build new empirical relations, to update the existing ones and to mitigate the effects of future earthquakes recognizing liquefaction prone areas for a correct land use planning, as for seismic microzonation studies.

Una ricostruzione geologica tridimensionale dell'area colpita dal terremoto dell'Italia centrale Sono stati pubblicati online e resi fruibili liberamente a tutta la comunità scientifica i risultati del Progetto RETRACE-3D (centRal italy EarThquakes integRAted Crustal model), che ha investigato il volume crostale in cui si è generato il terremoto di Amatrice del 24 agosto 2016 e la successiva sequenza sismica; un'area di circa 2.500 km2, a cavallo tra Lazio, Abruzzo, Umbria e Marche, interessata da oltre 118.000 eventi sismici tra il 2016 e il 2019. Il modello geologico 3D e il Rapporto finale sono il prodotto della collaborazione di più di 60 ricercatori ed esperti del Dipartimento della Protezione Civile e degli istituti di ricerca CNR-IGAG, CNR-IREA, INGV e ISPRA-Servizio Geologico d'Italia, Centri di Competenza che operano nel Servizio Nazionale della Protezione Civile. Al Progetto hanno inoltre collaborato Eni e Total, soggetti privati anch'essi parte del Servizio Nazionale della Protezione Civile, che hanno reso disponibili dati di sottosuolo e know-how. Il modello geologico 3D consente di mettere in evidenza la distribuzione e la geometria, nonché la complessità delle relazioni reciproche, delle più rilevanti unità geologiche e delle faglie principali, comprese quelle sismogeniche, fino a una profondità superiore a 10 km, riconciliando gli elementi osservati in superficie con le strutture in profondità. Il Rapporto finale e il modello geologico 3D sono rilasciati con licenza CC-BY 4.0, garantendo la possibilità di essere liberamente utilizzati per studi alla scala regionale con il riconoscimento della corretta attribuzione. Download dei risultati (http://www.retrace3d.it/contenuti.html) RETRACE-3D Central Italy Geological Model (2021). doi: 10.13127/retrace-3d/geomod.2021 RETRACE-3D Working Group (2021). RETRACE-3D: centRal italy EarThquakes integRAted Crustal modEl. Rapporto finale. Eds. INGV, ISPRA, CNR-IGAG, DPC. Roma, pp. 100. doi: 10.5281/zenodo.4604940

The Vettore–Bove normal fault system in central Italy ruptured during the 2016 MW 6.5 Norcia earthquake causing extensive surface faulting. At the Pian Grande di Castelluccio hanging wall basin, along the southern section of the fault ruptured during the MW 6.5 mainshock, we performed a high-resolution seismic reflection/refraction experiment aimed at (a) imaging the shallow pattern of the fault system, and (b) reconstructing the architecture of the continental infill. We collected three profiles for a total length of ∼8 km. We used a reflection processing flow and non-linear refraction tomography to obtain migrated stack sections and P-wave velocity images resolved down to the depth of the pre-Quaternary substratum. The main profile in the northern part of the basin crosses the westernmost splays of the ruptured fault zone striking N150°–170°. Seismic imaging unravels a ∼1 km-wide fault zone comprising three W-throwing splays and subsidiary faults, which affect the continental infill and produce a minimum aggregate Quaternary throw of ∼400 ± 100 m. Recent deformation is localized in this part of the surveyed fault section, attesting active displacement accumulation of the Vettore–Bove fault system. The other profiles in the central-southern part of the basin show additional faults, likely striking N20°–40° and which concurred to generate a ∼500 m-deep depocenter. These faults were mostly active during an early extensional phase; however, one of them likely displaces shallow layers with a throw close to the resolution limit of seismic data (

Active deformation in the Central Apennines is mostly accommodated by NW-SE normal faults systems that produce moderate to large earthquakes at shallow depth. Recent examples include the 1915 Mw≈7 Avezzano earthquake (Fucino basin) and the 2009 Mw=6.1 L’Aquila earthquake (Aterno basin) which were both associated with major loss of life and massive damage to buildings and infrastructure. Here, we study the 40-km-long Ovindoli – Piano di Pezza – Campo Felice – Monti d’Ocre (OPCM) fault system, a major NNW-SSE system that potentially links the Fucino and the Aterno fault-systems. The OPCM exhibits linear and arcuate sections with four main segments and borders the eastern margin of the Aterno basin. Paleo-earthquake rupture data on the Piano di Pezza (PPF) and Campo Felice (CFF) faults exhibit remarkable synchronicity with the Fucino fault system, with the most recent surface-rupturing earthquake likely occurring in the XIVth century. In order to better understand the relationships between and earthquake rupture scenarios, we focus on the basin geometry and fault surface expression of the Piano di Pezza fault combining geomorphology and subsurface geophysics. We map the fault trace with unprecedented detail using terrestrial laser scanner surveys and quantify surface deformation affecting alluvial fans as well as glacial moraines. We obtain a mean vertical offset of 2.5 m +/- 0.3 m for the most recent features, well in agreement with paleoseismological data. Furthermore, we document slip distributions at different time scales along strike with a maximum value at the connection between the PPF and the OF. Beneath the scarp, geophysical data reveal a complex faulting geometry with several parallel strands and two minor blind splays. Electrical resistivity tomography images show a cumulative vertical offset of ~ 15 m affecting an interface attributed to the Last Glacial Maximum and confirm the high vertical slip across the fault zone. Gravimetric anomalies across the basin also indicate the sedimentary fill has recorded a maximum finite cumulative throw of the PdP fault system of 110-140 m. This suggests a maximum vertical slip rate of 0.2-0.3 mm/year since the Pleistocene, which contrasts with the high post-LGM slip rate estimated from trenches. Overall, our observations suggest that the arcuate PPF originally formed as a reverse fault during the Mio-Pliocene compressive stage and is now reactivated as an extensional horsetail-like feature by ruptures along a major strike-slip fault (OF). This finding points to the PPF as mostly built through ruptures along the OF leaking onto an inherited structure. The time-varying slip rates may also denote an episodic behavior marked by short periods of high seismic activity (a few centuries) and long intervals of seismic quiescence (a few millennia). Furthermore, possible earthquake rupture scenarios along the OPCM may encompass the whole OPCM fault system (cumulative length ca. 40 km) or rupture termination along the PPF (cumulative length ca. 15-20 km) with significantly different impacts over the populated Fucino and Aterno basins.

The Castelluccio basin in the central Apennines (Italy) is a ~20-km2-wide intramontane Quaternary depression located in the hangingwall of the NW-trending and SW-dipping Vettore-Bove normal fault system (VBFS). This system is responsible for the 2016-2017 seismic sequence, culminated with the 30 October 2016 Mw 6.5 Norcia earthquake that caused widespread surface faulting affecting also the northern part of the Castelluccio basin. Available borehole and geophysical data are not enough to constrain the basin structure, infill architecture and their relations with the long-term activity of the VBFS. Therefore, we carried out an extensive 3D survey using the innovative Fullwaver (FW) technology, conceived to perform deep electrical resistivity tomography (DERT). We aimed at: a) mapping the geometry of the pre-Quaternary limestone basement and the basin infill thickness down to a depth of ~1 km; b) mapping the subsurface structure of known faults and their extent underneath the alluvial cover; c) mapping possible blind faults splays.The 3D survey covered a 23 km2 area and it was designed with the aim to map the region as regularly as possible, taking into account the rugged topography and logistic issues. We used a series of independent 2-channels receivers connected each to three grounded steel electrodes, 200 m spaced, to record the electrical field generated by a five kilowatt current regulated Time Domain Induced Polarization transmitter. Data were modelled with ViewLab software via a regularized inversion with smoothness constraints to cope with the expected subsurface strong resistivity changes, and to obtain a robust 3D resistivity model.The FW technology allowed us to constrain the geometry of the basin. The infill material is imaged as a wide, N-trending moderately resistive (< 300 Ωm) to conductive (< 100 Ωm) region, likely made of silty sands and gravels, deepening down to 500 m b.g.l. in the southern sector, suggesting the occurrence of two main depocenters. All over the basin, we identify paired high-resistivity (> 500-1000 Ωm) and low-resistivity (< 400 Ωm) belts related to the limestone basement and to the basin infill, respectively. They display NNE and NNW dominant trends. We interpret the sharp boundaries of NNE-trending belts as related to early extensional faults promoting the basin inception. The NNW-trending belts suggest the occurrence of faults that locally cross-cut the previous ones, and that we interpret as splays of the VBFS buried under the basin sedimentary cover. The recognition of different systems of extensional faults is coherent with results of high-resolution seismic profiling carried out recently in the basin. A high-resolution 2D transect with 15 m-spaced electrodes across the 2016 surface ruptures shows details of the active VBFS splay down to 300 m depth. Moreover, in the eastern sector of the survey area, low-resistivity round-shaped anomalies in the Mesozoic substratum hints for deep Miocene compressional structures. Therefore, our DERT imaging suggests a complex tectonics in the subsurface of the Norcia earthquake fault. In particular, the currently active NNW-trending faults seem to overprint a pre-existing structural framework, promoting fault segmentation at different spatial scales

We provide a database of the surface ruptures produced by the 26 December 2018 Mw 4.9 earthquake that struck the eastern flank of Mt. Etna volcano in Sicily (southern Italy). Despite its relatively small magnitude, this shallow earthquake caused about 8 km of surface faulting, along the trace of the NNW-trending active Fiandaca Fault. Detailed field surveys have been performed in the epicentral area to map the ruptures and to characterize their kinematics. The surface ruptures show a dominant right-oblique sense of displacement with an average slip of about 0.09 m and a maximum value of 0.35 m. We have parsed and organized all observations in a concise database, with 932 homogeneous georeferenced records. The Fiandaca Fault is part of the complex active Timpe faults system affecting the eastern flank of Etna, and its seismic history indicates a prominent surface-faulting potential. Therefore, this database is essential for unravelling the seismotectonics of shallow earthquakes in volcanic areas, and contributes updating empirical scaling regressions that relate magnitude and extent of surface faulting., Measurement(s)coseismic surface rupture • surface rupture kinematics • surface rupture displacement • surface rupture locationTechnology Type(s)field survey • GPS navigation systemFactor Type(s)offset • strike • angle • length • latitude • longitude • elevationSample Characteristic - Environmentvolcanic fieldSample Characteristic - LocationIsland of Sicily • Mount Etna Machine-accessible metadata file describing the reported data: 10.6084/m9.figshare.11673027

The Mediterranean area has been inhabited since long time and the development of human settlements along its coastlines allowed to “record” the occurrence of natural disasters as tsunamis in different coastal settings. Nowadays, the same coastal areas host large cities and critical infrastructures and the availability of historical and geological data on the tsunamis of the past, on their characteristics and effects can significantly contribute to the development of adequate mitigation measures, defense actions and better awareness of people living along the coastlines. The knowledge of past tsunamis is one of the keys to better forecast future events: our territory is a natural laboratory where past natural events can replicate in the future with similar characteristics. All the coastlines may be considered as the geological and geomorphological archive of past inundations. In the Mediterranean area, valuable historical evidence of past tsunamis and several paleotsunami traces exist, these were recently well organized in homogeneous and coherent geographical databases; these represent the core of this manuscript. This type of information, recently improved and enriched, allows to increase the knowledge of tsunami histories of the Mediterranean sites most prone to tsunamis, providing information on the intensity and frequency of past events.

In the engineering geology field increased attention has been posed in recent years to potential liquefaction mitigation interventions in susceptible sand formations. In silty sands this is a major challenge because, as the fines content increases, vibratory methods for densification become progressively less effective. An alternative mitigation technique can be the installation of Rammed Aggregate Pier® (RAP) columns that can increase the resistance of the soil, accounting for its lateral stress increase and for the stiffness increase from soil and RAP composite response. To investigate the influence of these factors on liquefaction resistance, full-scale blast tests were performed at a silty sand site in Bondeno (Ferrara, Italy) where liquefaction was observed after the 2012 Emilia-Romagna earthquake. A multidisciplinary team of forty researchers carried out devoted experimental activities aimed at better understanding the liquefaction process at the field scale and the effectiveness of the treatment using inter-related methods. Both natural and improved areas were investigated by in-situ tests and later subjected to controlled blasting. The blast tests were monitored with geotechnical and geophysical instrumentation, topographical surveying and geological analyses on the sand boils. Results showed the RAP effectiveness due to the improvement of soil properties within the liquefiable layer and a consequent reduction of the blast-induced liquefaction settlements, likely due to soil densification and increased lateral stress. The applied multidisciplinary approach adopted for the study allowed better understanding of the mechanism involved in the liquefaction mitigation intervention and provided a better overall evaluation of mitigation effectiveness.

Among the geologic data, trenching records of paleoearthquakes represent an important input for the development of seismic hazard evaluations and, specifically, for the characterization of fault rupture behavior in time. Based on the available paleoseismological trenching data in the Central Apennines (Italy) we have compiled a new database of surface faulting earthquakes for 10 faults. The compilation contains 109 surface faulting events, occurred in the past ~28 kyr. Events younger than 7 kyr (59 dated events) are much better constrained and, therefore, are the basis for most of the analyses. Through a quantitative multistep method, we integrate paleoseismic trench data and treat them statistically, providing relevant improvement needed for trenching-based seismic hazard evaluation. Indeed, the combined analysis of trenching data from multiple sites on the same fault led to a reduction of the approximations and uncertainties in the rupture history of individual faults. The procedure was also applied on specific fault systems and indicated that the single faults may have occasionally ruptured simultaneously, or close in time, in the past. The whole set of age ranges was also statistically analyzed to produce regional earthquake scenarios for a period much longer than the millennial historical record. The built scenarios for the last 7 kyr define a regional mean inter event time (IET) of 230–240 yr, comparable with the average recurrence time of about 200 yr calculated for the region. We also identify the possible occurrence of earthquake storms, i.e., concentration of surface faulting earthquakes in the region within time periods of 200 yr, suggesting a variability in time of the seismic behavior of the faults, with alternation of peaks of activity with more “quiet” periods. Even though the paleoseismic data from the Apennines contains uncertainties, the results appear quite stable and promising for future applications in earthquake-hazard assessment.

We investigated the late Quaternary throw distribution of the main normal fault that ruptured during the Mw 6.5 2016 earthquake in central Italy by means of a high-resolution structure-from-motion (SfM)-derived Digital Surface Model (DSM). We focused on a key area along the Cordone del Vettore fault (CDV), which is part of the Vettore-Bove fault system (VBFS). The CDV displays a prominent compound post-glacial scarp that allowed the reconstruction of the along-strike cumulative throw distribution. We propose a geometric approach to calculate the CDV fault throw distribution from the reconstruction of a displaced glaciation-related erosional surface, used as a geomorphic marker, and a series of closely spaced cross profiles. The proposed calculation accounts for both the slip vector direction and the degraded scarp top, including field data on fault dip angles. Following this approach, we recognized two scarps with a minimum average fault throw of ~21 m and ~35 m for this section of the investigated fault strand. The correlation with the possible post-LGM (Last Glacial Maximum) deglaciation phases of the erosional surface suggests a minimum scarp age of 25–27 ka cal BP. Such an age provides a reasonable CDV fault throw rate of ~0.8 mm/a, comparable with known long-term throw rates of the VBFS and active Apennines normal faults. By comparing the reconstructed long-term Cordone del Vettore throw distribution with the 2016 coseismic one, ~24 2016-like surface faulting events are required to generate the main cumulative scarp, under the assumption of constant slip per event. This, along with the age of the scarp, yields an average earthquake recurrence time interval of ~1100 a. These results suggest the presence of multiple regional markers that correlates with different LGM (if not pre-LGM) major glacial phases, whose erosional processes allow the preservation of pre-existing bedrock fault scarp remnants.

On 24 August 2016, a Mw 6.0 normal faulting earthquake struck central Italy, causing about 300 fatalities and heavy damage. A geological survey collected the coseismic effects observed at the surface in order to evaluate two competing hypotheses about their nature: surface faulting versus gravitational deformation. We find that the most significant geological effect is a 5.2 km-long alignment of ground ruptures along the Mt. Vettore Fault-System. These ruptures are independent from lithology, topography, morphology and change in slope and exhibit an average dip slip displacement of ~13 cm. Geometry, kinematics and dimensional properties of this zone of deformation strongly lead us to favor the primary surface faulting hypothesis that fits well the predicted estimates from experimental scaling-law relationships. Our study provides relevant hints for surface faulting in extensional domains, contributing to implement the worldwide database of the moderate earthquakes.

An Mw 6.1, devastating earthquake, on April 6, 2009, struck the Middle Aterno Valley (Abruzzi Apennines, Italy) due to the activation of a poorly known normal fault system. Structural analysis of the fault population and investigation of the relationships with the Quaternary continental deposits through integrated field and laboratory techniques were conducted in order to reconstruct the long-term, tectono-sedimentary evolution of the basin and hypothesize the size of the fault segment. A polyphasic evolution of the Middle Aterno Valley is characterized by a conjugate, ∼E-W and ∼N S-striking fault system, during the early stage of basin development, and by a dip-slip, NW-striking fault system in a later phase. The old conjugate fault system controlled the generation of the largest sedimentary traps in the area and is responsible for the horst and graben structures within the basin. During the Early Pleistocene the E-W and N S system reactivated with dip-slip kinematics. This gave rise to intra-basin bedrock highs and a significant syn-tectonic deposition, causing variable thickness and hiatuses of the continental infill. Subsequently, since the end of the Early Pleistocene, with the inception of the NW-striking fault system, several NW-strands linked into longer splays and their activity migrated toward a leading segment affecting the Paganica-San Demetrio basin: the Paganica-San Demetrio fault alignment. The findings from this work constrain and are consistent with the subsurface basin geometry inferred from previous geophysical investigations. Notably, two major elements of the ∼E-W and ∼N S-striking faults likely act as transfer to the nearby stepping active fault systems or form the boundaries, as geometric complexities, that limit the Paganica-San Demetrio fault segment overall length to 19 ± 3 km. The resulting size of the leading fault segment is coherent with the extent of the 6 April 2009 L'Aquila earthquake causative fault. The positive match between the geologic long-term and coseismic images of the 2009 seismogenic fault highlights that the comprehensive reconstruction of the deformation history offers a unique contribution to the understanding faults seismic potential.

International audience; We integrate paleoseismic data sets along the Mt. Vettore‐Mt. Bove normal fault system rupturing at the surface in the 30 October 2016 Norcia earthquake. Through the analysis of new trenches from this work and a review of the preexisting data, we correlate events among trench sites along antithetic and synthetic fault splays. We recognize seven M 6.5, 2016 Norcia‐type (or larger) surface‐faulting events in the last ~22 kyr, including 2016. Before 2016, one event occurred in the past two millennia (260–575 CE) and possibly corresponds to the event damaging Rome in 443 or 484/508 CE. Three previous events occurred between 10590 and 415 BCE, whereas the two oldest ones date between 19820 and 16540 BCE. The average recurrence time is 3,360–3,640 years for the last ~22 kyr and 1,220–1,970 years for the last ~4 kyr. We infer a minimum dip‐slip rate of 0.26–0.38 mm/year on the master fault in the central portion of the Mt. Vettore–Mt. Bove normal fault system and a dip‐slip rate of at least 0.10 mm/year on the southernmost portion. We infer a Middle–Late Pleistocene inception of the long‐term scarp of the investigated splays. The along‐strike variation of slip rates well reproduces the trend of the 2016 surface slip; thus, the time window exposed in the trenches is representative for the present fault activity. Based on trenching data, different earthquake rupture scenarios should be also considered for local hazard assessment.

We present a 1:20,000 scale map of the traces of the active Capitignano and San Giovanni faults in the area of the Montereale basin (central Apennines, Italy) covering an area of about 80 km2. Detailed fault mapping is based on high-resolution topography from airborne LiDAR imagery validated by extensive ground truthing and geophysical prospecting. Our analysis allowed the recognition of several features related to fault activity, even in scarcely accessible areas characterized by dense vegetation cover and rugged terrain. The identified fault traces run at the base of the NW-SE striking Montereale basin-bounding mountain front and along the base of the southwestern slope of the Monte Mozzano ridge, and have a length of about 12 and 8 km, respectively. Improving the knowledge of fault geometry is a critical issue not only for the recognition of seismogenic sources but also for surface fault hazard assessment and for local urban planning. The knowledge of the exact location of the fault traces is also crucial for the seismogenic characterization of the active faults by means of paleoseismological trenching.

We report a case study from the Po River plain region (northern Italy), where significant liquefaction-related land and property damage occurred during the 2012 Emilia seismic sequence. We took advantage of a 1 m pixel lidar digital terrain model (DTM) and of the 2012 Emilia coseismic liquefaction data set to (a) perform a detailed geomorphological study of the Po River plain area and (b) quantitatively define the liquefaction susceptibility of the geomorphologic features that experienced different abundance of liquefaction. One main finding is that linear topographic highs of fluvial origin – together with crevasse splays, abandoned riverbeds and very young land reclamation areas – acted as a preferential location for the occurrence of liquefaction phenomena. Moreover, we quantitatively defined a hierarchy in terms of liquefaction susceptibility for an ideal fluvial environment. We observed that a very high liquefaction susceptibility is found in coincidence with fluvial landforms, a high-to-moderate liquefaction susceptibility within a buffer distance of 100 and 200 m from mapped fluvial landforms and a low liquefaction susceptibility outside fluvial landforms and relative buffer areas. Lidar data allowed a significant improvement in mapping with respect to conventionally available topographic data and/or aerial imagery. These results have significant implications for accurate hazard and risk assessment as well as for land-use planning. We propose a simple geomorphological approach for liquefaction susceptibility estimation. Our findings can be applied to areas beyond Emilia that are characterized by similar fluvial-dominated environments and prone to significant seismic hazard.

We present geological evidence for paleotsunamis along the ~230 km-long coast of eastern Sicily (Italy); combining this information with historical data, we reconstruct a unique history of tsunami inundations. We integrate data on 38 paleotsunami deposits (from fine sand layers to boulders) collected at 11 sites (one offshore). The geological data record traces of large tsunamis which have occurred during the past 4 millennia. Chronological constrains include 14C, 210Pb and 137Cs, OSL and tephrochronology. When compatible, the age of the paleotsunami deposits is associated to historical events, but it is also used to highlight unknown tsunamis. Average tsunami recurrence interval (between 320 and 840 yr) and minimum inland tsunami ingressions (often greater than the historical ones) were estimated at several sites. On the basis of this work, the tsunami catalogue is implemented by two unknown tsunamis which occurred during the first millennium BC and by one unknown regional tsunami, which occurred in 650–770 AD. By including this latter event in the eastern Sicily catalogue, we estimate an average recurrence interval for strong tsunamis of ca. 385 yr. Comparison and merging of historical and geological data can definitely contribute to a better understanding of regional and local tsunami potential and provides robust parameters to be used in tsunami hazard estimates.

Two cores were sampled in the Fucino Basin (central Apennines, Italy), which represents an extensional intramountain basin filled by Pliocene to Quaternary continental alluvial and lacustrine deposits. The cores were investigated for tephra content and five visible tephras with thickness ranging from 1 to 8 cm were identified. Six additional cryptotephra were identified during the inspection of significant peaks of the magnetic susceptibility curve. Texture and mineralogy of five tephra and six cryptotephra layers were analyzed by means of scanning electron microscope coupled with energy-dispersive X-ray spectrometry system (SEM-EDS) and geochemical measurements were performed by an electron microprobe (EPMA) equipped with five wavelength-dispersive spectrometers (WDS) and using a laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) system on single glass shards. The results allowed us to assign tephra and cryptotephra to ten known volcanic eruptions that occurred over the last ca. 60 ka in the Campanian Province (Phlegrean Fields and Ischia Island), the Alban Hills volcanic complex, and Lipari island. In particular, we recognized the deposits of the Monte Epomeo Green Tuff and the Piroclastiti di Catavola eruptions of Ischia, the pre-Campanian Ignimbrite Tlc, the Campanian Ignimbrite and the Neapolitan Yellow Tuff eruptions of the Phlegrean Fields, the Gabellotto-Fiume Bianco eruption of Lipari, and all the four explosive events belonging to the last cycle of volcanic activity of Albano maar (Albano 4-7). Deposits from five of these identified events (i.e., Piroclastiti di Catavola, Gabellotto-Fiume Bianco, Albano 5 and 6 eruptions, and Campanian Ignimbrite) were previously un-reported in the Fucino basin. These findings add new tephra layers to the list of possible tephrochronologic markers in the region and highlight that a comprehensive tephra record may be constructed when the study of cryptotephra layers is included. Moreover, results provide insights into the most recent volcanic activity of Albano maar, allowing us to date the onset of activity at the maar system at ca. 40 ka and to estimate the ages of all four eruptions that made up this eruptive sequence at ca. 37.5 ka (Albano 5), ca. 36.5 ka (Albano 6) and ca. 36 ka (Albano 7), respectively. Our work extends the known dispersal of several major explosive events, suggesting the intensity and magnitude appraisals, and attended risk scenario's need to be revised using improved records of distal fall out.

In this paper we present a morphotectonic study of the Paganica–San Demetrio fault system (PSDFS) responsible for the Mw6.1 April 6, 2009 earthquake (L'Aquila, Central Italy). The discrepancy observed between the length of the seismologic–geodetic modeled fault, the limited size of the primary coseismic surface ruptures and the significant morphological expression of the PSDFS stimulated a debate about the maximum rupture length of the PSDFS and its capability to generate larger magnitude events. To image the PSDFS long-term morphological expression and define its surface geometrical arrangement (length, number of fault splays and boundaries), we took advantage of a high-resolution airborne LiDAR dataset. LiDAR topography substantially increased our confidence in detecting even subtle tectonic-controlled morphologies. We define the PSDFS as a ~ 19 km-long fault system that displays a complex structural setting characterized by two different sectors: 1) the Paganica sector to the NW, with a narrow deformation zone, and 2) the San Demetrio sector to SE, where the strain is accommodated by several fault-splays dissecting a wider Quaternary basin. We also defined a first-order hierarchy among the numerous fault splays across the PSDFS. The long-term geomorphic expression of the PSDFS suggests that it ruptured also involving the whole 19 km-long structure besides rupturing only small sections, as it occurred in 2009. This suggests a variable slip behavior. Empirical relations applied to this hypothesis allow up to M 6.6 earthquakes along the PSDFS. These results have a critical impact on the seismic hazard assessment of the area when compared with a M 6.1 event as the 2009.

FIRB Abruzzo: High-resolution analyses for assessing the seismic hazard and risk of the areas affected by the 2009 April 6 earthquake’, No RBAP10ZC8K_005

We present a 2-D subsurface image of the Paganica Fault from a high-resolution refraction tomography and detailed geological investigation carried out across part of the northwestern segment of the 20-km-long Paganica-San Demetrio fault-system, and which was responsible of the 2009 April 6 Mw 6.1 L'Aquila earthquake (central Italy). We acquired two seismic profiles crossing the Paganica basin with a dense-wide aperture configuration. More than 30 000 P wave first-arrival traveltimes were input to a non-linear tomographic inversion. The obtained 250-300 m deep 2-D Vp images illuminate the shallow portion of the Paganica Fault, and depict additional unreported splays defining a complex half-graben structure. We interpret local thickening of low-Vp (

We applied a joint survey approach integrating time-domain electromagnetic soundings and single-station ambient vibration surveys in the Middle Aterno Valley (MAV), an intermontane basin in central Italy and the locus of the 2009 L'Aquila earthquake. By imaging the buried interface between the infilling deposits and the top of the pre-Quaternary bedrock, we reveal the 3-D basin geometry and gain insights into the long-term basin evolution. We reconstruct a complex subsurface architecture, characterized by three main depocenters separated by thresholds. Basin infill thickness varies from ~200-300 m in the north to more than 450 m to the southeast. Our subsurface model suggests a strong structural control on the architecture of the basin and highlights that the MAV experienced considerable modifications in its architecture over time. The buried shape of the MAV suggests a recent and still ongoing predominant tectonic control by the NW-SE trending Paganica - San Demetrio Fault System (PSDFS), which cross-cuts older ca. ENE- and NNE- trending extensional faults. Furthermore, we postulate that the present-day arrangement of the PSDFS is the result of the linkage of two previously isolated fault segments. We provide constraints on the location of the southeastern boundary of the PSDFS, defining an overall ~19 km-long fault system. This allows for earthquakes larger than M6.5. This study emphasizes the benefit of combining two easily-deployable geophysical methods for reconstructing the 3-D geometry of a tectonically-controlled basin. Our joint approach provided us with a consistent match between these two independent estimations of the basin substratum depth within 15%.

In this paper we present the geological effects induced by the 2012 Emilia seismic sequence in the Po Plain. Extensive liquefaction phenomena were observed over an area of ~ 1200 km2 following the 20 May, ML 5.9 and 29 May, ML 5.8 mainshocks; both occurred on about E–W trending, S dipping blind thrust faults. We collected the coseismic geological evidence through field and aerial surveys, reports from local people and Web-based survey. On the basis of their morphologic and structural characteristics, we grouped the 1362 effects surveyed into three main categories: liquefaction (485), fractures with liquefaction (768), and fractures (109). We show that the quite uneven distribution of liquefaction effects, which appear concentrated and aligned, is mostly controlled by the presence of paleo-riverbeds, out-flow channels and fans of the main rivers crossing the area; these terrains are characterised by the pervasive presence of sandy layers in the uppermost 5 m, a local feature that, along with the presence of a high water table, greatly favours liquefaction. We also find that the maximum distance of observed liquefaction from the earthquake epicentre is ~ 30 km, in agreement with the regional empirical relations available for the Italian Peninsula. Finally, we observe that the contour of the liquefaction observations has an elongated shape almost coinciding with the aftershock area, the InSAR deformation area, and the I ≥ 6 EMS area. This observation confirms the control of the earthquake source on the liquefaction distribution, and provides useful hints in the characterisation of the seismogenic source responsible for historical and pre-historical liquefactions.

A sedimentary sequence that was highly sensitive to fault rupture–driven changes in water level and sediment supply has been used to extract a continuous record of 22 large earthquakes on the Alpine fault, the fastest-slipping fault in New Zealand. At Hokuri Creek, in South Westland, an 18 m thickness of Holocene sediments accumulated against the Alpine fault scarp from ca. A.D. 800 to 6000 B.C. We used geomorphological mapping, sedimentology, and paleoenvironmental reconstruction to investigate the relationship between these sediments and Alpine fault rupture. We found that repeated fault rupture is the most convincing mechanism for explaining all the features of the alternating peat and silt sedimentary sequence. Climate has contributed to sedimentation but is unlikely to be the driver of these cyclical changes in sediment type and paleoenvironment. Other nontectonic causes for the sedimentary alternations do not produce the incremental increase in basin accommodation space necessary to maintain the shallow-water environment for 6800 yr. Our detailed documentation of this near-fault sedimentary basin sequence highlights the advantages of extracting paleoearthquake records from such sites—the continuity of sedimentation, abundance of dateable material, and pristine preservation of older events.

Pucci et. al., Three 2-D Deep Electrical Resistivity Tomography (ERT) transects, up to 6.36 km long, were obtained across the Paganica-San Demetrio Basin, bounded by the 2009 L'Aquila M-w 6.1 normal-faulting earthquake causative fault (central Italy). The investigations allowed defining for the first time the shallow subsurface basin structure. The resistivity images, and their geological interpretation, show a dissected Mesozoic-Tertiary substratum buried under continental infill of mainly Quaternary age due to the long-term activity of the Paganica-San Demetrio normal faults system (PSDFS), ruling the most recent deformational phase. Our results indicate that the basin bottom deepens up to 600 m moving to the south, with the continental infill largely exceeding the known thickness of the Quaternary sequence. The causes of this increasing thickness can be: (1) the onset of the continental deposition in the southern sector took place before the Quaternary, (2) there was an early stage of the basin development driven by different fault systems that produced a depocentre in the southern sector not related to the present-day basin shape, or (3) the fault system slip rate in the southern sector was faster than in the northern sector. We were able to gain sights into the long-term PSDFS behaviour and evolution, by comparing throw rates at different timescales and discriminating the splays that lead deformation. Some fault splays exhibit large cumulative throws (> 300 m) in coincidence with large displacement of the continental deposits sequence (> 100 m), thus testifying a general persistence in time of their activity as leading splays of the fault system. We evaluate the long-term (3-2.5 Myr) cumulative and Quaternary throw rates of most of the leading splays to be 0.08-0.17 mm yr(-1), indicating a substantial stability of the faults activity. Among them, an individual leading fault splay extends from Paganica to San Demetrio ne' Vestini as a result of a post-Early Pleistocene linkage of two smaller splays. This 15 km long fault splay can explain the Holocene surface ruptures observed to be larger than those occurred during the 2009 L'Aquila earthquake, such as revealed by palaeoseismological investigations. Finally, the architecture of the basin at depth suggests that the PSDFS can also rupture a longer structure at the surface, allowing earthquakes larger than M 6.5, besides rupturing only small sections, as it occurred in 2009., This work was financially supported by Project ‘Fondo per gli Investimenti della Ricerca di Base (FIRB) Abruzzo: High-resolution analyses for assessing the seismic hazard and risk of the areas affected by the 2009 April 6 earthquake’, No. RBAP10ZC8K_005.

Offshore and inland geological evidence for multiple tsunami inundations was found in the Augusta Bay area: (1) the main local historical tsunamis (1908, 1693, 1169), (2) two far-generated tsunamis (i.e. Crete 365 AD and Santorini, 3600 BP), and (3) seven unknown tsunamis). Average tsunami recurrence intervals from inland and offshore investigations of about 550 and 320 yr, respectively were obtained for the past 4 ka. The history of paleotsunamis from the marine record appears to be as complete as the historical one for the past millennium, yielding an average tsunami recurrence interval of 250 yr for the Augusta Bay. Geological data allow also estimating a minimum tsunami inundation distance of 530 m and a minimum run-up of 5 m. The marine record contains evidence for more paleotsunamis with respect to the inland one because of continuous sedimentation and better preservation of stratigraphy in the offshore with respect to coastal areas, which are commonly affected by intermittent-erosion and sedimentation and anthropic activities. This work shows that the integration of geological and historical data can provide critical information regarding the extent and age of tsunamis of the past (e.g. inundation distance, age, and frequency), which is of immediate relevance for tsunami hazard assessment.

Twelve anomalous layers, marked by a high concentration of displaced epiphytic foraminifera (species growing in vegetated substrates like the Posidonia oceanica ) and subtle grain-size changes were found in a 6.7 m long, fine sediment core (MS-06), sampled 2 km off the shore of the Augusta Harbor (Eastern Sicily) at a depth of 72 m, recording the past 4500 yrs of deposition. Because concentrations of epiphytic foraminifera are quite common in infralittoral zones, but not expected at − 72 m, we believe that these anomalous layers might be related to the occurrence of tsunamis causing substantial uprooting and seaward displacement of P. oceanica blades with their benthic biota. Our approach involved the study of geophysical data (morphobathymetry, seismic reflection, and seafloor reflectivity) and sediment samples, including X-ray imaging, physical properties, isotopic dating, tephrochronology, grain-size and micropaleontology. Correlations between anomalous layers and tsunami events have been supported by a multivariate analysis on benthic foraminifera assemblage and dates of historical tsunami records. We found that four out of the eleven layers were embedded in age intervals encompassing the dates of major tsunamis that hit eastern Sicily (1908, 1693, and 1169) and the broader Eastern Mediterranean (Santorini at about BP 3600). One more layer, even if less distinct than the others, was also defined and may be the evidence for the AD 365 Crete tsunami.

SUMMARY We present a high-resolution palaeomagnetic and rock magnetic study of two cores, MS06 and MS06-SW (6.7 and 1.1 m long, respectively), collected at 72 m below sea level in the Augusta Bay shelf (Eastern Sicily, Ionian Sea, Italy) about 2.3 km from the coastline. Geophysical surveying carried out in the sampling area highlighted the presence of a homogeneous sedimentary sequence that most likely was deposited after the Last Glacial Maximum and was not affected by anthropogenic disturbances. The two cores penetrated a monotonous mud sedimentary sequence, interrupted at ∼3 m depth by a 3–4-cm-thick volcanic sandy layer that is correlated with the tephra fallout deposit produced by the 122 BC plinian eruption of Mt Etna. This tephra, along with radiocarbon dating of nine marine shells and with radioactive tracers for the uppermost 0.3 m (210Pb and 137Cs), provide the chronological constraints for the stratigraphic sequence that resulted younger than 4500 yr BP. Palaeomagnetic and rock magnetic data show that the sample sequence is magnetically homogeneous. A single peak of high magnetic mineral concentration is present and corresponds to the volcanic sandy layer. Palaeomagnetic data allowed the identification of a well-defined characteristic remanent magnetization that provides a high-resolution record of palaeosecular variation (PSV) at the sampling site. The reconstructed PSV curve is in good agreement with the available regional reference PSV curves and with the prediction from recent PSV modelling for Europe. The palaeomagnetic data obtained in this study on the one hand support and refine the age model for the cores, derived from other independent constraints, and on the other hand provide an original high-resolution PSV curve that can serve as a reference for the central Mediterranean over the last 4 ka.

We present the geological evidence for a 4000 year long record of multiple tsunami inundations along the coast of the Augusta Bay (eastern Sicily) and discuss its implications. The research was carried out through a multi-theme approach which benefited from an extraordinarily long historical record that we used to guide detailed geomorphologic and geologic surveys, coring campaigns and laboratory analyses. Two sites, named the Augusta Hospital and Priolo Reserve, were selected and investigated in detail along the 25 km-long coastline of Augusta Bay. We found evidence for six (possibly seven) tsunami deposits; three of them may be tentatively associated with the 1693 and 365 AD Ionian Sea historical tsunamis and the ~ 3600 BP Santorini event. The other three (possibly four) deposits are evidence for unknown paleo-inundations dated at about 650–770 AD, 600–400 BC and 975–800 BC (at Augusta Hospital site), and 800–600 BC (at Priolo Reserve site). We use these ages to extend further back the historical record of tsunamis available for this coastal area. The exceptional number of tsunami deposits found with this study allowed us to derive an average geologic tsunami recurrence interval in the Augusta Bay of about 600 years for the past 4 ka. Conversely, the historical tsunami data for the past millennium suggest an average tsunami recurrence interval of about 250 years. This difference in the average recurrence intervals suggests that only the strongest inundations may leave recognizable geological signatures at the investigated sites (i.e. the evidence for the 1908 and 1169 tsunamis is missing) but also that the geomorphological setting of the site and its erosional/depositional history are critical aspects for the data recording. Thus, an average recurrence interval derived from the geological record should be considered as a minimum figure. The identification and age estimation of tsunami deposits represent a new and independent contribution to tsunami scenarios and modeling for coastal hazard assessment in Civil Protection applications. Furthermore, our study cases provide new elements on tsunami deposit recognition related to exceptionally large events that occurred in the Aegean Sea.

Eastern Sicily is an area where some of the most catastrophic earthquakes in Italian history occurred. As reported by historical sources, these earthquakes induced liquefaction phenomena. In two areas along the Ionian coast of Sicily we found liquefaction evidence in Holocene deposits. In the Minissale site (east of Mount Etna), the observed liquefaction features can be related to the 1169 and 1693 earthquakes, while in the Agnone site (south of Catania), the deformational structures can be tentatively associated to the 1542 and 1693 events. Both sites locate in areas where historical liquefaction has been observed in the past, thus confirming the actual liquefaction susceptibility/potential in this region. The evidence of repeated liquefaction events at the same site, highlights the relevance of the paleoseismic approach for modeling the recurrence time and preparing scenarios of seismic effects in eastern Sicily, where seismogenic sources are scantly defined.

The source mechanisms responsible for large historical tsunamis that have struck eastern Sicily and southern Calabria are a topic of robust debate. We have compiled a database of historical coeval descriptions of three large tsunamis: 11 January 1693, 6 February 1783, and 28 December 1908. By using accounts of run-up and inundation and employing an approach proposed by Okal and Synolakis in 2004, we can provide discriminants to define the nature of the near-field tsunami sources (fault dislocation or landslide). Historical reports for the 1908 event describe affected localities, maximum run-ups, and inundation areas. However, for the 1693 and 1783 tsunamis, reports are limited to inundation and occasional run-up estimates. We calculate run-up values for these events using available relations between inundation and run-up. We employed the model of Okal and Synolakis to the obtained profiles of tsunami run-up along the inundated shorelines. The 1908 run-up data distribution confirms that the tsunami is compatible with a seismic dislocation source, whereas the 1783 data supports contemporary observations and recent offshore investigations suggesting that the tsunami was produced by an earthquake-triggered submarine landslide. Analysis of the 1693 event data suggests that the tsunami was generated during a tectonic event and thus a seismogenic source should be found offshore.