Forecasting Large-Scale Habitat Suitability of European Bustards under Climate Change: The Role of Environmental and Geographic Variables

We modelled the distribution of two vulnerable steppe birds, Otis tarda and Tetrax tetrax, in the Western Palearctic and projected their suitability up to the year 2080. We performed two types of models for each species: one that included environmental and geographic variables (space-included model) and a second one that only included environmental variables (space-excluded model). Our assumption was that ignoring geographic variables in the modelling procedure may result in inaccurate forecasting of species distributions. On the other hand, the inclusion of geographic variables may generate an artificial constraint on future projections. Our results show that space-included models performed better than space-excluded models. While distribution of suitable areas for T. tetrax in the future was approximately the same as at present in the space-included model, the space-excluded model predicted a pronounced geographic change of suitable areas for this species. In the case of O. tarda, the space-included model showed that many areas of current presence shifted to low or medium suitability in the future, whereas a northward expansion of intermediate suitable areas was predicted by the space-excluded one. According to the best models, current distribution of these species can restrict future distribution, probably due to dispersal constraints and site fidelity. Species ranges would be expected to shift gradually over the studied time period and, therefore, we consider it unlikely that most of the current distribution of these species in southern Europe will disappear in less than one hundred years. Therefore, populations currently occupying suitable areas should be a priority for conservation policies. Our results also show that climate-only models may have low explanatory power, and could benefit from adjustments using information on other environmental variables and biological traits; if the latter are not available, including the geographic predictor may improve the reliability of predicted results.
AE has a contract funded by the project 1098/2014 (Organismo Autónomo Parques Nacionales, Spain). This paper is a contribution to CGL2009-13029/BOS of the Spanish Ministry of Science, as well as to the REMEDINAL 3 (S2013/MAE-2719) network of the CAM. This work was partly supported by the Project CGL2009/11316/BOS (Spanish Ministry of Science and FEDER).