This study aims to estimate the responses of Cedrus libani, Abies cilicica, and Juniperus drupacea to the Last Glacial Maximum (LGM), present and future climate changes in Anatolia. For this purpose, the temporal and spatial distributions of these conifer species were modeled, based on global climate change scenarios. Accordingly, the temporal and spatial distributions of the studied species are predicted and back-projected using fossil pollen data, occurrence data, and 19 bioclimatic variables collected from the WorldClim database. The principal component analysis (PCA) method was applied to these variables, resulting in the identification of 8 variables for the species distribution model. Past and future climate information is based on CCSM4 and RCP 8.5 scenario was assumed for future projections. Furthermore, both MaxEnt 3.4.1 and ArcGIS 10.5 were utilized to develop the models.