Rupture Kinematics of the 11 January 2021 Mw 6.7 Hovsgol, Mongolia, Earthquake and Implications in the Western Baikal Rift Zone

The Mongolia plateau is the farthest intracontinental region of the India–Eurasia collision and is a transition zone between north–south convergence to the south in the Tien Shan and northwest–southeast extension to the north in the Baikal rift. Mongolia has experienced four M 8 earthquakes since 1905, but due to limited observations, the mechanism of these strong earthquakes and regional tectonics are poorly understood. The 11 January 2021 Mw 6.7 Hovsgol, Mongolia, earthquake is the largest event that has occurred in the Hovsgol graben, which is noted for being the northernmost convergence region of the India–Eurasia collision and the youngest extension region of the Baikal rift. In this article, the coseismic displacements are retrieved by space geodesy for the first time in this region, providing good constraints for the deformation pattern. We use a finite-fault inversion of InSAR and teleseismic data, and a backprojection analysis to reveal the rupture kinematics of this event. The geometry of the Hovsgol fault is determined as east-dipping with a dip of 45°. The rupture process is characterized by a northwestward propagation with a moderate average rupture velocity of ∼2.0 km/s and a complex slip pattern composed of two major slip patches with dimensions of 40 km×20 km. The oblique slip, illustrated by predominate extension and significant dextral striking, confirms the right-lateral-striking faulting in the Hovsgol rift, which indicates that the eastwardly north–south convergence across the southwest segment of the Baikal rift has decreased.

Consecutive ruptures on a complex conjugate fault system during the 2018 Gulf of Alaska earthquake

We developed a flexible finite-fault inversion method for teleseismic P waveforms to obtain a detailed rupture process of a complex multiple-fault earthquake. We estimate the distribution of potency-rate density tensors on an assumed model plane to clarify rupture evolution processes, including variations of fault geometry. We applied our method to the 23 January 2018 Gulf of Alaska earthquake by representing slip on a projected horizontal model plane at a depth of 33.6 km to fit the distribution of aftershocks occurring within one week of the mainshock. The obtained source model, which successfully explained the complex teleseismic P waveforms, shows that the 2018 earthquake ruptured a conjugate system of N-S and E-W faults. The spatiotemporal rupture evolution indicates irregular rupture behavior involving a multiple-shock sequence, which is likely associated with discontinuities in the fault geometry that originated from E-W sea-floor fracture zones and N-S plate-bending faults.

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

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.

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

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.