Potential habitat evaluation and studying of geographical species distribution is a key subjects in most of the ecological environment protection and wildlife as well as variation trend studies in different scales. Species distribution prediction models have been becoming one of the important tool to prepare potential conservation map in recent years. Typically, such models define relationship between environmental variables and species events (present and absence) to evaluate environmental condition as mathematical models. Species distribution models can be used to assess climate change scenarios and their effect on species geographical distribution. The current study is aimed to predict potential distribution of Daphne mucronata Royle and Kelossia odoratissima Mozaf and its distribution under different climate scenarios in years 2030 and 2080. For this, 19 climate variables and 3 physiographic variables with resolution 1 km was prepared using IPCC. For species Kelossia odoratissima Mozaf 30 sites and for Daphne mucronata Royle 75 sites of species Key words: Species distribution model, Geographic Information System, Generalized Additive Models, Bayesian Belief Networks, Climate change, Fereydoonshahr, Isfahan province