Presence-absence versus presence-only modelling methods for predicting bird habitat suitability

Habitat suitability models can be generated using methods requiring information on species presence or species presence and absence. Knowledge of the predictive performance of such methods becomes a critical issue to establish their optimal scope of application for mapping current species distributions under different constraints. Here, we use breeding bird atlas data in Catalonia as a working example and attempt to analyse the relative performance of two methods: the Ecological Niche factor Analysis (ENFA) using presence data only and Generalised Linear Models (GLM) using presence/absence data. Models were run on a set of forest species with similar habitat requirements, but with varying occurrence rates (prevalence) and niche positions (marginality). Our results support the idea that GLM predictions are more accurate than those obtained with ENFA. This was particularly true when species were using available habitats proportionally to their suitability, making absence data reliable and useful to enhance model calibration. Species marginality in niche space was also correlated to predictive accuracy, i.e. species with less restricted ecological requirements were modelled less accurately than species with more restricted requirements. This pattern was irrespective of the method employed. Models for wide-ranging and tolerant species were more sensitive to absence data, suggesting that presence/absence methods may be particularly important for predicting distributions of this type of species. We conclude that modellers should consider that species ecological characteristics are critical in determining the accuracy of models and that it is difficult to predict generalist species distributions accurately and this is independent of the method used. Being based on distinct approaches regarding adjustment to data and data quality, habitat distribution modelling methods cover different application areas, making it difficult to identify one that should be universally applicable. Our results suggest however, that if absence data is available, methods using this information should be preferably used in most situations

Predicting freshwater habitat integrity using land-use surrogates

Freshwater biodiversity is globally threatened due to human disturbances, but freshwater ecosystems have been accorded lessprotection than their terrestrial and marine counterparts. Few criteria exist for assessing the habitat integrity of rivers and data used for such assessments are generally of limited geographical coverage. Here, we use a fine-scale dataset describing river integrity in north-western South Africa to explore the extent to which measures of freshwater habitat integrity can be predicted from remotely sensed data, which are readily available in many parts of the world. A spatial statistical model was built using broad land-cover variables to predict the habitat integrity (subdivided into riparian and instream integrity) of rivers.We also explored the importance of the spatial scale. Results showed that riparian and, to a lesser degree, instream habitat integrity of river systems could be predicted with reasonable accuracy. The total area under natural vegetation was the most significant predictor of riparian integrity, which is best predicted by land-use activities at catchment level, rather than more locally. Our GIS-based model thus provides a fine-scale approach to assessing river habitat integrity as a supplement to landscape-level conservation plans for river systems, and represents a significant contribution towards the monitoring componentof the River Health Programme (RHP), which reports on the state of rivers in South Africa

How to best threshold and validate stacked species assemblages? Community optimisation might hold the answer

1. The popularity of species distribution models (SDMs) and the associated stacked species distribution models (S-SDMs), as tools for community ecologists, largely increased in recent years. However, while some consensus was reached about the best methods to threshold and evaluate individual SDMs, little agreement exists on how to best assemble individual SDMs into communities, i.e. how to build and assess S-SDM predictions. 2. Here, we used published data of insects and plants collected within the same study region to test (1) if the most established thresholding methods to optimize single species prediction are also the best choice for predicting species assemblage composition, or if community-based thresholding can be a better alternative, and (2) whether the optimal thresholding method depends on taxa, prevalence distribution and/or species richness. Based on a comparison of different evaluation approaches we provide guidelines for a robust community cross-validation framework, to use if spatial or temporal independent data are unavailable. 3. Our results showed that the selection of the “optimal” assembly strategy mostly depends on the evaluation approach rather than taxa, prevalence distribution, regional species pool or species richness. If evaluated with independent data or reliable cross-validation, community-based thresholding seems superior compared to single species optimisation. However, many published studies did not evaluate community projections with independent data, often leading to overoptimistic community evaluation metrics based on single species optimisation. 4. The fact that most of the reviewed S-SDM studies reported over-fitted community evaluation metrics highlights the importance of developing clear evaluation guidelines for community models. Here, we move a first step in this direction, providing a framework for cross-validation at the community level

Improving spatial predictions of taxonomic, functional and phylogenetic diversity

In this study, we compare two community modelling approaches to determine their ability to predict the taxonomic, functional and phylogenetic properties of plant assemblages along a broad elevation gradient and at a fine resolution. The first method is the standard stacking individual species distribution modelling (SSDM) approach, which applies a simple environmental filter to predict species assemblages. The second method couples the SSDM and macroecological modelling (MEMSSDM-MEM) approaches to impose a limit on the number of species co-occurring at each site. Because the detection of diversity patterns can be influenced by different levels of phylogenetic or functional trees, we also examine whether performing our analyses from broad to more exact structures in the trees influences the performance of the two modelling approaches when calculating diversity indices. We found that coupling the SSDM with the MEM improves the overall predictions for the three diversity facets compared with those of the SSDM alone. The accuracy of the SSDM predictions for the diversity indices varied greatly along the elevation gradient, and when considering broad to more exact structure in the functional and phylogenetic trees, the SSDM-MEM predictions were more stable. SSDM-MEM moderately but significantly improved the prediction of taxonomic diversity, which was mainly driven by the corrected number of predicted species. The performance of both modelling frameworks increased when predicting the functional and phylogenetic diversity indices. In particular, fair predictions of the taxonomic composition by SSDM-MEM led to increasingly accurate predictions of the functional and phylogenetic indices, suggesting that the compositional errors were associated with species that were functionally or phylogenetically close to the correct ones; however, this did not always hold for the SSDM predictions.Synthesis. In this study, we tested the use of a recently published approach that couples species distribution and macroecological models to provide the first predictions of the distribution of multiple facets of plant diversity: taxonomic, functional and phylogenetic. Moderate but significant improvements were obtained; thus, our results open promising avenues for improving our ability to predict the different facets of biodiversity in space and time across broad environmental gradients when functional and phylogenetic information is available

Evaluating the ecological realism of plant species distribution models with ecological indicator values

Species distribution models (SDMs) are routinely applied to assess current as well as future species distributions, for example to assess impacts of future environmental change on biodiversity or to underpin conservation planning. It has been repeatedly emphasized that SDMs should be evaluated based not only on their goodness of fit to the data, but also on the realism of the modelled ecological responses. However, possibilities for the latter are hampered by limited knowledge on the true responses as well as a lack of quantitative evaluation methods. Here we compared modelled niche optima obtained from European-scale SDMs of 1,476 terrestrial vascular plant species with empirical ecological indicator values indicating the preferences of plant species for key environmental conditions. For each plant species we first fitted an ensemble SDM including three modeling techniques (GLM, GAM and BRT) and extracted niche optima for climate, soil, land use and nitrogen deposition variables with a large explanatory power for the occurrence of that species. We then compared these SDM-derived niche optima with the ecological indicator values by means of bivariate correlation analysis. We found weak to moderate correlations in the expected direction between the SDM-derived niche optima and ecological indicator values. The strongest correlation occurred between the modelled optima for growing degree days and the ecological indicator values for temperature. Correlations were weaker for SDM-derived niche optima with a more distal relationship to ecological indicator values (notably precipitation and soil moisture). Further, correlations were consistently highest for BRT, followed by GLM and GAM. Our method gives insight into the ecological realism of modelled niche optima and projected core habitats and can be used to improve SDMs by making a more informed selection of environmental variables and modeling techniques

Back from a Predicted Climatic Extinction of an Island Endemic: A Future for the Corsican Nuthatch

The Corsican Nuthatch (Sitta whiteheadi) is red-listed as vulnerable to extinction by the IUCN because of its endemism, reduced population size, and recent decline. A further cause is the fragmentation and loss of its spatially-restricted favourite habitat, the Corsican pine (Pinus nigra laricio) forest. In this study, we aimed at estimating the potential impact of climate change on the distribution of the Corsican Nuthatch using species distribution models. Because this species has a strong trophic association with the Corsican and Maritime pines (P. nigra laricio and P. pinaster), we first modelled the current and future potential distribution of both pine species in order to use them as habitat variables when modelling the nuthatch distribution. However, the Corsican pine has suffered large distribution losses in the past centuries due to the development of anthropogenic activities, and is now restricted to mountainous woodland. As a consequence, its realized niche is likely significantly smaller than its fundamental niche, so that a projection of the current distribution under future climatic conditions would produce misleading results. To obtain a predicted pine distribution at closest to the geographic projection of the fundamental niche, we used available information on the current pine distribution associated to information on the persistence of isolated natural pine coppices. While common thresholds (maximizing the sum of sensitivity and specificity) predicted a potential large loss of the Corsican Nuthatch distribution by 2100, the use of more appropriate thresholds aiming at getting closer to the fundamental distribution of the Corsican pine predicted that 98% of the current presence points should remain potentially suitable for the nuthatch and its range could be 10% larger in the future. The habitat of the endemic Corsican Nuthatch is therefore more likely threatened by an increasing frequency and intensity of wildfires or anthropogenic activities than by climate change

Evaluation of the potential index model to predict habitat suitability of forest species: the potential distribution of mountain pine (Pinus uncinata) in the Iberian peninsula

Characterization of the suitability or potentiality of a territory for forest tree species is an important source of information for forest planning and managing. In this study, we compared a relatively simple methodology to generate potential habitat distribution areas that has been traditionally used in Spain (the potential index model) with a statistical modelling approach (generalized linear model). We modelled the potential distribution of mountain pine (Pinus uncinata) in the Iberian peninsula as a working example. The potential index model generated a map of habitat suitability according to the values of an index of potentiality, whose distribution has usually divided into four categories based on quartiles (from optimum to low suitability). Considering all values of the index of potentiality as presences of mountain pine resulted in a low to moderate degree of agreement between the potential index model and the generalized linear model according to the kappa coefficient. Using the cut-off value of the index of potentiality that maximized the degree of agreement between both modelling approaches resulted in a substantial similarity between the maps of the predicted distribution of mountain pine. This cut-off value did lie in the upper-third quartile of the potential index distribution (high suitability category), and roughly coincided with the upper 30th percentile. The use of statistical techniques, which have proved to be powerful and versatile for species distribution modelling, is recommended. However, the potential index model, together with the adjustments proposed here, could be a reasonably simple methodology to predict the potential distribution of forest tree species that forest managers should take into account when evaluating forestation and afforestation projects

Tree migration-rates : narrowing the gap between inferred post-glacial rates and projected rates

Faster-than-expected post-glacial migration rates of trees have puzzled ecologists for a long time. In Europe, post-glacial migration is assumed to have started from the three southern European peninsulas (southern refugia), where large areas remained free of permafrost and ice at the peak of the last glaciation. However, increasing palaeobotanical evidence for the presence of isolated tree populations in more northerly microrefugia has started to change this perception. Here we use the Northern Eurasian Plant Macrofossil Database and palaeoecological literature to show that post-glacial migration rates for trees may have been substantially lower (60–260 m yr–1) than those estimated by assuming migration from southern refugia only (115–550 m yr–1), and that early-successional trees migrated faster than mid- and late-successional trees. Post-glacial migration rates are in good agreement with those recently projected for the future with a population dynamical forest succession and dispersal model, mainly for early-successional trees and under optimal conditions. Although migration estimates presented here may be conservative because of our assumption of uniform dispersal, tree migration-rates clearly need reconsideration. We suggest that small outlier populations may be a key factor in understanding past migration rates and in predicting potential future range-shifts. The importance of outlier populations in the past may have an analogy in the future, as many tree species have been planted beyond their natural ranges, with a more beneficial microclimate than their regional surroundings. Therefore, climate-change-induced range-shifts in the future might well be influenced by such microrefugia

Equilibrium of Global Amphibian Species Distributions with Climate

A common assumption in bioclimatic envelope modeling is that species distributions are in equilibrium with contemporary climate. A number of studies have measured departures from equilibrium in species distributions in particular regions, but such investigations were never carried out for a complete lineage across its entire distribution. We measure departures of equilibrium with contemporary climate for the distributions of the world amphibian species. Specifically, we fitted bioclimatic envelopes for 5544 species using three presence-only models. We then measured the proportion of the modeled envelope that is currently occupied by the species, as a metric of equilibrium of species distributions with climate. The assumption was that the greater the difference between modeled bioclimatic envelope and the occupied distribution, the greater the likelihood that species distribution would not be at equilibrium with contemporary climate. On average, amphibians occupied 30% to 57% of their potential distributions. Although patterns differed across regions, there were no significant differences among lineages. Species in the Neotropic, Afrotropics, Indo-Malay, and Palaearctic occupied a smaller proportion of their potential distributions than species in the Nearctic, Madagascar, and Australasia. We acknowledge that our models underestimate non equilibrium, and discuss potential reasons for the observed patterns. From a modeling perspective our results support the view that at global scale bioclimatic envelope models might perform similarly across lineages but differently across regions

Predicting the Impact of Climate Change on Threatened Species in UK Waters

Global climate change is affecting the distribution of marine species and is thought to represent a threat to biodiversity. Previous studies project expansion of species range for some species and local extinction elsewhere under climate change. Such range shifts raise concern for species whose long-term persistence is already threatened by other human disturbances such as fishing. However, few studies have attempted to assess the effects of future climate change on threatened vertebrate marine species using a multi-model approach. There has also been a recent surge of interest in climate change impacts on protected areas. This study applies three species distribution models and two sets of climate model projections to explore the potential impacts of climate change on marine species by 2050. A set of species in the North Sea, including seven threatened and ten major commercial species were used as a case study. Changes in habitat suitability in selected candidate protected areas around the UK under future climatic scenarios were assessed for these species. Moreover, change in the degree of overlap between commercial and threatened species ranges was calculated as a proxy of the potential threat posed by overfishing through bycatch. The ensemble projections suggest northward shifts in species at an average rate of 27 km per decade, resulting in small average changes in range overlap between threatened and commercially exploited species. Furthermore, the adverse consequences of climate change on the habitat suitability of protected areas were projected to be small. Although the models show large variation in the predicted consequences of climate change, the multi-model approach helps identify the potential risk of increased exposure to human stressors of critically endangered species such as common skate (Dipturus batis) and angelshark (Squatina squatina)