Post-seismic deformation related to the 2016 central Italy seismic sequence from GPS displacement time-series

Author(s):  
Eugenio Mandler ◽  
Maria Elina Belardinelli ◽  
Enrico Serpelloni ◽  
Letizia Anderlini ◽  
Adriano Gualandi ◽  
...  

<p>The 2016-2017 Central Italy earthquake sequence was characterized by three main events striking the central Apennines between August 2016 and October 2016 with a Mw ∈ [5.9 to 6.5], plus four earthquakes occurring in January 2017 with a Mw ∈ [5.0; 5.5]. Here we study 85 Global Positioning System (GPS) stations active during the post-seismic phase in a region within a radius of 100 km around the epicentral area, including near and far-field domains. We separate the post-seismic deformation from other, mainly seasonal, deformation signals present in ground displacement time-series via a variational Bayesian Independent Component Analysis (vbICA) technique. Excluding the postseismic transient signal, we found that all the other components are due to hydrological processes, and found no evidence of pre-seismic deformation signals with a spatial and temporal pattern that can be ascribed to a precursory deformation. We study the role played by afterslip on the main structures activated during the co-seismic phase, and we infer the activation during the post-seismic phase of the Paganica fault, which is located further south of the 2016-2017 epicenters and did not rupture during the co-seismic phase. We investigate an aseismic activation of the ∼ 2 − 3 km thick subhorizontal layer of seismicity, which bounds at depth the SW-dipping normal faults where the mainshocks nucleated, and which has been interpreted as a shear zone. Moreover we consider the possibility that the shear zone marks the brittle-ductile transition including the viscoelastic relaxation of the lower crust and upper mantle as a driving mechanism of the post-seismic displacement. However, neither afterslip nor viscoelasticity can fully explain the observations alone: the former is capable of satisfactorily explaining only the data in the epicentral area but it generally underestimates the displacement in the far-field domain; the latter cannot simultaneously explain the displacement observed in the near-field and far-field domains. Hence we infer a mixed contribution of these two mechanisms. </p>

2021 ◽  
Vol 11 (11) ◽  
pp. 4817
Author(s):  
Filippos Vallianatos ◽  
Vassilis Sakkas

In the present work, a multiscale post-seismic relaxation mechanism, based on the existence of a distribution in relaxation time, is presented. Assuming an Arrhenius dependence of the relaxation time with uniform distributed activation energy in a mesoscopic scale, a generic logarithmic-type relaxation in a macroscopic scale results. The model was applied in the case of the strong 2015 Lefkas Mw6.5 (W. Greece) earthquake, where continuous GNSS (cGNSS) time series were recorded in a station located in the near vicinity of the epicentral area. The application of the present approach to the Lefkas event fits the observed displacements implied by a distribution of relaxation times in the range τmin ≈3.5 days to τmax ≈350 days.


2021 ◽  
Author(s):  
Eugenio Mandler ◽  
Letizia ANDERLINI ◽  
Adriano Gualandi ◽  
Francesco Pintori ◽  
Enrico Serpelloni ◽  
...  

Author(s):  
D Spallarossa ◽  
M Cattaneo ◽  
D Scafidi ◽  
M Michele ◽  
L Chiaraluce ◽  
...  

Summary The 2016–17 central Italy earthquake sequence began with the first mainshock near the town of Amatrice on August 24 (MW 6.0), and was followed by two subsequent large events near Visso on October 26 (MW 5.9) and Norcia on October 30 (MW 6.5), plus a cluster of 4 events with MW > 5.0 within few hours on January 18, 2017. The affected area had been monitored before the sequence started by the permanent Italian National Seismic Network (RSNC), and was enhanced during the sequence by temporary stations deployed by the National Institute of Geophysics and Volcanology and the British Geological Survey. By the middle of September, there was a dense network of 155 stations, with a mean separation in the epicentral area of 6–10 km, comparable to the most likely earthquake depth range in the region. This network configuration was kept stable for an entire year, producing 2.5 TB of continuous waveform recordings. Here we describe how this data was used to develop a large and comprehensive earthquake catalogue using the Complete Automatic Seismic Processor (CASP) procedure. This procedure detected more than 450,000 events in the year following the first mainshock, and determined their phase arrival times through an advanced picker engine (RSNI-Picker2), producing a set of about 7 million P- and 10 million S-wave arrival times. These were then used to locate the events using a non-linear location (NLL) algorithm, a 1D velocity model calibrated for the area, and station corrections and then to compute their local magnitudes (ML). The procedure was validated by comparison of the derived data for phase picks and earthquake parameters with a handpicked reference catalogue (hereinafter referred to as ‘RefCat’). The automated procedure takes less than 12 hours on an Intel Core-i7 workstation to analyse the primary waveform data and to detect and locate 3000 events on the most seismically active day of the sequence. This proves the concept that the CASP algorithm can provide effectively real-time data for input into daily operational earthquake forecasts, The results show that there have been significant improvements compared to RefCat obtained in the same period using manual phase picks. The number of detected and located events is higher (from 84,401 to 450,000), the magnitude of completeness is lower (from ML 1.4 to 0.6), and also the number of phase picks is greater with an average number of 72 picked arrival for a ML = 1.4 compared with 30 phases for RefCat using manual phase picking. These propagate into formal uncertainties of ± 0.9km in epicentral location and ± 1.5km in depth for the enhanced catalogue for the vast majority of the events. Together, these provide a significant improvement in the resolution of fine structures such as local planar structures and clusters, in particular the identification of shallow events occurring in parts of the crust previously thought to be inactive. The lower completeness magnitude provides a rich data set for development and testing of analysis techniques of seismic sequences evolution, including real-time, operational monitoring of b-value, time-dependent hazard evaluation and aftershock forecasting.


2014 ◽  
Vol 409 (1) ◽  
pp. 149-166 ◽  
Author(s):  
G. Lena ◽  
M. R. Barchi ◽  
W. Alvarez ◽  
F. Felici ◽  
G. Minelli

2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Sijia Wang ◽  
Yongzhi Zhang ◽  
Yipeng Wang ◽  
Jiashuang Jiao ◽  
Zongtong Ji ◽  
...  

Abstract On 3 July 2015, the Mw 6.5 Pishan earthquake occurred at the junction of the southwestern margin of the Tarim Basin and the northwestern margin of the Tibetan Plateau. To understand the seismogenic mechanism and the post-seismic deformation behavior, we investigated the characteristics of the post-seismic deformation fields in the seismic area, using 9 Sentinel-1A TOPS synthetic aperture radar (SAR) images acquired from 18 July 2015 to 22 September 2016 with the Small Baseline Subset Interferometric SAR (SBAS-InSAR) technique. Postseismic LOS deformation displayed logarithmic behavior, and the temporal evolution of the post-seismic deformation is consistent with the aftershock sequence. The main driving mechanism of near-field post-seismic displacement was most likely to be afterslip on the fault and the entire creep process consists of three creeping stages. Afterward, we used the steepest descent method to invert the afterslip evolution process and analyzed the relationship between post-seismic afterslip and co-seismic slip. The results witness that 447 days after the mainshock (22 September 2016), the afterslip was concentrated within one principal slip center. It was located 5–25 km along the fault strike, 0–10 km along with the fault dip, with a cumulative peak slip of 0.18 m. The 447 days afterslip seismic moment was approximately 2.65 × 1017 N m, accounting for approximately 4.1% of the co-seismic geodetic moment. The deep afterslip revealed that a creeping process from steady-state “secondary” creeping to accelerating “tertiary” creep in the deep of fault. The future seismic hazard deserves further attention and research.


2003 ◽  
Vol 40 (9) ◽  
pp. 1203-1218 ◽  
Author(s):  
David W Eaton ◽  
Jacqueline Hope

The Great Slave Lake shear zone (GSLsz) exposes lower crustal rocks analogous to deep-seated segments of modern strike-slip fault zones, such as the San Andreas fault. Extending for 1300 km beneath the Western Canada Sedimentary Basin to the southern margin of the Slave Province, the GSLsz produces one of the most prominent linear magnetic anomalies in Canada. From May to October 1999, 13 three-component portable broadband seismograph stations were deployed in a 150-km profile across a buried segment of the shear zone to investigate its lithospheric structure. Splitting analysis of core-refracted teleseismic shear waves reveals an average fast-polarization direction (N49°E ± 19°) that is approximately parallel to the shear zone. Individual stations near the axis of the shear zone show more northerly splitting directions, which we attribute to interference between regional anisotropy in the upper mantle (fast axis ~N60°E) and crustal anisotropy within the shear zone (fast axis ~N30°E). At the location of our profile, the shear zone is characterized by a 10-mGal axial gravity high with a wavelength of 30 km, superimposed on a longer wavelength 12-mGal low. This gravity signature is consistent with the basic features of the crustal model derived from receiver-function analysis: a Moho that dips inward toward the shear-zone axis and a mid-crustal zone with high S-wave velocity (ΔVs = 0.6 ± 0.2 km/s). The axial gravity high may be related to uplift of deeper crustal material within the shear zone, or protolith-dependent compositional differences between the shear zone and surrounding wall rocks.


2017 ◽  
Vol 17 (11) ◽  
pp. 1885-1892 ◽  
Author(s):  
Giorgio De Guidi ◽  
Alessia Vecchio ◽  
Fabio Brighenti ◽  
Riccardo Caputo ◽  
Francesco Carnemolla ◽  
...  

Abstract. On 24 August 2016 a strong earthquake (Mw = 6.0) affected central Italy and an intense seismic sequence started. Field observations, DInSAR (Differential INterferometry Synthetic-Aperture Radar) analyses and preliminary focal mechanisms, as well as the distribution of aftershocks, suggested the reactivation of the northern sector of the Laga fault, the southern part of which was already rebooted during the 2009 L'Aquila sequence, and of the southern segment of the Mt Vettore fault system (MVFS). Based on this preliminary information and following the stress-triggering concept (Stein, 1999; Steacy et al., 2005), we tentatively identified a potential fault zone that is very vulnerable to future seismic events just north of the earlier epicentral area. Accordingly, we planned a local geodetic network consisting of five new GNSS (Global Navigation Satellite System) stations located a few kilometres away from both sides of the MVFS. This network was devoted to working out, at least partially but in some detail, the possible northward propagation of the crustal network ruptures. The building of the stations and a first set of measurements were carried out during a first campaign (30 September and 2 October 2016). On 26 October 2016, immediately north of the epicentral area of the 24 August event, another earthquake (Mw = 5.9) occurred, followed 4 days later (30 October) by the main shock (Mw = 6.5) of the whole 2016 summer–autumn seismic sequence. Our local geodetic network was fully affected by the new events and therefore we performed a second campaign soon after (11–13 November 2016). In this brief note, we provide the results of our geodetic measurements that registered the co-seismic and immediately post-seismic deformation of the two major October shocks, documenting in some detail the surface deformation close to the fault trace. We also compare our results with the available surface deformation field of the broader area, obtained on the basis of the DInSAR technique, and show an overall good fit.


2020 ◽  
Vol 178 (1) ◽  
pp. jgs2020-109
Author(s):  
Paulo Castellan ◽  
Gustavo Viegas ◽  
Frederico M. Faleiros

Fabrics of the East Pernambuco shear zone (EPSZ) were studied via microstructural analysis, mineral chemistry and isochemical phase diagram modelling to constrain the pressure and temperature conditions of deformation during shearing. Granitic mylonites show fractured feldspar porphyroclasts embedded in a fine-grained, recrystallized quartzo-feldspathic matrix. These mylonites grade laterally into banded ultramylonites characterized by stretched feldspar clasts alternated with recrystallized quartz bands. Fractures in these ultramylonites are filled by phyllosilicates. The mineral chemistry of the feldspars points to systematic changes between porphyroclasts, grains within fractures and fine-grained mixtures. Quartz crystallographic fabrics in the mylonites suggest activation of prism slip, while the ultramylonites show the activation of both rhomb and basal slip systems. Thermodynamic modelling suggests that the mylonites were formed at 4.75 ± 0.25 kbar and 526 ± 9°C, while the ultramylonites yield conditions of 5.9 ± 1 kbar and 437 ± 17°C. These observations suggest that the EPSZ records a heterogeneous path of strain accommodation, marked by decreasing temperature from its western sector to its eastern termination. The differences in metamorphic conditions are consistent with a transitional, brittle–ductile strain regime. Such characteristics indicate that the EPSZ is a Neoproterozoic shear belt nucleated and heterogeneously exhumed at the brittle–ductile transition, possibly in an intracontinental setting.Supplementary Material: EPMA analysis of feldspars in Caruaru and Gravatá domains and T-X(O2) pseudosections are available at https://doi.org/10.6084/m9.figshare.c.5125957


Sign in / Sign up

Export Citation Format

Share Document