Integrating seismological data, DInSAR measurements and numerical modelling to analyse seismic events: the Mw 6.5 Norcia earthquake case-study

This Ph.D. thesis is focused on the detailed analysis of seismic sequences and on the application of a multidisciplinary approach based on the integration of several geophysical and geodetic techniques. In the context of a more general study of seismic sequences, we focus this work on the analysis and comparison of five extensional and five compressional seismic sequences to understand the differences between these two tectonic settings in terms of aftershocks duration. In fact, aftershocks number decay through time, depending on several parameters peculiar to each seismogenic regions, including mainshock magnitude, crustal rheology, and stress changes along the fault. However, the exact role of these parameters in controlling the duration of the aftershock sequence is still unknown. Here, by using two methodologies, we show that the tectonic setting primarily controls the duration of aftershocks. On average and for a given mainshock magnitude, (i) aftershock sequences are longer and (ii) the number of earthquakes is greater in extensional than in compressional tectonic settings. We suggest as possible explanation that this difference can be related to the different type of energy dissipated during earthquakes; in detail, (i) a joint effect of gravitational forces and pure elastic stress release governs extensional earthquakes, whereas (ii) pure elastic stress release controls compressional earthquakes. Accordingly, normal faults operate in favour of gravity, preserving inertia for a longer period and seismicity lasts until gravitational equilibrium is reached. Vice versa, thrusts act against gravity, exhaust their inertia faster and the elastic energy dissipation is buffered by the gravitational force. Hence, for seismic sequences of comparable magnitude and rheological parameters, aftershocks last longer in extensional settings because gravity favours the collapse of the hangingwall volumes. The occurrence of the 2016 Central Italy seismic sequence furnishes a test-bed for a detailed analysis of a normal fault earthquake. Therefore, we analyse also the Mw 6.5 Norcia (Central Italy) earthquake to add another extensional seismic sequence to the previously examined case-studies. The results of this analysis show that, with respect to the other considered extensional seismic sequences, also the Mw 6.5 Norcia seismic sequence present the same behaviour about the aftershocks temporal and spatial evolution. Moreover, we decide to take into account the Mw 6.5 Norcia mainshock as case-study for the application of a multidisciplinary approach, in order to understand the kinematics and the role of gravity during nucleation processes of extensional events. In particular, we investigate the Mw 6.5 Norcia earthquake by exploiting seismological data, DInSAR measurements and a numerical modelling approach. In particular, we first take into consideration the relocated hypocentres with 0.1≤Mw≤ 6.5 that occurred between August 24th and November 29th, 2016, recorded by the INGV seismometric network; the projection onto sections and the subsequent analysis of the considered hypocentres allow us to identify the geological structures that were involved during earthquake nucleation process. Then, we retrieve the vertical component (uplift and subsidence) of the displacements affecting the hangingwall and the footwall blocks of the seismogenic faults identified, at depth, through the hypocentres distribution analysis; to do this, we combine the DInSAR measurements obtained from coseismic SAR data pairs collected by the ALOS-2 sensor from ascending and descending orbits. The achieved vertical deformation map displays three main deformation patterns: (i) a major subsidence that reaches the maximum value of about 98 cm near the epicentral zones nearby the town of Norcia; (ii) two smaller uplift lobes that affect both the hangingwall (reaching maximum values of about 14 cm) and the footwall blocks (reaching maximum values of about 10 cm). Also GPS measurements were used to compare the displacements recorded next to the epicentral area. Starting from this evidence, we compute the rock volumes affected by uplift and subsidence phenomena, highlighting that those involved by the retrieved subsidence are characterized by significantly higher deformation values than those affected by uplift (about 14 times). In order to provide a possible interpretation of this volumetric asymmetry, we extend the analysis by running a 2D numerical model based on the finite element method, implemented in a structural-mechanic framework and exploiting the available geological and seismological data. Modelling results are compared with the ground deformation measurements retrieved from the multi-orbit ALOS-2 DInSAR analysis. In the modelling approach, we test the effects of different geometries, by considering two different scenarios: the first is based on including only a single SW-dipping fault, the second includes a main SW-dipping fault and an antithetic zone. In this context, the model characterized by the occurrence of an antithetic zone presents the retrieved best fit coseismic surface deformation pattern. This result allows us to interpret the subsidence and uplift phenomena caused by the Mw 6.5 Norcia earthquake as the result of the gravitational sliding of the hangingwall along the main fault plane and of the frictional force acting in the opposite direction, consistently with the double couple fault plane mechanism.

Questa tesi di dottorato è incentrata sull’analisi dettagliata di sequenze sismiche e sull’applicazione di un approccio multidisciplinare basato sull’integrazione di diverse tecniche geofisiche e geodetiche. Nel contesto di uno studio più generale delle sequenze sismiche, abbiamo concentrato il lavoro sull’analisi ed il confronto di dieci sequenze sismiche (cinque estensionali e cinque compressive), al fine di comprendere le differenze tra questi due ambienti tettonici in termini di durata degli aftershocks. Infatti, il numero di aftershocks decade nel tempo in funzione di vari parametri che risultano essere peculiari di ogni area sismogenetica; tra questi si possono annoverare la magnitudo del mainshock, la reologia crostale e le variazioni dello stress lungo la faglia. Tuttavia, il ruolo esatto svolto da questi parametri nel controllo della durata delle sequenze di aftershocks non è ancora noto. Utilizzando due diverse metodologie, abbiamo evidenziato che l’ambiente tettonico gioca un ruolo primario nell’influenzare la durata degli aftershocks. In media e per una data magnitudo del mainshock, (i) le sequenze di aftershocks sono più lunghe e (ii) il numero di terremoti è maggiore negli ambienti tettonici estensionali rispetto a quelli compressivi. Una possibile spiegazione consiste nel fatto che questa differenza possa essere correlata al diverso tipo di energia dissipata durante i terremoti; in dettaglio, (i) un effetto congiunto di forza gravitazionale e di energia elastica governerebbe i terremoti estensionali, mentre (ii) il rilascio di pura energia elastica controllerebbe i terremoti compressivi. Infatti, le faglie normali operano a favore della gravità, preservando così l'inerzia per un periodo più lungo, e la sismicità dura fino a quando l'equilibrio gravitazionale non viene nuovamente raggiunto dal sistema. Viceversa, i thrusts agiscono contro la gravità, esauriscono la loro inerzia più velocemente e la dissipazione di energia elastica viene controbilanciata dalla forza gravitazionale. Quindi, per sequenze sismiche con magnitudo e parametri reologici paragonabili, gli aftershocks durano più a lungo negli ambienti estensionali poiché la gravità favorisce il collasso dei volumi di hangingwall. Il verificarsi della sequenza sismica del Centro Italia nel 2016 ha fornito un banco di prova per un'analisi dettagliata di un altro terremoto estensionale. Per questo motivo, abbiamo analizzato il terremoto di Norcia (Mw 6.5; Italia Centrale) per aggiungere un’altra sequenza sismica estensionale ai casi di studio precedentemente esaminati. I risultati di questa analisi mostrano che anche la sequenza sismica di Norcia presenta lo stesso comportamento delle altre sequenze estensionali in termini di evoluzione temporale e spaziale degli aftershocks. Inoltre, abbiamo deciso di prendere in considerazione il terremoto di Norcia come caso di studio per l’applicazione di un approccio multidisciplinare, al fine di cercare di comprendere la possibile cinematica e il ruolo della gravità durante i processi di enucleazione degli eventi estensionali. In particolare, abbiamo investigato il terremoto di Norcia, ricorrendo all’utilizzo di dati sismologici, di misure DInSAR e della modellazione numerica. In particolare, abbiamo prima di tutto preso in considerazione gli ipocentri rilocalizzati con 0.1≤Mw≤ 6.5, verificatisi tra il 24 agosto e il 29 novembre 2016 e registrati dalla rete sismometrica INGV; la proiezione su sezioni e la successiva analisi degli ipocentri considerati hanno consentito di comprendere quali strutture geologiche siano state coinvolte durante il processo di enucleazione del terremoto. In seguito, abbiamo analizzato la componente verticale (sollevamento e subsidenza) dei displacements che hanno interessato i blocchi di hangingwall e di footwall della faglia sismogenetica, precedentemente identificata in profondità mediante l’analisi della distribuzione ipocentrale; per fare ciò, abbiamo utilizzato le misure DInSAR ottenute dalla combinazione delle coppie di dati SAR cosismici acquisite dal sensore ALOS-2 lungo orbite ascendenti e discendenti. La mappa di deformazione verticale ottenuta mostra tre aree di deformazione principali: (i) una maggiore subsidenza che raggiunge il valore massimo di circa 98 cm in prossimità delle zone epicentrali vicine alla città di Norcia; (ii) due piccoli lobi di sollevamento che interessano sia il blocco di hangingwall (dove raggiunge valori massimi di circa 14 cm) sia quello di footwall (dove raggiunge valori massimi di circa 10 cm). Partendo da queste evidenze, abbiamo calcolato i volumi di roccia interessati dai fenomeni di sollevamento e subsidenza, evidenziando che quelli coinvolti dal fenomeno di subsidenza sono caratterizzati da valori di deformazione significativamente più alti di quelli affetti da sollevamento (circa 14 volte). Al fine di fornire una possibile interpretazione di questa asimmetria volumetrica, abbiamo esteso l'analisi elaborando un modello numerico 2D basato sul metodo degli elementi finiti, implementandolo in un quadro strutturale-meccanico e sfruttando i dati geologici e sismologici disponibili. I risultati della modellazione sono stati poi confrontati con le misure della deformazione del suolo ottenute dall'analisi DInSAR. Nel corso della realizzazione del modello numerico, abbiamo collaudato gli effetti di geometrie diverse, considerando in particolare due scenari: il primo si basa su una singola faglia immergente a sud-ovest, il secondo su una faglia principale immergente a sud-ovest e una fascia antitetica. In questo contesto, il modello caratterizzato dalla presenza della fascia antitetica fornisce il miglior fit quando confrontato con il pattern cosismico di deformazione superficiale. Questo risultato consente di interpretare i fenomeni di subsidenza e sollevamento causati dal terremoto di Norcia come il risultato di un collasso gravitazionale del blocco di hangigwall lungo la faglia principale e della forza frizionale che agisce in direzione opposta, consistentemente con il meccanismo di doppia coppia lungo il piano di faglia.