11 research outputs found

    Ground deformation and source geometry of the 30 October 2016 Mw 6.5 Norcia earthquake (Central Italy) investigated through seismological data, DInSAR measurements, and numerical modelling

    Get PDF
    We investigate the Mw 6.5 Norcia (Central Italy) earthquake by exploiting seismological data, DInSAR measurements, and a numerical modelling approach. In particular, we first 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 hypocenters 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). 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 our analysis by applying a 2D numerical modelling approach based on the finite element method, implemented in a structural-mechanic framework, and exploiting the available geological and seismological data, and the ground deformation measurements retrieved from the multi-orbit ALOS-2 DInSAR analysis. In this case, we consider two different scenarios: the first one based on a single SW-dipping fault, the latter on 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 the frictional force acting in the opposite direction, consistently with the double couple fault plane mechanism

    Dual control of fault intersections on stop-start rupture in the 2016 Central Italy seismic sequence

    Get PDF
    Large continental earthquakes necessarily involve failure of multiple faults or segments. But these same critically-stressed systems sometimes fail in drawn-out sequences of smaller earthquakes over days or years instead. These two modes of failure have vastly different implications for seismic hazard and it is not known why fault systems sometimes fail in one mode or the other, or what controls the termination and reinitiation of slip in protracted seismic sequences. A paucity of modern observations of seismic sequences has hampered our understanding to-date, but a series of three Mw>6 earthquakes from August to November 2016 in Central Italy represents a uniquely well-observed example. Here we exploit a wealth of geodetic, seismological and field data to understand the spatio-temporal evolution of the sequence. Our results suggest that pre-existing fault structures controlled the extent and termination of rupture in each event in the sequence, and that fluid diffusion, channelled along these same structures, may have also determined the timing of rupture reinitiation. This dual control of subsurface structure on the stop-start rupture in seismic sequences may be common; future efforts should focus on investigating its prevalence

    Testing Different Tectonic Models for the Source of the M w 6.5, 30 October 2016, Norcia Earthquake (Central Italy): A Youthful Normal Fault, or Negative Inversion of an Old Thrust?

    Get PDF
    We adopted a multidisciplinary approach to investigate the seismotectonic scenario of the 30 October 2016, Mw 6.5, Norcia earthquake, the largest shock of the 2016\u20132017 central Italy earthquake sequence. First, we used seismological and geodetic data to infer the dip of the main slip patch of the seismogenic fault that turned out to be rather low\u2010angle (~37\ub0). To evaluate whether this is an acceptable dip for the main seismogenic source, we modeled earthquake deformation using single\u2010 and multiple\u2010fault models deduced from aftershock pattern analyses. These models show that the coseismic deformation generated by the Norcia earthquake is coherent with slip along a rather shallow\u2010dipping plane. To understand the geological significance of this solution, we reconstructed the subsurface architecture of the epicentral area. As the available data are not robust enough to converge on a single fault model, we built three different models encompassing all major geological evidence and the associated uncertainties, including the tectonic style and the location of major d\ue9collement levels. In all models the structures derived from the contractional phase play a significant role: from controlling segmentation to partially reusing inherited faults, to fully reactivating in extension a regional thrust, geometrically compatible with the source of the Norcia earthquake. Based on our conclusions, some additional seismogenic sources falling in the eastern, external portions of the Apennines may coincide with inherited structures. This may be a common occurrence in this region of the chain, where the inception of extension is as recent as Middle\u2010Upper Pleistocene

    Source Parameters of the 2016–2017 Central Italy Earthquake Sequence from the Sentinel-1, ALOS-2 and GPS Data

    No full text
    In this study, joint inversions of Synthetic Aperture Radar (SAR) and Global Position System (GPS) measurements are used to investigate the source parameters of four Mw > 5 events of the 2016–2017 Central Italy earthquake sequence. The results show that the four events are all associated with a normal fault striking northwest–southeast and dipping southwest. The observations, in all cases, are consistent with slip on a rupture plane, with strike in the range of 157° to 164° and dip in the range of 39° to 44° that penetrates the uppermost crust to a depth of 0 to 8 km. The primary characteristics of these four events are that the 24 August 2016 Mw 6.2 Amatrice earthquake had pronounced heterogeneity of the slip distribution marked by two main slip patches, the 26 October 2016 Mw 6.1 Visso earthquake had a concentrated slip at 3–6 km, and the predominant slip of the 30 October 2016 Mw 6.6 Norcia earthquake occurred on the fault with a peak magnitude of 2.5 m at a depth of 0–6 km, suggesting that the rupture may have reached the surface, and the 18 January 2017 Mw 5.7 Campotosto earthquake had a large area of sliding at depth 3–9 km. The positive static stress changes on the fault planes of the latter three events demonstrate that the 24 August 2016 Amatrice earthquake may have triggered a cascading failure of earthquakes along the complex normal fault system in Central Italy

    Measuring Coseismic Deformation With Spaceborne Synthetic Aperture Radar: A Review

    Get PDF
    In the past 25 years, space-borne Synthetic Aperture Radar imagery has become an increasingly available data source for the study of crustal deformation associated with moderate to large earthquakes (M > 4.0). Coseismic surface deformation can be measured with several well-established techniques, the applicability of which depends on the ground displacement pattern, on several radar parameters, and on the surface properties at the time of the radar acquisitions. The state-of-the-art concerning the measurement techniques is reviewed, and their application to over 100 case-studies since the launch of the Sentinel-1a satellite is discussed, including the performance of the different methods and the data processing aspects, which still constitute topics of ongoing research

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

    Get PDF
    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.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

    Multiple Lines of Evidence for a Potentially Seismogenic Fault Along the Central-Apennine (Italy) Active Extensional Belt–An Unexpected Outcome of the MW6.5 Norcia 2016 Earthquake

    Get PDF
    The Apenninic chain, in central Italy, has been recently struck by the Norcia 2016 seismic sequence. Three mainshocks, in 2016, occurred on August 24 (MW6.0), October 26 (MW 5.9) and October 30 (MW6.5) along well-known late Quaternary active WSW-dipping normal faults. Coseismic fractures and hypocentral seismicity distribution are mostly associated with failure along the Mt Vettore-Mt Bove (VBF) fault. Nevertheless, following the October 26 shock, the aftershock spatial distribution suggests the activation of a source not previously mapped beyond the northern tip of the VBF system. In this area, a remarkable seismicity rate was observed also during 2017 and 2018, the most energetic event being the April 10, 2018 (MW4.6) normal fault earthquake. In this paper, we advance the hypothesis that the Norcia seismic sequence activated a previously unknown seismogenic source. We constrain its geometry and seismogenic behavior by exploiting: 1) morphometric analysis of high-resolution topographic data; 2) field geologic- and morphotectonic evidence within the context of long-term deformation constraints; 3) 3D seismological validation of fault activity, and 4) Coulomb stress transfer modeling. Our results support the existence of distributed and subtle deformation along normal fault segments related to an immature structure, the Pievebovigliana fault (PBF). The fault strikes in NNW-SSE direction, dips to SW and is in right-lateral en echelon setting with the VBF system. Its activation has been highlighted by most of the seismicity observed in the sector. The geometry and location are compatible with volumes of enhanced stress identified by Coulomb stress-transfer computations. Its reconstructed length (at least 13 km) is compatible with the occurrence of MW≥6.0 earthquakes in a sector heretofore characterized by low seismic activity. The evidence for PBF is a new observation associated with the Norcia 2016 seismic sequence and is consistent with the overall tectonic setting of the area. Its existence implies a northward extent of the intra-Apennine extensional domain and should be considered to address seismic hazard assessments in central Italy

    Quantifying co-seismic and post-seismic slip on fault scarps and their erosional modification using high-resolution Pleiades optical satellite data and repeat Terrestrial Laser Scanning: the 2016 Mw 6.6 Norcia earthquake (Central Italy)

    Get PDF
    Fault scarps are a topographical expression of sharp gradients in ground movements in an active tectonic region. However, inferences of slip history and causative earthquake sizes may be biased by co-seismic slip gradients, near-fault deformation, afterslip and erosional processes. To address these biases, I investigate co-seismic and continuing post-seismic deformation of near-fault areas and degradation of fault scarps, using high resolution Pleiades optical satellite images and repeat Terrestrial Laser Scans (TLS). The study area is Monte Vettore in the Apennines, Cental Italy, which has extensive surface ruptures associated with the Mw 6.6 30th October 2016 Norcia earthquake, part of the Central Italy Earthquake Sequence. I combine image correlation techniques with novel median-based filtering to effectively de-noise the Pleiades data, creating Digital Elevation Models (DEMs) from before and after the Norcia earthquake. Those DEMs are then differenced horizontally and vertically. The results identify detail of near-fault co-seismic surface deformation. I jointly invert those data with far-field InSAR (Interferometric Synthetic Aperture Radar) and GNSS (Global Navigation Satellite System) datasets to model co-seismic slip at depth. My model reveals detail of slip transfer from the Monte Vettore fault at shallow depth. This provides insights into the distribution of near-fault co-seismic slip in an area of complex faulting by slip being partitioned onto minor near-surface hanging wall structures, with slip vectors diverging from those at greater depth. The causes of post-seismic alteration or degradation of fault scarps are expected to be tectonic-related after-slip and/or erosion. Combining careful alignment of repeat TLS, use of an ICP (Iterative Closest Point) algorithm, filtering and detrending techniques, I characterise post-seismic deformation at 6 individual sites at ~centimetre scale. This provides insights into how individual factors (e.g. underlying geology, topography, and co-seismic slip gradients and distribution) influence which causes dominate and how degradation develops spatially and temporally. I show that fault scarps are highly variable records of a fault’s slip history. Any assessment of previous slip history using fault scarps as evidence needs to have regard to all those factors.

    Variations in slip-rate and earthquake occurrence across 3D structural complexities on active normal faults

    Get PDF
    This PhD thesis provides a series of studies on the relationship between the non-planar geometry and the seismic behaviour of active normal faults. Herein, several examples show that throw and fault dip increase within along-strike fault bends in order to preserve the horizontal strain-rate within the bend and along the fault. This has been demonstrated for a variety of normal faults (a) located in different geodynamic domains and (b) for measurements of throw taken over different time periods. Furthermore, throw enhancement within fault bends has been observed on (1) immature faults, where fault bends are still propagating up to the surface and are not yet fully established, (2) well-established single fault segments, where fault bends affect one continuous fault segment, and (3) densely-spaced fault systems arranged across strike (with fault spacing < 5 km), where a change in strike across several fault segments creates an overall bend in the system. The results presented in this thesis suggest that the relationship between the non-planar geometry and the distribution of throw is scale-independent, and can act across systems of faults if they are closely spaced across strike. Moreover, 36Cl-cosmogenic dating of tectonically-exposed fault planes on faults spaced > 5 km across-strike shows that these faults are clustered, with a non-systematic alternance of periods of rapid slip accumulation (i.e. earthquake clustering) and periods of quiescence (i.e. earthquake anti-clustering); this suggests that parallel faults interact in terms of sharing the regional strain-rate, with switching activity that affects the slip-rate on a single fault. The results have implications for numerous and diverse aspects of the earthquake geology, such as interpretation of palaeoseismology studies including at trench sites and cosmogenic dating of fault planes, empirical scaling relationships, PSHA, and off-fault deformation
    corecore