Crustal Density Structure of the Jiuzhaigou Ms7.0 Earthquake Area Revealed by the Barkam–Jiuzhaigou–Wuqi Gravity Profile

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.

Inversion of gravity data with isostatic constraints

We have developed a simple iterative gravity-inversion approach to map the basement and Moho surfaces of a rift basin simultaneously. Gravity anomalies in rift basins commonly consist of interfering broad, positive crustal-thinning anomalies and narrow, negative sedimentary-basin anomalies. In our model, we assumed that the Moho and basement surfaces are in Airy isostatic equilibrium. An initial plane-layered model was iterated to fit the gravity data. We applied the process to a model in which the inverted basement and Moho surfaces matched the model surfaces well and to a gravity profile across the Kosti Basin in Sudan.