SPRD6 (Ishii & Tromp, 1999, 2001)

SPRD6 contains 3D models of shear-wave velocity, compressional-wave velocity, and density within the mantle and 2D models of topography on the free-surface (dynamic), the 660-km discontinuity, and the core-mantle boundary. The 3D models are parameterised in terms of spherical harmonics laterally (up to degree 6) and Chebyshev polynomials radially (up to order 13).
There are two spherical harmonic conventions used. One is the fully normalised "complex" spherical harmonics (Edmonds, 1960) and the other is the "real" spherical harmonics in which most Harvard 3D model coefficients are given (e.g., Su & Dziewonski, 1997). Model coefficients are available in both conventions.
The model coefficients are given in relative perturbations (delta x/x).
For convenience, there is also a tar file (plot_sprd6.tar) which contains files of model coefficients, FORTRAN programs and a README file.

• Shear Velocity (Real Spherical Harmonics, Complex Spherical Harmonics)
• Compressional Velocity (Real Spherical Harmonics, Complex Spherical Harmonics)
• Density (Real Spherical Harmonics, Complex Spherical Harmonics)
• Free-Surface Dynamic Topography (Real Spherical Harmonics, Complex Spherical Harmonics)
• 660-km Discontinuity Topography (Real Spherical Harmonics, Complex Spherical Harmonics)
• Core-Mantle Boundary Topography (Real Spherical Harmonics, Complex Spherical Harmonics)

The files (except for those containing real spherical harmonic coefficients of 3-D models) contain

There are 5 FORTRAN source codes (which have to be compiled once they are transferred to your machine).

• plot.f
  This program reads a model (either 3D or 2D) and outputs a file (with specified name) which contains longitude, latitude and value of the model. For 3D models, the values are in percentage, and for 2D models, the values are in kilometres. The increment of latitude/longitude can be changed by changing the value of parameter "block_size".
  In order to run this program, one needs all the following subroutines and also 3D models (in real spherical harmonics - sprd6_s, sprd6_p, and sprd6_r) and 2D models (in complex spherical harmonics - sprd6_fs.XCS, sprd6_d660.XCS, and sprd6_dcmb.XCS).
  
• sub_chebyshev.f
  This subroutine calculates the values of Chebyshev polynomials at a given depth (within the mantle) with appropriate normalisation as discussed in Su (1992).
• sub_kernel.f
  This subroutine builds a matrix that contains values of (real) spherical harmonics at each latitude and longitude.
• sub_readModel.f
  This subroutine reads a file containing model coefficients and converts them (if necessary) into appropriate spherical harmonics convention to be consistent with spherical harmonic convention used to create the matrix by the subprogram sub_kernel.f.
• sub_switchLatLon.f
  This subroutine switches latitude and longitude of the data vector so that the output file (from plot.f) will contain values of longitude, latitude and model.
• Makefile
  If one prefers, download this file and type
  make plot
  to compile plot.f with its subroutines (all the FORTRAN codes must be in the same directory). Then, typing
  plot.ex
  will run the program plot.f.

• Supplementary Table 1: Statistical results for various inversions.
• Supplementary Figure 1: Comparison of density models obtained from different starting models.
• Supplementary Figure 2: Comparison of density and P velocity models obtained from different starting models.
• Supplementary Figure 3: Comparison of density models with various parameterisations.
• Supplementary Figure 4: Comparison of density models with and without boundary topography.
• Supplementary Figure 5: Resolution tests with input peak at a depth.
• Supplementary Figure 6: Backus-Gilbert resolution kernels.
• Supplementary Figure 7: Checker-board resolution tests.
• Supplementary Figure 8: Resolution tests for the 660 and the CMB.
• Supplementary Figure 9: Mode-by-mode improvements in fit.
• Supplementary Figure 10: Effects of random noise in the data on the inversion.
• Supplementary Figure 11: Resolution matrix for the real part of degree 2 order 1.
• Supplementary Figure 12: Dependence on data sets (Density models at 2850 km depth).
• Supplementary Figure 13: Dependence on data sets (Comparison of models).
• Supplementary Figure 14: Density model using different radial basis function.

• Edmonds, A.R., 1960.
  Angular Momentum in Quantum Mechanics, Princeton Univ. Press, Princeton, N.J.
• Ishii, M., and Tromp, J., 1999.
  Normal-mode and free-air gravity constraints on lateral variations in velocity and density of Earth's mantle. Science 285, 1231-1236.
• Ishii, M., and Tromp, J., 2001.
  Even-degree lateral variations in the Earth's mantle constrained by free oscillations and the free-air gravity anomaly. Geophys. J. Int., 145, 77-96.
• Su, W.-J., 1992.
  The Three-Dimensional Shear-Wave Velocity Structure of the Earth's Mantle. Ph.D. Thesis, Harvard Univ.
• Su, W.-J., and Dziewoński, A.M., 1997.
  Simultaneous inversion for 3-D variations in shear and bulk velocity in the mantle. Phys. Earth Plnet. Inter. 100, 135-156.