The pattern evolution and dynamic mechanism of the dynamic changes of regional gravity fields occurring before and after the Wenchuan Ms8.0 earthquake are analyzed, based on five epochs of 1998–2007 mobile gravity data from the middle-south section of the north–south seismic belt, and two epochs of field research data collected after the 2008 Wenchuan earthquake in combination with GPS data, leveling observations, and geotectonic environment data. The regional dynamic gravity changes demonstrate the effects of the eastward flow of solid matter in the Qinghai-Tibetan plateau and the preparation of the 2008 Wenchuan earthquake (2–10yr). The two most meaningful gravity indicators of the Wenchuan earthquake preparation are the positive (increasing) gravity changes occurring over many years in the southwest epicenter and the large-scale gradient zone of gravity variation, with the cumulative difference between the two sides of the gradient zone of gravity exceeding 200 μGal. The positive gravity changes may facilitate a constant energy accumulation and the gradient belt may support seismic shear breakage. Overall, the gravity changes associated with the earthquake preparation indicate a pattern of accelerating increase-decelerating increase-earthquake occurrence. The Songpan-Ganzi block generally displays a negative gravity change, providing evidence for a local upwarping of the deep crust-mantle and an interior expansion of the deep crust attributable to high temperatures. The viewpoint is consistent with the dilatant mechanism for earthquake preparation

By using GPS and gravity data before and after the Wenchuan Ms8.0 earthquake and combining data from geological surveys and geophysical inversion studies, an initial coseismic fault model is constructed. The dip angle changes of the fault slip distribution on the fault plane are inversed, and the inversion results show that the shape of the fault resembles a double-shovel. The Yingxiu–Beichuan Fault is approximately 330 km long, the surface fault dip angle is 65.1°, which gradually reduces with increasing depth to 0° at the detachment layer at a depth of 19.62 km. The Guanxian–Jiangyou Fault is approximately 90 km long, and its dip angle at the surface is 55.3°, which gradually reduces with increasing depth; the fault joins the Yingxiu–Beichuan Fault at 13.75 km. Coseismic slip mainly occurs above a depth of 19 km. There are five concentrated rupture areas, Yingxiu, Wenchuan, Hanwang, Beichuan, and Pingwu, which are consistent with geological survey results and analyses of the aftershock distribution. The rupture mainly has a thrust component with a small dextral strike–slip component. The maximum slip was more than 10 m, which occurred near Beichuan and Hanwang. The seismic moment is 7.84 × 1020 Nm (Mw7.9), which is consistent with the seismological results

The Barkam–Jiuzhaigou–Wuqi gravity profile extends across the Jiuzhaigou Ms7.0 earthquake (in 2017) zone and passes through several historical big earthquakes’ zones. We have obtained Bouguer gravity anomalies along the profile composed of 365 gravity observation stations with Global Positioning System (GPS) coordinates, analyzed the observed data and inverted subsurface density structure. The results show that the Moho depth has a big lateral variation from southwest to northeast, which shallows from 57 km to 43 km with maximum variation up to 14 km within 800 km. The most acute depth change of the Moho is in the boundary region between the Bayan Har block and West Qinling–Qilian block. According to our analysis, it is related to the eastward movement of the Bayan Har block. There are three main pieces of evidence that support it: (1) Density is higher in the east of the Bayan Har block and smaller in the west, which is the same as seismic activity; (2) Two thin low-density layers exist in the upper and middle crust of the Bayan Har block, which may promote inter-layer slip and the Jiuzhaigou Ms7.0 earthquake occurred in the boundary area of the two low-density layers, where the crustal density and Moho surface fluctuate sharply; (3) the GPS velocity field in the southwestern part gravity profile is significantly larger than that of the northeastern part, which is consistent with the density structure. Our studies also suggest that the large undulation of the Moho prevents the movement of the Bayan Har block, and strain is prone to accumulate here. The dynamic background analysis of the crust in this area indicates that the Moho surface uplifts in the West Qinling–Qilian block, which decelerates the eastern migration of material on the Qinghai–Tibet Plateau, and leads to the weak tectonic activity of the north part of the Bayan Har block.

Solid Earth is affected by tidal cycles triggered by the gravity attraction of the celestial bodies. However, about 70% the Earth is covered with seawater which is also affected by the tidal forces. In the coastal areas, the ocean tide loading (OTL) can reach up to 10% of the earth tide, 90% for tilt, and 25% for strain (Farrell, 1972). Since 2007, a high-precision continuous gravity observation network in China has been established with 78 stations. The long-term high-precision tidal data of the network can be used to validate, verifying and even improve the ocean tide model (OTM).In this paper, tidal parameters of each station were extracted using the harmonic analysis method after a careful editing of the data. 8 OTMs were used for calculating the OTL. The results show that the Root-Mean-Square of the tidal residuals (M0) vary between 0.078-1.77 μgal, and the average errors as function of the distance from the sea for near(0-60km), middle(60-1000km) and far(>1000km) stations are 0.76, 0.30 and 0.21 μgal. The total final gravity residuals (Tx) of the 8 major constituents (M2, S2, N2, K2, K1, O1, P1, Q1) for the best OTM has amplitude ranging from 0.14 to 3.45 μgal. The average efficiency for O1 is 77.0%, while 73.1%, 59.6% and 62.6% for K1, M2 and Tx. FES2014b provides the best corrections for O1 at 12 stations, while SCHW provides the best for K1 ,M2and Tx at 12,8and 9 stations. For the 11 costal stations, there is not an obvious best OTM. The models of DTU10, EOT11a and TPXO8 look a litter better than FES2014b, HAMTIDE and SCHW. For the 17 middle distance stations, SCHW is the best OTM obviously. For the 7 far distance stations, FES2014b and SCHW model are the best models. But the correction efficiency is worse than the near and middle stations’.The outcome is mixed: none of the recent OTMs performs the best for all tidal waves at all stations. Surprisingly, the Schwiderski’s model although is 40 years old with a coarse resolution of 1° x 1° is performing relative well with respect to the more recent OTM. Similar results are obtained in Southeast Asia (Francis and van Dam, 2014). It could be due to systematic errors in the surroundings seas affecting all the ocean tides models. It's difficult to detect, but invert the gravity attraction and loading effect to map the ocean tides in the vicinity of China would be one way.