Dataset and Software for An anisotropic shear velocity model of the Earth's mantle using normal modes, body waves, surface waves and long-period waveforms

What is the nature of flow in the mantle?
How fast do waves travel anywhere on Earth?
Where can radial anisotropy be robustly detected?
Can we reconcile a broad spectrum of seismic data? What are the benefits?

We use normal-mode splitting functions in addition to surface-wave phase anomalies, body-wave travel times and long-period waveforms to construct a three-dimensional model of anisotropic shear-wave velocity in the Earth's mantle. This is the first tomographic study to exploit the sensitivity of mode-splitting data to constrain radial anisotropy in the Earth's mantle jointly with several other types of data. Our modeling approach inverts for mantle velocity and anisotropy as well as transition-zone discontinuity topographies, and incorporates new crustal corrections for the splitting functions that are consistent with the nonlinear corrections we employ for the waveforms. Our preferred anisotropic model, S362ANI+M, is an update to the earlier model S362ANI, which did not include normal-mode splitting functions in its derivation.

Feedback/Questions? Please contact Raj Moulik (rajmoulik.com) at moulik@caa.columbia.edu 

Reference:

Please cite the following work if you use this data or software.

• Moulik, P. & Ekström, G., 2014. An anisotropic shear velocity model of the Earth's mantle using normal modes, body waves, surface waves and long-period waveforms, Geophys. J. Int., 199(3), 1713-1738, doi: 10.1093/gji/ggu356. pdf

You can also cite the dataset and software from this Zenodo page (Optional).

• Moulik, P. & Ekström, G., 2014. Dataset and Software for An anisotropic shear velocity model of the Earth's mantle using normal modes, body waves, surface waves and long-period waveforms. In Geophys. J. Int. (v1.0, Vol. 199, pp. 1713–1738). Zenodo. doi: 10.5281/zenodo.8357379

Data Products:

• S362ANIplusM_Figures.tar.gz - contains all figures from the paper in .png format
• S362ANI+M_MapViewsCrossSections.pdf - Some cross sections and map views at various depths in the mantle
• S362ANI+M - Coefficients of the spline basis functions for each parameter. Refer c_ij in equation 11.  This is our preferred global model of shear-wave velocity. In this model, radial anisotropy is confined to the uppermost mantle (that is, since the anisotropy is parameterized with only the four uppermost splines, it becomes very small below a depth of 250 km, and vanishes at 410 km). This is an updated version of S362ANI (Kustowski et al., 2008) which did not include normal modes in its derivation. Please note the stronger isotropic shear velocity anomalies in the transition zone.
• STW105 - reference model used in S362ANI+M. Described in Kustowski et al. (2008)
• setup.cfg -  Some configuration metadata relevant to this model for reproducibility.
• epix.tar.gz - Perturbations in horizontally (vsh) and vertically polarized shear velocity (vsv), Voigt-average isotropic velocity (vs), ansotropy (as) and topography of the internal boundaries. This is calculated from the spline coefficients at every 1 by 1 degree cell-centered pixel and at every ~25 km depth region from Moho to the core-mantle boundary and stored in extended pixel format (.epix) ASCII files. 
• S362ANI+M.BOX25km_PIX1X1.avni.nc4 -  The perturbations in a standard AVNI format that utilizes the NETCDF4 container format. This file can be read in Python using either xarray or AVNI libraries. For example, to plot vs perturbations at 24.4-50 km depth range
  • import xarray as xr
  • ds = xr.open_dataset('S362ANI+M.BOX25km_PIX1X1.avni.nc4')
  • ds.vs[0].plot()
• Fortran Code
  • readme - contains a description of all files in the folder below
  • PROGRAMS.tar.gz - tools for obtaining model values at specific locations and some GMT plotting tools
• GRD.tar.gz - longitude-latitide-velocity files with vsh, vsv, and Voigt average in km/s evaluated on a grid of points at many depths in the mantle