Two endemic plant species in the Mediterranean island of Cyprus, Crocus cyprius and Ophrys kotschyi, were used as a case study. We have coupled climate change scenarios, and land use change models with species distribution models. Future land use scenarios were modelled by initially calculating the rate of current land use changes between two time snapshots (2000 and 2006) on the island, and based on these transition probabilities markov-chain cellular automata were used to generate future land use changes for 2050. Climate change scenarios A1B, A2, B1 and B2A were derived from the IPCC reports. Species climatic preferences were derived from their current distributions using classification trees while habitats preferences were derived from the Red Data Book of the Flora of Cyprus. A bioclimatic model for Crocus cyprius was built using mean temperature of wettest quarter, max temperature of warmest month and precipitation seasonality, while for Ophrys kotchyi the bioclimatic model was built using precipitation of wettest month, mean temperature of warmest quarter, isothermality, precipitation of coldest quarter, and annual precipitation. Sequentially, simulation scenarios were performed regarding future species distributions by accounting climate alone and both climate and land use changes. The distribution of the two species resulting from the bioclimatic models was then filtered by future land use changes, providing the species projected potential distribution. The species projected potential distribution varies depending on the type and scenario used, but many of both species current sites/locations are projected to be outside their future potential distribution. Our results demonstrate the importance of including both land use and climatic changes in predictive species modeling. ; Comment: to appear (in press in Ecological Informatics (2015))