Estimating gravity changes caused by crustal strain: application to the Tibetan Plateau

Temporal gravimetry is an efficient tool for monitoring mass transfers, but distinguishing the contribution of each process to the measured signals is challenging. Few effective methods have been developed to estimate the changes in gravity caused by crustal strain for large-scale geophysical problems. To fill this research gap, we proposed a formula that describes a negative linear correlation between changes in gravity and crustal dilatational strain. Surface observations of gravity changes and dilatational strains were simulated using PSGRN/PSCMP, which is a numerical code used to calculate the surface response to fault dislocations, and the accuracy of the formula was quantitatively verified. Four parameters are required for this formula: the crustal dilatational strain, the crustal density, the Moho depth, and a coefficient that characterizes the degree of crust-mantle coupling. To illustrate the application of this new method to a natural case study, including specifying the values of the necessary parameters, the crustal strain-caused gravity changes (CSGCs) were calculated at 1 degrees x 1 degrees grid nodes over the Tibetan Plateau (TP). The CSGC model shows that most of the crust of the TP is undergoing extension, which generates negative gravity signals. The magnitude of the Tibetan CSGC model is approximately 0.2 mu Gal yr(-1), which is similar to the results obtained from numerical modelling of the crustal tectonics of the Taiwanese Orogen. To evaluate the reliability of the Tibetan CSGC model, the uncertainties in the crustal dilatational strain, crustal density, Moho depth, and crust-mantle coupling factor were evaluated and then used to estimate the CSGC uncertainty by applying the error propagation law. The CSGC model was used to analyse the mass transfers of the TP. The results suggest that a significant mass accumulation process may be occurring beneath the crust of the northern TP.