Lithospheric 3-D density structure beneath the Tibetan plateau and adjacent areas derived from joint inversion of satellite gravity and gravity-gradient data

Li Hong-Lei,Fang Jian,Wang Xin-Sheng,Liu Jie,Cui Rong-Hua,Chen Ming

DOI: https://doi.org/10.6038/cjg20170634

2017-01-01

Abstract:A three-dimensional density model of the crust and uppermost mantle beneath the Tibetan plateau and adjacent areas is derived from joint inversion of satellite gravity and gravity gradient data. In this work, we choose five independent gravity gradients (T-xx, T-zz, T-xy, T-xz,T-yz) of GOCE mission L2 gradient products combined with vertical gravity (Gz) of EGM2008 to perform density inversion. The objective function is given based on the Tikhonov regularization theory. The signal of topographic masses, the effect of density interfaces and the influence of density changes below the uppermost mantle are reduced. Using standard deviation ratios of the residual anomalies, we calculate weights for different components. Seismic S-wave velocities provide initial constraints on the inversion based on a relationship between density and S-wave velocity. The damped least square method is used in the inversion. Final density results suggest some new insights to understanding the underlying geodynamic processes as follows: (1) The lowest density layers are distributed mainly on the boundary of the plateau. This phenomenon supports the existing of the high strain accumulation and locking regions caused by the combined effects of the subduction of the Indian plate and the obstruction of the surrounding rigid blocks along the borders of the plateau. (2)At the depth of 80 km, the top of upper mantle inside the Tibetan plateau shows low density. Significantly different density characteristics are present between eastern, central and western parts of the plateau. Low densities are symmetrically distributed on either side of the central line of 95 degrees E. Taking the Bangong-Nujiang suture zone as the center, "low-high-low" density characteristize the Lhasa and Qiangtang blocks from north to south. The density fluctuation is consistent with the undulation of the Moho derived from the satellite GRACE. According to tectonics, such an undulation also supports the bidirectional subduction model that the India plate and Qiangtang block underthrust to the Himalaya and Lhasa blocks, respectively. (3)Density anomalies are significantly lower than velocity anomalies in the Sichuan basin and the Ordos basin, implying they are colder and harder than originally thought. The high densities of the crust and upper mantle in the Qaidam and Tarim basins correspond to the upwelling of heavy materials from the upper mantle.