https://twitter.com/allixender/status/1325843393050710016
The DEM is a pre-clipped and reprojected (in QGIS) GeoTiff from the awesome MERIT-Hydro global hydrologically pre-processed DEM:
http://hydro.iis.u-tokyo.ac.jp/~yamadai/MERIT_Hydro/
It would be quite straight-forward to do the clipping and reprojecting in Rasterio, too (e.g. here).
For the scientific plot style I experimented with the HoloViews python package using a Surface3D and Scatter3D plot.
That’s all, folks :-)
The Jupyter notebook (view here) is provided in my GitHub repo.