
SpecFWAT: Illuminate Refined Subsurface Structures
An easy, fast, and powerful full-waveform adjoint tomography (FWAT) tool for multiple seismic data.
Features
SpecFWAT employ adjoint waveform tomography method and support multiple types of seismic data and is designed to be user-friendly and fast
SPECFEM3D
It is built on top of the open-source software package SPECFEM3D for simulating seismic wave propagation in 3D heterogeneous media.
Flexible data support
Teleseismic, ambient noise, and receiver function data are supported and their joint inversion can also be used.
Multi-parameter inversion
Both isotropic (Vp, Vs, and density) and anisotropic (Gc and Gs) elastic parameters can be inverted for.
Fast and efficient
SpecFWAT is designed to be fast and efficient, allowing users to perform FWI on supercomputer accelerated by GPUs and CPUs.
Easy
Simple format for model files (.h5), data files (.sac) and parameter files (.yml) that are easy to read and write.
Scalable
It is designed to be scalable under Object-Oriented Fortran by modularizing the code and using CMake for building.
Showcases
SpecFWAT has been used to invert for the subsurface structure of the Earth in various regions
Receiver Function Adjoint Tomography
We have developed a novel technique that allows for direct inversion of P-wave receiver functions for a detailed 3-D S-wave velocity model. By applying this technique to a dense seismic array located in the southeastern margin of the Tibetan Plateau, we can resolve detailed crustal and uppermost mantle structures with strong lateral variations. The results indicate that the boundary faults of the geological blocks play a critical role in controlling the lateral variations of the crustal and upper mantle structures. Our findings provide new insights into the dynamics of lithosphere-asthenosphere interactions and continental collision processes in this region.

Joint Full-waveform Adjoint Tomography
The genesis of the Cenozoic intraplate volcanism in Central Mongolia, characterized by sustained and low-volume eruptions remains debated due to the lack of a comprehensive model to interpret the Cenozoic volcanic activities. Here, we introduce a high-resolution 3D velocity model of the Hangay Dome, using a novel joint method which combines receiver function adjoint tomography and ambient noise adjoint tomography. The small-scale low-velocity zones in the crust and uppermost mantle reveal a crustal magma reservoir and partially molten subcontinental lithospheric mantle (SCLM). Melt fraction estimation indicates low-degree partially molten crust and SCLM. Combining previous geophysical and geochemical observations, we suggest that the volcanism in the Hangay Dome is driven by multilevel mechanisms. The remnant Mesozoic volatiles triggered upper mantle upwelling. This upwelling accumulated in the asthenosphere, heating the SCLM, and prompted its low-degree partial melting. The molten SCLM caused local lithospheric thinning and facilitated the magmatic underplating in the lower crust, eventually leading to the formation of the crustal magma reservoir.
