Abstract We present the first-generation full-waveform tomographic model (SinoScope 1.0) for the crust-mantle structure beneath China and adjacent regions. The three-component seismograms from 410 earthquakes recorded at 2,427 stations are employed in iterative gradient-based inversions for three successively broadened period bands of 70–120 s, 50–120 s, and 30–120 s. Synthetic seismograms were computed using GPU-accelerated spectral-element simulations of seismic wave propagation in 3-D anelastic models, and Fréchet derivatives were calculated based on an adjoint-state method facilitated by a checkpointing algorithm. The inversion involved 352 iterations, which required 18,600 wavefield simulations. SinoScope 1.0 is described in terms of isotropic P-wave (VP), horizontally and vertically polarized S-wave velocities (VSH and VSV), and mass density (ρ), which are independently constrained with the same data set coupled with a stochastic L-BFGS quasi-Newton optimization scheme. It systematically reduced differences between observed and synthetic full-length seismograms. We performed a detailed resolution analysis by repairing input random parametric perturbations, indicating that resolution lengths can approach the half propagated wavelength within the well-covered areas. SinoScope 1.0 reveals strong lateral heterogeneities in the lithosphere, and features correlate well with geological observations, such as sedimentary basins, Holocene volcanoes, Tibetan Plateau, Philippine Sea Plate, and various tectonic units. The asthenosphere lies below the lithosphere beneath East and Southeast Asia, bounded by subduction trenches and cratonic blocks. Furthermore, we observe an enhanced image of well-known slabs along strongly curved subduction zones, which do not exist in the initial model.

Full-waveform inversion (FWI) is considered as a high-resolution imaging technique to recover the geophysical parameters of the elastic subsurface from the entire content of the seismic signals. However, the subsurface material properties are less well estimated with elastic constraints, especially for the near-surface structure, which usually contains fluid contents. Since Biot theory has provided a framework to describe seismic wave propagations in the poroelastic media, in this work, we propose an algorithm for the 2-D time-domain (TD) poroelastic FWI (PFWI) when the fluid-saturated poroelastic equations are applied to carve the physical mechanism in the shallow subsurface. To detect the contribution of the poroelastic parameters to shallow seismic wavefields, the scattered P-SV\&SH wavefields corresponding to a single model parameter are derived explicitly by Born approximation and shown numerically afterward. The Fréchet kernels are also derived and exhibited in P-SV\&SH schemes to analyse the sensitivities of the objective function to different poroelastic parameters. Furthermore, we verify the accuracy of the derivations through model parameter reconstructions. We perform a series of numerical tests on gradients with respect to different model parameters to further evaluate inter-parameter trade-offs. PFWI holds potential possibilities to directly invert fluid-related physical parameters of the shallow subsurface.