Zoonotic cutaneous leishmaniasis (ZCL) is among the most endemic zoonotic diseases in Golestan Province of Iran. The aim of this study was to find the high risk areas of this infection by considering the distribution of reservoirs and human infection. A cross-sectional study was conducted, in which Rhombomys opimus (Gerbils) were captured from different collection sites across the Golestan Province, Iran. records about the occurrence of Rh. opimus in the Province was obtained from earlier studies and were gathered in a database. Furthermore, records about the disease existence were also obtained from the health system database of Golestan Province. Villages with at least three cases of ZcL were considered endemic foci and used as presence sites for Leishmania major. ArcGIS and MaxEnt model were used to map and predict the best ecological niches for both reservoir and parasite. According to the MaxEnt model, the area under Roc curve for Rh. opimus and L. major was 0.92 and 0.89, respectively. The probability of presence for both species in the northeastern part of Golestan Province was more than the other parts. The Jackknife test indicated that factors like temperature and altitude plays significant role in predicting the environmental suitability for ZcL reservoir and parasite, respectively. this modeling approach predicted the areas suitable for reservoir host and circulation of parasite to human. These findings can be used in proper mapping, surveillance and control of the CL.