3-D joint geodetic and strong-motion finite fault inversion of the 2008 May 12, Wenchuan, China Earthquake

2020 ◽  
Vol 222 (2) ◽  
pp. 1390-1404
Author(s):  
Leonardo Ramirez-Guzman ◽  
Stephen Hartzell
Keyword(s):  
Time Window ◽  
Joint Inversion ◽  
Strong Motion ◽  
Multiple Time ◽  
Source Inversion ◽  
Geodetic Inversion ◽  
Finite Fault Inversion ◽  
Finite Fault ◽  
China Earthquake ◽  
The Moment

SUMMARY We present a source inversion of the 2008 Wenchuan, China earthquake, using strong-motion waveforms and geodetic offsets together with 3-D synthetic ground motions. We applied the linear multiple time window technique considering geodetic and dynamic Green's functions computed with the finite-element method and the reciprocity and Strain Green's Tensor formalism. All ground motion estimates, valid up to 1 Hz, accounted for 3-D effects, including the topography and the geometry of the Beichuan and Pengguan faults. Our joint inversion has a higher moment (M0) than a purely geodetic inversion and the slip distribution presents differences when compared to 1-D model source inversions. The moment is estimated to be M0 = 1.2 × 1021 N·m, slightly larger than other works. Our results show that considering a complex 3-D structure reduces the size of large areas of 10 m slip or greater by distributing it in wider zones, with reduced slips, in the central portion of the Beichuan and the Pengguan faults. Finally, we compare our source with a relocated aftershock catalogue and conclude that the 4–5 m slip contours approximately bound the absence or presence of aftershocks.

Comparison of local kinematic rupture joint inversion approaches for tsunami early warning: Examples of the 2017-2020 Mw > 6.3 East Aegean earthquakes

2021 ◽  
Author(s):  
Malte Metz ◽  
Marius Isken ◽  
Rongjiang Wang ◽  
Torsten Dahm ◽  
Haluk Özener ◽  
...  
Keyword(s):  
Early Warning ◽  
Joint Inversion ◽  
Moment Tensor ◽  
Strong Motion ◽  
Aegean Sea ◽  
High Rate ◽  
Source Inversion ◽  
Dynamic Rupture ◽  
Centroid Moment Tensor ◽  
Tsunami Early Warning

<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>

scholarly journals RETRACTED ARTICLE: Estimation of the source process and forward simulation of long-period ground motion of the 2018 Hokkaido Eastern Iburi, Japan, earthquake

2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Hisahiko Kubo ◽  
Asako Iwaki ◽  
Wataru Suzuki ◽  
Shin Aoi ◽  
Haruko Sekiguchi
Keyword(s):  
Ground Motion ◽  
Velocity Structure ◽  
Time Window ◽  
Waveform Inversion ◽  
Strong Motion ◽  
Source Model ◽  
Multiple Time ◽  
Velocity Models ◽  
Long Period ◽  
Slip Area

Abstract In this study, we investigate the source rupture process of the 2018 Hokkaido Eastern Iburi earthquake in Japan (MJMA 6.7) and how the ground motion can be reproduced using available source and velocity models. First, we conduct a multiple-time-window kinematic waveform inversion using strong-motion waveforms, which indicates that a large-slip area located at a depth of 25–30 km in the up-dip direction from the hypocenter was caused by a rupture propagating upward 6–12 s after its initiation. Moreover, the high-seismicity area of aftershocks did not overlap with the large-slip area. Subsequently, using the obtained source model and a three-dimensional velocity structure model, we conduct a forward long-period (< 0.5 Hz) ground-motion simulation. The simulation was able to reproduce the overall ground-motion characteristics in the sedimentary layers of the Ishikari Lowland.

scholarly journals Triggered slip on a back reverse fault in the Mw6.8 2013 Lushan, China earthquake revealed by joint inversion of local strong motion accelerograms and geodetic measurements

2016 ◽  
Vol 672-673 ◽  
pp. 24-33 ◽  
Author(s):  
Guohong Zhang ◽  
Eric A. Hetland ◽  
Xinjian Shan ◽  
Martin Vallée ◽  
Yunhua Liu ◽  
...  
Keyword(s):  
Joint Inversion ◽  
Strong Motion ◽  
Reverse Fault ◽  
Geodetic Measurements ◽  
China Earthquake

scholarly journals Development of fully Bayesian multiple-time-window source inversion

2016 ◽  
Vol 204 (3) ◽  
pp. 1601-1619 ◽  
Author(s):  
Hisahiko Kubo ◽  
Kimiyuki Asano ◽  
Tomotaka Iwata ◽  
Shin Aoi
Keyword(s):  
Time Window ◽  
Multiple Time ◽  
Source Inversion ◽  
Fully Bayesian

scholarly journals Estimation of the Wenchuan Earthquake Rupture Sequence Utilizing Teleseismic Records and Coseismic Displacements

2018 ◽  
Vol 2018 ◽  
pp. 1-13 ◽  
Author(s):  
Deyu Yin ◽  
Qifang Liu ◽  
Jingke Wu
Keyword(s):  
Wenchuan Earthquake ◽  
Time Window ◽  
Temporal Distribution ◽  
Field Investigation ◽  
Rupture Process ◽  
Least Square ◽  
Multiple Time ◽  
Finite Fault ◽  
Beichuan Fault ◽  
Dip Angle

For the 12 May 2008 Mw 7.9 Wenchuan earthquake, two imbricate faults, Beichuan fault and Pengguan fault, have ruptured simultaneously. Special attention should be paid to the point of 40 km northeast of the epicenter, in which the Xiaoyudong fault intersects the above two faults, creating a complex fault structure. Surface rupture data from field surveys and previous research of dynamics studies indicate that an important transformation may take place at the intersection. But, few studies about inversion of source rupture process have focused on this issue. We establish a multiple-segment, variable-slip, finite-fault model to reproduce the rupture process and distinguish rupture sequence. Based on the nonnegative least square method and multiple-time-window approach, the spatial and temporal distribution of slip for three rupture sequences are exhibited, using teleseismic records and coseismic displacements. The conformity between synthetic and observed teleseismic records as well as the slip value of the shallowest subfaults and the coseismic displacements is utilized to calibrate the model. The results are as follows: (1) The teleseismic records inversion alone could not distinguish different rupture sequences. However, in order to make the slip of the Hongkou and Yingxiu area coincide with the field investigation, only the Beichuan fault has a bilateral rupture on the point of intersection of Xiaoyudong fault. So the possible rupture sequence is that the earthquake started at the low dip angle part of southern Beichuan fault, and then it propagated to the Pengguan fault, which caused the rupture of Xiaoyudong fault. Then the southern part of Beichuan fault with high dip angle is triggered by the Xiaoyudong fault. (2) The coseismic displacements constraint can control the slip of subfaults near the surface and has little impact on the deeper subfaults. (3) The maximum slip on the fault is located near the Yingxiu and Beichuan area; moreover, the slip is mainly distributed at the shallow region rather than at the deep, which led to serious disasters. Meanwhile, majority of the aftershocks occur in the periphery of large slip.

scholarly journals Construction of fault geometry by finite-fault inversion of teleseismic data

2020 ◽  
Vol 224 (2) ◽  
pp. 1003-1014
Author(s):  
Kousuke Shimizu ◽  
Yuji Yagi ◽  
Ryo Okuwaki ◽  
Yukitoshi Fukahata
Keyword(s):  
Fault Plane ◽  
Seismic Source ◽  
Rupture Process ◽  
Inversion Method ◽  
Source Inversion ◽  
Fault Geometry ◽  
Waveform Data ◽  
Finite Fault Inversion ◽  
Finite Fault ◽  
Double Couple

SUMMARY Conventional seismic source inversion estimates the earthquake rupture process on an assumed fault plane that is determined a priori. It has been a difficult challenge to obtain the fault geometry together with the rupture process by seismic source inversion because of the nonlinearity of the inversion technique. In this study, we propose an inversion method to estimate the fault geometry and the rupture process of an earthquake from teleseismic P waveform data, through an elaboration of our previously published finite-fault inversion analysis (Shimizu et al. 2020). That method differs from conventional methods by representing slip on a fault plane with five basis double-couple components, expressed by potency density tensors, instead of two double-couple components compatible with the fault direction. Because the slip direction obtained from the potency density tensors should be compatible with the fault direction, we can obtain the fault geometry consistent with the rupture process. In practice we rely on an iterative process, first assuming a flat fault plane and then updating the fault geometry by using the information included in the obtained potency density tensors. In constructing a non-planar model-fault surface, we assume for simplicity that the fault direction changes only in either the strike or the dip direction. After checking the validity of the proposed method through synthetic tests, we applied it to the MW 7.7 2013 Balochistan, Pakistan, and MW 7.9 2015 Gorkha, Nepal, earthquakes, which occurred along geometrically complex fault systems. The modelled fault for the Balochistan earthquake is a curved strike-slip fault convex to the south-east, which is consistent with the observed surface ruptures. The modelled fault for the Gorkha earthquake is a reverse fault with a ramp-flat-ramp structure, which is also consistent with the fault geometry derived from geodetic and geological data. These results exhibit that the proposed method works well for constraining fault geometry of an earthquake.

Two main rupture stages during the 2018 magnitude 7.5 Sulawesi earthquake

2020 ◽  
Vol 221 (3) ◽  
pp. 1873-1882
Author(s):  
Qi Li ◽  
Bin Zhao ◽  
Kai Tan ◽  
Wenbin Xu
Keyword(s):  
Body Wave ◽  
Geographical Location ◽  
Strike Slip ◽  
Inversion Method ◽  
The North ◽  
Finite Fault Inversion ◽  
Strike Direction ◽  
Finite Fault ◽  
The Moment ◽  
North And South

SUMMARY On 28 September 2018, a Mw 7.5 strike-slip earthquake occurred in Sulawesi Island, Indonesia, and it unexpectedly triggered a tsunami. To clearly understand the spatiotemporal evolution process of source rupture, we collected the far-field body wave data and utilized the back-projection method together with finite fault inversion method to investigate the rupture kinematics of this earthquake. Results obtained with the two methods have good consistency and complementarity. We hold that the rupture expanded from the epicentre and propagated bilaterally towards the north and south along the strike direction during the first 24 s, and then to the south. Therefore, the whole rupture process consists of two main stages. For the second stage, the fault segment experienced most of the moment release between 0 and 15 km depth, while the fault plane tended to slip at greater depth (down to 20 km) in the first stage. The total length of the rupture was about 200 km and the seismic moment was ∼2.48 × 1020 Nm, which was equivalent to Mw 7.5. The surface rupture was evident and the maximum slip of 6.24 m was observed in the Palu basin, which was close to Palu city. The rupture was dominated by left-lateral strike-slip with both normal and thrust components as well. The normal slip exhibited in the shallow part of the fault on the north side of Palu bay together with the special geographical location of Palu bay likely favored tsunami genesis.

Assessing the role of selected constraints in Bayesian dynamic source inversion: Application to the 2017 Mw 6.3 Lesvos earthquake

2021 ◽  
Author(s):  
Filip Kostka ◽  
Jiří Zahradník ◽  
Efthimios Sokos ◽  
František Gallovič
Keyword(s):  
Posterior Distribution ◽  
Moment Tensor ◽  
Strong Motion ◽  
Monte Carlo Algorithm ◽  
Source Inversion ◽  
Dynamic Rupture ◽  
Centroid Moment Tensor ◽  
Trade Offs ◽  
Finite Fault ◽  
The Mean

Summary A dynamic finite-fault source inversion for stress and frictional parameters of the Mw 6.3 2017 Lesvos earthquake is carried out. The mainshock occurred on June 12, offshore the southeastern coast of the Greek island of Lesvos in the north Aegean Sea. It caused 1 fatality, 15 injuries, and extensive damage to the southern part of the island. Dynamic rupture evolution is modeled on an elliptic patch, using the linear slip-weakening friction law. The inversion is posed as a Bayesian problem and the Parallel Tempering Markov Chain Monte Carlo algorithm is used to obtain posterior probability distributions by updating the prior distribution with progressively more constraints. To calculate the first posterior distribution, only the constraint that the model should expand beyond the nucleation patch is used. Then, we add the constraint that the model should reach a moment magnitude similar to that obtained from our centroid moment tensor inversion. For the final posterior distribution, 15 acceleration records from Greek and Turkish strong motion networks at near regional distances ($\approx 30 - 150$ km) in the frequency range of 0.05–0.15 Hz are used. The three posterior distributions are compared to understand how much each constraint contributes to resolving different quantities. The most probable values and uncertainties of individual parameters are also calculated, along with their mutual trade-offs. The features best determined by seismograms in the final posterior distribution include the position of the nucleation region, the mean direction of rupture (towards WNW), the mean rupture speed (with 68 per cent of the distribution lying between 1.4–2.6 km/s), radiated energy (12–65 TJ), radiation efficiency (0.09–0.38), and the mean stress drop (2.2–6.5 MPa).

Slip history of the 1994 Sanriku-Haruka-Oki, Japan, earthquake deduced from strong-motion data

1997 ◽  
Vol 87 (4) ◽  
pp. 918-931 ◽  
Author(s):  
Wataru Nakayama ◽  
Minoru Takeo
Keyword(s):  
Seismic Waves ◽  
Time Window ◽  
Strong Motion ◽  
Rupture Process ◽  
Multiple Time ◽  
Strong Motion Data ◽  
Motion Data ◽  
The Third ◽  
Slip History ◽  
Double Couple

Abstract We analyzed the seismic waves of the 1994 Sanriku-Haruka-Oki earthquake (Mw = 7.7), which occurred in the aftershock area of the 1968 Tokachi-Oki earthquake (Mw = 8.2). Applying a multiple-time window inversion scheme to near-source strong-motion data, we obtained a detailed spatiotemporal rupture process and compared it with that of the Tokachi-Oki earthquake. The fault geometry is constructed based on the aftershock distribution. The obtained rupture model is consistent with the CMT solution even for a non-double-couple component. The total seismic moment is 4.0 × 1020 N-m. Large slips are concentrated in three asperities: the first asperity centers about 40 km south and 50 km west from the hypocenter with a maximum slip of 4.4 m, the second one centers about 60 km west from the hypocenter with a maximum slip of 2.2 m, and the third one lies about 110 km west from the hypocenter with a maximum slip of 2.6 m. The obtained moment rate and the duration on the first and second asperities are lower and much longer than those on the third asperity, respectively. The first asperity does not overlap with an area of large slip during the Tokachi-Oki earthquake, but the second or third seem to overlap with or be adjacent to the asperity of the Tokachi-Oki earthquake. Our inversion result also shows an abrupt change of the rupture velocity (from 1.8 to 3.0 km/sec) at the central part of the fault plane. A difference of the seismic coupling between the oceanic and the continental lithospheres at the trenchward side and at the landward side of the 143° E meridian seems to affect the rupture process of this earthquake.

Sign in / Sign up

Export Citation Format

Share Document