No abstract

Determining high‐resolution three‐dimensional (3‐D) crustal structures of the Qinling Orogenic Belt (QOB) can reveal significant evidences to enhance our understanding of its tectonic evolution. By jointly inverting group and phase velocity dispersion curves, we obtain a high‐resolution 3‐D shear wave velocity model for the QOB. According to the obviously layered features, which include an E‐W trending high‐velocity structure in the upper crust (0–10 km), a transitional zone in the middle crust (10–30 km), and an approximately N‐S trending low‐velocity zone in the middle‐lower crust (20–40 km), we elucidate the crustal structure as a flyover model to explain how the crust beneath the QOB was deformed under the perpendicular force systems from the NE‐SW and E‐W directions during the Meso‐Cenozoic. First, the apparent NE‐SW trending low‐velocity zone in the middle‐lower crust indicates that the low‐viscosity crustal materials beneath the northeastern Tibetan Plateau did not flow eastward into the eastern Qinling terrane during the Cenozoic. Second, the extension in the northern Qinling Mountains and Weihe Basin was not largely affected by the lateral crustal growth of the NE Tibetan Plateau, which might be the result of deep mantle upwelling and/or residual effects from the subduction of the western Pacific plate. Third, during the Mesozoic, the high‐velocity structures beneath the Hannan‐Micang and Shennongjia‐Huangling domes served as anchors resisting the northward subduction, rotation, and continuous collision of the Yangtze Block. The eastward extruding Hannan‐Micang dome served as an indenter that eventually shaped the present‐day asymmetric style of the Dabashan Orocline during the Cenozoic.

Abstract. The development of an efficient algorithm for teleseismic wave field modeling is valuable for calculating the gradients of the misfit function (termed "misfit gradients") or Fréchet derivatives when the teleseismic waveform is used for adjoint tomography. Here, we introduce an element-byelement parallel spectral-element method (EBE-SEM) for the efficient modeling of teleseismic wave field propagation in a reduced geology model. Under the plane-wave assumption, the frequency-wavenumber (FK) technique is implemented to compute the boundary wave field used to construct the boundary condition of the teleseismic wave incidence. To reduce the memory required for the storage of the boundary wave field for the incidence boundary condition, a strategy is introduced to efficiently store the boundary wave field on the model boundary. The perfectly matched layers absorbing boundary condition (PML ABC) is formulated using the EBE-SEM to absorb the scattered wave field from the model interior. The misfit gradient can easily be constructed in each time step during the calculation of the adjoint wave field. Three synthetic examples demonstrate the validity of the EBE-SEM for use in teleseismic wave field modeling and the misfit gradient calculation.

Full waveform tomography is an effective method to obtain high‐resolution subsurface velocity structures. It is, however, difficult for the conventional L2‐norm‐based inversion to reach the global minimum because of cycle‐skipping problem. In this study, we investigated the feasibility of using the quadratic Wasserstein metric distance (W2 norm) as the misfit function for passive adjoint tomography. We first derived equations of the Fréchet gradient and adjoint source under W2‐norm misfit function. We then conducted numerical experiments to illustrate the effectiveness of the proposed method in avoiding cycle skipping. We finally applied the adjoint tomography to eastern Tibet and obtained a 3‐D velocity model of the lithosphere. The adjoint tomography revealed two low‐velocity channels beneath the NE and SE margins of the Tibetan plateau, which were also observed by previous studies. Our results are consistent with the lower crust flow model that was proposed to explain the deformation occurring at the two margins.

The North China craton (NCC) was dominated by tectonic extension from late Cretaceous to Cenozoic, yet seismic studies on the relationship between crust extension and lithospheric mantle deformation are scarce. Here we present a three dimensional radially anisotropic model of NCC derived from adjoint traveltime tomography to address this issue. We find a prominent low S-wave velocity anomaly at lithospheric mantle depths beneath the Taihang Mountains, which extends eastward with a gradually decreasing amplitude. The horizontally elongated low-velocity anomaly is also featured by a distinctive positive radial anisotropy (VSH > VSV). Combining geodetic and other seismic measurements, we speculate the presence of a horizontal mantle flow beneath central and eastern NCC, which led to the extension of the overlying crust. We suggest that the rollback of Western Pacific slab likely played a pivotal role in generating the horizontal mantle flow at lithospheric depth beneath the central and eastern NCC.