Global climate change poses a new challenge for species and can even push some species toward an extinction vortex. The most affected organisms are those with narrow tolerance to the climatic factors but many large mammals such as ungulates with a wider ecological niche are also being affected indirectly. Our research mainly used wild sheep in central Iran as a model species to explore how the suitable habitats will change under different climatic scenarios and to determine if current borders of protected areas will adequately protect habitat requirements. To create habitat models we used animal-vehicle collision points as an input for species presence data. We ran habitat models using MaxEnt modeling approach under ifferent climatic scenarios of the past, present and future (under the climatic scenarios for minimum (RCP2.6) and maximum (RCP8.5) CO2 concentration trajectories). We tried to estimate the overlap and the width of the ecological niche using relevant metrics. In order to analyze the effectiveness of the protected areas, suitable maps were concerted to binary maps using True Skill Statistic (TSS) threshold and measured the similarity of the binary maps for each scenario using Kappa index. In order to assess the competence of the present protected areas boundary in covering the distribution of species, two different scenarios were employed, which are ensemble scenario 1: an ensemble of the binary maps of the species distribution in Mid-Holocene, present, and RCP2.6; and ensemble scenario 2: an ensemble of binary suitability maps in Mid-Holocene, present, and RCP8.5. Then, the borders of modeled habitats with the boundaries of 23 existing protected areas in two central provinces in Iran were compared. The predicted species distribution under scenario 1 (RCP2.6) was mostly similar to its current distribution (Kappa = 0.53) while the output model under scenario 2 (RCP8.5) indicated a decline in the species distribution range. Under the first ensemble scenario
