Global Topography Data

Local Scale Data Sets

This pages deals with global topography data sets on Earth up to a maximum spatial resolution of 1.8 km. If you need even higher resolution data please use for example the AppEEARS service from the USGS.

A key feature of Earth and other planets are its external (e.g. the surface) and internal (e.g. the Moho) boundaries. Salvus requires access to to this data in order to create a computational model of the planet.

Mondaic maintains an open-source tool SalvusTopo to create these files. Please have a look at the repository for more information and also if you are interested in creating your own files. It teaches you everything that is required. The files themselves are stored in NetCDF files following the CF convention.

Data Sources

Most of the downloadable files combine a number of freely available data sets. If you use it in a scientific context, please remember to cite the original data source, where appropriate. Files we offer on this page specifically utilize three external data sets:

Earth Gravitational Model 2008 (EGM2008)

Pavlis, Nikolaos K., Simon A. Holmes, Steve C. Kenyon, and John K. Factor. (2012).
The Development and Evaluation of the Earth Gravitational Model 2008 (EGM2008)
Journal of Geophysical Research: Solid Earth, 117 (4): 1–38

Specifically we use the Geoid undulations with respect to WGS84. In the following we'll refer to this as Geoid.

Earth2014 global topography (relief) model

Hirt, C. and Rexer, M. (2015)
Earth2014: 1 arc-min shape, topography, bedrock and ice-sheet models - available as gridded data and degree-10,800 spherical harmonics
International Journal of Applied Earth Observation and Geoinformation, 39, 103-112

Height in these is always relative to mean sea level.

Specifically we use two pieces of data: The spherical harmonics coefficients for the Earth surface (referred to as Earth2014_SUR) and the for the Topo bedrock, bathymetry, and ice (referred to as Earth2014_TBI).


Laske, G., Masters., G., Ma, Z. and Pasyanos, M. (2013)
Update on CRUST1.0 - A 1-degree Global Model of Earth's Crust
Geophysical Research Abstracts, 15, EGU2013-2658, EGU General Assembly 2013

We use CRUST1.0 to compute the position of the Moho boundary on Earth.

Specifically we use three pieces of data: The topography (referred to as CRUST1.0_TOPO), the ocean (referred to as CRUST1.0_OCEAN), and the crustal thickness (referred to as CRUST1.0_CRUSTAL_THICKNESS).

Preprocessed Data Sets for Download

The following data sets can be directly used by SalvusMesh to generate realistic Earth meshes. Please let us know if you'd like to see other data sets hosted. Aside from file size and data transfer concerns there is no penalty for using the higher resolution models as every file always contains the lower angular orders as well. What we mean by angular order is explained at the end of this page.

Global Topography Data Set

This data set combines the EGM2008 Geoid with surface data from Earth2014 to yield the absolute surface elevation relative to the WGS84 ellipsoid.

Topography = Geoid + Earth2014_TBI

lmaxl_{max} orderFile SizeDownload URL

Global Bathymetry Data Set

This data set represents the bathymetry, e.g. the depth of the oceans. It is important that this data set is consistent with the used topography data set. You only need this if you intend to simulation an ocean layer either directly or via an effective ocean load.

Bathymetry = Earth2014_SUR - Earth2014_TBI

lmaxl_{max} orderFile SizeDownload URL

Global Crustal Data Set

This data set represents the position of the Moho boundary on Earth computed as:

Moho Topography = Geoid + CRUST1.0_TOPO - CRUST_1.0_OCEAN - CRUST1.0_CRUSTAL_THICKNESS

lmaxl_{max} orderFile SizeDownload URL

Angular Order

To avoid spatial aliasing it is of crucial importance to use an appropriately filtered version for the chosen mesh. SalvusTopo internally represents models as spherical harmonics and can thus easily produce differently filtered versions depending on the maximum degree. In the following we'll refer to this as lmaxl_{max}. Some interfaces in Salvus choose this automatically, some require you to make a choice. The maximum angular order lmaxl_{max} approximately corresponds to spatial resolutions on Earth as given in the following table:

lmaxl_{max}Approximate equivalent spatial resolution on Earth's surface
16~ 1250 km
32~ 625 km
64~ 310 km
128~ 155 km
256~ 80 km
512~ 40 km
1024~ 20 km
2048~ 10 km
4096~ 5 km
10800~ 1.8 km

A single output file of SalvusTopo will typically contain a range of maximum angular orders. An advantage of using the CF convention is that the models can be readily visualized in Paraview:

This images demonstrates the difference between a global Earth surface topography model of degree 64 vs degree 1024.