The occurrence of strong earthquake is closely related to the distribution of crustal velocity anomalies. Some studies have shown that strong earthquakes occur in the transition zone between high-velocity anomalies and low-velocity anomalies in the middle and upper crust or inside the low-velocity anomaly zone. Thus, high-resolution imaging of the velocity structure in the seismic source area and accurate earthquake location can assist the evaluation of seismogenic settings of strong earthquakes. On May 21, 2021, an M_S6.4 earthquake occurred in Yangbi, Yunnan with casualties and property losses. The epicenter region of the Yangbi earthquake is in the western Yunnan area of the Sichuan-Yunnan block, which is located on the southeastern edge of the Qinghai-Tibet Plateau and characterized with intensive tectonic activity. Previous studies in this area are mostly on regional scales, and lacking on the three-dimensional fine crustal velocity structure in the Yangbi earthquakes area. To investigate the seismogenic environment and source characteristics of the 2021 Yangbi M_S6.4 sequence in Yunnan, we used the P-wave and S-wave arrival data of 12 652 earthquakes recorded by both the Yunnan regional digital network and the mobile observation arrays over a 10-year period(May 1, 2011, to May 31, 2021) and obtained the average V_P/V_S ratio of 1.79 via fitting the P-wave and S-wave arrival-time curves with the Wadati method. The magnitude ranges from M_S0.0 to M_S6.4, and the original focal depth ranges from 0 to 35km. To ensure the reliability of the calculation results, at least 4 stations records are required, and the maximum station azimuth gap allowed is 120°. Furthermore, the event-station distance is restricted to 400km and only earthquakes with travel time residuals<0.5s are retained. Our final velocity model is further refined via gridding(i.e., nodes)with an optimal horizontal grid of 0.25°×0.25° and a range between 0~65km vertically. A checkerboard test is also conduced to validate our inversion results. The test results showed that the recovery degree is high except for the depths of 0 and 65km, which were impacted by the uneven seismic distribution and rays. The high degree of recovery of 5~45km suggests high-resolution and robust imaging at these depths. Finally, the double-difference tomography method(TomoDD)was used to invert the three-dimensional P-wave and S-wave velocity structures in the Yangbi and its surrounding areas(24.5°~26.5°N, 99°~101°E). According to the result of precise location, the M_S6.4 main shock is located at 99.89°E, 25.70°N with a focal depth of 7.9km. The Yangbi M_S6.4 earthquake sequence is mainly distributed along the NW direction. Least-squares fitting prefers a~20km long axis with a strike of 312°, and the hypocenter depths are 5~20km. In general, the studied sequence is shallow and located within the upper crust, consistent with the depth distribution characteristics of historical earthquakes in this area. According to the spatio-temporal evolution characteristics of the aftershock sequence, the aftershocks of the M_S6.4 earthquake mainly spread unilaterally toward SE direction. Thus, we speculate that the overall medium in the NW of the mainshock is rigid and hinders aftershocks evolution. On the north side of the M_S6.4 mainshock epicenter, a group of earthquakes spread along the NNE direction and extended to the Weixi-Qiaohou Fault that hosted the M_S4.1 earthquake on May 27, 2021. Considering the geological and structural background, we believe this earthquake occurred on a parallel but unmapped fault on the SE side of the Weixi-Qiaohou Fault. In contrast, the earthquakes spreading in the NNE direction on the north side of the main shock maybe occurred on an unknown fault in the NNE direction. Therefore, the two faults form a conjugate structure. From the imaging results, the upper crustal velocity structure in the study area is consistent with the geological structure changes and the active faults, where the velocities are low. At 0km depth, the extremely low P-wave and S-wave velocities may reflect impacts from surface sediments. A velocity contrast is observed at a depth of 5km near the mainshock. In addition, a high-velocity anomaly was observed to the southeast side of the mainshock at 10-km depth, with a length of about 0.6°(EW)and a width of about 0.2°(SN). Within the depth range of 10~20km, the distribution of earthquakes near the mainshock shows a clear strip-like distribution, delineating the geometry of the fault. The velocity structure and seismic relocation results at 10-km depth suggest that majority of the events locate around the high-velocity anomaly on the west side of the Weixi-Qiaohou Fault. From the AA' profile, both P- and S-wave velocities suggest high-velocity anomalies in the SE direction of the mainshock. Combining with the distribution characteristics of aftershocks, the non-uniform variations of velocity structure are probably the major factor controlling the distribution of aftershocks, leading to the aftershock distribution extending along the SE direction.

A M_W6.4(M_S7.0)earthquake occurred in Gonghe, Qinghai on 26 April 1990. The Gonghe area is located on the northeastern margin of the Qinghai-Tibet Plateau. The geological tectonic movement in this area is mainly affected by the uplift of the Qinghai-Tibet Plateau. There has been no earthquakes larger than moderate strength in the Gonghe Basin since the historical records, and there are no large-scale active faults on the surface of the epicenter area, so the earthquake has aroused great concern. No major earthquakes have occurred in the Gonghe area since 1995, but the data of small earthquakes is very rich, which ensures the completion of the research. The TomoDD method combines the double-difference relocation method with seismic tomography, and solves two problems at the same time, one is the problem of fine positioning of the earthquake, and the other is the calculation of the 3D velocity structure of the earth’s crust. In this paper, we collected 63872 P and S wave arrival time data in Gonghe and surrounding area recorded by Qinghai, Gansu seismic networks and temporary seismic array from January 2009 to January 2019. The 3D crustal velocity structure and source position parameters of the region are inversed. The relationship between the geological structure setting of the main shock and the velocity structure and seismicity of the region was analyzed. The results show that the crustal velocity structure in the Gonghe area shows lateral inhomogeneity. The Gonghe mainshock is located in the low-velocity anomaly directly below the Gonghe Basin, close to the high-low-velocity anomaly boundary. There is an obvious high-speed anomaly in the southwest of the mainshock, which thrusts from underground to near-surface in the northeast direction. It is estimated that the Wayuxiang-Lagan concealed fault is located at 35.95°N, the dip of the fault is about 45° at the deep part. It is inferred that the occurrence of the Gonghe main shock is caused by the sliding of the Wayuxiangka-Lagan Fault whose strike is NWW and dip is SSW under the action of horizontal tectonic stress. The high-velocity anomaly is about 5~40km deep underground in the northeast direction of the Riyueshan Fault, and a large number of small earthquakes occurred around the high- and low-velocity transition zone. It is presumed that under the action of the near-horizontal NE-directed tectonic stress, the high- and low-velocity zones were further interacted to generate faults and ground folds, and a large number of small earthquakes occurred during the fusion process.