The DEM problem

One of the main goals of rgitools is to provide not only aggregated statistics (hypsometry), but also a local topography map for each single glacier in the RGI.

Unfortunately, there is no gap-free and freely available global DEM to date. For most regions we can rely on a number of sources, which all have their own issues. Here is the list currently supported by OGGM/rgitools:

In theory, this should be enough data. In practice, none of the solutions above is 100% satisfying. We didn’t do a thorough assesment yet, but you can have a look at the following examples:

TODO: examples in Antarctica

Summary

We have to keep in mind that any kind of glacier bed estimate inversion or glacier simulation based on ice-dynamics cannot deal with artefacts. Therefore, robust and gap-free datasets are much preferred over more accurate but incomplete DEMS. Furthermore, the concurrent timing of the glacier outline with the DEM is another important criterion, as shown by the Columbia example.

Altogether, we are confident in: - SRTM for all latitudes below 60° N and S - GIMP for Greenland - RAMP for Antarctica (but the peripheral Islands are an issue)

For everything else, more investigation is needed. DEM3 offers the stability and timeliness that TanDEM-X and ArcticDEM cannot offer (yet), and is the current default in OGGM.