CryoGrid 3 is a land-surface scheme dedicated to modeling of ground temperatures in permafrost environments. Its excess ice module Xice is capable of simulating ground subsidence and thermokarst lake formation due to melting of excess ground ice. The basic version of CryoGrid 3 and Xice is described in the following article:

Westermann, S., Langer, M., Boike, J., Heikenfeld, M., Peter, M., Etzelmüller, B., & Krinner, G. (2016). Simulating the thermal regime and thaw processes of ice-rich permafrost ground with the land-surface model CryoGrid 3. Geosci. Model Dev., 9(2), 523–546. https://doi.org/10.5194/gmd-9-523-2016

This version v1.2.0 has been extended by a multi-scale tiling scheme. It has been used for simulations presented in a research article submitted to The Cryosphere and published as a preprint in The Cryosphere Discussions:

Nitzbon, J., Langer, M., Martin, L. C. P., Westermann, S., Schneider von Deimling, T., & Boike, J. (2020). Effects of multi-scale heterogeneity on the simulated evolution of ice-rich permafrost lowlands under a warming climate. The Cryosphere Discussions, 1–29. https://doi.org/10.5194/tc-2020-137