This virtual globe is created using the data from the paper below.
Juan Carlos Afonso, Farshad Salajegheh, Wolfgang Szwillus, Jorg Ebbing, Carmen Gaina, A global reference model of the lithosphere and upper mantle from joint inversion and analysis of multiple data sets, Geophysical Journal International, Volume 217, Issue 3, June 2019, Pages 1602–1628, https://doi.org/10.1093/gji/ggz094
The original data is in CSV text format( Afonso_etal_GJI_2019_LithoRef18.xyz).
The NetCDF files being created by GMT Surface and scipy.interpolate.griddata() method can be downloaded using the links below.
The "gmt surface" command can create smoother grids than scipy.interpolate.griddata() using nearest neighbour interpolation algorithm. But it also creates artifacts. The "Nearest Neighbour" interpolation algorithm was used in scipy.interpolate.griddata() method. Although it is less smoother, the data in the grids are more accurate. The crustal thickness was calculated from MOHO - ELEVATION(Column 4 - Column 3). The lithospheric thickness was calculated from LAB - ELEVATION(Column 5 - Column 3).
