SphModel is a subpackage of SphGLLTools, a toolbox for visualization, processing, sharing, and spherical harmonics analysis of tomographic models defined on GLL meshes.
Author: Caio Ciardelli
If you use SphModel, please, cite the following paper:
Ciardelli, C., Bozdağ, E., Peter, D., and van der Lee, S., 2021. SphGLLTools: A toolbox for visualization of large seismic model files based on 3D spectral-element meshes. Computer & Geosciences, v. 159, 105007, doi: https://doi.org/10.1016/j.cageo.2021.105007.
Please, also make sure to cite the authors of GLAD-M15:
Bozdağ, E., Peter, D., Lefebvre, M., Komatitsch, D., Tromp, J., Hill, J., Podhorszki, N., Pugmire, D., 2016. Global adjoint tomography: first-generation model. Geophysical Journal International 207, 1739–1766. URL: https://doi.org/10.1093/gji/ggw356.
SphModel requires no installation. Just take the following steps:
- Make sure you have:
- GCC 5.5.0 or greater;
- Python 2.7.12 or greater;
- GMT 6.0.0 or greater;
-
Download the latest version of GLAD-M15 at SphModel:GLAD-M15.
-
Untar the crustal block model using:
$ for file in crust/*tar.xz; do tar -xvJf "$file" -C crust/; done- Run the Makefile with:
$ makeAll routines have a help menu that shows up whenever you run them with no or with a wrong number of command-line parameters. The
same menu is also at the beginning of each source code. For more details, please, refer to the manual of SphGLLTools.
Development is hosted on GitHub at SphGLLTools.
This parametrization includes Vpv (km/s), Vph (km/s), Vsv (km/s), Vsh (km/s), Eta, and Rho (g/cm^3). The perturbations for the default parameters, as well as some derived parameters such as the isotropic velocities (Vp and Vs), the bulk sound speed, the Vp/Vs ratios, and the transverse isotropy, can also be computed using these routines.
The model is composed of four layers:
Zone 1: Represents the crust using a block model. The horizontal resolution is 0.5 degrees, both in latitude and longitude, and the vertical resolution is 1 km. This zone stretches from 5 km above sea level to 80 km depth.
Zone 2: Represents the upper mantle using a combination of spherical harmonics of various degrees and six cubic B-splines. This zone stretches from 80 km to 410 km depth.
Zone 3: Represents the transition zone using a combination of spherical harmonics of various degrees and five cubic B-splines. This zone stretches from 410 km to 650 km depth.
Zone 4: Represents the lower mantle using a combination of spherical harmonics of various degrees and fourteen cubic B-splines. This zone stretches from 650 km to 2891 km depth.
The plate boundaries are from http://geoscience.wisc.edu/~chuck/MORVEL/citation.html (access on June 20, 2020)
DeMets, C., Gordon, R. G., and Argus, D. F., 2010. Geologically current plate motions, Geophysical Journal International, v. 181, no. 1, p. 1-80, doi: 10.1111/j.1365-246X.2009.04491.x
The hot spot locations are from the table compiled by Don L. Anderson using multiple sources, available at http://www.mantleplumes.org/CompleateHotspot.html (access on November 27, 2020)
If you have any questions, suggestions, and bug reports, you can email caio.ciardelli@gmail.com