Modeling plant species distributions under future climates: how fine scale do climate projections need to be?

Recent studies suggest that species distribution models (SDMs) based on fine‐scale climate data may provide markedly different estimates of climate‐change impacts than coarse‐scale models. However, these studies disagree in their conclusions of how scale influences projected species distributions. In rugged terrain, coarse‐scale climate grids may not capture topographically controlled climate variation at the scale that constitutes microhabitat or refugia for some species. Although finer scale data are therefore considered to better reflect climatic conditions experienced by species, there have been few formal analyses of how modeled distributions differ with scale. We modeled distributions for 52 plant species endemic to the California Floristic Province of different life forms and range sizes under recent and future climate across a 2000‐fold range of spatial scales (0.008–16 km²). We produced unique current and future climate datasets by separately downscaling 4 km climate models to three finer resolutions based on 800, 270, and 90 m digital elevation models and deriving bioclimatic predictors from them. As climate‐data resolution became coarser, SDMs predicted larger habitat area with diminishing spatial congruence between fine‐ and coarse‐scale predictions. These trends were most pronounced at the coarsest resolutions and depended on climate scenario and species' range size. On average, SDMs projected onto 4 km climate data predicted 42% more stable habitat (the amount of spatial overlap between predicted current and future climatically suitable habitat) compared with 800 m data. We found only modest agreement between areas predicted to be stable by 90 m models generalized to 4 km grids compared with areas classified as stable based on 4 km models, suggesting that some climate refugia captured at finer scales may be missed using coarser scale data. These differences in projected locations of habitat change may have more serious implications than net habitat area when predictive maps form the basis of conservation decision making.     

Coarse climate change projections for species living in a fine‐scaled world

Accurately predicting biological impacts of climate change is necessary to guide policy. However, the resolution of climate data could be affecting the accuracy of climate change impact assessments. Here, we review the spatial and temporal resolution of climate data used in impact assessments and demonstrate that these resolutions are often too coarse relative to biologically relevant scales. We then develop a framework that partitions climate into three important components: trend, variance, and autocorrelation. We apply this framework to map different global climate regimes and identify where coarse climate data is most and least likely to reduce the accuracy of impact assessments. We show that impact assessments for many large mammals and birds use climate data with a spatial resolution similar to the biologically relevant area encompassing population dynamics. Conversely, impact assessments for many small mammals, herpetofauna, and plants use climate data with a spatial resolution that is orders of magnitude larger than the area encompassing population dynamics. Most impact assessments also use climate data with a coarse temporal resolution. We suggest that climate data with a coarse spatial resolution is likely to reduce the accuracy of impact assessments the most in climates with high spatial trend and variance (e.g., much of western North and South America) and the least in climates with low spatial trend and variance (e.g., the Great Plains of the USA). Climate data with a coarse temporal resolution is likely to reduce the accuracy of impact assessments the most in the northern half of the northern hemisphere where temporal climatic variance is high. Our framework provides one way to identify where improving the resolution of climate data will have the largest impact on the accuracy of biological predictions under climate change.     

How to survive a glaciation: the challenge of estimating biologically realistic potential distributions under freezing conditions

Correlative ecological niche models are increasingly used to estimate potential distributions during the Last Glacial Maximum (LGM) for biogeographical research. In the case of presence‐background/pseudoabsences techniques, cold environments that are poorly represented in existing geography can complicate the process of model calibration and transfer into more extreme cold environments that were very common during the LGM (non‐analog conditions). This may lead to biologically unrealistic estimations. Using one cold‐adapted North American mammal, we explore a real scenario to better understand the effect of restricting the range of environmental conditions over which niche models are calibrated and then transferred to LGM conditions. We performed two sets of experiments in Maxent: 1) we calibrated models in the context of only present‐day climate conditions, which is the most common practice, and compared predictions under LGM conditions based on two extrapolation methods (clamping versus unconstrained); 2) we calibrated single models using both present‐day and LGM conditions as part of the same background in order to include more extreme environments in the model calibration. Our experiments led to dramatically different estimates of species’ potential distributions, showing notable differences with respect to latitudinal and elevational shifts during the LGM. Models calibrated using present‐day climates yielded biologically unrealistic estimations, suggesting that species survived in the glaciers during the LGM. Even more unrealistic estimations were achieved when clamping was enforced as the method to extrapolate. Models calibrated in the context of both modern and past climates reduced the required degree of extrapolation and allowed more realistic potential distributions, suggesting that the species avoided extremely cold conditions during the LGM. This study alerts to the possibility of obtaining implausible potential distributions during the LGM due to restricted background datasets and offers recommendations that should promote better strategies to estimate distributional changes during glaciations.     

How will climate novelty influence ecological forecasts? Using the Quaternary to assess future reliability

Future climates are projected to be highly novel relative to recent climates. Climate novelty challenges models that correlate ecological patterns to climate variables and then use these relationships to forecast ecological responses to future climate change. Here, we quantify the magnitude and ecological significance of future climate novelty by comparing it to novel climates over the past 21,000 years in North America. We then use relationships between model performance and climate novelty derived from the fossil pollen record from eastern North America to estimate the expected decrease in predictive skill of ecological forecasting models as future climate novelty increases. We show that, in the high emissions scenario (RCP 8.5) and by late 21st century, future climate novelty is similar to or higher than peak levels of climate novelty over the last 21,000 years. The accuracy of ecological forecasting models is projected to decline steadily over the coming decades in response to increasing climate novelty, although models that incorporate co‐occurrences among species may retain somewhat higher predictive skill. In addition to quantifying future climate novelty in the context of late Quaternary climate change, this work underscores the challenges of making reliable forecasts to an increasingly novel future, while highlighting the need to assess potential avenues for improvement, such as increased reliance on geological analogs for future novel climates and improving existing models by pooling data through time and incorporating assemblage‐level information.     

The role of land cover in bioclimatic models depends on spatial resolution

Aim We explored the importance of climate and land cover in bird species distribution models on multiple spatial scales. In particular, we tested whether the integration of land cover data improves the performance of pure bioclimatic models. Location Finland, northern Europe. Methods The data of the bird atlas survey carried out in 1986–89 using a 10 × 10 km uniform grid system in Finland were employed in the analyses. Land cover and climatic variables were compiled using the same grid system. The dependent and explanatory variables were resampled to 20‐km, 40‐km and 80‐km resolutions. Generalized additive models (GAM) were constructed for each of the 88 land bird species studied in order to estimate the probability of occurrence as a function of (1) climate and (2) climate and land cover variables. Model accuracy was measured by a cross‐validation approach using the area under the curve (AUC) of a receiver operating characteristic (ROC) plot. Results In general, the accuracies of the 88 bird–climate models were good at all studied resolutions. However, the inclusion of land cover increased the performance of 79 and 78 of the 88 bioclimatic models at 10‐km and 20‐km resolutions, respectively. There was no significant improvement at the 40‐km resolution. In contrast to the finer resolutions, the inclusion of land cover variables decreased the modelling accuracy at 80km resolution. Main conclusions Our results suggest that the determinants of bird species distributions are hierarchically structured: climatic variables are large‐scale determinants, followed by land cover at finer resolutions. The majority of the land bird species in Finland are rather clearly correlated with climate, and bioclimate envelope models can provide useful tools for identifying the relationships between these species and the environment at resolutions ranging from 10 km to 80 km. However, the notable contribution of land cover to the accuracy of bioclimatic models at 10–20‐km resolutions indicates that the integration of climate and land cover information can improve our understanding and model predictions of biogeographical patterns under global change.     

Time of emergence of novel climates for North American migratory bird populations

To better understand the ecological implications of global climate change for species that display geographically and seasonally dynamic life‐history strategies, we need to determine where and when novel climates are projected to first emerge. Here, we use a multivariate approach to estimate time of emergence (ToE) of novel climates based on three climate variables (precipitation, minimum and maximum temperature) at a weekly temporal resolution within the Western Hemisphere over a 280‐yr period (2021–2300) under a high emissions scenario (RCP8.5). We intersect ToE estimates with weekly estimates of relative abundance for 77 passerine bird species that migrate between temperate breeding grounds in North America and southern tropical and subtropical wintering grounds using observations from the eBird citizen‐science database. During the non‐breeding season, migrants that winter within the tropics are projected to encounter novel climates during the second half of this century. Migrants that winter in the subtropics are projected to encounter novel climates during the first half of the next century. During the beginning of the breeding season, migrants on their temperate breeding grounds are projected to encounter novel climates during the first half of the next century. During the end of the breeding season, migrants are projected to encounter novel climates during the second half of this century. Thus, novel climates will first emerge ca 40–50 yr earlier during the second half of the breeding season. These results emphasize the large seasonal and spatial variation in the formation of novel climates, and the pronounced challenges migratory birds are likely to encounter during this century, especially on their tropical wintering grounds and during the transition from breeding to migration. When assessing the ecological implications of climate change, our findings emphasize the value of applying a full annual cycle perspective using standardized metrics that promote comparisons across space and time.     

Spatial scale,topography and thermoregulatory behaviour interact when modelling species’ thermal niches

The spatial scale at which climate and species’ occupancy data are gathered, and the resolution at which ecological models are run, can strongly influence predictions of species performance and distributions. Running model simulations at coarse rather than fine spatial resolutions, for example, can determine if a model accurately predicts the distribution of a species. The impacts of spatial scale on a model's accuracy are particularly pronounced across mountainous terrain. Understanding how these discrepancies arise requires a modelling approach in which the underlying processes that determine a species’ distribution are explicitly described. Here we use a process‐based model to explore how spatial resolution, topography and behaviour alter predictions of a species thermal niche, which in turn constrains its survival and geographic distribution. The model incorporates biophysical equations to predict the operative temperature (T_e), thermal‐dependent performance and survival of a typical insect, with a complex life‐cycle, in its microclimate. We run this model with geographic data from a mountainous terrain in South Africa using climate data at three spatial resolutions. We also explore how behavioural thermoregulation affects predictions of a species performance and survival by allowing the animal to select the optimum thermal location within each coarse‐grid cell. At the regional level, coarse‐resolution models predicted lower T_e at low elevations and higher T_e at high elevations than models run at fine‐resolutions. These differences were more prominent on steep, north‐facing slopes. The discrepancies in T_e in turn affected estimates of the species thermal niche. The modelling framework revealed how spatial resolution and topography influence predictions of species distribution models, including the potential impacts of climate change. These systematic biases must be accounted for when interpreting the outputs of future modelling studies, particularly when species distributions are predicted to shift from uniform to topographically heterogeneous landscapes.     

Geographic patterns and dynamics of Alaskan climate interpolated from a sparse station record

Data from a sparse network of climate stations in Alaska were interpolated to provide 1‐km resolution maps of mean monthly temperature and precipitation–‐variables that are required at high spatial resolution for input into regional models of ecological processes and resource management. The interpolation model is based on thin‐plate smoothing splines, which uses the spatial data along with a digital elevation model to incorporate local topography. The model provides maps that are consistent with regional climatology and with patterns recognized by experienced weather forecasters. The broad patterns of Alaskan climate are well represented and include latitudinal and altitudinal trends in temperature and precipitation and gradients in continentality. Variations within these broad patterns reflect both the weakening and reduction in frequency of low‐pressure centres in their eastward movement across southern Alaska during the summer, and the shift of the storm tracks into central and northern Alaska in late summer. Not surprisingly, apparent artifacts of the interpolated climate occur primarily in regions with few or no stations. The interpolation model did not accurately represent low‐level winter temperature inversions that occur within large valleys and basins. Along with well‐recognized climate patterns, the model captures local topographic effects that would not be depicted using standard interpolation techniques. This suggests that similar procedures could be used to generate high‐ resolution maps for other high‐latitude regions with a sparse density of data.     

A closer look at novel climates: new methods and insights at continental to landscape scales

Novel climates – emerging conditions with no analog in the observational record – are an open problem in ecological modeling. Detecting extrapolation into novel conditions is a critical step in evaluating bioclimatic projections of how species and ecosystems will respond to climate change. However, biologically informed novelty detection methods remain elusive for many modeling algorithms. To assist with bioclimatic model design and evaluation, we present a first‐approximation assessment of general novelty based on a simple and consistent characterization of climate. We build on the seminal global analysis of Williams et al. (2007 PNAS, 104, 5738) by assessing of end‐of‐21st‐century novelty for North America at high spatial resolution and by refining their standardized Euclidean distance into an intuitive Mahalanobian metric called sigma dissimilarity. Like this previous study, we found extensive novelty in end‐of‐21st‐century projections for the warm southern margin of the continent as well as the western Arctic. In addition, we detected localized novelty in lower topographic positions at all latitudes: By the end of the 21st century, novel climates are projected to emerge at low elevations in 80% and 99% of ecoregions in the RCP4.5 and RCP8.5 emissions scenarios, respectively. Novel climates are limited to 7% of the continent's area in RCP4.5, but are much more extensive in RCP8.5 (40% of area). These three risk factors for novel climates – regional susceptibility, topographic position, and the magnitude of projected climate change – represent a priori evaluation criteria for the credibility of bioclimatic projections. Our findings indicate that novel climates can emerge in any landscape. Interpreting climatic novelty in the context of nonlinear biological responses to climate is an important challenge for future research.     

Fine‐grain,large‐domain climate models based on climate station and comprehensive topographic information improve microrefugia detection

Large‐domain species distribution models (SDMs) fail to identify microrefugia, as they are based on climate estimates that are either too coarse or that ignore relevant topographic climate‐forcing factors. Climate station data are considered inadequate to produce such estimates, a viewpoint we challenge here. Using climate stations and topographic data, we developed three sets of large‐domain (450 000 km²), fine‐grain (50 m) temperature grids accounting for different levels of topographic complexity. Using these fine‐grain grids and the Worldclim data, we fitted SDMs for 78 alpine species over Sweden, and assessed over‐ versus underestimations of local extinction and area of microrefugia by comparing modelled distributions at species' rear edges. Accounting for well‐known topographic climate‐forcing factors improved our ability to model fine‐scale climate, despite using only climate station data. This approach captured the effect of cool air pooling, distance to sea, and relative humidity on local‐scale temperature, but the effect of solar radiation could not be accurately accounted for. Predicted extinction rate decreased with increasing spatial resolution of the climate models and with increasing number of topographic climate‐forcing factors accounted for. About half of the microrefugia detected in the most topographically complete models were not detected in the coarser SDMs and in the models calibrated from climate variables extracted from elevation only. Although major limitations remain, climate station data can potentially be used to produce fine‐grain topoclimate grids, opening up the opportunity to model local‐scale ecological processes over large domains. Accounting for the topographic complexity encountered within landscapes permits the detection of microrefugia that would otherwise remain undetected. Topographic heterogeneity is likely to have a massive impact on species persistence, and should be included in studies on the effects of climate change.     

Impact of landscape predictors on climate change modelling of species distributions: a case study with Eucalyptus fastigata in southern New South Wales,Australia

Aim We consider three questions. (1) How different are the predicted distribution maps when climate‐only and climate‐plus‐terrain models are developed from high‐resolution data? (2) What are the implications of differences between the models when predicting future distributions under climate change scenarios, particularly for climate‐only models at coarse resolution? (3) Does the use of high‐resolution data and climate‐plus‐terrain models predict an increase in the number of local refugia? Location South‐eastern New South Wales, Australia. Methods We developed two species distribution models for Eucalyptus fastigata under current climate conditions using generalized additive modelling. One used only climate variables as predictors (mean annual temperature, mean annual rainfall, mean summer rainfall); the other used both climate and landscape (June daily radiation, topographic position, lithology, nutrients) variables as predictors. Predictions of the distribution under current climate and climate change were then made for both models at a pixel resolution of 100 m. Results The model using climate and landscape variables as predictors explained a significantly greater proportion of the deviance than the climate‐only model. Inclusion of landscape variables resulted in the prediction of much larger areas of existing optimal habitat. An overlay of predicted future climate on the current climate space indicated that extrapolation of the statistical models was not occurring and models were therefore more robust. Under climate change, landscape‐defined refugia persisted in areas where the climate‐only model predicted major declines. In areas where expansion was predicted, the increase in optimal habitat was always greater with landscape predictors. Recognition of extensive optimal habitat conditions and potential refugia was dependent on the use of high‐resolution landscape data. Main conclusions Using only climate variables as predictors for assessing species responses to climate change ignores the accepted conceptual model of plant species distribution. Explicit statements justifying the selection of predictors based on ecological principles are needed. Models using only climate variables overestimate range reduction under climate change and fail to predict potential refugia. Fine‐scale‐resolution data are required to capture important climate/landscape interactions. Extrapolation of statistical models to regions in climate space outside the region where they were fitted is risky.     

No‐analog climates and shifting realized niches during the late quaternary: implications for 21st‐century predictions by species distribution models

Empirically derived species distributions models (SDMs) are increasingly relied upon to forecast species vulnerabilities to future climate change. However, many of the assumptions of SDMs may be violated when they are used to project species distributions across significant climate change events. In particular, SDM's in theory assume stable fundamental niches, but in practice, they assume stable realized niches. The assumption of a fixed realized niche relative to climate variables remains unlikely for various reasons, particularly if novel future climates open up currently unavailable portions of species’ fundamental niches. To demonstrate this effect, we compare the climate distributions for fossil‐pollen data from 21 to 15 ka bp (relying on paleoclimate simulations) when communities and climates with no modern analog were common across North America to observed modern pollen assemblages. We test how well SDMs are able to project 20th century pollen‐based taxon distributions with models calibrated using data from 21 to 15 ka. We find that taxa which were abundant in areas with no‐analog late glacial climates, such as Fraxinus, Ostrya/Carpinus and Ulmus, substantially shifted their realized niches from the late glacial period to present. SDMs for these taxa had low predictive accuracy when projected to modern climates despite demonstrating high predictive accuracy for late glacial pollen distributions. For other taxa, e.g. Quercus, Picea, Pinus strobus, had relatively stable realized niches and models for these taxa tended to have higher predictive accuracy when projected to present. Our findings reinforce the point that a realized niche at any one time often represents only a subset of the climate conditions in which a taxon can persist. Projections from SDMs into future climate conditions that are based solely on contemporary realized distributions are potentially misleading for assessing the vulnerability of species to future climate change.     

Projected Future Vegetation Changes for the Northwest United States and Southwest Canada at a Fine Spatial Resolution Using a Dynamic Global Vegetation Model

Future climate change may significantly alter the distributions of many plant taxa. The effects of climate change may be particularly large in mountainous regions where climate can vary significantly with elevation. Understanding potential future vegetation changes in these regions requires methods that can resolve vegetation responses to climate change at fine spatial resolutions. We used LPJ, a dynamic global vegetation model, to assess potential future vegetation changes for a large topographically complex area of the northwest United States and southwest Canada (38.0–58.0°N latitude by 136.6–103.0°W longitude). LPJ is a process-based vegetation model that mechanistically simulates the effect of changing climate and atmospheric CO₂ concentrations on vegetation. It was developed and has been mostly applied at spatial resolutions of 10-minutes or coarser. In this study, we used LPJ at a 30-second (~1-km) spatial resolution to simulate potential vegetation changes for 2070–2099. LPJ was run using downscaled future climate simulations from five coupled atmosphere-ocean general circulation models (CCSM3, CGCM3.1(T47), GISS-ER, MIROC3.2(medres), UKMO-HadCM3) produced using the A2 greenhouse gases emissions scenario. Under projected future climate and atmospheric CO₂ concentrations, the simulated vegetation changes result in the contraction of alpine, shrub-steppe, and xeric shrub vegetation across the study area and the expansion of woodland and forest vegetation. Large areas of maritime cool forest and cold forest are simulated to persist under projected future conditions. The fine spatial-scale vegetation simulations resolve patterns of vegetation change that are not visible at coarser resolutions and these fine-scale patterns are particularly important for understanding potential future vegetation changes in topographically complex areas.     

Back to the future: using historical climate variation to project near‐term shifts in habitat suitable for coast redwood

Studies that model the effect of climate change on terrestrial ecosystems often use climate projections from downscaled global climate models (GCMs). These simulations are generally too coarse to capture patterns of fine‐scale climate variation, such as the sharp coastal energy and moisture gradients associated with wind‐driven upwelling of cold water. Coastal upwelling may limit future increases in coastal temperatures, compromising GCMs’ ability to provide realistic scenarios of future climate in these coastal ecosystems. Taking advantage of naturally occurring variability in the high‐resolution historic climatic record, we developed multiple fine‐scale scenarios of California climate that maintain coherent relationships between regional climate and coastal upwelling. We compared these scenarios against coarse resolution GCM projections at a regional scale to evaluate their temporal equivalency. We used these historically based scenarios to estimate potential suitable habitat for coast redwood (Sequoia sempervirens D. Don) under ‘normal’ combinations of temperature and precipitation, and under anomalous combinations representative of potential future climates. We found that a scenario of warmer temperature with historically normal precipitation is equivalent to climate projected by GCMs for California by 2020–2030 and that under these conditions, climatically suitable habitat for coast redwood significantly contracts at the southern end of its current range. Our results suggest that historical climate data provide a high‐resolution alternative to downscaled GCM outputs for near‐term ecological forecasts. This method may be particularly useful in other regions where local climate is strongly influenced by ocean–atmosphere dynamics that are not represented by coarse‐scale GCMs.     

Designing ecological climate change impact assessments to reflect key climatic drivers

Identifying the climatic drivers of an ecological system is a key step in assessing its vulnerability to climate change. The climatic dimensions to which a species or system is most sensitive – such as means or extremes – can guide methodological decisions for projections of ecological impacts and vulnerabilities. However, scientific workflows for combining climate projections with ecological models have received little explicit attention. We review Global Climate Model (GCM) performance along different dimensions of change and compare frameworks for integrating GCM output into ecological models. In systems sensitive to climatological means, it is straightforward to base ecological impact assessments on mean projected changes from several GCMs. Ecological systems sensitive to climatic extremes may benefit from what we term the ‘model space’ approach: a comparison of ecological projections based on simulated climate from historical and future time periods. This approach leverages the experimental framework used in climate modeling, in which historical climate simulations serve as controls for future projections. Moreover, it can capture projected changes in the intensity and frequency of climatic extremes, rather than assuming that future means will determine future extremes. Given the recent emphasis on the ecological impacts of climatic extremes, the strategies we describe will be applicable across species and systems. We also highlight practical considerations for the selection of climate models and data products, emphasizing that the spatial resolution of the climate change signal is generally coarser than the grid cell size of downscaled climate model output. Our review illustrates how an understanding of how climate model outputs are derived and downscaled can improve the selection and application of climatic data used in ecological modeling.     

Global patterns of nonanalogous climates in the past and future derived from thermal and hydraulic factors

Nonanalogous climates (NACs), climates without modern analogs on Earth, challenge our understanding of eco‐evolutionary processes that shape global biodiversity, particularly because of their propensity to promote novel ecosystems. However, NAC studies are generally inadequate and partial. Specifically, systematic comparisons between the future and the past are generally lacking, and hydraulic NACs tend to be underemphasized. In the present study, by adopting a frequency distribution‐based method that facilitates the procedures of contributions parsing and conducting multiple comparisons, we provide a global overview of multidimensional NACs for both the past and the future within a unified framework. We show that NACs are globally prevalent, covering roughly half of the land area across the time‐periods under investigation, and have a high degree of spatial structure. Patterns of NACs differ dramatically between the past and the future. Hydraulic NACs are more complex both in spatial patterns and in major contributions of variables than are thermal NACs. However, hydraulic NACs are more predictable than originally thought. Generally, hydraulic NACs in the future (2100 AD) exhibit comparable predictability to thermal NACs in the last glacial maximum (LGM) (21k BP). Identifying these NAC patterns has potential implications on climate‐adaptive managements and preparing in advance to possibly frequent novel ecosystems. However, a learning‐from‐the‐past strategy might be of limited utility for management under present circumstances.     

Ice age legacies in the geographical distribution of tree species richness in Europe

Aim This study uses a high‐resolution simulation of the Last Glacial Maximum (LGM) climate to assess: (1) whether LGM climate still affects the geographical species richness patterns in the European tree flora and (2) the relative importance of modern and LGM climate as controls of tree species richness in Europe. Location The parts of Europe that were unglaciated during the LGM. Methods Atlas data on the distributions of 55 tree species were linked with data on modern and LGM climate and climatic heterogeneity in a geographical information system with a 60‐km grid. Four measures of species richness were computed: total richness, and richness of the 18 most restricted species, 19 species of medium incidence (intermediate species) and 18 most widespread species. We used ordinary least‐squares regression and spatial autoregressive modelling to test and estimate the richness–climate relationships. Results LGM climate constituted the best single set of explanatory variables for richness of restricted species, while modern climate and climatic heterogeneity was best for total and widespread species richness and richness of intermediate species, respectively. The autoregressive model with all climatic predictors was supported for all richness measures using an information‐theoretic approach, albeit only weakly so for total species richness. Among the strongest relationships were increases in total and intermediate richness with climatic heterogeneity and in restricted richness with LGM growing‐degree‐days. Partial regression showed that climatic heterogeneity accounted for the largest unique variation fraction for intermediate richness, while LGM climate was particularly important for restricted richness. Main conclusions LGM climate appears to still affect geographical patterns of tree species richness in Europe, albeit the relative importance of modern and LGM climate depends on range size. Notably, LGM climate is a strong richness control for species with a restricted range, which appear to still be associated with their glacial refugia.     

Do we need land‐cover data to model species distributions in Europe?

Aim To assess the influence of land cover and climate on species distributions across Europe. To quantify the importance of land cover to describe and predict species distributions after using climate as the main driver. Location The study area is Europe. Methods (1) A multivariate analysis was applied to describe land‐cover distribution across Europe and assess if the land cover is determined by climate at large spatial scales. (2) To evaluate the importance of land cover to predict species distributions, we implemented a spatially explicit iterative procedure to predict species distributions of plants (2603 species), mammals (186 species), breeding birds (440 species), amphibian and reptiles (143 species). First, we ran bioclimatic models using stepwise generalized additive models using bioclimatic variables. Secondly, we carried out a regression of land cover (LC) variables against residuals from the bioclimatic models to select the most relevant LC variables. Finally, we produced mixed models including climatic variables and those LC variables selected as decreasing the residual of bioclimatic models. Then we compared the explanatory and predictive power of the pure bioclimatic against the mixed model. Results (1) At the European coarse resolution, land cover is mainly driven by climate. Two bioclimatic axes representing a gradient of temperature and a gradient of precipitation explained most variation of land‐cover distribution. (2) The inclusion of land cover improved significantly the explanatory power of bioclimatic models and the most relevant variables across groups were those not explained or poorly explained by climate. However, the predictive power of bioclimatic model was not improved by the inclusion of LC variables in the iterative model selection process. Main conclusion Climate is the major driver of both species and land‐cover distributions over Europe. Yet, LC variables that are not explained or weakly associated with climate (inland water, sea or arable land) are interesting to describe particular mammal, bird and tree distributions. However, the addition of LC variables to pure bioclimatic models does not improve their predictive accuracy.     

Species richness, environmental correlates, and spatial scale: a test using South African birds

Energy and habitat heterogeneity are important correlates of spatial variation in species richness, though few investigations have sought to determine simultaneously their relative influences. Here we use the South African avifauna to examine the extent to which species richness is related to these variables and how these relationships depend on spatial grain. Taking spatial autocorrelation and area effects into account, we find that primary productivity, precipitation, absolute minimum temperature, and, at coarser resolutions, habitat heterogeneity account for most of the variation in species richness. Species richness and productivity are positively related, whereas the relationship between potential evapotranspiration (PET) and richness is unimodal. This is largely because of the constraining effects of low rainfall on productivity in high-PET areas. The increase in the importance of vegetation heterogeneity as an explanatory variable is caused largely by an increase in the range of vegetation heterogeneity included at coarse resolutions and is probably also a result of the positive effects of environmental heterogeneity on species richness. Our findings indicate that species richness is correlated with, and hence likely a function of, several variables, that spatial resolution and extent must be taken into account during investigations of these relationships, and that surrogate measures for productivity should be interpreted cautiously.