scholarly journals Interplay of seismic and a-seismic deformation during the 2020 sequence of Atacama, Chile.

2021 ◽  
Author(s):  
Emilied Klein ◽  
Bertrand Potin ◽  
Francisco Pasten-Araya ◽  
Roxane Tissandier ◽  
Kellen Azua ◽  
...  

An earthquake sequence occurred in the Atacama region of Chile throughout September 2020. The sequence initiated by a mainshock of magnitude Mw6.9, followed 17 hours later by a Mw6.4 aftershock. The sequence lasted several weeks, during which more than a thousand events larger than Ml 1 occurred, including several larger earthquakes of magnitudes between 5.5 and 6.4. Using a dense network that includes broad-band, strong motion and GPS sites, we study in details the seismic sources of the mainshock and its largest aftershock, the afterslip they generate and their aftershock, shedding light of the spatial temporal evolution of seismic and aseismic slip during the sequence. Dynamic inversions show that the two largest earthquakes are located on the subduction interface and have a stress drop and rupture times which are characteristics of subduction earthquakes. The mainshock and the aftershocks, localised in a 3D velocity model, occur in a narrow region of interseismic coupling (ranging 40%-80%), i.e. between two large highly coupled areas, North and South of the sequence, both ruptured by the great Mw~8.5 1922 megathrust earthquake. High rate GPS data (1 Hz) allow to determine instantaneous coseismic displacements and to infer coseismic slip models, not contaminated by early afterslip. We find that the total slip over 24 hours inferred from precise daily solutions is larger than the sum of the two instantaneous coseismic slip models. Differencing the two models indicates that rapid aseismic slip developed up-dip the mainshock rupture area and down-dip of the largest aftershock. During the 17 hours separating the two earthquakes, micro-seismicity migrated from the mainshock rupture area up-dip towards the epicenter of the Mw6.4 aftershocks and continued to propagate upwards at ~0.7 km/day. The bulk of the afterslip is located up-dip the mainshock and down-dip the largest aftershock, and is accompanied with the migration of seismicity, from the mainshock rupture to the aftershock area, suggesting that this aseismic slip triggered the Mw6.4 aftershock. Unusually large post-seismic slip, equivalent to Mw6.8 developed during three weeks to the North, in low coupling areas located both up-dip and downdip the narrow strip of higher coupling, and possibly connecting to the area of the deep Slow Sleep Event detected in the Copiapo area in 2014. The sequence highlights how seismic and aseismic slip interacted and witness short scale lateral variations of friction properties at the megathrust.

Author(s):  
Yong Zhang ◽  
Wanpeng Feng ◽  
Xingxing Li ◽  
Yajing Liu ◽  
Jieyuan Ning ◽  
...  

Abstract The 8 August 2017 Mw 6.5 Jiuzhaigou earthquake occurred in a tectonically fractured region in southwest China. We investigate the multifault coseismic rupture process by jointly analyzing teleseismic, strong-motion, high-rate Global Positioning System, and Interferometric Synthetic Aperture Radar (InSAR) datasets. We clearly identify two right-stepping fault segments and a compressional stepover based on variations in focal mechanisms constrained by coseismic InSAR deformation data. The average geometric parameters of the northwest and southeast segments are strike = 130°/dip = 57° and strike = 151°/dip = 70°, respectively. The rupture model estimated from a joint inversion of the seismic and geodetic datasets indicates that the rupture initiated on the southeastern segment and jumped to the northwestern segment, resulting in distinctive slip patches on the two segments. A 4-km-long coseismic slip gap was identified around the stepover, consistent with the aftershock locations and mechanisms. The right-stepping segmentation and coseismic rupture across the compressional stepover exhibited by the 2017 Jiuzhaigou earthquake are reminiscent of the multifault rupture pattern during the 1976 Songpan earthquake sequence farther south along the Huya fault system in three successive Ms∼7 events. Although the common features of fault geometry and stepover may control the similarity in event locations and focal mechanisms of the 2017 and 1976 sequences, the significantly wider (~15 km) stepover in the 1976 sequence likely prohibited coseismic rupture jumping and hence reduced seismic hazard.


2019 ◽  
Vol 219 (3) ◽  
pp. 1757-1772 ◽  
Author(s):  
Jianfei Zang ◽  
Caijun Xu ◽  
Guanxu Chen ◽  
Qiang Wen ◽  
Shijie Fan

SUMMARY In traditional tight integration of high-rate GNSS and strong motion sensors, an appropriate process variance is crucial for obtaining accurate broad-band coseismic deformations. In this paper, instead of using a subjectively empirical value, we present an approach for determining the process variance adaptively based on the adaptive Kalman filter for real-time use. The performance of the approach was validated by the colocated stations collected during the 2010 Mw 7.2 earthquake in El-Mayor, 2016 Mw 7.8 earthquake in New Zealand and 2016 Mw 6.5 earthquake in central Italy. The results show that this method complements the advantages of GNSS and strong motion accelerometers and can provide more accurate coseismic waveforms especially during the strong shaking period, due to the ability of the method to adjust the process variance in real time according to the actual status of the station. In addition, this method is also free from the influence of the baseline shift. Testing of the new method for the integration of strong motion and multi-GNSS indicates that multi-GNSS has an obvious improvement in the precision while single GPS has a poor observation condition.


Geophysics ◽  
2019 ◽  
Vol 84 (1) ◽  
pp. B41-B57 ◽  
Author(s):  
Himanshu Barthwal ◽  
Mirko van der Baan

Microseismicity is recorded during an underground mine development by a network of seven boreholes. After an initial preprocessing, 488 events are identified with a minimum of 12 P-wave arrival-time picks per event. We have developed a three-step approach for P-wave passive seismic tomography: (1) a probabilistic grid search algorithm for locating the events, (2) joint inversion for a 1D velocity model and event locations using absolute arrival times, and (3) double-difference tomography using reliable differential arrival times obtained from waveform crosscorrelation. The originally diffusive microseismic-event cloud tightens after tomography between depths of 0.45 and 0.5 km toward the center of the tunnel network. The geometry of the event clusters suggests occurrence on a planar geologic fault. The best-fitting plane has a strike of 164.7° north and dip angle of 55.0° toward the west. The study region has known faults striking in the north-northwest–south-southeast direction with a dip angle of 60°, but the relocated event clusters do not fall along any mapped fault. Based on the cluster geometry and the waveform similarity, we hypothesize that the microseismic events occur due to slips along an unmapped fault facilitated by the mining activity. The 3D velocity model we obtained from double-difference tomography indicates lateral velocity contrasts between depths of 0.4 and 0.5 km. We interpret the lateral velocity contrasts in terms of the altered rock types due to ore deposition. The known geotechnical zones in the mine indicate a good correlation with the inverted velocities. Thus, we conclude that passive seismic tomography using microseismic data could provide information beyond the excavation damaged zones and can act as an effective tool to complement geotechnical evaluations.


2017 ◽  
Vol 50 (3) ◽  
pp. 1583
Author(s):  
V. Saltogianni ◽  
M. Gianniou ◽  
T. Taymaz ◽  
S. Yolsal-Çevikbilen ◽  
S. Stiros

A strong earthquake (Mw 6.9) on 24 May 2014 ruptured the North Aegean Trough (NAT) in Greece, west of the North Anatolian Fault Zone (NAFZ). In order to provide unbiased constrains of the rupture process and fault geometry of the earthquake, seismological and geodetic data were analyzed independently. First, based on teleseismic long-period P- and SH- waveforms a point-source solution yielded dominantly right-lateral strike-slip faulting mechanism. Furthermore, finite fault inversion of broad-band data revealed the slip history of the earthquake. Second, GPS slip vectors derived from 11 permanent GPS stations uniformly distributed around the meizoseismal area of the earthquake indicated significant horizontal coseismic slip. Inversion of GPS-derived displacements on the basis of Okada model and using the new TOPological INVersion (TOPINV) algorithm permitted to model a vertical strike slip fault, consistent with that derived from seismological data. Obtained results are consistent with the NAT structure and constrain well the fault geometry and the dynamics of the 2014 earthquake. The latter seems to fill a gap in seismicity along the NAT in the last 50 years, but seems not to have a direct relationship with the sequence of recent faulting farther east, along the NAFZ.


1995 ◽  
Vol 38 (2) ◽  
Author(s):  
M. Di Bona ◽  
M. Cocco ◽  
A. Rovelli ◽  
R. Berardi ◽  
E. Boschi

The strong motion accelerograms recorded during the 1990 Eastern Sicily earthquake have been analyzed to investigate source and attenuation parameters. Peak ground motions (peak acceleration, velocity and displacement) overestimate the values predicted by the empirical scaling law proposed for other Italian earthquakes, suggesting that local site response and propagation path effects play an important role in interpreting the observed time histories. The local magnitude, computed from the strong motion accelerograms by synthesizing the Wood-Anderson response, is ML = 5.9, that is sensibly larger than the local magnitude estimated at regional distances from broad-band seismograms (ML = 5.4). The standard omega-square source spectral model seems to be inadequate to describe the observed spectra over the entire frequency band from 0.2 to 20 Hz. The seismic moment estimated from the strong motion accelerogram recorded at the closest rock site (Sortino) is Mo = 0.8 x 1024 dyne.cm, that is roughly 4.5 times lower than the value estimated at regional distances (Mo = 3.7 x 1024 dyne.cm) from broad-band seismograms. The corner frequency estimated from the accelera- tion spectra i.5 J; = 1.3 Hz, that is close to the inverse of the dUl.ation of displacement pulses at the two closest recording sites. This value of corner tì.equency and the two values of seismic moment yield a Brune stress drop larger than 500 bars. However, a corner frequency value off; = 0.6 Hz and the seismic moment resulting from regional data allows the acceleration spectra to be reproduced on the entire available frequency band yielding to a Brune stress drop of 210 bars. The ambiguity on the corner frequency value associated to this earthquake is due to the limited frequency bandwidth available on the strong motion recordil1gs. Assuming the seismic moment estimated at regional distances from broad-band data, the moment magnitude for this earthquake is 5.7. The higher local magnitude (5.9) compared with the moment magnitude (5.7) is due to the weak regional attenuation. Beside this, site amplifications due to surface geology have produced the highest peak ground motions among those observed at the strong motion sites.


2021 ◽  
Author(s):  
Malte Metz ◽  
Marius Isken ◽  
Rongjiang Wang ◽  
Torsten Dahm ◽  
Haluk Özener ◽  
...  

<p>The fast inversion of reliable centroid moment tensor and kinematic rupture parameters of earthquakes occurring near coastal margins is a key for the assessment of the tsunamigenic potential and early tsunami warning (TEW). In recent years, more and more multi-channel seismic and geodetic online station networks have been built-up to improve the TEW, for instance the GNSS and strong motion networks in Italy, Greece, and Turkey, additionally to the broadband seismological monitoring. Inclusion of such data for the fast kinematic source inversion can improve the resolution and robustness of its’ solutions. However, methods have to be further developed and tested to fully exploit the potential of such rich joint dataset.</p><p>In this frame, we compare and test two in-house developed, kinematic / dynamic rupture inversion methods which are based on completely different approaches. The IDS (Iterative Deconvolution and Stacking, Zhang et al., 2014) combines an iterative seismic network inversion with back projection techniques to retrieve subfault source time functions. The pseudo dynamic rupture model (Dahm et al., in review) links the rupture front propagation estimate based on the Eikonal equation with the dislocation derived from a boundary element method to model dislocation snapshots. We used the latter in both a fast rupture estimate and a fully probabilistic source inversion.</p><p>We use some Mw > 6.3 earthquakes that occurred in the coastal range of the Aegean Sea as an example for comparison: the Mw 6.3 Lesbos earthquake (12 June 2017), the Mw 6.6 Bodrum earthquake (20 July 2017), and the recent Mw 7.0 earthquake which occurred at Samos on 30 October 2020. The latter earthquake and the resulting tsunami caused fatalities and severe damage at the shorelines of Samos and around the city of Izmir, Turkey.<br>The fast estimates are based on only little data and/or prior information obtained from the regional seismicity catalogue and available active fault information. The large number of seismic (broadband, strong motion) and geodetic (high-rate GNSS) stations in local and regional distance from the earthquake with good azimuthal coverage jointly inverted with InSAR data allows for robust inversion results. These, and other solutions, are used as a reference for the comparison of our fast source estimates.<br>Preliminary results of the slip distribution and the source time function are in good agreement with modelling results from other authors.</p><p>We present our insights into the kinematics of the chosen earthquakes investigated by means of joint inversions. Finally, the accuracy of our first fast source estimates, which could be of potential use in tsunami early warning, will be discussed.</p>


Sign in / Sign up

Export Citation Format

Share Document