5.5 KiB
TomoATT
TomoATT is a library which implements an eikonal equation solver based Adjoint-state Traveltime Tomography for a very large-scale computation, which implements the methods described in the publications:
-
The TomoATT software package
- Chen, J., Nagaso, M., Xu, M., & Tong, P. (2025). TomoATT: An open-source package for Eikonal equation-based adjoint-state traveltime tomography for seismic velocity and azimuthal anisotropy. Computers & Geosciences, 105995. DOI.
-
Regional tomography in Cartesian coordinates
-
Tong, P. (2021). Adjoint‐state traveltime tomography: Eikonal equation‐based methods and application to the Anza area in southern California. Journal of Geophysical Research: Solid Earth, 126(5), e2021JB021818. DOI.
-
Tong, P. (2021). Adjoint‐state traveltime tomography for azimuthally anisotropic media and insight into the crustal structure of central California near Parkfield. Journal of Geophysical Research: Solid Earth, 126(10), e2021JB022365. DOI.
-
-
Regional tomography in Spherical coordinates
- Chen, J., Chen, G., Nagaso, M., & Tong, P. (2023). Adjoint-state traveltime tomography for azimuthally anisotropic media in spherical coordinates. Geophysical Journal International, 234(1), 712-736. DOI.
-
Teleseismic tomography in Spherical coordinates
- Chen, J., Wu, S., Xu, M., Nagaso, M., Yao, J., Wang, K., ... & Tong, P. (2023). Adjoint‐state teleseismic traveltime tomography: Method and application to Thailand in Indochina Peninsula. Journal of Geophysical Research: Solid Earth, 128(12), e2023JB027348. DOI.
Thanks to the efficiency of an eikonal equation solver, the computation of the travel-time is very fast and requires less amount of computational resources. As an input data for TomoATT is travel times at seismic stations, we can easily prepare a great amount of input data for the computation.
This library is developped to be used for modeling a very-large domain. For this purpose, 3-layer parallelization is applied, which are:
- layer 1: simulutaneous run parallelization (travel times for multiple seismic sources may be calculated simultaneously)
- layer 2: subdomain decomposition (If the number of computational nodes requires too large memory, we can separate the domain into subdomains and run each subdomain in a separate compute node)
- layer 3: sweeping parallelization (in each subdomain, sweeping layers are also parallelized)
The details of the parallelization method applied in this library are described in the paper Miles Detrixhe and Frédéric Gibou (2016).
Regional events (sources within the global domain) and teleseismic events (sources outside the global domain) may be used for inversion.
Quick installation
This library is available in the conda-forge channel. You can install it on personal computer by the following command:
conda install -c conda-forge tomoatt pytomoatt
TomoATT is also capable of running on high-performance computing (HPC) systems and. Detailed installation instructions are described in the installation manual.
to run an example
mpirun -n 4 ./TOMOATT -i ./input_params.yml
Please check the user manual and examples directory for the details.