Zoonotic cutaneous leishmaniasis (ZCL) constitutes a serious public health problem in many parts of the world including Iran. This study was carried out to assess the risk of the disease in an endemic province by developing spatial environmentally based models in yearly intervals. To fill the gap of underestimated true burden of ZCL and short study period, analytical hierarchy process (AHP) and fuzzy AHP decision-making methods were used to determine the ZCL risk zones in a Geographic Information System platform. Generated risk maps showed that high-risk areas were predominantly located at the northern and northeastern parts in each of the three study years. Comparison of the generated risk maps with geocoded ZCL cases at the village level demonstrated that in both methods more than 90%, 70% and 80% of the cases occurred in high and very high risk areas for the years 2010, 2011, and 2012, respectively. Moreover, comparison of the risk categories with spatially averaged normalized difference vegetation index (NDVI) images and a digital elevation model of the study region indicated persistent strong negative relationships between these environmental variables and ZCL risk degrees. These findings identified more susceptible areas of ZCL and will help the monitoring of this zoonosis to be more targeted.