Random weighting error estimation for the inversion result of finite-fault rupture history

Since the non-unique solution exists in the inversion for finite-fault rupture history, the random weighting method hasbeen used to estimate error of the inversion results in this paper. The resolution distributions of slip amplitude, rake,rupture time and rise time on the finite fault were deduced quantitatively by model calculation. By using the randomweighting method, the inversion results of Taiwan Strait earthquake and Myanmar-China boundal earthquake showthat the parameters related to the rupture centers of two events have the highest resolution, and the solutinn are the mostreliable(otherwise the resolution of the slip amplitudes and rise time on the finite-fault boundary is low.

A new method to calculate near-field radiation excited by heterogeneous fault rupture

This paper derives from the representation theory the formula for calculating the radiation excited by heterogeneous fault rupture based on box-like discretization scheme. Preliminary validation indicates that our algorithm has very high computation precision and efficiency; therefore, it is a very practical tool to investigate strong ground motion problems. Additionally, the equations given in this study can also be used to invert the fault rupture process.

Rapid rupture characterization for the 2023 M_(S)6.2 Jishishan earthquake

On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending In SAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array) or 1.01 km/s(AU array).

Source model of the 11th July 2004 Zhongba earthquake revealed from the joint inversion of InSAR and seismological data

We use interferometric synthetic aperture radar （InSAR） and broadband seismic waveform data to estimate a source model of the 11th July, 2004 M W 6.2 Zhongba earthquake, Tibet of China. This event occurred within the seismically active zone of southwestern Tibetan Plateau where the east-west extension of the upper crust is observed. Because of limitations in one pair of InSAR data available, there are trade-offs among centroid depth, rupture area and amount of slip. Available seismic data tightly constrain the focal mechanism and centroid depth of the earthquake but not the horizontal location. Together, two complementary data sets can be used to identify the actual fault plane, better constrain the slip model and event location. We first use regional seismic waveform to estimate point source mechanism, then InSAR data is used to obtain better location. Finally, a joint inversion of teleseismic P-waves and InSAR data is performed to obtain a distributed model. Our preferred point source mechanism indicates a seismic moment of ～2.2×10 18 N·m （～M W 6.2）, a fault plane solution of 171° （342 ° ）/42 ° （48 ° ）/-83 ° （-97 ° ）, corresponding to strike/dip/rake, and a depth of 11 km. The fault plane with strike of 171 ？ and dip of 42° is identified as the ruptured fault with the aid of InSAR data. The preferred source model features compact area of slips between depth of 5–11 km and 10 km along strike with maximum slip amplitude of about 1.5 m.

Stress triggering of the Lushan M7. 0 earthquake by the Wenchuan Ms8. 0 earthquake

The Wenchuan Ms8.0 earthquake and the Lushan M7.0 earthquake occurred in the north and south segments of the Longmenshan nappe tectonic belt, respectively. Based on the focal mechanism and finite fault model of the Wenchuan Ms8.0 earthquake, we calculated the coulomb failure stress change. The inverted coulomb stress changes based on the Nishimura and Chenji models both show that the Lushan MT. 0 earth- quake occurred in the increased area of coulomb failure stress induced by the Wenchuan Ms8. 0 earthquake. The coulomb failure stress increased by approximately 0. 135 - 0. 152 bar in the source of the Lushan M7.0 earthquake, which is far more than the stress triggering threshold. Therefore, the Lushan M7.0 earthquake was most likely triggered by the coulomb failure stress change.

Simulation of surface displacement and strain field of the 2011 Japan Mw9. 0 earthquake

Based on dislocation theory of 0kada, we adopted a finite-element fault model inverted by Gavin Hayes from seismic data for the 2011 Japan Mw9.0 earthquake, and obtained the corresponding surface displacement and strain fields. The calculated displacement field is consistent with the observed GPS results in the trend of changes. Also the surface displacement and strain fields both show large variations in space.

Source parameter and rupture process of the M_W6.3 early strong aftershock immediately following the 2016 M_W7.8 Kaikoura earthquake(New Zealand)

The 2016 A/w7.8 Kaikoura(New Zealand)earthquake was the most complex event ever instrumentally recorded and geologically investigated,as it ruptured on more than 12 fault segments of various geometries.To study the mainshock rupture characteristics,geodetic methods like InSAR and GPS play an essential role in providing satisfactory spatial resolution.However,early strong aftershocks may cause extra ground deformation which bias the mainshock rupture inversion result.In this paper,we will focus on studying the Mw 6.3 aftershock,which is the only A/6+thrust slip aftershock that occurred only 30 minutes after the Kaikoura mainshock.We will relocate the hypocenter of this event using the hypo 2000 method,make the finite fault model(FFM)inversion for the detailed rupture processes and calculate the synthetic surface displacement to compare with the observed GPS data and figure out its influence on the mainshock study.Although we are not able to resolve the real ruptured fault of this event because of limited observation data,we infer that it is a west-ward dipping event of oblique slip mechanism,consistent with the subfault geometries of the Kaikoura mainshock.According to the inverted FFM,this event can generate 10-20 cm ground surface displacement and affect the ground displacement observation at nearby GPS stations.

The May 23 and 24, 1994 Taiwan earth quakes: comparison of source process between the two events

Inversion for the seismic fault rupture history is an important way to study the nature of the earthquake source. Inthis paper, we have selected two Taiwan earthquakes that occurred closely in time and located in the same region,inversed the distribution of the slip amplitudes, rakes, risetimes and the rupture times on the fault planes by usingGDSN broad-band and long-period records and the adaptive hybrid global search algorithm, and compared the twoevents. The slip rate of every subfault calculated provides information about the distribution of tectonic stress andfault strength. To the former event (Ms=6.0), the maximum slip amplitude 2.4 m and the minimum risetime 1.2 sare both located at the hypocentre. The latter earthquake (Ms=6.6) consisted of two subevents and the second source has 4 s delay. The maximum slip amplitUde 0.9 m located near hypocentre is corresponding to the minimumrisetime l.4 s, and the corresponding maximum slip rate 0.7 m.s~-1 is similar to the peak value of other large sliprate areas. We consider that the latter event has more complicated temporal-spatial distribution than the former.

Hybrid Slip Model for Near-Field Ground Motion Estimation Based on Uncertainty of Source Parameters

The hybrid slip model used to generate a finite fault model for near-field ground motion estimation and seismic hazard assessment was improved to express the uncertainty of the source form of a future earthquake.In this process, source parameters were treated as normal random variables, and the Fortran code of hybrid slip model was modified by adding a random number generator so that the code could generate many finite fault models with different dimensions and slip distributions for a given magnitude.Furth...

Source rupture model of the 2018 M_(W)6.7 Iburi,Hokkaido earthquake from joint inversion of strong motion and InSAR observations

The 2018 M_(W)6.7 Iburi earthquake shocked the eastern Iburi region to the west of the Hidaka Collision Zone in Hokkaido,which is a destructive inland earthquake.We resolved the kinematic rupture process of the event by combining strong motions(SM)and synthetic aperture radar(SAR)images in a joint inversion.The results reveal that the duration of the whole rupture is about 17s,yielding a total seismic moment of 1.4×10^(19)N·m(M_(W)=6.7).The main slip area is located at a depth of approximately 24 km with a peak slip of~0.8m above the hypocenter.The comparison with the regional velocity model shows the earthquake was initiated in the upper mantle,while the majority of slips are located in the lower crust,which is an“aseismic”domain in the typical sandwich model.The location of the major slip area is consistent with a high-conductivity volume.We proposed a mechanism of low frictional property(<0.3)produced by high pore pressure to explain the abnormal high dip angle and centroid depth located in the ductile lower-crust.Aftershocks are distributed in areas where the Coulomb frictional stress increases due to co-seismic displacement with a mechanism conjugating to the mainshock.

Rupture process of the M_s 7.0 Lushan earthquake, 2013

On April 20, 2013 at 8：02 am, a magnitude 7.0 earthquake occurred in Lushan County, Sichuan Province, China, which induces massive landslides, causes great losses to life and property. Based on the locations of aftershocks provided by the China Earthquake Network Center and the characteristic of Longmenshan active faults system, combined with the current preliminary focal mechanism solution, the fault rupture direction is determined. With the finite fault inversion method, we invert the rupture process of the Lusban Ms7.0 earthquake by teleseismic waveforms data. The inversion results indicate that the main shock is dominated by thrust fault component and the rupture initiated at depth of 15 km, and most of slip ruptured around the hypocenter with the peak slip of about 1.5 m. Most of rupture slips released at the first 20 s and the main rupture occurred at the first 10 s after the onsets of the mainshock. Most of seismic energy released near the hypocenter with a length of 28 km, especially on both sides of the hypocenter with the range of 20 km, and the seismic energy released relatively smaller in other areas. There is a large area with weak slip between the main rupture and another two asperities on both sides of the hypocenter; it may imply that the accumulated strain on the rupture fault has not been completely released. Therefore, there is a significant possibility of having strong aftershocks in the areas where energy is not fully released. This is also the main reason why there are a lot of moderate to strong aftershocks in the Lushan aftershock sequence. In addition, there is an earthquake vacant zone with a length of about 50 km between the Wenchuan Mw7.9 earthquake and this event, which is of high earthquake risk and is deserved to be paid close attention to.