Slip-rate variations over multiple seismic cycles play a fundamental role in controlling the behaviour of active fault systems, as they are linked to spatio-temporal earthquake clustering and can influence the recurrence patterns of adjacent faults. However, processes that produce slip-rate fluctuations are yet to be fully defined. Despite their importance, the physical mechanisms responsible for such slip-rate fluctuations remain only partially understood. In this study, we investigate whether interactions between neighbouring along-strike brittle faults and their underlying viscous shear zones can generate slip-rate variability associated with synchronous earthquake clustering and fault system synchronization. We focus on nine normal faults and related shear zones within the Central Apennines fault system (Italy), arranged in six along-strike fault pairs characterized by different fault spacings and strike geometries. We integrate cosmogenic 36Cl dating of tectonically exhumed fault scarps with numerical modelling of differential stress transfer between interacting fault–shear-zone pairs. The results identify a mechanism capable of producing simultaneous earthquake clusters, driven by the synchronization of high driving stresses within the viscous shear zones beneath the brittle faults. This behaviour is strongly modulated by along-strike fault spacing and strike variations. In settings with closely spaced fault pairs and limited strike variations, earthquake clusters induce positive differential stress variations on neighbouring shear-zones of sufficient magnitude to induce positive slip-rate variations on their overlying brittle faults. This produces positive feedback mechanism that sustains the occurrence of earthquake clusters that will continue to positively load the neighbouring shear zones. These findings provide new insights into fault system dynamics across multiple timescales and have important implications for seismic hazard evaluation.

Paleoseismology is a vital tool for the study of earthquake hazard and active tectonics. Its application in the context of Late Quaternary basaltic volcanoes encounters considerable limitations due to the inherent highly dynamic nature of such systems. Etna volcano, however, provides an ideal setting for such studies. In particular, the densely populated Mt. Etna eastern flank record frequent surface faulting earthquakes and aseismic fault creep, which result in significant offsets of well-dated historical landforms and stratigraphy, including lava flows, interlayered pyroclastic deposits, and anthropic structures. This allows for the analysis of fault slip rates across various time scales.We present the first paleoseismological results along the Fiandaca Fault, the source of the 26 December 2018, Mw 4.9 Fleri earthquake. We excavated two exploratory trenches along the coseismic surface ruptures at the Collegio Fiandaca site. Analysis of trench walls allow identifying, besides the 2018 event, two historical surface faulting events. The youngest one occurred in the period 1281-1926 CE, and most likely during the 8 August 1894 Fiandaca earthquake. The oldest one, previously unknown, occurred in the Early Middle Ages (757-894 CE). This paleoseismic evidence strongly suggest increased seismic activity along the Fiandaca Fault in the last centuries. In order to verify this hypothesis, we conducted detailed morphotectonic analyses and throw rate measurements along the Fiandaca and other capable normal faults in the Mt. Etna eastern flank. Throw rates mean values show an increase from 1.4 mm/yr during the 15-3.9 ka time interval to 3.4 mm/yr between 3.9 ka and the Greek-Roman period, with a further increase since the late Middle Ages, reaching 10 mm/yr. This trend suggests a very recent growth in flank instability, in agreement with current geodetic data but also with historical eruptive activity.These findings highlight an increase of the associated geological hazards along the inhabited eastern flank, emphasizing the need for further research and a multi-hazard approach to risk assessment and land planning for Mt. Etna and similar volcanic regions.

The 1688 Sannio–Matese earthquake, with a macroseismically derived magnitude of Mw = 7 and an epicentral intensity of IMCS = XI, had a deep impact on Southern Italy, causing thousands of casualties, extensive damage and significant environmental effects (EEEs) in the epicentral area. Despite a comprehensive knowledge of its economic and social impacts, information regarding the earthquake’s environmental effects remains poorly studied and far from complete, hindering accurate intensity calculations by the Environmental Seismic Intensity Scale (ESI-07). This study aims to address this knowledge gap by compiling a thorough dataset of the EEEs induced by the earthquake. By consulting over one hundred historical, geological and scientific reports, we have collected and classified, using the ESI-07 scale, its primary and secondary EEEs, most of which were previously undocumented in the literature. We verified the historical sources regarding some of these effects through reconnaissance field mapping. Analysis of the obtained dataset reveals some primary effects (surface faulting) and extensive secondary effects, such as slope movements, ground cracks, hydrological anomalies, liquefaction and gas exhalation, which affected numerous towns. These findings enabled us to reassess the Sannio earthquake intensity, considering its environmental impact and comparing traditional macroseismic scales with the ESI-07. Our analysis allowed us to provide an epicentral intensity ESI of I = X, one degree lower than the published IMCS = XI. This study highlights the importance of combining traditional scales with the ESI-07 for more accurate hazard assessments. The macroseismic revision provides valuable insights for seismic hazard evaluation and land-use planning in the Sannio–Matese region, especially considering the distribution of the secondary effects.

We present slip versus time histories derived from in situ 36 Cl cosmogenic dating for three active normal faults in the southern Apennines, Italy. In this region the total extensional strain is accommodated by either a small number of faults located across strike from each other or, in places, a single fault where no other active faults exist across strike. We investigate how strain‐rates on individual faults vary through time in the context of the overall geometry of the fault system. The 36 Cl results confirm that the San Gregorio Magno, Auletta, and Vallo di Diano faults were active in the Holocene, with each fault exhibiting alternating periods of relatively rapid and slow, or even absence of, slip. During periods of rapid slip, lasting a few millennia, the faults accumulate up to ∼5 m of slip, which we interpret as earthquake clusters. At other times, the faults exhibit no slip for time periods lasting multiple millennia. The fluctuations in slip‐rates reveal the migration of activity between faults and out‐of‐phase behavior. Such fluctuations have important consequences for tectonic evolution and crustal rheology, and in particular for hazard estimation because they introduce considerable variability and hence uncertainty in earthquake probability calculations.

We present an in situ 36Cl dataset recording the exhumation of 27 active normal fault planes by earthquake slip for the central Apennines, Italy. We do this to constrain the characteristics of earthquake clustering and anticlustering across the entire extending orogen, and in an attempt to constrain the reasons why clustering and anticlustering occurs. We show that duration and magnitude of clustering and anticlustering, and their characteristics, can be explained by a model where the transfer of differential stress between faults and their underlying shear-zones, and between neighbouring fault/shear-zone structures, produces changes in strain-rates on underlying viscous shear zones which drive periods of rapid or reduced slip-rate on their overlying faults. We suggest that stress increase on an underlying shear zone produced by coseismic slip on its overlying fault could be the mechanism that initiates an earthquake cluster. We suggest that stress reductions on shear-zones from coseismic slip located across strike could be the mechanism that initiates an earthquake anticluster. The durations of anticlusters are controlled by the summed stress decreases through time on shear zones, because although these shear zones are slipping relatively slowly, eventually they will load their overlying fault to failure initiating a new cluster, with anticlusters induced across strike. Thus, there is dynamic feedback both up and down dip between faults and their underlying shear zones and crucially across strike between neighbouring fault/shear-zone structures. If the dynamics producing clustering and anticlustering can be constrained, it may be that observations of these phenomena should be included in probabilistic seismic hazard assessments (PSHA) and also interpretations of regional deformation rates and crustal rheologies based on geodetic data. Multi-millennial clustering and anticlustering should become a subject for discussion in these scientific communities.

We present a high-resolution topographic map at 1:5000 scale of the northern slope of Mt Epomeo (Ischia Island, Italy), where the active Casamicciola fault system is exposed, and where coseismic surface ruptures associated with the 21 August 2017 earthquake (Md 4.0) were documented. This fault system is also considered the likely source of the destructive 28 July 1883 earthquake (Io = XI MCS), which caused extensive damage and approximately 2300 fatalities. The fault mapping was conducted by integrating UAV-based LiDAR-derived topography with extensive field investigations, 2-D Deep Electrical Resistivity Tomography, and seismic ambient noise measurements. Despite logistical challenges due to urban density, vegetation cover, and rugged terrain, fault scarps were identified and mapped in detail within the Casamicciola Holocene Graben. The map reveals a nested graben structure defined by two principal faults (Monte Cito and Sentinella) and two subsidiary ones (Purgatorio and Nizzola). Improved understanding of fault geometry is essential for seismic source characterisation and surface faulting hazard assessment in this densely populated, high-tourism region.

Arsenic (As) mobility in aquifer systems is mainly governed by its adsorption and desorption behaviour at the sediment-water interface, directly influencing its environmental availability and risks to water quality. This study explores the adsorption-desorption behaviour of inorganic As species through batch experiments on environmental sediments collected from three representative depths, selected to reflect local contrasting geochemical, mineralogical, and granulometric characteristics of the Como basin aquifer (Northern Italy). This setting was selected as a case study owing to its notable gradient in As concentration in groundwater: the shallow aquifers host concentrations typically below 10 µg/L, while the deep aquifer reaches concentrations of about 250 µg/L. Statistical analyses (ANOVA and simple linear regression) identified Mn- and Al-(hydr)oxide content, grain size, and mineralogy as strong predictors of As(V) retention, whereas As(III) showed no significant correlation with individual sediment properties within the tested conditions. Shallow, Mn- and Al-rich sediments exhibited higher adsorption capacity and corresponded to lower dissolved As in groundwater, while deeper, finer-grained sediments with lower oxide content coincided with elevated groundwater As concentrations. Desorption experiments indicated that As(III) dominated the released fraction, reflecting its greater mobility under variable pH and redox aquifer conditions. These results provide mechanistic insight into sediment-water interactions controlling As distribution in multilayer aquifers, supporting improved risk assessment and management of As in complex groundwater systems.

In this study, we present scaling relationships between fault lengths, fault slip-rates and historical seismicity for an active normal fault system, seismically accommodating crustal extension within the upper plate of the Ionian subduction zone (southern Italy). This crustal extension is confirmed by historical seismicity and instrumental geodesy, with GNSS-derived values of horizonal deformation within a range of 2-3 mm/yr throughout Calabria and the Messina Strait region. We collated data for fault slip-rates, fault lengths and historical earthquakes for a given fault to explore whether fault slip-rates are correlated with fault size and their geometric moment. We present new results showing a robust correlation between fault lengths and fault slip-rates, which supports the idea of a relationship for a given fault between fault slip-rates and the geometric moment. We discuss our results in terms of how these correlations should be used if regional deformation is accommodated by localised strain on faults mostly arranged along strike rather than distributed strain on multiple faults across-strike. For instance, we compare our empirical correlation between fault lengths and fault throw-rates over the Middle-Late Pleistocene in Calabria and the Messina Strait with those from Central and Southern Apennines over the Holocene, characterized by strain distributed on multiple faults across-strike and strain localised on faults mostly arranged along-strike, respectively. Tectonic and seismic hazard implications are discussed for future investigations based on fault slip-rates, fault size and historical seismicity.

Uncertainty concerning the processes responsible for slip-rate fluctuations associated with temporal clustering of surface faulting earthquakes is a fundamental, unresolved issue in tectonics, because strain-rates accommodated by fault/shear-zone structures are the key to understanding the viscosity structure of the crust and seismic hazard. We constrain the timing and amplitude of slip-rate fluctuations that occurred on three active normal faults in central Italy over a time period of 20–30 kyrs, using in situ 36Cl cosmogenic dating of fault planes. We identify five periods of rapid slip on individual faults lasting a few millennia, separated time periods of up to 10 millennia with low or zero slip-rate. The rapid slip pulses migrated across the strike between the faults in two waves from SW to NE. We replicate this migration with a model where rapid slip induces changes in differential stress that drive changes in strain-rate on viscous shear zones that drive slip-rate variability on overlying brittle faults. Earthquakes increase the differential stress and strain-rate on underlying shear zones, which in turn accumulate strain, re-loading stress onto the overlying brittle fault. This positive feedback produces high strain-rate episodes containing several large magnitude surface faulting earthquakes (earthquake clusters), but also reduce the differential stress on the viscous portions of neighbouring fault/shear-zones slowing the occurrence of large-magnitude surface faulting earthquakes (earthquake anticlusters). Shear-zones on faults experiencing anticlusters continue to accumulate viscous strain at a lowered rate, and eventually this loads the overlying brittle fault to failure, initiating a period of rapid slip through the positive feedback process described above, and inducing lowered strain-rates onto neighbouring fault/shear-zones. We show that these patterns of differential stress change can replicate the measured earthquake clustering implied by the 36Cl data. The stress changes are related to the fault geometry in terms of distance and azimuth from the slipping structure, implying that (a) strain-rate and viscosity fluctuations for studies of continental rheology, and (b) slip-rates for seismic hazard purposes are to an extent predictable given knowledge of the fault system geometry.

In the period of October–December 2019, the Cotabato–Davao del Sur region (Philippines) was hit by a seismic sequence comprising four earthquakes with magnitude M W > 6.0 (EQ1-4; max magnitude M W 6.8). The earthquakes triggered widespread environmental effects, including landslides and liquefaction features. We documented such effects by means of field surveys, which we supplemented with landslide mapping from satellite images. Field surveys allowed us to gather information on 43 points after EQ1, 202 points after EQs2–3 and 87 points after EQ4. Additionally, we built a multi-temporal inventory of landslides from remote sensing, comprising 190 slope movements triggered by EQ1, 4737 after EQs2–3, and 5666 at the end of the sequence. We assigned an intensity value to each environmental effect using the environmental seismic intensity (ESI-07) scale. Our preferred estimates of ESI-07 epicentral intensity are VIII for the first earthquake and IX at the end of the sequence, which is in broad agreement with other events of similar magnitude globally. This study, which is the first case of the application of the ESI-07 scale to a seismic sequence in the Philippines, shows that repeated documentation of environmental damage and the evaluation of the progression through time may be useful for providing input data for derivative products, such as susceptibility assessment, evaluation of residual risk or investigation of the role played by ground shaking and by other mechanisms able to trigger environmental effects.

Most seismic risk assessments primarily focus on the impact of ground acceleration on infrastructure. However, in Mexico, along with numerous countries in Latin America and beyond, a significant portion of earthquake risk stems from secondary environmental effects such as tsunamis, landslides, and liquefaction processes. These secondary effects can often prove more lethal than the earthquake shaking itself. We used the Environmental Seismic Intensity scale (ESI-07) to assess the Earthquake Environmental Effects (EEEs) of the M-w 7.7 Michoacan earthquake on September 19th, 2022. This made it possible to comprehensively characterize the EEEs and their potential social consequences in the epicentral area along the Mexican subduction zone. Our study draws on data gathered from extensive field surveys, satellite imagery analysis, social media monitoring, and online resources, totalling data from over 8000 sites with EEEs. Through rigorous analysis, ESI-07 intensity values illuminate the vulnerability of coastal areas and coastal ranges to specific secondary effects triggered by large seismic events. We meticulously map the spatial distribution of these triggered EEEs, presenting A) an ESI-07 isoseismal map, and B) comparative analysis with other significant earthquakes, particularly subduction megathrust interface events. Our findings underscore the critical importance of considering environmental coseismic effects and using the ESI-07 scale in territorial planning and seismic risk preparedness. Furthermore, we highlight the unique characteristics of the Mexican Pacific coastal region and, more broadly, the distinct risks posed by coastal areas along subduction zones in terms of secondary seismic environmental hazards. This work emphasizes the urgent need to acknowledge the vulnerability of these regions to secondary effects and the imperative for resilience-building measures to safeguard human well-being and mitigate economic repercussions in the face of future seismic events.

Over the past two decades, the airborne Light Detection and Ranging (LiDAR) system has become a useful tool for acquiring high-resolution topographic data, especially in active tectonics studies. Analyzing Digital Terrain Models (DTMs) from LiDAR exposes morpho-structural elements, aiding in the understanding of fault zones, among other applications. Despite its effectiveness, challenges persist in regions with rapid deformation, dense vegetation, and human impact. We propose an adapted workflow transitioning from the conventional airborne LiDAR system to the usage of drone-based LiDAR technology for higher-resolution data acquisition. Additionally, drones offer a more cost-effective solution, both in an initial investment and ongoing operational expenses. Our goal is to demonstrate how drone-based LiDAR enhances the identification of active deformation features, particularly for earthquake-induced surface faulting. To evaluate the potential of our technique, we conducted a drone-based LiDAR survey in the Casamicciola Terme area, north of Ischia Island, Italy, known for the occurrence of destructive shallow earthquakes, including the 2017 Md = 4 event. We assessed the quality of our acquired DTM by comparing it with existing elevation datasets for the same area. We discuss the advantages and limitations of each DTM product in relation to our results, particularly when applied to fault mapping. By analyzing derivative DTM products, we identified the fault scarps within the Casamicciola Holocene Graben (CHG) and mapped its structural geometry in detail. The analysis of both linear and areal geomorphic features allowed us to identify the primary factors influencing the current morphological arrangement of the CHG area. Our detailed map depicts a nested graben formed by two main structures (the Maio and Sentinella faults) and minor internal faults (the Purgatorio and Nizzola faults). High-resolution DEMs acquired by drone-based LiDAR facilitated detailed studies of the geomorphology and fault activity. A similar approach can be applied in regions where the evidence of high slip-rate faults is difficult to identify due to vegetation cover and inaccessibility.

Earthquake environmental effects may significantly contribute to the damage caused by seismic events; similar to ground motion, the environmental effects are globally stronger in the vicinity and decrease moving away from the epicenter or seismogenic source. To date, a single intensity prediction equation (IPE) has been proposed in the Italian Apennines for intensity scale dealings with environmental effects: the Environmental Seismic Intensity (ESI-07). Here, we evaluate the sensitivity of the IPE with respect to input data and methodological choices and we propose IPEs with global validity for crustal normal faults. We show the strong influence of input data on the obtained attenuation investigating the 1980 Irpinia–Basilicata (Southern Italy) earthquake. We exploit a dataset of 26 earthquakes to build an IPE considering the epicentral distance. We also propose an IPE considering the distance from the fault rupture, which is derived from a dataset of 10 earthquakes. The proposed equations are valid for normal faults up to 40 km from the epicenter/fault and may flank other models predicting ground motion or damage to the built environment. Our work thus contributes to the use of the ESI-07 scale for hazard purposes.

Surface faulting earthquakes are known to cluster in time from historical and palaeoseismic studies in multiple active tectonic settings, including central Greece, southern California and central Italy. However, the mechanism(s) responsible for clustering, such as fault interaction, strain-storage, and evolving dynamic topography, are poorly quantified and hence not well understood. We combine surface dating of active normal fault scarps in central Italy with stress modelling and quartz flow laws, to produce a quantified replication of observed earthquake clustering.We study six active normal faults (including the Mt Vettore fault which ruptured during the 2016 central Italy earthquake sequence) using 36Cl cosmogenic dating. This reveals periods of high and low slip rate, which we interpret to be earthquake clusters/anti-clusters. Interestingly, these changes in slip rate (or clustering) are out-of-phase between neighbouring faults, i.e. when one fault slows down, nearby faults speed up at the same time. To explore the underlying processes driving this out-of-phase clustering behaviour, we link stress transfer caused by slip over clusters/anti-clusters on coupled fault/shear-zone structures with viscous quartz flow laws derived from laboratory experiments.We show that differential stress fluctuates due to fault/shear-zone interactions, and that the magnitude of these fluctuations are sufficient to induce changes in strain-rate and associated slip-rate on neighbouring faults and shear zones. Our results suggest that fault/shear-zone interactions are a plausible and quantifiable explanation for earthquake clustering, thus opening possibilities for process-led and time-dependent seismic hazard assessments.

On December 26, 2018 (2:19 UTC), during a volcanic eruption on the Mt. Etna eastern flank (Sicily, southern Italy), the largest instrumental earthquake ever recorded in the volcano ruptured the Fiandaca Fault, with epicenter between Fleri and Pennisi villages (hypocenter at ca. 300 m a. s. l., Mw 4.9). This was the mainshock of an earthquake swarm and it was accompanied by widespread surface faulting and extensive damage along a narrow belt near the fault trace. Few hours after the mainshock, an episodic aseismic creep event occurred along the Aci Platani Fault, a SE extension of the Fiandaca Fault, which caused several damages in the Aci Platani village. We surveyed and mapped the coseismic and aseismic ground ruptures, and collected structural data on their geometry, displacement, and fault zone fabric. We compared the mapped surface ruptures with topography, lithology, and morphology of the buried top of the sedimentary basement. We conclude that the geometry of the volcanic pile influenced the surface expression of faulting during the December 26, 2018 event. The top surface of the marly clay basement should be considered as a detachment surface for shallow sliding blocks. The earthquake occurred on top of a depression of the sedimentary basement forcing the sliding eastward, causing at surface the re-arrangement of the fault strand pattern and deformation style, switching from shear faulting to a tensile failure. The Fleri earthquake therefore provides an unprecedented dataset for 1) understanding active faulting in the European largest onshore volcano, 2) modeling its complex dynamics, and 3) contributing to a more refined surface faulting hazard assessment at Mt. Etna. Results from this investigation might be useful for characterizing capable faulting in similar volcano-tectonic settings worldwide.

The macroseismic intensity of the February 6, 2018, Mw 6.4, Hualien earthquake, which caused extensive damage around the Hualien area of eastern Taiwan is reassessed using the Environmental Seismic Intensity (ESI-07) scale. We compiled data on earthquake environmental effects (EEEs) caused by the 2018 Hualien earthquake, which includes surface ruptures, ground cracks, liquefaction, and occasional landslides, and estimated the epicentral intensity (I0) as well as site-specific intensities. We found that the ESI-07 epicentral intensity of the Hualien quake in 2018 is IX. We note that the epicentral area of the 2018 Hualien earthquake was the mesoseismal area of the October 22, 1951, (Mw 6.6) Hualien earthquake, as reported in primary contemporary sources and historical earthquake catalogs. The 1951 Hualien earthquakes also caused prominent surface ruptures, liquefaction, and ground cracks. Consequently, we reassess the macroseismic intensities of this historical seismic event and compare it to the Hualien earthquake in 2018. The comparison suggests similar epicentral intensities for the two earthquakes (IX and X ESI-07). Moreover, we conducted a systematic comparison between intensity obtained using different scales which revealed the differences of two to three degrees between the ESI-07 and traditional intensity scales. This result reconfirms the significance of documentation and recording of earthquake environmental effects to make intensity assessments for modern seismic events consistent with the historical earthquake records. Moreover, a re-evaluation of historical earthquake intensity in eastern Taiwan could be performed in order to update the seismic hazard map. Application of the ESI-07 intensity scale of recent and historical earthquakes will be helpful in post-earthquake recovery efforts for a future earthquake. The prepared ShakeMaps from the ESI-07 values suggests completely different shapes to the previously generated ShakeMaps considering the peak ground acceleration or peak ground velocity. It suggests that the ShakeMaps prepared from the earthquake environmental effects can be complemented with the instrumental based intensity map to have a better seismic hazard prediction and future land use planning for the region.

Abstract. Low-deforming regions are characterized by long earthquake recurrence intervals; thus, it is fundamental to extend back as much as possible the record of past events. Evidence from single sites or proxies may be not compelling, whereas a more substantial picture may be obtained from the integration of paleo- and archaeo-seismic evidence at multiple sites, eventually supplemented with historical chronicles. Here, we document deformations observed in a stratigraphic sequence (i.e., Via Manzoni Site) and in an archaeological site (Roman Baths) in the city of Como (N. Italy). We perform stratigraphic and sedimentological analyses on the sedimentary sequences at via Manzoni and we document Earthquake Archaeological Effects at the Roman Baths by means of Structure from Motion and field surveys. We interpret the observed deformations as due to earthquake ground shaking and provide constraints on the lower threshold for the triggering of such evidence. Radiocarbon datings and chronological constraints from the archaeological site allow to bracket the time of occurrence of the deformations at the VI century AD. We move toward a more regional view to infer possible seismogenic sources by exploiting a dataset of published paleoseismic evidence in Swiss and N. Italy lakes. We perform an inverse grid search to identify magnitude and location of an earthquake that can explain all the positive and negative evidence consistent with the time interval of the event dated at Como. Our results show that a so far undocumented earthquake (minimum Mw 6.32) with epicenter located at the border between Italy and Switzerland may account for all the observed effects. Our study calls for the need to refine the characterization of the local seismic hazard, especially considering that this region seems unprepared to face the effects of a potential earthquake similar to the VI century AD one.

Fault creep along the lower eastern flank of Mt. Etna volcano has been documented since the end of the 19th century and significantly contributes to the surface faulting hazard in the area. On 29 October 2002, during a seismic swarm related to dyke intrusions, two earthquakes caused extensive damage and surface faulting in an area between the Santa Venerina and Santa Tecla villages. On the same day after the two earthquakes, an episodic aseismic creep occurred along the Scalo Pennisi Fault close to the Santa Tecla coastline. On 8 February 2022, during another aseismic creep event along the Scalo Pennisi Fault, we observed the reopening of the preexisting 2002 ground ruptures mostly as pure dilational fractures. We mapped the 2002 and 2022 surface ruptures, and collected data on displacement, length, and pattern of ground breaks. Ground ruptures affected structures located along the activated fault segments, including roads, walls and buildings. The 2002 surface faulting propagation can be ascribed to a sliding of the Mt. Etna eastern flank toward the SE, as also suggested by the related shallow seismicity, and InSAR and geodetic data between 2002 and 2005. For the 2022 event, differential InSAR data, acquired in both descending and ascending views, allowed us to decompose Line of Sight (LOS) displacement into horizontal and vertical components. We detect a similar to 700 m long and similar to 500 m wide deformation zone with a downward and eastward motion (max displacement similar to 1,5 cm) consistent with a normal fault. We inverted the InSAR-detected surface deformation using a uniform-slip fault model and obtained a shallow detachment for the causative fault, located at similar to 300 m depth, within the volcanic pile. This is the first indepth study along the Scalo Pennisi Fault to suggest a shallow faulting that accommodates Mt. Etna E flank gravitational sliding.

This dataset includes documentation of Earthquake Environmental Effects (EEEs) triggered by the Mw 7.7, September 19, 2022, Michoacán (Mexico) earthquake. Data derive from original field surveys, published reports and papers, and from a search for EEEs posted online in social media and other websites. The files include: Shapefile of the sites where each EEE has been documented Shapefile of earthquake-triggered landslides, mapped as polygons from the interpretation of satellite images Shapefile of the area investigated for landslide mapping Pdf document providing a description of each EEE. Shapefiles are in coordinates WGS84 UTM Zone 13 N Credits: EEE documentation from field surveys: Velázquez Bucio M. Magdalena, Pizza Marco, Andres David Nuñez-Meneses EEE documentation from social media: Muccignato Eliana, Velázquez Bucio M. Magdalena, Pizza Marco Landslide mapping from interpretation of satellite images: Ferrario Francesca, Sridharan Aadityan ESI intensity assessment: Velázquez Bucio M. Magdalena, Ferrario Francesca, Lacan Pierre, Muccignato Eliana, Pizza Marco, Porfido Sabina, Michetti Alessandro M. Supervision: Velázquez Bucio M. Magdalena, Lacan Pierre, Michetti Alessandro M.