Crustal and uppermost mantle S velocity structure under Hi-CLIMB seismic array in central Tibetan Plateau from joint inversion of surface wave dispersion and receiver function data

Trade-off between the depth of an interface and the average velocity above the interface is well known in receiver function inversion. To resolve model ambiguity, surface wave dispersion data, which are sensitive to vertical shear velocity average, are commonly combined with P receiver functions in the inversion. In this paper, we design a joint inversion scheme using the neighborhood searching algorithm to solve shear velocity structure of the crust and uppermost mantle with receiver function and surface wave dispersion data. Two sets of stretched splines are implemented in the inversion to obtain smooth velocity profiles, separated by a specifically defined Moho discontinuity. We apply this method to the dense Hi-CLIMB seismic array in the central Tibetan Plateau (TP). The Indian crust is observed to subduct and underplate the Tibetan crust between the Yarlung-Zangbo suture and the Bangong-Nujiang suture, at the latitude of about 31.5°N. The crust thickens from about 50. km below the Indian foreland to over 75. km south of Bangong-Nujiang suture and turns shallower at around 65. km under the Qiangtang block. We identify two zones with complicated Moho structure beneath the northern Lhasa block and southern Qiangtang block, indicating highly deformed Moho morphology, and beneath the Himalaya block associated with the underthrusting of the Indian Plate, respectively. A low velocity zone is present in the mid-crust in most parts of the profile under the TP with up to more than 10% velocity reduction. The low velocity structure is not continuous, forming discontinuous patches rather than channels, calling for a revision of the channel flow model.