Moho modeling in spatial domain: A case study under Tibet

The functional relation between the (known) refined gravity field and the (unknown and sought) Moho depth is defined by means of solving the Fredholm integral equation of the first kind. Linearization is applied to define the Moho depth corrections with respect to the mean Moho depth. The system of linearized observation equations is solved to find the Moho depth corrections, and Tikhonov's regularization is applied to stabilize this ill-posed inverse problem. Developed algorithm is applied to model the Moho depth regionally under the Tibetan Plateau and Himalayas, while adopting the uniform and variable Moho density contrast models, and gravimetric results are validated using the CRUST1.0 seismic model. Our results show a relatively good agreement between gravimetric and seismic models, without presence of a significant systematic bias. Our result, however, indicates that for a regional study the variable Moho density contrast might not improve the results especially when density structure of the lower crust and uppermost mantle is not know accurately. In that case, the use of the uniform Moho density contrast for a regional Moho recovery is more appropriate.