Using 35 presence-only data samples and five uncorrelated bioclimatic variables, we made species distribution models (SDMs) for 4 species of critically endangered (CR) liverworts from genus Jungermanniales and Marchantiales (Cephaloziella elegans, Leiocolea heterocolpos, Lophozia wenzelii and Riccia papillosa) using the maximum entropy modelling method (MaxEnt). Since we were modelling CR species, only one model proved to be strong enough to be used in the field. However, SDMs can serve as effective and fast tools for acceleration of the discovery of the rare and endangered species. The final model presented in this study can serve as a guide to future survey expeditions, the conservation of the target species and also to help understand their ecology.