Moment tensor inversion for focal mechanism of the Beibuwan earthquakes

Two earthquakes of Ms=6.0 and Ms=6. 1 consecutively occurred on December 31, 1994 and January 10, 1995 in Beibuwan region, China. By using the generalized reflection-transmission coefficient matrix and the discrete slowness integration method in the calculation of Green's functions, we obtained the focal mechanisms of these earthquakes using long-period waveforms of regional body waves recorded by the China Digital Seismograph Network (CDSN) by means of moment tensor inversion method in frequency domain. The results inverted indicate that the focal mechanisms of these two earthquakes were similar to each other. Their principal compressional stresses are in NW-SE direction and principal tensional stresses are in NE-SW direction. It turns out that the occurrence of the two earthquakes was controlled by the same tectonic environment related to the collision of the Philippine Plate and the Eurasian Plates. On the other hand, the results imply that the stress field in the seismogenic region has a significant change after the Ms=6.0 earthquake. It may be proposed that the occurrence of the Ms=6. 1 earthquake could be related to the stress field adjustment caused by the Ms=6.0 earthquake.

Analysis and Inverse Substructuring Computation on Dynamic Quality of Mechanical Assembly

Mechanical assembly has its own dynamic quality directly affecting the dynamic quality of whole product and should be considered in quality inspection and estimation of mechanical assembly. Based on functional relations between dynamic characteristics involved in mechanical assembly, the effects of assembling process on dynamic characteristics of substructural components of an assembly system are investigated by substructuring analysis. Assembly-coupling dynamic stiffness is clarified as the dominant factor of the effects and can be used as a quantitative measure of assembly dynamic quality. Two computational schemes using frequency response functions（FRFs） to determine the stiffness are provided and discussed by inverse substructuring analysis, including their applicable conditions and implementation procedure in application. Eigenvalue analysis on matrix-ratios of FRFs before and after assembling is employed and well validates the analytical outcomes and the schemes via both a lumped-parameter model and its analogic experimental counterpart. Applying the two schemes to inspect the dynamic quality provides the message of dynamic performance of the assembly system, and therefore improves conventional quality inspection and estimation of mechanical assembly in completeness.

Inversion of focal mechanism and events identification using an improved relative amplitude method

The relative amplitude method （RAM） is more suitable for source inversion of low magnitude earthquakes because it avoids the modeling of short-period waveforms. We introduced an improved relative amplitude method （IRAM） which is more robust in practical cases. The IRAM uses a certain function to quantify the fitness between the observed and the predicted relative amplitudes among direct P wave, surface reflected pP and sP waves for a given focal mechanism. Using the IRAM, we got the fault-plane solutions of two earthquakes of mb4.9 and mb3.8, occurred in Issyk-Kul lake, Kyrgyzstan. For the larger event, its fault-plane solutions are consistent with the Harvard＇s CMT solutions. As to the smaller one, the strikes of the solution are consistent with those of the main faults near the epicenter. The synthetic long period waveforms and the predicted P wave first motions of the solutions are consistent with observations at some of regional stations. Finally, we demonstrated that fault-solutions cannot interpret the characteristics of teleseismic P waveforms of the underground nuclear explosion detonated in Democratic People＇s Republic of Korea （DPRK） on October 9, 2006.

Location and moment tensor inversion of small earthquakes using 3D Green＇s functions in models with rugged topography： application to the Longmenshan fault zone

With dense seismic arrays and advanced imaging methods, regional three-dimensional （3D） Earth models have become more accurate. It is now increasingly feasible and advantageous to use a 3D Earth model to better locate earthquakes and invert their source mechanisms by fitting synthetics to observed waveforms. In this study, we develop an approach to determine both the earthquake location and source mechanism from waveform information. The observed waveforms are filtered in different frequency bands and separated into windows for the individual phases. Instead of picking the arrival times, the traveltime differences are measured by cross-correlation between synthetic waveforms based on the 3D Earth model and observed waveforms. The earthquake location is determined by minimizing the cross-correlation traveltime differences. We then fix the horizontal location of the earthquake and perform a grid search in depth to determine the source mechanism at each point by fitting the synthetic and observed waveforms. This new method is verified by a synthetic test with noise added to the synthetic waveforms and a realistic station distribution. We apply this method to a series of Mw3.4-5.6 earthquakes in the Longmenshan fault （LMSF） zone, a region with rugged topography between the eastern margin of the Tibetan plateau and the western part of the Sichuan basin. The results show that our solutions result in improved waveform fits compared to the source parameters from the catalogs we used and the location can be better constrained than the amplitude-only approach. Furthermore, the source solutions with realistic topography provide a better fit to the observed waveforms than those without the topography, indicating the need to take the topography into account in regions with rugged topography.

Source model of 2016 Mw6.6 Aketao earthquake,Xinjiang derived from Sentinel-1 InSAR observation

On 25 November 2016（14:24:30 UTC）, an Mw6.6 earthquake occurred in Aketao county of Xinjiang. We derived the coseismic deformation field of the earthquake from ESA’s Sentinel-1 B images and inverted the fault geometry and slip distribution by using a homogenous half-space elastic model and the Steepest Decent Method（SDM） program. The slip model shows that the rupture plane is about 68 km in length and 40 km in width and has an average strike of 105.2°SE and an average dip of 76.0°S. The average amount of slip is-0.16 m and the maximum amount of slip is 0.59 m. The estimated seismic moment is 1.31 × 1019 N·m and the corresponding moment magnitude is Mw6.68. The dislocation is mainly distributed in the depth of 5-22 km, where the rupture center is located at 39.170 E, 74.38°N, and a depth of-8.8 km, so this earthquake is a shallow earthquake. In terms of the tectonic settings, the earthquake was triggered by the release of residual stresses on the seismogenic Muji fault. The ruptures were stimulated simultaneously by the main shock, which made the slip distribution pattern of a strike-slip fault with a double fracture. Due to the tectonic setting and the seismic gap between the ruptures, the Muji fault would have a risk of rupture to some extent and trigger an earthquake in future.