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 "Nearest Neighbour" interpolation algorithm was used in scipy.interpolate.griddata() method. The crustal thickness was calculated from MOHO - ELEVATION(Column 4 - Column 3). The lithospheric thickness was calculated from LAB - ELEVATION(Column 5 - Column 3).
